Changes of Tumor Infiltrating Lymphocytes after Core Needle Biopsy and the Prognostic Implications in Early Stage Breast Cancer: A Retrospective Study.

PURPOSE: The purpose of this study was to investigate the changes of tumor infiltrating lymphocytes (TILs) between core needle biopsy (CNB) and surgery removed sample (SRS) in early stage breast cancer patients and to identify the correlating factors and prognostic significance of TILs changes.
MATERIALS AND METHODS: A retrospective study was carried out on 255 patients who received CNB and underwent surgical resection for invasive breast cancer. Stromal TILs levels of CNB and SRS were evaluated respectively. Tumors with ≥50% stromal TILs were defined as lymphocyte-predominant breast cancer (LPBC). Clinicopathological variables were analyzed to determine whether there were factors associated with TILs changes. Log-rank tests and Cox proportional hazards models were used to analyze the influences of TILs and TILs changes on survival.
RESULTS: SRS-TILs (median, 10.0%) were significant higher than CNB-TILs (median, 5.0%; p<0.001). Younger age (<60 years, p=0.016) and long surgery time interval (STI, ≥4 days; p=0.003) were independent factors correlating with higher TILs changes. CNB-LPBC patients showed better breast cancer-free interval (BCFI, p=0.021) than CNB-non-LPBC (CNB-nLPBC) patients. Patients were categorized into four groups according to the LPBC change pattern from CNB to SRS: LPBC→LPBC, LPBC→nLPBC, nLPBC→LPBC, and nLPBC→nLPBC, with estimated 5-year BCFI 100%, 100%, 69.7%, and 86.0% (p=0.016). nLPBC→LPBC pattern was an independent prognostic factor of worse BCFI (hazard ratio, 2.19; 95% confidence interval, 1.06 to 4.53; p=0.035) compared with other patterns.
CONCLUSION: TILs were significantly higher in SRS than in CNB. Higher TILs changes were associated with younger age and long STI. Changing from nLPBC to LPBC after CNB indicated a worse BCFI, which needs further validation.

Introduction

Tumor immune microenvironment plays an important role in tumor development, progression, and control. Effective immune surveillance by the host immune system protects the body against cancer. In breast cancer, lymphocyte infiltration was detectable in more than 80% of tumors [1-3], and accumulating evidence from several recent large clinical studies indicates that high levels of tumor infiltrating lymphocytes (TILs) are predictive of better response to neoadjuvant chemotherapy, and increased survival [2-9]. These results suggest that tumors with different TILs levels may have different biologic behaviors and different treatment responses, and there is a role to evaluate TILs in clinical trial cohorts as well as in daily histopathological practice. However, as a quantitative histological biomarker, the utilization of TILs depends on the standardization and reproducibility of TILs evaluation. A guideline for evaluation of TILs in breast cancer has been published in 2014 [10], and two multicenter international ring studies have been conducted to make TILs a reproducible marker [11].

Neoadjuvant therapy (NAT) is an increasingly used therapeutic strategy for initially large, inoperable, or locally advanced breast tumors, as well as operable cancers susceptible to specific treatments. NAT not only induce clinical responses that allow breast-conserving surgery and better survival, but also allows observation of tumor response to specific regimen. Sequential biopsies of the same tumor to be taken at different time points during NAT, allows assessment of biologic changes in the tumor long before clinical evidence of response can be determined. The changes of biomarkers like Ki67 after a short exposure to pre-surgery therapy has been found to be related with survival [12,13]. Post-treatment infiltration of lymphocytes has also been found to be associated with breast cancer outcome [14-18]. Early detection of changes of TILs during NAT may play a potential role in treatment response evaluation for breast cancer patients.

