{'fraction_non_alphanumeric': 0.03197776821010607, 'fraction_numerical': 0.015577909115816666, 'mean_word_length': 5.177466150870406, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 1, '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}}
Joint Tumor Bud–MMP/TIMP Count at the Invasive Front Improves the Prognostic Evaluation of Invasive Breast Carcinoma

Background: Tumor budding is a histological phenomenon consisting of the formation of small clusters of one to five undifferentiated malignant cells detached from the main tumor mass which are observed in the tumor stroma. In the present study, we investigated the prognostic significance of tumor budding in breast cancer and its relationship with the expressions of matrix metalloproteases (MMPs) and their tissue inhibitors (TIMPs). Methods: The number of buds was counted in whole-tissue sections from 153 patients with invasive ductal carcinomas who underwent a long follow-up period. In addition, an immunohistochemical study of MMP-9, -11, and -14 TIMP-1 and -2 expression by cell types at the invasive tumor front was carried out. Results: There was a wide variability in the number of buds among tumors, ranging from 0 to 28 (median = 5). Tumor budding count ≥ 4 was the optimal cut-off to predict both relapse-free and overall survival. High-grade tumor budding was associated with MMP/TIMP expression by cancer-associated fibroblasts. In addition, we found that the combination of tumor budding grade with MMP/TIMP expression by stromal cells, and especially with MMP-11 expression by mononuclear inflammatory cells, significantly improved the prognostic evaluation. Conclusion: High-grade tumor budding is associated with a more aggressive tumor phenotype, which, combined with MMP/TIMP expression by stromal cells at the invasive front of the tumor, identifies patients with poor prognosis.


Introduction
From a histological and molecular point of view, breast cancer is a heterogeneous disease, with around 30% of patients developing metastasis [1]. Although molecular tests are being used to assess the likelihood of treatment response and/or recurrence, additional prognosticators are needed to enhance personalized treatment and, especially, to overcome the over-and undertreatment of patients. Tumor budding is a morphological phenomenon found in various cancers. This histological finding consists of single or small clusters of one to five undifferentiated malignant cells detached from the main tumor mass, which are seen in the tumor stroma in close proximity ahead of the invasive front of a tumor (reviewed in by the authors of [2,3]). The method of scoring tumor budding, based in hematoxylin and/or eosin-stained sections and/or pancytokeratin staining, is a single and low-cost method with high reproductivity [4].
Colorectal cancer was the first cancer type in which tumor budding was addressed systematically. In fact, tumor budding is an additional prognostic factor for this tumor type according to the Union for International Cancer Control (UICC) [5], and a potential prognostic factor in its early stage according to the European Society for Medical Oncology consensus guidelines [6]. More recently, some of these associations were also found in esophageal [7], laryngeal [8], tongue squamous cell carcinoma [9], gingival buccal complex squamous cell carcinoma [10], and lung [11], pancreatic [12], bladder [13], and gastric carcinomas [14]. In invasive ductal breast cancer, high budding has been correlated with clinical-pathological parameters such as larger tumor size [15], lymph vessel invasion [15][16][17][18], and nodal metastasis [16][17][18], as well as shorter survival [15,17,18]. However, the definitive implementation of tumor budding into clinical practice of breast cancer is currently limited by the vast heterogeneity in its exact definition, methodology of assessment, and patient stratification. Therefore, in order to better characterized this histological finding, it may be relevant to integrate tumor budding in its biological context as part of the invasive tumor front. Interestingly, it has been suggested that cells associated with tumor budding share similar properties with malignant stem cells [19], as well as a manifestation of one hybrid epithelial/mesenchymal phenotype displaying collective cell migration [3]. Apart from this last event, disruption of the barriers is a prerequisite for invasion in breast cancer, and the matrix metalloproteases (MMPs) play a key role in this process [20]. In fact, they have been clinically associated with metastasis development in many tumors, including breast cancer.
In the present study, we investigated the prognostic signification of tumor budding grade in patients with breast cancer who underwent a long follow-up period. In addition, we evaluated the relationship of this histological finding with the expression of MMPs and their tissue inhibitors (TIMPs) at the invasive front of the tumors in order to improve the prognostic evaluation of breast cancer based just on the morphological context.

Patients
This study comprised 153 women with a histologically confirmed diagnosis of invasive ductal breast cancer treated between 1990 and 2001, some of which were previously included in our preliminary studies on the expression of MMPs and TIMPs in breast cancer [21][22][23][24]. We selected women with the following inclusion criteria: Invasive ductal carcinoma, at least 6 histopathologically assessed axillary lymph nodes, and a minimum of 10 years of follow-up without tumor recurrence. The exclusion criteria were the following: Metastatic disease at presentation, prior history of any type of malignant tumor, bilateral breast cancer at presentation, having received any type of neoadjuvant therapy, the development of locoregional recurrence during the follow-up period, or the development of a second primary cancer. We randomly selected a sample size of 153 patients, in accordance with 4 different groups of similar size, stratified with regard to nodal status and the development or not of metastatic disease. This was in order to include a sufficient number of patients with recurrence for securing the statistical power of the survival analysis. Note that approximately half of the cases with distant metastasis during the follow-up period occurred in each of the node-negative and node-positive subgroups. Patients' characteristics included in the 2 main groups, with or without distant metastases (recurrence), are listed in Table 1. Patients underwent either modified radical mastectomy or wide resection with axillary lymphadenectomy. The median follow-up period was about 187 months in patients without metastasis and 52 months in patients with metastatic disease. The study adhered to national regulations and was approved by our Institution's Ethics and Research Committee. Women were treated according to the guidelines used in our Institution (Fundación Hospital de Jove). Written informed consent was obtained from all patients and controls.

Definition of Tumor Budding
Tumor buds are defined as a small number of cells (up to 5) which have detached from the bulk of the tumor and are observed as isolated cells or small clusters of cells in histologic sections. In this study, we analyzed peritumoral buds defined as those observed in areas near the margin of the tumor at the invasive tumor front [25]. The number of buds was counted in whole-tissue sections of surgical resection specimens. Prior to the selection of the areas where the count was carried out, a scan of the entire peritumoral area was made, and the areas with the highest budding density were chosen. The number of buds was counted in the hematoxylin and eosin-stained sections with the maximal invasive region and in areas with the highest concentration of buds ("hotspots"), as considered by other authors [15][16][17][18]. The areas with the highest budding density were chosen, but if 1 case had fewer high-power fields (HPFs: 0.55 mm) with any buds, the HPFs without buds (n = 0) were chosen to achieve 10 evaluated HPFs. The number of buds was highly variable along the invasive front, and consequently, the final number obtained was the average of the areas with the highest budding density [26]. Pan-cytokeratin immunostaining (antibody: AE1/AE3) was carried out in exceptional cases when it was difficult to distinguish budded tumor cells from fibroblasts or inflammatory cells on morphological criteria. We used predefined criteria for the assessment of tumor buds in accordance with recent publications describing the scoring of tumor buds in colon cancer [26,27]. A breast pathologist (LOG) scored 10 high-power fields (0.55 mm) blindly without knowledge of tumor characteristics or the results of the present studies. In a recent report [15], tumors were considered to have a high tumor budding if the average number of tumor buds in 10 HPF was > 4. In contrast, tumors were considered to have a low tumor budding if the average number of buds in 10 HPF was ≤ 4. Figure 1 shows representative examples of either low or high tumor budding.

Tissue Arrays
Breast carcinoma tissue samples were obtained at the time of surgery. Tumor tissue array (TA) blocks were obtained by punching a tissue cylinder (core) with a diameter of 1.5 mm through a histologically representative area of each 'donor' tumor block (routinely fixed (overnight in 10% buffered formalin), paraffin-embedded tumor samples), using a manual tissue arrayer (Beecker Instruments, Sun Praerie, WI, USA) as described elsewhere [22]. A total of 2 cores was employed for each case, which corresponded to the invasive front and have been shown to correlate well with conventional immunohistochemical staining [22]. The invasive front was defined as the area within 2 mm surrounding the tumors which contained cancerous cells.

Immunohistochemistry
Histopathologically representative tumor areas were defined in hematoxylin-and eosin-stained sections and marked on the slide. TA blocks were obtained from areas of nonnecrotic cancerous tissues. Serial 5 µm sections of the high-density TA blocks were consecutively cut with a microtome (Leica Microsystems GmbH, Wetzlar, Germany) and transferred to adhesive-coated slides. One section of each tissue array block was stained with H&E, and these slides were then reviewed to confirm that the sample was representative of the original tumor area. Immunohistochemical staining was carried out on these sections using a TechMate TM50 autostainer (Dako, Glostrup, Denmark). Antibodies for MMPs and TIMPs were obtained from Neomarker (Lab Vision Corporation, Fremont, CA, USA). The dilution for each antibody was established based on negative and positive controls (1/50 for MMP-14 and TIMP-2, 1/100 for MMP-9 and TIMP-1, and 1/200 for MMP-11). The negative control was DakoCytomation mouse serum diluted to the same mouse IgG concentration as the primary antibody. All the dilutions were made in Antibody Diluent (Dako) and incubated for 30 min at room temperature. In a prior report, we confirmed the presence of the evaluated proteins by Western blot analysis of breast tumor cytosol samples. A single band of the expected molecular mass was observed for each protein [28]. We also used other antibodies for several factors, such as MMP-11 (clone SC3-05, 1/100, Calbiochem MERCK KgaA Darmstadt Germany). On the other hand, we also used antibodies against cytokeratins (AE1-AE3, DAKO 1/1) and vimentin (DAKO 1/100) to distinguish fibroblasts from tumor cells and CD45 (leukocyte common antigen, LCA) for macrophages. Tissue sections were deparaffinized in xylene and then rehydrated in graded concentrations of ethyl alcohol (100%, 96%, 80%, 70%, then water). To enhance antigen retrieval only for some antibodies, TA sections were microwave-treated in a H2800 Microwave Processor (EBSciences, East Granby, CT, USA) in citrate buffer (Target Retrieval Solution; Dako) at 99 • C for 16 min. Endogenous peroxidase activity was blocked by incubating the slides in peroxidase-blocking solution (Dako) for 5 min. The EnVision Detection Kit (Dako) was used as the reactivity detection system. Sections were counterstained with hematoxylin, dehydrated with ethanol, and permanently cover-slipped. We evaluated the immunohistochemical staining in the whole-tissue array section for each main cell type: Cancer cells; mononuclear inflammatory cells, including lymphocytes, plasma cells, and monocytes/macrophages (MICs); and cancer-associated fibroblasts (CAFs). We distinguished stromal cells from cancer cells because these latter cells are larger in size. In addition, fibroblasts are spindle cells, whereas mononuclear inflammatory cells are round cells. On the other hand, while cancer cells are arranged to form either an acinar or a trabecular pattern, stromal cells are isolated. Moreover, we used 2 markers to distinguish fibroblasts from tumor cells, cytokeratins and vimentin, as described above.

Statistical Analysis
The Chi-Square test was used to calculate significant differences between categorical variables. p < 0.05 was considered significant. For metastasis-free survival analysis, we used Cox's univariate method. Cox's regression model was used to examine interactions of different prognostic factors in a multivariate analysis, and all factors with p < 0.05 in the univariate analysis were included in the multivariate analysis. The PASW statistics 18 program (SPSS Inc., Chicago, IL, USA) was used for all calculations.

Bud Quantification and Its Relationship with Clinical Outcome and Clinicopathological Features
As can be seen in Figure 2, there was a wide variability of bud number in all tumors, ranging from 0 to 28, with 5 as the median value. In general, we found that areas with high budding had a more infiltrative pattern and peritumoral stroma response. The potential relationships between tumor budding and relapse-free survival and overall survival were evaluated in all patients included in the present study who were metastasis-free at the time of initial diagnosis. We examined all possible values obtained by counting tumor budding as cut-off points for predicting relapse-free survival. This analysis led us to define a tumor budding count of four as the optimal cut-off (X 2 = 19.1; p < 0.0001) ( Figure 3A). Differences in both relapse-free survival and overall survival curves were significant taking this cut-off value of four (p < 0.001, for both) ( Figure 3B Table 2 shows the relationship between tumor budding grade and clinicopathological features. Budding grade correlated significantly with patient age (p = 0.011). Thus, tumors from older women had more tumor buds that tumors from younger patients. However, our results did not show significant associations between this histological finding and other clinicopathological features, including menopausal status, tumor size, nodal status, histological grade, and ER (estrogen receptors), PgR (progesterone receptors), and HER-2 (human epidermal growth factor receptor 2) status.

Expression of MMPs and TIMPs at the Invasive Front and Their Relationship with Tumor Budding Grade
Considering the importance of MMPs/TIMPs in tumor progression, especially in invasion and metastasis, we also investigated their possible relationship with tumor budding grade at the invasive front of breast carcinomas in all cases included in the present study. In addition, we investigated the expressions of these factors according to each cell type. In tumors positive for cells expressing either MMPs or TIMPs, at least 70% of these cells showed a positive immunostaining of each evaluated field. The majority of MMPs and TIMPs were mainly expressed in cancer cells at the invasive front (MMP-9: 97.3% of the tumors, MMP-11: 98.7%, MMP-14: 94.0%, TIMP-1: 86.7%, TIMP-2: 96.7%). Remarkably, as shown in Figure 4A, we found a similar MMP/TIMP expression both in cancer cells from tumor buds and those from the tumor mass at the invasive front, with differences for MMPs and TIMPs immunostaining of less than 5%. These proteins were also expressed by stromal cells but in a lower percentage of tumors. We found no significant differences for MMPs and TIMPs immunostaining between tissue array sections of each tumor or between the different areas of each section. Figure 4 shows representative examples of CAFs and MICs expressing MMPs and TIMPs localized at the invasive front in breast carcinomas. Immunostaining for these proteins revealed a cytoplasmic location in cancer cells and tumor-associated stromal cells, including both CAFs and MICs. In neoplasms positive for CAFs and MICs expressing either MMPs or TIMPs, at least 70% of these cells showed a positive immunostaining in each evaluated field. Immunostaining for these proteins had a cytoplasmic location in all positive cases. With regard to MMP-14 expression, the immunostaining showed both the cytoplasmic and membrane location. As can be seen in Table 3

Single Combination of Tumor Budding Grade and MMP/TIMP Expression Improves Prognostic Evaluation
To assess the contribution of the joint evaluation of tumor budding and MMP/TIMP expression to the prediction of clinical outcome, we explored all possible combinations.

Discussion
We selected invasive ductal breast carcinomas because they are the most prevalent histological subtype. In addition, they represent a heterogeneous group with varying tumor features and clinical outcomes in which it would be of great clinical value to find additional prognostic factors. Our results are in agreement with previous reports demonstrating that tumor budding is a morphological finding related to worse prognosis in several malignant tumors [7][8][9][10][11][12][13][14][29][30][31], including breast cancer [15][16][17][18]. In the present study, we found that tumor budding count was associated with MMP/TIMP expression, and that the combination between both factors identified a population of patients differing in their clinical outcome.
Since tumor budding at the invasive front has been postulated as the first step of invasion and metastasis [32], this histological feature could be a sensitive indicator of tumor aggressiveness. However, it has been reported that tumor cells from buds do not differ from those belonging to the tumor mass with regard to their both biological and therapeutic target factors in breast cancer, such as ER or HER-2 status [33], and even show a lower Ki67 index compared to those from the tumor mass [15]. In the present study, we found a similar MMP/TIMP expression in tumor cells from buds and those from the invasive front mass. Nevertheless, although analyzed as a histological finding represented by a tiny number of tumor cells compared with those of the tumor body [34], tumor budding is considered to reflect tumor migration in the context of the epithelialmesenchymal transition (EMT), or at least constitutes a hybrid epithelial/mesenchymal phenotype displaying collective cell migration [3] This migration consists of epithelial cells, which lose their cell-cell adhesion and gain migratory and invasive traits that are typical of mesenchymal cells. In accordance with this biological event, it has been shown that tumor cells from buds from breast cancer show phenotypic features, such as the lack of E-cadherin expression and expression of vimentin [15,35]. However, it is also known that tumor invasion depends not only on the biological aggressiveness of tumor cells but also on the tumor microenvironment. In fact, it has been previously reported that tumor budding correlates significantly and positively with the tumor stroma percentage [17]. It has been proposed that tumor stroma has an important role in facilitating tumor cell dedifferentiation and dissemination, perhaps providing suitable energy substrate and reducing the buildup of metabolic waste [36]. A functional variability of tumor stroma cells which influences the tumor progression can also be hypothesized. Two well-studied cellular components of the tumor stroma are CAFs and MICs, both of which may be present within the invasive front, contributing to tumor progression through several biological mechanisms such as the secretion of several molecules [24,37]. One of these stromal derived molecular factors which may influence tumor aggressiveness are MMPs. It is known that these enzymes play a key role in the degradation of the extracellular matrix and basal membrane, which are prerequisites for invasion and metastasis [20]. In addition, it has also been shown that MMPs display other biological activities related to tumor progression, such as the induction of proliferation [38,39], apoptosis-resistant cells [40], neoangiogenesis [41], and the facilitation of EMT by increasing cell mobility [42,43] or by repressing the expression of E-cadherin and other cell surface adhesion molecules [44,45]. On the other hand, TIMPs do not only have an inhibitory effect against MMP actions but also are involved in protumoral activities such as the induction of proliferation and the inhibition of apoptosis [46,47].
Taken together, in the present work, we integrated a morphological feature, tumor budding grading, in the context of the biological phenotype of the stromal cells from the invasive tumor front. Our own data show that the combination between budding count and MMP/TIMP expression by stromal cells strongly improved the prognostic evaluation. Thus, whereas the concurrence of low budding count and the absence of MMP/TIMP expression by stromal cells were significantly associated with an excellent prognosis, conversely, the concurrence of high budding count with the expression of MMPs/TIMPs by these stromal cells were associated with the occurrence of distant metastasis and shortened survival. Especially relevant was our finding that the combination of tumor budding grade with MMP-11 expression by stromal MICs was significantly and independently associated with the highest prognostic accuracy. This finding points to the importance of the functional type of the tumor inflammatory infiltrate, as well as MMP-11 expression in the tumor progression of breast carcinomas. MMP-11 (also named stromalysin 3) was found to increase cancer cell survival and implantation during the early steps of the adjacent connective tissue invasion [48]. Likewise, our group previously reported that tumors with MMP-11 expression by MICs showed a high cell immune ratio ((CD68 + ) macrophages/(CD3 + ) T cells + (CD20 + ) B cells) at the invasive front [49], an upregulation of inflammatory-related genes [50,51], and a poor prognosis in breast carcinomas [24,28]. Interestingly, in a recent report by our team, we showed that the coculture of CAFs, and especially those derived from breast carcinomas showing MMP-11 expression by stromal MICs, significantly increased the invasive capability of tumor cells with the cell line MDA-MB-231 [52]. Therefore, all these data led us to consider the importance of integrating the morphological evaluation of tumor budding grade at the invasive tumor front with biological aspects of stromal cells in this relevant scenario from breast carcinomas.
This study presents some limitations. First, the study had the difficulty of evaluating the heterogeneity of counted tumor budding along the tumor sample. We suggest that further studies be carried out to explore a more accurate and rapid evaluation of this morphologic finding, such as the use of artificial intelligence tools. Second, it would be interesting to integrate an evaluation of intratumor budding, which is quantified within the main tumor mass and has been associated with clinical outcome in colorectal cancers [53,54] and observed in breast cancer [16]. Third, further studies may be of interest in other breast cancer patient populations in order to evaluate the optimal cut-off point of tumor budding count for prognosis when we combine it with MMP/TIMP expression. In this sense, further studies with patients who underwent neoadjuvant treatment may be of interest due to the fact that is the standard of care in patients with breast cancers, especially in these ones with stage III tumors. In the present study, in a population of women with primary breast cancer randomized according to distant metastasis development and a long follow-up period, we found that a value of tumor budding count of 4 (out of a range from 0 to 28) was the optimal cut-off point to predict distant relapse-free survival. Interestingly, this cut-off for tumor budding at the invasive tumor front was also associated with metastatic properties in breast cancer by other authors [16]. However, a cut-off in a similar range of seven buds per field was also considered as an optimal cut-off value to predict clinical outcome in breast cancer by other authors [15].
In conclusion, our data show that a high tumor budding grade is a histological feature in breast carcinomas associated with a more aggressive tumor phenotype and which may increase metastatic potential. In addition, our results enhance the clinical relevance of tumor budding at the invasive front when this morphological finding is evaluated in the context of the biological phenotype of stromal cells which contribute to tumor progression. Therefore, our data may contribute to a better prognosis interpretation of tumor budding, avoiding several issues related with reproducibility among pathologies or quantification criteria.  Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author

Conflicts of Interest:
The authors declare that they have no conflicts of interest.