Expression Patterns of Microenvironmental Factors and Tenascin-C at the Invasive Front of Stage II and III Colorectal Cancer: Novel Tumor Prognostic Markers

Background: Biological markers expressed in cancer cells and the surrounding cancer-associated fibroblasts (CAF) can be used for prediction of patient prognosis in colorectal cancer (CRC). Here, we used immunohistochemical techniques to evaluate cancer cells’ expression of specific biomarkers that are closely associated with neoplastic progression.

Methods: Immunohistochemical markers included Ki-67, p53, β-catenin, MMP7, E-cadherin and HIF1-α. We also characterized microenvironmental markers expressed by CAF, including expression of α-smooth muscle actin, CD10, podoplanin, fibroblast specific protein 1, platelet derived growth factor β, fibroblast association protein, tenascin-C (TNC), ZEB1 and TWIST1. The study population consisted of 286 CRC patients with stage II and III disease. Stage II and III CRC were divided into a first and a second cohort (for validation). The CRCs were stratified using cluster analysis. To identify the utility of prognostic markers in stage II and III CRC, univariate and multivariate analyses were performed in both cohorts.

Results: Stage II and III CRCs were stratified into 3 subgroups. Specific subgroups were significantly correlated to disease-free survival using univariate and multivariate analyses in the first cohort. High expression of TNC was identified as a single prognostic marker in both cohorts by univariate and multivariate analyses.

Conclusions: We suggest that the presence of a specific subgroup defined by multiple markers can be used for prediction of CRC outcome in stages II and III. In addition, we showed that high expression of TNC was correlated with a poorer prognosis in stages II and III of CRC.

Introduction

Colorectal cancer (CRC) is the third most commonly diagnosed cancer and the third leading cause of cancer death in both men and women in the United States (1). These trends for incidence and mortality are common worldwide (1). Remarkable progress has been made in the diagnosis and treatment of CRC. In spite of such advances, CRC is often discovered at an advanced stage at which point achieving a cure is very difficult (2). Therefore, the development of effective markers to predict patient prognosis of CRC is greatly needed.

The outcome of patients with CRC can be predicted by prognostic factors, such as the TNM staging system proposed by the UICC and AJCC (3, 4). Additionally, novel and promising prognostic biomarkers are listed in the WHO classification 2019 (5). There are 2 histological processes that are present within the tumor microenvironment at the invasive front of CRC: tumor budding and the desmoplastic reaction (DR) (6–9). Tumor budding is defined as single cells or clusters of up to four tumor cells at the invasion front of CRC (6–8). It is closely associated with both local and distant metastases and is therefore a histological biomarker of tumor progression and a poor prognosis (6–8). The classification of the DR was recently proposed by Ueno et al. as a prognostic histological marker (9). A pronounced desmoplastic stromal reaction in the microenvironment involves complex cellular interactions at the invasive front (10). This theory posits that cooperation between cancer cells and cancer associated fibroblasts (CAFs) present within the tumor microenvironment is necessary to support tumor growth and progression (10, 11). In addition, the microenvironment itself plays an important role in neoplastic progression and metastasis in CRC (10, 11). Whereas such histological findings are widely used as markers for establishing a patient’s prognosis, they do not explain the underlying cellular processes that promote tumor growth and metastasis (12, 13). Therefore, the discovery of additional markers would be very beneficial. We propose that identification of protein expression patterns in cancer cells and CAFs could provide new biological insights and guide the development of new therapies for CRC (12, 13).

In this study, we analyzed immunohistochemical data to identify possible protein expression patterns in stages II and III of CRC that predict patient outcome. We focused on markers that are closely associated with tumor growth and progression within the microenvironment.

Materials And Methods

Patients

CRC patients who underwent curative surgery at stages II or III at Iwate Medical University Hospital from January 2009 to December 2015 were included in the present study. In total, 286 patients were included the first cohort (148 cases) and in a second cohort for validation (138 cases), which were evaluated through a retrospective analysis. We used a block randomization method in the research design to select and divide participants into different groups or conditions in order to avoid bias in the selection of two cohorts. Paraffin embedded tissues were well preserved, medical records were complete and patient status had been followed up, including overall survival and disease-free survival data that were confirmed through telephone interviews and by the mail. In addition, cases with invasion beyond the proper muscular layer were included for determination of the desmoplastic reaction (9). Finally, patients who underwent preoperative chemoradiotherapy and emergency surgery were excluded. In addition, patients who had evidence of hereditary non-polyposis colorectal cancer or familial adenomatous polyposis were not enrolled. The clinicopathological variables characterizing the patients included tumor location, stage and t stage, histological type, lymphatic/venous invasion and tumor budding. The variables were recorded according to the General Rules for Management of the Japanese Colorectal Cancer Association (Table 1) (14). In addition, DR classification was determined based on Ueno’s classification (9).