Core needle biopsy (CNB) is routinely used in breast cancer biopsy. Discordance of biomarkers, including estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), and Ki67 exist between CNB and surgery removed samples (SRS) [19,20]. Tumor heterogeneity, the process of specimen fixation and the procedure of CNB can cause discordance of biomarkers evaluation between CNB and SRS [21-23]. If we observed changes of biomarkers in sequential biopsies, we need to consider whether the variations have occurred as a result of treatment, but not for other reasons. We need to confirm that the information we got from the core biopsies can represent the whole tumor. However, there is limited data on whether TILs status in CNBs can represent those in the whole tumor section. To solve this question, the present study was designed to compare the TILs status between sequential CNB samples and SRS of breast cancer patients who received surgery soon after CNB without NAT, to find which factors are related with TILs changes, and to further explore the prognostic significance of TILs and TILs changes.

Materials and Methods

1. Patients

We retrospectively analyzed patients who received CNB and underwent the following surgical resection for breast cancer between January 2009 and December 2011 in Comprehensive Breast Health Center, Shanghai Ruijin Hospital. The CNB procedure was done under the guide of ultrasound and at least four 14-gauge CNB samples were collected for pathological examination.

Patient who met all the following criteria were included in this study: (1) female sex; (2) found invasive carcinoma in both CNB and SRS; (3) no preoperative anti-cancer therapy; and (4) paired CNB and SRS sections available for TILs evaluation. Patients with stage IV breast cancer, with bilateral breast cancer or with other malignant diseases at the time of breast cancer diagnosis were excluded from this study.

The baseline clinicopathological characteristics of the patients were achieved from Shanghai Jiaotong University Breast Cancer Database, including age, menopausal status, performing date of CNB and definite excision surgery, surgery type, and histopathological information. Patients received adjuvant treatments according to physicians' choices. Adjuvant treatment information and survival status were obtained through outpatient medical history and/or phone calls.

2. Pathologic assessment

The common histopathological information, including breast cancer pathological type, histological grade, tumor size, lymph node status, immunohistochemical (IHC) results of ER, PR, HER2, and Ki67, were achieved from patients' medical records. ER positivity (ER+) and PR positivity (PR+) were defined as 1% or more positive invasive tumor cells with nuclear staining. Hormonal receptor positivity (HR+) was defined as ER+ or PR+. HER2 status was firstly determined by IHC staining and scored as 0, 1+, 2+, 3+ according to the 2007 American Society of Clinical Oncology/College of American Pathologists (ASCO/CAP) guideline. Samples with HER2 IHC 2+ were further examined by fluorescence in situ hybridization (FISH). Tumors with HER2 IHC 3+ or FISH positive were regarded as HER2 positive. Five breast cancer molecular subtypes were classified according to the 2013 St. Gallen breast cancer consensus [24]: luminal A (ER+/HER2–, Ki67 < 14% and PR ≥ 20%), luminal B-HER2– (ER+/HER2–, Ki67 ≥ 14%, or ER+/HER2–, PR < 20%, or ER–/PR+/HER2–), luminal B-HER2+ (HR+/HER2+), triple-negative (TN, HR–/HER2–), and HER2+ (HR–/HER2+).

Evaluations of stromal TILs were performed on readymade full-face hematoxylin and eosin stained sections of paraffin-embedded CNB and SRS according to the 2014 recommendation of international TILs working group [10] by an experienced pathologist (X.C. Fei). TILs were reported for the stromal compartment (=% stromal TILs).

TILs levels were evaluated separately for CNB and SRS, recorded as CNB-TILs and SRS-TILs. Tumors with ≥ 10% stromal TILs were regarded as high-TILs tumors, while tumors with < 10% TILs were regarded as low-TILs tumors. Tumors with ≥ 50% stromal TILs were defined as lymphocyte-predominant breast cancer (LPBC) [2-5,10], while tumors with < 50% TILs were defined as non-lymphocytepredominant breast cancer (nLPBC). Tumors were categorized into CNB-High-TILs, CNB-Low-TILs, SRS-High-TILs, SRS-Low-TILs, CNB-LPBC, SRS-LPBC, CNB-nLPBC, and SRS-nLPBC according to CNB-TILs and SRS-TILs status.

3. Statistical analysis

Change of stromal tumor infiltrating lymphocytes (ΔTILs) was defined as the value of SRS-TILs minus CNB-TILs of each patient. Surgery time interval (STI) was defined as the day interval between CNB and definite surgery. Breast cancer-free interval (BCFI) was defined as the interval from the first day after surgery to the date of any of the following events happened: breast cancer recurrence at local or regional sites, metastasis to distant sites or newly diagnosed breast cancer in the contralateral breast. Overall survival (OS) was defined as time from the date of definite surgery to the date of death from any cause.

Kappa test was used to calculate the concordance rate of ER, PR, HER2, Ki67, and TILs between CNB and SRS. Values of κ > 0.6 were correlated with good agreement, values between 0.4 and 0.6 considered moderate agreements, values < 0.4 corresponded to fair, and values < 0.2 reflected poor agreement. TIL change after CNB was compared by using the Wilcoxon signed-rank test. Mann-Whitney U test and Kruskal-Wallis test were used to calculate the association between CNB-TILs, SRS-TILs, and ΔTILs with clinicopathological factors. Multivariate ANOVA analysis was then done to find the association and interaction between the rank of ΔTILs and the tested factors. Kaplan-Meier method was used to estimate survival curves and log-rank test was used to compare survival between groups. Multivariate Cox proportional hazards models were used to estimate the hazard ratio, 95% confidence interval (CI), and the effects of clinicopathological variables. All statistical tests were two-sided and p < 0.05 was considered as statistically significant. All the analyses were done by SPSS ver. 23.0 for Mac (IBM Corp., Armonk, NY).

4. Ethical statement

All procedures performed in studies involving human participants were in accordance with the ethical standards of the independent Ethical Committee/Institutional Review Board of Ruijin Hospital and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Informed consent was exempted due to the retrospective format of this study.

Results

1. Patient characteristics

The clinicopathological characteristics of the 255 patients included in this study are detailed in Table 1. The median patients' age was 56 years (range, 28 to 86 years). Mastectomy and breast-conserving surgery were performed in 223 (87.5%) and 32 (12.5%) patients, respectively. Tumor size was pT1 (≤ 2 cm) in 41.6% and pT2 (2-5 cm) in 54.5% patients. The majority (94.1%) of the patients had an invasive ductal carcinoma and 43.1% of the patients had a grade III tumor. Axillary lymph node status was pN0 in 51.4% and pN1 (1-3 positive nodes) in 29.4% of the cases; the remaining 19.2% of the patients had four or more metastatic lymph nodes. Most (74.1%) of the patients had an ER-positive tumor and 20% had a HER2-positive tumor. Median Ki67 was 20%. There were 14.1%, 51.0%, 9.8%, 10.2%, and 14.9% patients classified as luminal A, luminal B-HER2–, luminal B-HER2+, HER2+, and TN subtype, respectively. Median STI was 4 days, ranging from 1 to 32 days. Forty point eight percent of patients received surgery within 1-3 days after CNB. Seventy-six point one percent of the patients received chemotherapy; 40.4% received radiotherapy; 74.5% received endocrine therapy; 11.8% received trastuzumab as targeted therapy, accounting for 58.8% (30/51) of all HER2-positive patients.

Table 1.

Baseline patient characteristics