TABLE 1

Table 1 Clinicopathological findings in stage II and III colorectal cancer.

This study was approved by the local ethics committee of Iwate Medical University (approval number MH2020-070), and all patients provided informed consent.

Determination of Disease-Free Survival

We determined the duration of disease-free survival at which metastasis was discovered during the follow-up period (2 times/year to 3 times/year) using computed tomography.

Chemotherapeutic Treatment After Surgery for Stage II or III CRC

Following surgery, Capecitabine or UFT/UZEL (Tegafur Uracil + Calcium Folinate) were administered in stage II CRC (20/140 cases), whereas FOLFOX, including the drugs leucovorin calcium (folinic acid), fluorouracil and oxaliplatin were used in stage III CRC (85/146 cases). The other 181 patients, including 120 cases in stage II and 61 cases in stage III did not receive additional chemotherapy following surgery.

Determination of Sample Size

The sample size required to identify differences in overall and disease-free survival between cohorts was determined using JMP Pro 13.0 software (SAS, Tokyo, Japan). From the calculation, at least 120 cases were required. The statistical power (detection power) was set to 0.8, which is commonly used in medical studies.

Tissue Microarray Construction (TMA)

The TMAs were assembled using a manual tissue array (Azumaya Co, Tokyo, Japan). Five mm tissue cores were taken from each targeted lesion and placed into a recipient block containing 12 cores including 10 cancer tissues and 2 cores for control tissues (normal colon; CRC). After construction, 3-micron sections were cut and stained with hematoxylin and eosin on the initial slides to verify the histologic diagnosis. Serial sections were cut from the TMA block for immunohistochemical staining.

Immunohistochemistry

Tumors were routinely fixed in 20% neutral-buffered formalin and embedded in paraffin wax. Three-micron-thick paraffin sections were cut, dewaxed, and rehydrated. Microarray slides were incubated in 3% hydrogen peroxide to block endogenous peroxidase. Antigen retrieval was performed using an autoclave-based method, followed by incubation with the primary antibody overnight at 4°C in a high humidity cabinet. Slides were processed using the Dako Autostainer Universal Staining System (Dako, Glostrup, Denmark) (12). The specimens were treated with citrate buffer (pH 6.0) using a microwave [three times for 5 min, 750 W; cat. no. H2500; Microwave Processor (Bio-Rad Laboratories, CA, USA)] and then reacted with antibodies, as previously described. Antibodies used in this study were classified into 2 subgroups: epithelial (cancer cells) and interstitial (cancer associated fibroblasts, CAF) markers. Antibodies targeting CAFs included the following: α-smooth muscle actin (α-SMA, Dako 1A4), CD10 (Dako, 56C6), podoplanin (Dako, D2-40), fibroblast specific protein 1 (FSP1; S100A4, Dako, polyclonal), platelet derived growth factor receptor (PDGFR-β; 28E1, Cell Signaling Technology), fibroblast association protein (FAP, Abcam, EPR20021) and tenascin-C (IBL, 4F10TT). For EMT, we utilized zinc finger E-box binding homeobox 1 (ZEB1, Sigma-Aldrich, polyclonal) and Twist-related protein 1 (TWIST1, Abcam, Twist2C1a). CAFs were recognized as “spindle-shaped cells” by experienced pathologists (T.S. and N.U.). Cytoplasmic staining of tumor cells was conducted with antibodies against α-SMA, CD10, podoplanin, FSP1, PDGFR-β, FAP and tenascin-C. Nuclear staining of fibroblasts was based on positivity for ZEB1 and TWIST1 expression. Furthermore, antibodies targeting cancer cells in this study included Ki-67 (Dako, MIB1) for proliferative activity, p53 (Dako, Do7) for p53 mutation, β-catenin (Dako, β-catenin-1) for activation of Wnt signaling, a central signal transducer in CRC, MMP7 (Daiichi Fine Chemical, 141-7B2) for cancer progression, E-cadherin (Dako, NCH-38) for cellular adhesion and HIF1-α (Novus Biologicals, polyclonal) for cancer-specific metabolic marker which may be associated with tumor progression. Detailed information of antibodies is summarized in Supplementary Table 1.

Assessment of Scoring of Immunohistochemical Expression