Characteristic	No. (%) (n=255)
Age (yr)
< 40	16 (6.3)
40-49	60 (23.5)
50-59	78 (30.6)
60-69	61 (23.9)
≥ 70	40 (15.7)
Menstrual status
Premenopausal	86 (33.7)
Postmenopausal	169 (66.3)
Breast surgery type
Mastectomy	223 (87.5)
Breast-conserving surgery	32 (12.5)
Pathological type
Invasive ductal carcinoma	241 (94.9)
Invasive lobular carcinoma	9 (3.5)
Others	5 (1.6)
Histological grade
I-II	132 (51.8)
III	110 (43.1)
NA	13 (5.1)
Tumor size (cm)
≤ 2	106 (41.6)
2-5	139 (54.5)
> 5	10 (3.9)
Axillary lymph node status
Negative	131 (51.4)
1-3 Positive	75 (29.4)
4-9 Positive	24 (9.4)
> 10 Positive	25 (9.8)
ER
Negative	66 (25.9)
Positive	189 (74.1)
PR
Negative	93 (36.5)
Positive	162 (63.5)
HER2
Negative	204 (80.0)
Positive	51 (20.0)
Ki67 (%)
< 20	86 (33.7)
≥ 20	169 (66.3)
Molecular subtype
Luminal A	36 (14.1)
Luminal B-HER2–	130 (51.0)
Luminal B-HER2+	25 (9.8)
HER2 enriched	26 (10.2)
Triple negative	38 (14.9)
Surgery time interval (day)
< 4	104 (40.8)
≥ 4	151 (59.2)
Chemotherapy
Yes	194 (76.1)
No	61 (23.9)
Radiotherapy
Yes	103 (40.4)
No	152 (59.6)
Targeted therapy (trastuzumab)
Yes	30 (11.8)
No	225 (88.2)
Endocrine therapy
Yes	190 (74.5)
No	65 (25.5)

Pathological data is based on surgery removed sample.

ER, estrogen receptor; PR, progesterone receptor; HER2, human epidermal growth factor receptor 2.

2. Association of TILs with clinicopathological characteristics

The distribution of CNB-TILs and SRS-TILs by different subtypes is shown in Fig. 1A. Overall, the median value of CNB-TILs was 5.0% (interquartile range [IQR], 1.0% to 20.0%), and 10.2% of the patients had CNB-LPBC (TILs ≥ 50% in CNB samples). CNB-TILs were correlated with molecular subtypes (p < 0.001), high histological grade (p < 0.001), ER negativity (p < 0.001), and high Ki67 expression (≥ 20%, p=0.003). The median values of CNB-TILs of luminal A, luminal B-HER2–, luminal B-HER2+, HER2+, and TN were 1.0% (IQR, 0.0% to 5.0%), 5.0% (IQR, 1.0% to 16.3%), 8.0% (IQR, 1.0% to 25.0%), 8.0% (IQR, 2.8% to 16.3%), and 20.0% (IQR 5.0 to 50.0), respectively (p < 0.001). CNB-LPBC phenotype comprised 2.8%, 6.9%, 12.0%, 11.5%, and 26.3% in luminal A, luminal B-HER2–, luminal B-HER2+, HER2+, and TN subtypes, respectively (p=0.004). CNB-TILs were not statistically associated with age, menopausal status, lymph node involvement, tumor size, PR status or HER2 status (p > 0.05).

Fig. 1.

Distribution of TILs and ΔTILs. (A) Distribution of CNB-TILs% and SRS-TILs% in the global population and per breast cancer subtype. The box shows the limits of the middle half of the data (1st to 3rd quartile), and the line inside represents the median. Whiskers are drawn to represent values within 1.5 times the length of the box from either end of the box (1.5 times the interquartile range). Open circles correspond to mild outliers. Asterisks correspond to extreme outliers. (B) Distribution of CNB-TILs and SRS-TILs per patient. Each open circle represents one patient; (C) Distribution of ΔTILs (per 10%). TILs, tumor infiltrating lymphocytes; ΔTILs, change of TILs; CNB, core needle biopsy; SRS, surgery removed sample, Max, maximum; Min, minimum; Q, quartile.

The median value of SRS-TILs was 10.0% (IQR, 5.0% to 35.0%), and 22.0% of the patients had SRS-LPBC (TILs ≥ 50% in SRS). SRS-TILs were associated with molecular subtype (p < 0.001), high histological grade (p < 0.001), ER negativity (p < 0.001), PR negativity (p=0.036), high Ki67 expression (≥ 20%, p < 0.001), younger age (< 60 years, p=0.014), and long STI (≥ 4 days, p=0.006). The median values of SRS-TILs of luminal A, luminal B-HER2–, luminal B-HER2+, HER2+, and TN were 5.0% (IQR, 3.0% to 15.0%), 10.0% (IQR, 5.0% to 30.0%), 10.0% (IQR, 4.0% to 55.0%), 20.0% (IQR, 7.3% to 56.3%), and 35.0% (IQR, 18.8 to 60.0), respectively (p < 0.001). SRS-LPBC phenotype comprised 5.6%, 14.6%, 28.0%, 38.5%, and 47.4% in the luminal A, luminal B-HER2–, luminal B-HER2+, HER2+, and TN subtypes, respectively (p < 0.001). SRS-TILs were not statistically associated with menopausal status, lymph node involvement, tumor size, or HER2 status (p > 0.05).

3. Comparison of ER, PR, HER2, Ki67, and TILs between CNB and SRS

The concordance rates of ER, PR, and HER2 between CNB and SRS were 94.1% (κ=0.841), 87.1% (κ=0.724), and 97.3% (κ=0.911), respectively. CNB and SRS had a good concordance in ER, PR, and HER2.

Ki67 value was much higher in SRS comparing with CNB samples by using two-paired samples t test, with mean values of 29.5% and 26.4%, respectively (p < 0.001). Take ≥ 20% as the cutoff value of high Ki67, the concordance rate of Ki67 was 82.4% (κ=0.635).

The distribution of CNB-TIL and SRS-TIL of each patient was shown in Fig. 1B. The median and mean value of CNB-TILs was 5.0% and 13.8%, while in SRS-TILs was 10.0% and 22.3%. SRS-TILs were significantly higher than CNB-TILs by using the Wilcoxon test (p < 0.001). Respectively use ≥ 10%, ≥ 20%, ≥ 30%, ≥ 40%, ≥ 50%, ≥ 60%, and ≥ 70% as the cutoff value of high TILs and low TILs, the concordance rates between CNB-TILs and SRS-TILs were 75.3% (κ=0.522), 82.4% (κ=0.607), 82.4% (κ=0.554), 85.5% (κ=0.525), 86.3% (κ=0.500), 91.0% (κ=0.579), and 98.4% (κ=0.593). The agreement of TILs between CNB and SRS was just moderate.

ΔTILs ranged from –45.0% to +79.0% (Fig. 1C), with median value of 1.0%, and mean value of 8.6%. High ΔTILs was correlated with younger age (< 60 years, p=0.016), ER negativity (p=0.006), and long STI (≥ 4 days, p=0.003) (Table 2). Multivariate ANOVA analysis showed that younger age (< 60 years, p=0.007), and long STI (≥ 4 days, p=0.004) were still significantly associated with ΔTILs.

Table 2.

ΔTILs among patients with different clinicopathological characteristics

Characteristic	ΔTILs (%)		p-value
	Median (IQR)	Mean±SD
Age (yr)
< 60	2.0 (0.0-17.0)	11.1±19.6	0.016
≥ 60	0.0 (0.0-5.0)	4.9±13.3
Menstrual status
Premenopausal	0.0 (0.0-15.0)	11.2±20.5	0.888
Postmenopausal	1.0 (0.0-10.0)	7.24±15.8
Tumor size (cm)
≤ 2	0.0 (0.0-10.0)	8.0±17.8	0.517
> 2	1.0 (0.0-12.0)	9.0±17.5
Axillary lymph node status
Negative	0.0 (0.0-10.0)	8.7±17.3	0.973
Positive	1.0 (0.0-12.0)	8.4±18.0
Histological grade
I-II	0.0 (0.0-10.0)	7.2±15.3	0.360
III	2.0 (0.0-20.0)	10.4±20.5
NA	2.0 (0.0-12.5)	7.2±11.7
ER
Negative	6.0 (0.0-30.0)	15.3±24.5	0.006
Positive	0.0 (0.0-7.0)	6.2±13.7
PR
Negative	2.0 (0.0-17.0)	10.5±20.0	0.110
Positive	0.0 (0.0-10.0)	7.4±16.0
HER2
Negative	1.0 (0.0-10.0)	7.5±15.2	0.602
Positive	0.0 (0.0-27.5)	13.0±24.7
Ki67 (%)
< 20	0.0 (0.0-5.0)	4.8±11.5	0.062
≥ 20	2.0 (0.0-15.0)	10.5±19.8
Molecular subtype
Luminal A	2.0 (0.0-13.0)	6.8±11.7	0.074
Luminal B-HER2–	0.0 (0.0-7.0)	6.1±13.4
Luminal B-HER2+	0.0 (0.0-14.5)	8.3±21.2
HER2 positive	8.5 (0.0-38.5)	17.5±27.4
Triple negative	5.0 (0.0-20.5)	13.0±21.7
Surgery time interval (day)
< 4	0.0 (0.0-5.0)	4.9±15.4	0.003
≥ 4	2.0 (0.0-16.0)	11.1±18.6