The expression of the markers was scored for both the intensity and extent of immunopositivity, as described in a previous report with slight modification (15). The immunostaining intensity of the cancer cells and CAFs in the CRCs was classified into 4 categories as follows: negative, weak, moderate and strong. The immunostaining extent was semi-quantified as follows: 0%, 1-25%, 26-50%, 51-100%. The combination of intensity and extent was scored. Scores 2–3 were defined as a positive staining pattern, as shown in Supplementary Table 2. In addition, the score was also sub-classified into low (score 0-1) and high expression (score 2-3). Assessment of scoring was performed by two pathologists. If agreement was not obtained between the pathologists, we asked an additional pathologist regarding the assessment. Finally, the score was determined by agreement of more than two pathologists.

In the present study, a wide range of expression levels was observed for all the markers. Thus, we selected the deepest invasive region as a target area to measure the expression levels of markers.

Hierarchical Analysis of the Expression of CAF and EMT Markers

Hierarchical cluster analysis was performed for clustering of the samples according to the expression level in order to achieve maximal homogeneity for each group and the greatest differences between the groups using open-access clustering software (Cluster 3.0 software; bonsai.hgc.jp/~mdehoon/software/cluster/software.htm). The clustering algorithm was set to centroid linkage clustering, which is the standard hierarchical clustering method used in biological studies.

Statistical Analysis

Data were analyzed using JMP Pro 13.0 software (SAS, Tokyo, Japan). Data obtained for clinicopathological features (sex, location, pT, stage, histological type, lymphatic invasion, venous invasion, tumor budding, desmoplastic reaction, overall survival, disease-free survival) and subgroup (subgroups 1, 2 and 3) were analyzed using Fisher’s exact test. In addition, the comparison of the age distributions within each subgroup was performed using the Kruskal-Wallis test. If multigroup comparisons were needed for statistical analysis, we used Bonferroni corrections.

Kaplan-Meier analyses were performed using a log-rank test for survival analyses. Univariate and multivariate analyses were conducted with Cox proportional hazards model to identify statistical differences for prediction of overall and disease-free survival. The level of significance was p < 0.05, and the confidence interval (CI) was determined at the 95% level.

Results

A representative figure is shown in Figure 1. In addition, the cancer invasive front is depicted in Supplementary Figure 1.

FIGURE 1

Figure 1 Representative features of immunohistochemical staining of biological markers we examined based on expression level (low and high). (A) α-SMA. (B) CD10. (C) Podoplanin. (D) FSP1. (E) FAP. (F) Tenascin-C. (G) PDGFR-β. (H) ZEB1. (I) TWIST1. (J) Ki-67. (K) p53. (L) MMP7. (M) β-catenin. (N) E-cadherin. (O) HIF1-α.

Analyses of Clinicopathological Variables and Biological Markers in the First Cohort

Hierarchical Clustering Based on Marker Scores in First Cohort

We performed hierarchical clustering based on marker scores to evaluate differences in expression patterns of cancer cell-, CAF- and EMT-related markers in stage II and III CRC. Three distinct subgroups were stratified, as shown in Figure 2. The vertical line shows the expression of each marker in cancer cells and fibroblasts and the horizontal lines denote “relatedness” between samples. There was no statistical difference in the frequency of clinicopathological variables among subgroups 1, 2 and 3. Although immature desmoplastic reaction present in subgroup 1 showed a high frequency among the 3 subgroups, such association between the 3 subgroups did not quite reach a statistically significant level (p = 0.0508). However, the frequency of disease-free survival was significantly higher in subgroup 1 than in subgroup 2 (p<0.0001). Detailed data are shown in Table 2.

TABLE 2

Table 2 Clinicopathological variables according to each subgroup in the first cohort.

FIGURE 2

Figure 2 Hierarchical cluster analysis of colorectal cancer patients with stage II or III disease based on the expression patterns of cancer cells and cancer-associated fibroblast (CAF) proteins in the first cohort. The examined CRCs were subclassified into 3 subgroups.

Survival Analyses of Each Subgroup in the First Cohort

Kaplan–Meier analyses were performed to determine the association between the disease-free survival frequencies and the subgroups. Subgroup 1 had a poorer disease-free survival, compared to subgroup 2 (p < 0.0001). However, overall survival did not differ among the subgroups (Supplementary Figure 2).

The Association of Clinicopathological Variables and Subgroups With Survival of Stage II and III CRC Patients: Univariate and Multivariate Analyses of the First Cohort Using a Cox Proportional Hazards Model

The univariate analysis of stage II and III CRC patients (Table 3a) identified 5 factors: histologic type (mucinous carcinoma vs. well differentiated adenocarcinoma), stage (II vs III), desmoplastic reaction (mature vs immature) and subgroup (1 vs 2; 1 vs 3). Table 3b reveals that 3 factors (mucinous carcinoma vs well differentiated adenocarcinoma, mature DR versus immature DR, subgroup 1 versus 2) were retained in the multivariate analysis using a Cox proportional hazards model