Pathological data is based on surgery removed sample. p-value calculated by Mann-Whitney U test. TILs, tumor infiltrating lymphocytes; IQR, interquartile range; SD, standard deviation; ER, estrogen receptor; PR, progesterone receptor; HER2, human epidermal growth factor receptor 2.

4. Association between TILs and prognosis

After a median follow-up of 68 months, there were 44 BCFI events being recorded. Twenty-two deaths were recorded in our study until the last follow-up. Fifteen patients died of breast cancer and the other seven died from other reasons.

At univariate survival analysis, CNB-LPBC was a significant predictor of better BCFI compared with CNB-nLPBC patients (estimated 5-year BCFIrate, 100% vs. 83.6%; p=0.021), while no significant BCFI difference was observed by CNB-High TILs or CNB-Low TILs (p=0.215), SRS-High TILs or SRS-Low TILs (p=0.457), SRS-LPBC or SRS-nLPBC (p=0.752).

Patients were categorized into four groups according to the LPBC change pattern from CNB to SRS: LPBC→LPBC (n=23), LPBC→nLPBC (n=2), nLPBC→LPBC (n=33), and nLPBC→nLPBC (n=197). Details of TILs value of CNB and SRS according to LPBC change pattern was listed in Table 3. The estimated 5-year BCFI of patients with LPBC→LPBC, LPBC→nLPBC, nLPBC→LPBC, and nLPBC→nLPBC pattern was 100%, 100%, 69.7%, and 86.0%, respectively (p=0.016) (Fig. 2A). nLPBC→LPBC patients had a worse BCFI than others, p=0.013 (Fig. 2B). We also group patients according to ΔTILs values (ΔTILs < 0% vs. ≥ 0%, < 5% vs. ≥ 5%, < 10% vs. ≥ 10%, < 20% vs. ≥ 20%, and < 30% vs. ≥ 30%) and found no statistical difference of BCFI between groups (p > 0.05, respectively).

Table 3.

CNB-TILs, SRS-TILs and ΔTILs according to LPBC change pattern

LPBC change pattern	CNB-TILs%		SRS-TILs%		ΔTILs (%)
	Median (IQR)	Mean±SD	Median (IQR)	Mean±SD	Median (IQR)	Mean±SD
nLPBC→nLPBC (n=197)	5.0 (1.0 to 10.0)	7.5±9.5	8.0 (4.0 to 15.0)	11.3±10.7	0.0 (0.0 to 5.0)	3.8±8.7
nLPBC→LPBC (n=33)	10.0 (2.5 to 27.5)	1.4±12.9	60.0 (50.0 to 65.0)	59.2±7.9	45.0 (30.0 to 59.0)	44.7±16.5
LPBC→LPBC (n=23)	60.0 (60.0 to 70.0)	63.3±8.3	65.0 (60.0 to 70.0)	65.4±8.8	0.0 (0.0 to 0.0)	2.2±5.8
LPBC→nLPBC (n=2)	55.0 (55.0 to 55.0)	55.0± 7.1	10.0 (5.0 to 10.0)	10.0±7.1	–45.0 (–45.0 to –45.0)	–45.0±0.0
Total	5.0 (1.0 to 20.0)	13.8±19.1	10.0 (5.0 to 35.0)	22.3±23.3	1.0 (0.0 to 10.0)	8.6±17.6