TABLE 3

Table 3 Association of clinicopathological variables and subgroups with disease-free survival and overall survival in the first cohort in univariate and multivariate analyses.

Using a similar method, we performed univariate analysis for screening of overall survival of stage II and III CRC patients. As a result, 3 factors, including stage (II vs III), desmoplastic reaction (mature vs immature), and subgroup (1 vs 2) were identified in univariate analysis (Table 3c). However, no factors were retained in multivariate analysis (Table 3d).

Association of Individual Markers With Individual Subgroups in the First Cohort

The frequency of positive scores (score 2 or 3) of SMA was higher in subgroup 2 than in subgroup 1. There were statistically significant differences in the frequencies of positive scores among subgroups 1, 2 and 3 (subgroup 1, 2 > 3). In addition, significant differences in the frequencies of positive scores for tenascin-C between subgroups 1 and 2, and 3 were found (subgroup 1 > 2, 3). The frequency of the positive score for ZEB1 was statistically higher in subgroup 2 than in subgroup 3. Next, there was a statistically significant difference in the frequencies of positive scores for TWIST1 between subgroup 3 and subgroup 1 (subgroup 1 > 3). The positive score for p53 was significantly greater in subgroup 2 than in subgroups 1 and 3. Furthermore, there was a significant difference in the frequencies of positive scores for p53 between subgroups 1 and 3. Finally, we observed statistically significant differences in the frequencies of positive MMP7 scores among subgroups 1 and 2, and 3 (subgroup 1, 2 > 3). Detailed data are shown in Figure 3.

FIGURE 3

Figure 3 Marker expression levels. (A) Expression level of each marker in the first cohort. (A) α-SMA, (B) CD10, (C) Podoplanin, (D) FSP1, (E) FAP, (F) Tenascin- C, (G) PDGFR-β, (H) ZEB1, (I) TWIST, (J) Ki-67, (K) p53, (L) MMP7, (M) β-catenin, (N) E-cadherin and (O) HIF1-α.

The Association of Clinicopathological Variables and Individual Markers With the Survival of Stage II and III CRC Patients: Univariate and Multivariate Analyses of the First Cohort

With regard to disease-free survival, 3 variables (stage II vs III; mature vs immature; mucinous carcinoma vs well differentiated adenocarcinoma) and one marker (tenascin-C) were identified in univariate analysis (Table 4a). Among those 4 parameters, 2 variables, including desmoplastic reaction and histological type and one marker, tenascin-C, were retained in multivariate analysis (Table 4b). In overall survival, stages (II vs III) and desmoplastic reaction (mature vs immature) were identified in univariate analysis (Table 4c). Desmoplastic reaction (mature vs immature) was retained in multivariate analysis (Table 4d).

TABLE 4

Table 4 Association of clinicopathological variables and individual marker with disease-free survival and overall survival in the first cohort in univariate and multivariate analyses.

Analyses of Clinicopathological Variables and Individual Markers in the Second Cohort (Validation)

The Association of Clinicopathological Variables and Individual Markers With the Survival of Stage II and III CRC Patients: Univariate and Multivariate Analyses of the Second Cohort