CNB, core needle biopsy; TILs, tumor infiltrating lymphocytes; SRS, surgery removed sample; LPBC, lymphocyte-predominant breast cancer; IQR, interquartile range; SD, standard deviation; nLPBC, non-lymphocyte-predominant breast cancer.

Fig. 2.

Breast cancer-free interval by LPBC change status in overall population (A, B), luminal A (C), luminal B-HER2– (D), luminal B-HER2+ (E), HER2+ (F), and TNBC (G) patients. LPBC, lymphocyte predominant breast cancer; nLPBC, non-lymphocyte-predominant breast cancer; HER2, human epidermal growth factor receptor 2; TNBC, triple-negative breast cancer.

In subgroup analysis by molecular subtypes, nLPBC→ LPBC pattern was associated with worse BCFI in HER2-enriched subtype (p=0.038) but not in luminal A, luminal B-HER2–, luminal B-HER2+, or TN subtypes (Fig. 2C-G). The difference of BCFI according to CNB-LPBC status, SRS-LPBC status, CNB-TILs status (high vs. low) and SRS-TILs (high vs. low) status was not seen in any of the five molecular subtypes.

We generate a multivariate Cox proportional hazards model including tumor size, age, nodal status, histological grade, STI, molecular subtype, CNB-LPBC status, and LPBC change pattern (Table 4) and found LPBC change pattern of nLPBC→LPBC (hazard ratio, 2.19; 95% CI, 1.06 to 4.53; p=0.035) and tumor size > 2 cm (hazard ratio, 2.85; 95% CI, 1.36 to 6.00; p=0.006) were independent prognostic factor of worse BCFI.

Table 4.

Multivariate survival analysis of BCFI

Characteristic	BCFI
	p-value	Hazard ratio (95% CI)
Age (yr)
< 60	0.156	1
≥ 60		1.58 (0.84-2.96)
Tumor size (cm)
< 2	0.006	1
≥ 2		2.85 (1.36-6.00)
Axillary lymph node status
Negative	0.230	1
Positive		1.45 (0.79-2.66)
Histological grade
I-II	0.072	1
III		1.83 (0.95-3.51)
NA		3.01 (0.99-9.09)
Molecular subtype
Luminal A	0.857
Luminal B-HER2–		1.59 (0.46-5.49)
Luminal B-HER2+		2.21 (0.55-8.91)
HER2 positive		1.83 (0.42-7.92)
Triple-negative		0.62 (0.33-1.17)
Surgery time interval (day)
< 4	0.211	1
≥ 4		0.67 (0.36-1.25)
CNB-LPBC status
nLPBC	0.969	1
LPBC		0.00 (0.00-8.85279)
LPBC change pattern
nLPBC→nLPBC, LPBC→LPBC, or LPBC→nLPBC	0.035	1
nLPBC→LPBC		2.19 (1.06-4.53)

Pathological data is based on surgery removed sample if not specified. BCFI, breast cancer-free interval; CI, confidence interval; NA, not available; HER2, human epidermal growth factor receptor 2; CNB, core needle biopsy; LPBC, lymphocytepredominant breast cancer; nLPBC, non-lymphocyte-predominant breast cancer.

There were limited OS events in our study. No significant OS difference was detected between groups by different clinicopathological variables including TILs.

Discussion

We have already known that variations of biomarkers exist between CNB and SRS, which is supported by multiple studies that investigated the concordance of ER, PR, HER2, and Ki67 status in CNB compared to SRS [19,20]. Our results showed a similar concordance of ER, PR, HER2, and Ki67 with our previous study, in which CNB was accurate in evaluating ER, PR, and HER2 while Ki67 value significantly increased after CNB [20]. Few studies pay attention to the variability of immune profiles between CNB and SRS. Jeselsohn et al. [25] have reported activation of immune-related genes after CNB in a retrospective study evaluating 21 paired core biopsies and excisions from breast cancer patients who did not receive NAT. However, animal experiments conducted by Mathenge et al. [26] showed that tumors from biopsied mice contained a lower frequency of CD4+ T cells, CD8+ T cells, and macrophages than that of un-biopsied mice, suggesting biopsy-mediated development of an increasingly immunosuppressive tumor microenvironment [26]. In this study, we found that TILs were higher in SRS than in CNB. The reasons why TILs changes after CNB were hypothetical: (1) with limited evidence, CNB was considered to associate with activation of immune response in breast cancer [25], which indicates the difference between CNB and SRS may be caused by the process of CNB. (2) Core biopsy sample is supposed to contain more tissue from the periphery of the tumor and fewer tissue from the center of the tumor, although knowledge of which area of the tumor the needle targets, is difficult to elucidate and might be regarded as random. Malignant tumors usually develop hypoxic and necrotic cores with higher microvascular density [27], which may result in more immune infiltration in tumor center than in the periphery area. This hypothesis indicates that it might be intra-tumor heterogeneity that caused the difference of TILs between CNB and SRS samples.

We found younger age and long STI were correlated with higher TILs changes after CNB. The level of TILs was more stable in elder patients. Aging is associated with the decrease of autoimmune function, due to the degeneration of thymus gland, the decreased export and replication of T cells, and the reduced immune cell function [28]. For the above reasons, elder patients may have insufficient immune activation after CNB compared with younger patients so that TILs may remain stable. We also assumed that the tumor biological behavior is always milder in elderly patients and tend to generate a more homogeneous tumor, which result in a similar TILs between CNB and SRS. Longer STI was also an indicator of higher chance of TILs increase after CNB in our study. It indicated that the change of TILs after CNB can be dynamic and contribute to the hypothesis that CNB process caused TILs changes. Further studies are required to verify he difference we found between CNB and SRS and to elucidate the underlying mechanism.

The prognostic value of TILs had been studied in a number of trials. Meta-analysis has indicated that TILs were not prognostic markers for disease-free survival (DFS) and OS in overall population, but related to improved DFS and OS in triple-negative breast cancer (TNBC) [29]. In our study, we failed to detect the association between CNB-TILs or SRS-TILs and survival, neither in the overall population nor in the TNBC subgroup. At this point, we cannot conclude whether CNB-TILs and SRS-TILs have equal value in survival prediction and further validation is needed. This is also the first study assessing the correlation between changes of TILs after CNB and survival. We found that patients changed from nLPBC to LPBC after CNB experienced a worse BCFI than patients with other LPBC change patterns, especially in HER2+ patients. This result was quite different from the common view that TILs increase is related with better survival [14,16,18], although the related studies were mainly conducted in TNBC. In HER2+ tumor, a recent study showed spatial heterogeneity of lymphocytic infiltrate was associated with poor prognosis [30], which is similar with our study. The difference of LPBC status between CNB and SRS may reflect strong regional segregation of immune cells and its association with poor prognosis may indicate dysfunctional immune infiltrate. These results indicated that when we do research about the clinical significance of TILs change during NAT, we need to take account of the difference of TILs between CNB sample and the whole tumor at the time of biopsy.

There were several limitations of this study. The sample size was too small especially when we analyze patients in different subgroups. In addition, it was a retrospective study, the treatments of patients were not standardized, which may affect the result of survival analysis. At last, we only analyzed total TILs without analyzing TILs subsets and were lack of mechanism study. Our next step is to enlarge the sample size and initiate a prospective cohort to validate our results, to figure out the detailed composition of TILs subsets and the changes of gene profiles to study the underlying mechanism.

In summary, our study demonstrated that for patients who received CNB before surgery, TILs were significantly higher in SRS than in CNB. ΔTILs were associated with younger age and long STI. Patients who change from nLPBC to LPBC after CNB tend to have the worst BCFI among the study population, which warrants further validation and further mechanism studies.