With regard to disease-free survival, 5 variables (pT3 vs. pT4; stage II vs. III; positive venous invasion vs. negative venous invasion; low grade budding vs. high grade budding; mature vs. immature) and 2 markers (tenascin-C and β-catenin) were identified in univariate analysis (Table 5a). However, only 1 factor (tenascin-C) was retained in multivariate analysis (Table 5b). In overall survival, 4 variables (pT3 vs. pT4; stage II vs. III and desmoplastic reactions (mature vs. immature; and intermediate vs. mature) and 2 individual markers (tenascin-C and Ki-67) were detected in univariate analysis (Table 5c). Only the positive expression of tenascin-C was retained in multivariate analysis (Table 5d).

TABLE 5

Table 5 Association of clinicopathological variables and individual markers with disease-free survival and overall survival in the second cohort in univariate and multivariate analyses.

Discussion

Certain proteins expressed by microenvironmental cells play crucial roles in neoplastic progression of CRC. Those proteins may be derived from cancer cells or from stromal cells (sometimes termed “cancer-associated fibroblasts” (CAFs) (12, 13). Proteins expressed by cancer cells and CAFs interact with one another and this interaction is likely important at the invasive front (12, 13). According to that theory, the combination of proteins from cancer cells and CAFs mediate tumor growth and progression (12, 13). In the present study, specific expression patterns could be correlated with the prognosis of stage II and III CRC patients. Therefore, the current results suggest that a specific subgroup (identified here by stratification) can be used to evaluate the role and significance of various proteins produced by microenvironmental cells. Finally, in the present study, subgroup 1 was correlated with disease-free survival. However, the presence in subgroup 1 did not correlate with overall survival. The reason remains unknown.

In the current study, we used 15 microenvironment-related markers (cancer cell markers and CAF markers) to identify associations of expression patterns with patient outcomes. Among the cancer cell-related markers, a high Ki-67-positive rate and overexpression of p53 were considered to reflect the characteristics of tumors. Intranuclear expression of β-catenin and high expression of MMP7, E-cadherin, and HIF1-α are closely associated with tumor budding, which is a key histological feature occurring in the cancer microenvironment (16–18). By contrast, stromal markers, including α-SMA, CD10, podoplanin, FSP1, PDGFR β, FAP, and TNC, were used as CAF markers. These markers are thought to be associated with enhanced progression of CAFs. Based on these findings, we suggest that the microenvironment-related markers used in the current study may be suitable for identification of the molecular mechanisms of neoplastic progression and cancer metastasis in the tumor microenvironment.

Tenascin-C (TNC) is an extracellular matrix molecule that drives the progression of many types of human cancer. The basis for its actions remains unclear (19). TNC is associated with organogenesis accompanying cell proliferation and migration, resulting in the epithelial-mesenchymal transition (EMT) that might result from interactions between cancer cells and stromal cells (20). EMT is the process by which polarized epithelial cells are converted into mesenchymal cells during cancer progression. As a result, carcinoma cells lose their epithelial polarity and intercellular connections, allowing them to escape the surrounding epithelium (20, 21). The expression of TNC facilitates such phenotypic changes, alterations that are enhanced by TGF-β, a promoter of EMT (19–21). Murakami et al. revealed that TNC in primary CRC stroma might be a novel biomarker that is predictive of postoperative prognosis (21). Finally, TNC may promote EMT-like change and proliferation, alterations that lead to poor prognosis in CRC patients (20).

TNC may be involved in cancer growth and metastatic processes via the Hedgehog (HH) signaling pathway, caused either by mutations in the pathway (ligand independent) or through HH overexpression (ligand dependent) (22). HH signaling starts with secretion of the HH ligand, followed by secretion of Patched (PTC), the transmembrane protein Smoothened (SMO) and three GLI (Glioma-associated oncogene) zinc finger transcription factors (23). The HH/GLI1 pathway promotes cancer growth, stem cell self-renewal and metastatic behavior in advanced CRC (24). Human CRC stem cells require active HH/GLI1 signaling for survival and self-renewal (25). Our finding suggests that activation of CAF at the invasive front is caused by high expression of TNC facilitated via HH signaling (26). In addition, accumulating evidence suggests that activated HH signaling plays an important role in neoplastic transformation as well as the development of drug resistance of human cancers (27). Thus, HH signaling during tumorigenesis and the development of chemo-resistance are closely associated. Those findings suggest that therapeutic strategies might target such signals in human cancers and their relapse (26, 27). For example, cyclopamine is an HH signal pathway antagonist and consequently is expected to improve the survival of patients with CRC by inhibiting the proliferation of colon cancer cells (28). Previous study showed that cyclopamine treatment results in decreased levels of mRNA coding for HH, SMO and PTCH, all of which were highly expressed in colon cancer cell lines (28). These findings may influence potential therapeutic strategies because TNC expression by CAF may be targeted in future molecular therapies.

High expression of TNC was reported to be a prognostic marker for CRC through induction of EMT and cell proliferative activity (20). According to that study, TNC may facilitate EMT-like changes and could be associated with a poor prognosis of CRC patients. This finding is consistent with other data showing that cancer cell-derived TNC promotes cancer cell invasion via EMT regulation. Thus, it is a novel indicator of poor prognosis (29). In the present study, we found that even in stages II and III, intermediate stages that account for the majority of surgically resected CRC, TNC was an independent prognostic marker. This result was validated by analysis of a second cohort. The present results showed that TNC in primary CRC stroma has the potential to be a novel biomarker that predicts postoperative prognosis.

There are some limitations to this study. First, the immunohistochemical markers we used in the present study may not yield consistent results. For clinical application, immunohistochemical reagents must be reliable and reproducible. In that regard, many immunohistochemical markers that are closely associated with the formation of the microenvironment have been analyzed (12, 13). In the current study, 15 microenvironment-related markers, including Ki-67, p53, β-catenin, MMP7, E-cadherin, and HIF1-α (for cancer cells) and CD10, podoplanin, FSP 1, PDGFR β, FAP, TNC, ZEB1, and TWIST1 (for CAFs) were used. Briefly, Ki-67 positivity and p53 overexpression have been widely used as characteristics of tumors. The remaining factors, including β-catenin, MMP7, E-cadherin, and HIF1-α, are closely associated with the formation of the cancer microenvironment. In addition, stromal factors could be classified as CAF or EMT markers. The two stromal markers used in this study were considered CAF markers given that all markers we used were expressed in CAFs. These CAF markers are suitable for identifying the functions of CAFs. Therefore, we concluded that the immunohistochemical markers examined in this study were all involved in generation of the tumor microenvironment at the invasive front. Finally, analysis of these immunohistochemical markers should yield reliable and reproducible results, as demonstrated in the current study. Second, the heterogeneous expression of the markers examined in this study may be problematic when determining marker expression levels (30). Although it may be difficult to avoid this problem, we suggest that the invasive front of cancer cells, which is critical for tumor progression, may be the best region for measuring the immunohistochemical expression levels of the chosen markers (10, 11). Finally, although there are many different reports regarding prognostic factors in CRC (31, 32), the different results may reflect the choice of markers, patient stage, heterogeneity of expression, staining platform, judging methods and cut-off value. In the present study, we suggest that the current results are reliable and reproducible under the conditions we employed.

Conclusions

Cancer cells and CAFs express many proteins that modulate neoplastic progression and metastasis. In the present study, we found that specific expression patterns may allow the prediction of patient outcome in CRC. In addition, the expression of TNC by CAFs might be a potential prognostic biomarker in stage II and III CRC patients. These results highlight a potential role for TNC in CRC tumor progression and provide novel mechanistic insights into the roles of HH, as it is associated with high expression of TNC in driving CRC progression. Our findings also suggest that TNC could be a critical target gene for the treatment of CRC. However, further study will be needed in the near future to confirm these results.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the ethics committee of Iwate Medical University Hospital (approval number MH2020-070). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

MH, who is the first author, constructed the figures and tables and performed the statistical analyses. MO assisted statistical analyses. YK and HU supported pathological interpretation of desmoplastic reactions. NY and NU helped in the interpretation of pathological findings. KO and AS provided clinical support during the preparation of the manuscript. TS, who is the corresponding author, contributed to the preparation of the manuscript, including all aspects of the data collection and analysis. All authors contributed to the article and approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We gratefully acknowledge the technical assistance of members of the Department of Molecular Diagnostic Pathology, Iwate Medical University for their support.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.690816/full#supplementary-material

Supplementary Figure 1 | Histological features of the invasive front. (A) Low-power view of the invasive area of CRC. (B) High-power view of the invasive area of CRC (mature type desmoplastic reaction). (C) Low-power view of the invasive area of CRC. (D) High-power view of the invasive area of CRC (intermediate type desmoplastic reaction). (E) Low-power view of the invasive area of CRC. (F) High-power view of the invasive area of CRC (immature type desmoplastic reaction).

Supplementary Figure 2 | Kaplan-Meier analyses of the disease-free survival (A) and overall survival (B) based on each subgroup of the first cohort.

Abbreviations

CAF, Cancer associated fibroblast; CRC, colorectal cancer; TMA, tissue microarray; MMP7, matrix metalloproteinase-7; FSP1, fibroblast specific protein 1; PDGFR-b, platelet derived growth factor receptor beta; FAP, fibroblast associated protein; ZEB1, zinc finger E-box binding homeobox 1; TWIST1, twistrelated protein 1.

References

1. Siegel RL, Miller KD, Goding Sauer A, Sauer A, Fedewa SA, Butterly LF, et al. Colorectal Cancer Statistics, 2020. CA Cancer J Clin (2020) 70:145–64. doi: 10.3322/caac.21601

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Chen FW, Sundaram V, Chew TA, Ladabaum U. Advanced-Stage Colorectal Cancer in Persons Younger Than 50 Years Not Associated With Longer Duration of Symptoms or Time to Diagnosis. Clin Gastroenterol Hepatol (2017) 15:728–37.e3. doi: 10.1016/j.cgh.2016.10.038

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Moccia F, Tolone S, Allaria A, Napolitano V, Rosa D, Ilaria F, et al. Lymph Node Ratio Versus TNM System As Prognostic Factor in Colorectal Cancer Staging. A Single Center Experience. Open Med (Wars) (2019) 14:523–31. doi: 10.1515/med-2019-0058

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Puppa G, Sonzogni A, Colombari R, Pelosi G. TNM Staging System of Colorectal Carcinoma: A Critical Appraisal of Challenging Issues. Arch Pathol Lab Med (2010) 134:837–52. doi: 10.1043/1543-2165-134.6.837

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Nagtegaal ID, Odze RD, Kimstra D, Paradis V, Rugge M, Schirmacher P, et al. The 2019 WHO Classification of Tumours of the Digestive System. Histopathology (2020) 76:182–8. doi: 10.1111/his.13975

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Lugli A, Kirsch R, Ajioka Y, Bosman F, Cathomas G, Dawson H, et al. Recommendations for Reporting Tumor Budding in Colorectal Cancer Based on the International Tumor Budding Consensus Conference (ITBCC) 2016. Mod Pathol (2017) 30:1299–311. doi: 10.1038/modpathol.2017.46

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Watanabe T, Muro K, Ajioka Y, Hashiguchi Y, Ito Y, Saito Y, et al. Japanese Society for Cancer of the Colon and Rectum. Japanese Society for Cancer of the Colon and Rectum (JSCCR) Guidelines 2016 for the Treatment of Colorectal Cancer. Int J Clin Oncol (2018) 23:1–34. doi: 10.1007/s10147-017-1101-6

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Mitrovic B, Schaeffer DF, Riddell RH, Kirsch R. Tumor Budding in Colorectal Carcinoma: Time to Take Notice. Mod Pathol (2012) 25:1315–25. doi: 10.1038/modpathol.2012.94

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Ueno H, Ishiguro M, Nakatani E, Ishikawa T, Uetake H, Murotani K, et al. Prognostic Value of Desmoplastic Reaction Characterization in Stage II Colon Cancer: Prospective Validation in a Phase 3 Study (SACURA Trial). Br J Cancer (2021) 124:1088–97. doi: 10.1038/s41416-020-01222-8

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zlobec I, Lugli A. Epithelial Mesenchymal Transition and Tumor Budding in Aggressive Colorectal Cancer: Tumor Budding. Oncotarget (2010) 1:651–61. doi: 10.18632/oncotarget.199

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Tommelein J, Verset L, Boterberg T, Demetter P, Bracke M, De Wever O. Cancer-Associated Fibroblasts Connect Metastasis-Promoting Communication in Colorectal Cancer. Front Oncol (2015) 5:63. doi: 10.3389/fonc.2015.00063

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sugai T, Yamada N, Eizuka M, Sugimoto R, Uesugi N, Osakabe M, et al. Vascular Invasion and Stromal S100A4 Expression at the Invasive Front of Colorectal Cancer Are Novel Determinants and Tumor Prognostic Markers. J Cancer (2017) 8:1552–61. doi: 10.7150/jca.18685

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sugai T, Uesugi N, Kitada Y, Yamada N, Osakabe M, Eizuka M, et al. Analysis of the Expression of Cancer-Associated Fibroblast- and EMT-Related Proteins in Submucosal Invasive Colorectal Cancer. J Cancer (2018) 9:2702–12. doi: 10.7150/jca.25646

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Hashiguchi Y, Muro K, Saito Y, Yamada N, Osakabe M, Eizuka M, et al. Japanese Society for Cancer of the Colon and Rectum. Japanese Society for Cancer of the Colon and Rectum (JSCCR) Guidelines 2019 for the Treatment of Colorectal Cancer. Int J Clin Oncol (2020) 25:1–42. doi: 10.1007/s10147-019-01485-z

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Sasaki K, Sugai T, Ishida K, Osakabe M, Amano H, Kimura H, et al. Analysis of Cancer-Associated Fibroblasts and the Epithelial-Mesenchymal Transition in Cutaneous Basal Cell Carcinoma, Squamous Cell Carcinoma, and Malignant Melanoma. Hum Pathol (2018) 79:1–8. doi: 10.1016/j.humpath.2018.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhang W, Shi X, Peng Y, Wu M, Zhang P, Xie R, et al. HIF-1α Promotes Epithelial-Mesenchymal Transition and Metastasis Through Direct Regulation of ZEB1 in Colorectal Cancer. PloS One (2015) 10:e0129603. doi: 10.1371/journal.pone.0129603

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Mansour RN, Enderami SE, Ardeshirylajimi A, Fooladsaz K, Fathi M, Ganji SM. Evaluation of Hypoxia Inducible Factor-1 Alpha Gene Expression in Colorectal Cancer Stages of Iranian Patients. J Cancer Res Ther (2016) 12:1313–7. doi: 10.4103/0973-1482.199542

PubMed Abstract | CrossRef Full Text | Google Scholar

18. van Wyk HC, Roseweir A, Alexander P, Park JH, Horgan PG, McMillan DC, et al. The Relationship Between Tumor Budding, Tumor Microenvironment, and Survival in Patients With Primary Operable Colorectal Cancer. Ann Surg Oncol (2019) 26:4397–404. doi: 10.1245/s10434-019-07931-6

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Sun Z, Schwenzer A, Rupp T, Murdamoothoo D, Vegliante R, Lefebvre O, et al. Tenascin-C Promotes Tumor Cell Migration and Metastasis Through Integrin Alpha9beta1-Mediated YAP Inhibition. Cancer Res (2018) 78:950–61. doi: 10.1158/0008-5472.CAN-17-1597

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Yang Z, Zhang C, Qi W, Cui C, Cui Y, Xuan Y. Tenascin-C as a Prognostic Determinant of Colorectal Cancer Through Induction of Epithelial-to-Mesenchymal Transition and Proliferation. Exp Mol Pathol (2018) 105:216–22. doi: 10.1016/j.yexmp.2018.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Murakami T, Kikuchi H, Ishimatsu H, Iino I, Hirotsu A, Matsumoto T, et al. Tenascin C in Colorectal Cancer Stroma Is a Predictive Marker for Liver Metastasis and Is a Potent Target of miR-198 as Identified by microRNA Analysis. Br J Cancer (2017) 117:1360–70. doi: 10.1038/bjc.2017.291

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Yang Z, Zhang C, Feng Y, Quan M, Cui Y, Xuan Y. Tenascin-C Predicts Poor Outcomes for Patients With Colorectal Cancer and Drives Cancer Stemness via Hedgehog Signaling Pathway. Cancer Cell Int (2020) 20:122. doi: 10.1186/s12935-020-01188-w

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ruiz i Altaba A. Hedgehog Signaling and the Gli Code in Stem Cells, Cancer, and Metastases. Sci Signal (2011) 4:Pt9. doi: 10.1126/scisignal.2002540

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Gulino A, Ferretti E, De Smaele E. Hedgehog Signalling in Colon Cancer and Stem Cells. EMBO Mol Med (2009) 1:300–2. doi: 10.1002/emmm.200900042

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Stecca B, Mas C, Clement V, Zbinden M, Correa R, Piguet V, et al. Melanomas Require HEDGEHOG-GLI Signaling Regulated by Interactions Between GLI1 and the RAS-MEK/AKT Pathways. Proc Natl Acad Sci USA (2007) 104:5895–900. doi: 10.1073/pnas.0700776104

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Orend G, Chiquet-Ehrismann R. Tenascin-C Induced Signaling in Cancer. Cancer Lett (2006) 244:143–63. doi: 10.1016/j.canlet.2006.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Sari IN, Phi LTH, Jun N, Wijaya YT, Lee S, Kwon HY. Hedgehog Signaling in Cancer: A Prospective Therapeutic Target for Eradicating Cancer Stem Cells. Cells (2018) 7:208. doi: 10.3390/cells7110208

CrossRef Full Text | Google Scholar

28. Wu JY, Xu XF, Xu L, Niu PQ, Wang F, Hu GY, et al. Cyclopamine Blocked the Growth of Colorectal Cancer SW116 Cells by Modulating Some Target Genes of Gli1 In Vitro. Hepatogastroenterology (2011) 58:1511–8. doi: 10.5754/hge10765

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Takahashi Y, Sawada G, Kurashige J, Matsumura T, Uchi R, Ueo H, et al. Tumor-Derived Tenascin-C Promotes the Epithelial-Mesenchymal Transition in Colorectal Cancer Cells. Anticancer Res (2013) 33:192.

Google Scholar

30. Sagaert X, Vanstapel A, Verbeek S. Tumor Heterogeneity in Colorectal Cancer: What Do We Know So Far? Pathobiology (2018) 85:72–84. doi: 10.1159/000486721

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Jang BG, Kim HS, Chang WY, Bae JM, Oh HJ, Wen X, et al. Prognostic Significance of Stromal GREM1 Expression in Colorectal Cancer. Hum Pathol (2017) 62:56–65. doi: 10.1016/j.humpath.2016.12.018

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Oh HJ, Bae JM, Wen XY, Cho NY, Kim JH, Kang GH. Overexpression of POSTN in Tumor Stroma Is a Poor Prognostic Indicator of Colorectal Cancer. J Pathol Transl Med (2017) 51:306–13. doi: 10.4132/jptm.2017.01.19

PubMed Abstract | CrossRef Full Text | Google Scholar