ORIGINAL RESEARCH article

Front. Oncol., 19 August 2021

Sec. Surgical Oncology

Volume 11 - 2021 | https://doi.org/10.3389/fonc.2021.690816

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

  • 1. Department of Molecular Diagnostic Pathology, School of Medicine, Iwate Medical University, Shiwagun’yahabachou, Japan

  • 2. Department of Surgery, School of Medicine, Iwate Medical University, Shiwagun’yahabachou, Japan

  • 3. Department of Surgery, National Defense Medical College, Tokorozawa, Japan

Article metrics

View details

19

Citations

3,1k

Views

1,8k

Downloads

Abstract

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) (69). Tumor budding is defined as single cells or clusters of up to four tumor cells at the invasion front of CRC (68). It is closely associated with both local and distant metastases and is therefore a histological biomarker of tumor progression and a poor prognosis (68). 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

Cohort 1 (%)Cohort 2 (%)
Total148138p value
Age, median (range) (y)67.5 (34–94)70.0 (41–88)0.1583
SexMan90 (60.8)81 (58.7)0.7193
Woman58 (39.2)57 (41.3)
LocationRight colon30 (20.3)32 (23.2)0.4271
Left colon64 (43.2)49 (35.5)
Rectum54 (36.5)57 (41.3)
pTpT3129 (86.5)111 (80.4)0.1474
pT419 (13.5)27 (19.6)
StageII71 (48.0)69 (50.0)0.8129
III77 (52.0)69 (50.0)
Histological typeWDA17 (11.5)26 (18.8)0.0681
MDA121 (81.8)109 (79.9)
PDA2 (1.4)1 (0.7)
PAP6 (4.1)1 (0.7)
MUC2 (1.4)1 (0.7)
Lymphatic invasionPositive130 (87.8)129 (93.5)0.4404
Negative18 (12.2)9 (6.5)
Venous invasionPositive129 (87.2)128 (92.8)0.4534
Negative19 (12.8)10 (7.2)
Tumor buddingLow117 (79.1)108 (78.3)0.8861
High31 (20.9)30 (21.7)
Desmoplastic reactionMature65 (43.9)61 (44.2)0.9877
Intermediate53 (35.8)49 (35.5)
Immature30 (20.3)28 (20.3)
Disease-free survival, median (range) (d)1857 (33–3196)1835 (93–3308)
Overall survival, median (range) (d)3077 (52–3196)2195 (93–3308)

Clinicopathological findings in stage II and III colorectal cancer.

WDA, well-differentiated adenocarcinoma; MDA, moderately differentiated adenocarcinoma; PDA, poorly differenced adenocarcinoma; PAP, papillary carcinoma; MUC, mucinous carcinoma.

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

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

Subgroup 1 (%)Subgroup 2 (%)Subgroup 3 (%)p value
Total327442
Age median (range) (y)69.0 (43–94)67.0 (34–92)67.5 (42–88)0.1971
SexMan19 (59.4)45 (60.8)26 (61.9)1.0000
Woman13 (41.6)29 (39.2)16 (38.1)
LocationRight colon7 (21.9)15 (20.3)8 (19.0)0.9516
Left colon15 (46.9)30 (40.5)19 (45.2)
Rectum10 (31.2)29 (39.2)15 (35.7)
pTpT328 (87.5)63 (85.1)38 (90.5)0.7678
pT44 (12.5)11 (14.9)4 (9.5)
StageII13 (40.6)37 (50.0)21 (50.0)0.6456
III19 (59.4)37 (50.0)21 (50.0)
Histological typeWDA1 (3.1)10 (13.5)6 (14.3)0.1125
MDA28 (87.5)60 (81.1)33 (78.6)
PDA 0 (0.0)0 (0.0)2 (4.8)
PAP3 (9.4)3 (4.1)0 (0.0)
MUC0 (0.0)1 (1.4)1 (2.4)
Lymphatic invasionPositive27 (84.4)65 (87.8)38 (90.5)0.7149
Negative5 (15.6)9 (12.2)4 (9.5)
Venous invasionPositive29 (90.6)61 (82.4)39 (92.9)0.2694
Negative3 (9.4)13 (17.6)3 (7.1)
Tumor buddingLow28 (87.5)58 (78.4)31 (73.8)0.3738
High4 (12.5)16 (21.6)11 (26.2)
Desmoplastic reactionMature8 (25.0)36 (48.6)21 (50.0)0.0508
Intermediate11 (34.4)27 (36.5)15 (35.7)
Immature13 (40.6)11 (14.9)6 (14.3)
Disease-free survivalPositive20 (62.5)*15 (20.3)*15 (35.7)0.0002
Negative12 (37.5)59 (79.7)27 (64.3)
Overall survivalDead10 (31.6)10 (13.5)7 (16.7)0.0999
Alive20 (62.5)64 (86.5)35 (83.3)

Clinicopathological variables according to each subgroup in the first cohort.

WDA, well-differentiated adenocarcinoma; MDA, moderately differentiated adenocarcinoma; PDA, poorly differenced adenocarcinoma; PAP, papillary carcinoma; MUC, mucinous carcinoma; *p < 0.0001.

Figure 2

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

Variablesa. Univariate analysisb. Multivariate analysisc. Univariate analysisd. Multivariate analysis
HR95%CIp valueHR95%CIp valueHR95%CIp valueHR95%CIp value
Disease-free survivalOverall survival
SexWoman vs man1.4200.807-2.4770.22071.4560.674-3.1260.3335
LocationRectum vs Left colon1.4880.784-2.8640.22381.3360.549-3.4290.5257
Right colon vs Left colon1.6710.778-3.4760.18211.4110.492-3.9500.5104
Right colon vs Rectum1.1230.534-2.2600.75081.0560.393-2.6270.9084
pTpT4 vs pT31.4880.646-3.0070.32612.4770.967-5.6300.0581
StageIII vs II2.3191.299-4.3220.00411.8160.981-3.5200.05772.5751.139-6.5700.02222.2640.972-5.9230.0586
Histological typeMUC vs WDA10.4431.857-58.7260.02539.5561.537-59.3960.03504.9940.448-55.6760.2455
Lymphatic invasionPositive vs Negative2.4440.895-10.0690.08653.5960.762-64.2380.1226
Venous invasionPositive vs Negative1.4730.643-4.2550.38781.6770.499-10.4310.4493
Tumor buddingHigh vs Low1.3860.707-2.5480.32751.5020.636-3.2990.3383
Desmoplastic reactionImmature vs Intermediate1.7660.894-3.4460.10031.7880.725-4.2570.2004
Immature vs Mature3.0231.486-6.1830.00252.1901.012-4.7930.04673.6601.316-10.9470.01322.4730.821-7.8700.107
Intermediate vs Mature1.7110.869-3.4350.12002.0480.793-5.8930.1407
SubgroupSubgroup 1 vs 23.6261.860-7.2240.00023.0121.512-6.1270.00182.6431.082-6.4590.03352.0820.825-5.2610.1185
Subgroup 1 vs 32.0241.036-4.0490.03921.7100.819-3.6720.15382.0580.789-5.6760.1393
Subgroup 3 vs 21.7910.868-3.6970.11351.2840.466-3.3470.6154

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

WDA, well-differentiated adenocarcinoma; MUC, mucinous carcinoma; HR, hazard ratio; 95%CI, 95% confidence interval.

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

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

Variablesa. Univariate analysisb. Multivariate analysisc. Univariate analysisd. Multivariate analysis
HR95%CIp valueHR95%CIp valueHR95%CIp valueHR95%CIp value
Disease-free survivalOverall survival
StageIII vs II2.3191.299-4.3220.00411.7600.943-3.4300.07622.5751.139-6.5700.02222.2370.959-5.8560.0629
Histological typeMUC vs WDA10.4431.857-58.7260.025315.0972.505-90.9730.01325.0010.929-33.6670.2451
Desmoplastic reactionImmature vs Intermediate1.7660.894-3.4460.10031.7880.725-4.2570.2004
Immature vs Mature3.0231.486-6.1830.00252.3581.115-5.0330.02503.6601.316-10.9470.01323.0761.085-9.3520.0347
Intermediate vs Mature1.7110.869-3.4350.12002.0480.793-5.8930.1407
α-SMAPositive vs Negative0.6840.250-2.8210.54680.6280.187-3.9050.5532
CD10Positive vs Negative1.6900.798-3.2540.16061.4120.472-3.4530.5025
PodoplaninPositive vs Negative1.0190.583-1.8210.94771.0390.486-2.3070.9215
FSP1Positive vs Negative0.6770.209-4.1500.61110.3570.072-6.4610.3876
FAPPositive vs Negative0.5010.120-1.3920.20670.2780.154-1.334 0.1270
Tenascin-CPositive vs Negative2.6941.293-6.5570.00652.3171.068-5.8700.03242.0520.840-6.1300.1201
PDGFR-βPositive vs Negative0.9880.550-1.8630.96800.7640.351-1.7860.5176
ZEB1Positive vs Negative1.2690.584-3.3250.57302.0940.623-13.0180.2638
TWIST1Positive vs Negative1.5120,851-2.6390.15541.1850.533-2.5400.6673
Ki-67Positive vs Negative1.1390.653-2.0110.64771.2450.579-2.7360.5754
p53Positive vs Negative0.6290.359-1.0980.10230.5040.231-1.0770.0767
MMP7Positive vs Negative1.6210.890-3.0690.10961.2050.560-2.7340.6383
β-cateninPositive vs Negative0.7700.401-1.3950.39740.9790.419-2.1270.9578
E-cadherin*
HIF1-αPositive vs Negative0.4610.101-8.1720.49610.2970.062-5.3440.3187

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

WDA, Well-differentiated adenocarcinoma; MUC, mucinous carcinoma; HR, Hazard ratio; 95%CI, 95%confidence interval. *could not analyze, why all cases were positive expression of the marker.

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

Variablesa. Univariate analysisb. Multivariate analysisc. Univariate analysisd. Multivariate analysis
HR95%CIp valueHR95%CIp valueHR95%CIp valueHR95%CIp value
Disease-free survivalOverall survival
SexWoman vs Man0.8190.430-1.5120.52830.9190.395-2.0480.8320
LocationRectum vs Left colon1.9820.976-4.3450.05881.6170.639-4.6030.3180
Right colon vs Left colon1.0940.417-2.7740.85051.0880.313-3.6130.8894
Right colon vs Rectum1.8130.850-4.3100.12781.4860.560-4.6320.4393
pTpT4 vs pT32.2791.142-4.3070.02081.5270.714-3.1120.26622.6301.065-5.9860.03691.9360.746-4.6700.1667
StageIII vs II3.1141.633-6.3340.00041.9910.977-4.3010.05814.5991.847-13.8850.00072.6190.980-8.3210.0552
Histological typeWDA vs MDA0.6270.264-1.4920.26570.7870.269-2.3030.6533
Lymphatic invasionPositive vs Negative4.3340.944-76.8330.06202.4470.516-43.7770.3125
Venous invasionPositive vs Negative4.8551.057-86.0760.04014.0730.838-73.4230.09042.7220.574-48.6900.2501
Tumor buddingHigh vs Low2.2911.170-4.2850.01691.6670.825-3.2370.14972.0900.847-4.7520.1052
Desmoplastic reactionImmature vs Intermediate1.5200.712-3.2040.27361.2270.451-3.1180.6754
Immature vs Mature2.3221.076-4.9760.03240.9930.415-2.3400.98673.2931.092-10.2450.03481.5960.493-5.3050.4312
Intermediate vs Mature1.5270.733-3.2020.25582.6841.021-7.7930.04521.8950.696-5.6860.2143
α-SMAPositive vs Negative1.4230.437-8.7360.60781.5800.333-28.2640.6303
CD10Positive vs Negative0.6570.267-1.3910.28830.7040.295-1.9410.4698
PodoplaminPositive vs Negative0.7520.399-1.4960.40200.5640.251-1.3430.1872
FSP1Positive vs Negative0.3550.108-2.1850.21810.2960.061-5.1390.3172
FAPPositive vs Negative0.8190.425-1.5200.53234.0450.853-72.3580.0866
Tenascin-CPositive vs Negative5.0252.164-14.6210.00013.9731.656-11.7950.00124.5271.559-19.1720.00363.1881.038-13.8820.0421
PDGFR-βPositive vs Negative1.2900.685-2.5660.43951.4160.611-3.6610.4288
ZEB1Positive vs Negative0.6830.248-2.8250.54630.5620.166-3.5100.4713
TWIST1Positive vs Negative1.0670.581-1.9710.83311.2370.553-2.8180.6028
Ki-67Positive vs Negative1.6670.885-3.3160.11572.9521.187-8.9060.01852.1140.834-6.4650.1191
p53Positive vs Negative1.3620.730-2.6630.33781.2550.552-3.0960.5955
MMP7Positive vs Negative0.6510.354-1.1990.16710.8230.368-1.8750.6358
β-cateninPositive vs Negative2.1191.143-4.0840.01691.5020.796-2.9400.21211.8650.830-4.4400.1323
E-cadherin*
HIF1-α*

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

WDA, Well-differentiated adenocarcinoma; MUC, mucinous carcinoma; HR, Hazard ratio; 95%CI, 95%confidence interval. *Could not analyze, why all cases were positive expression of the marker.

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 (1618). 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 (1921). 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.

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.

Statements

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.

Acknowledgments

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

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.

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

    SiegelRLMillerKDGoding SauerASauerAFedewaSAButterlyLFet al. Colorectal Cancer Statistics, 2020. CA Cancer J Clin (2020) 70:145–64. doi: 10.3322/caac.21601

  • 2

    ChenFWSundaramVChewTALadabaumU. 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:72837.e3. doi: 10.1016/j.cgh.2016.10.038

  • 3

    MocciaFToloneSAllariaANapolitanoVRosaDIlariaFet 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

  • 4

    PuppaGSonzogniAColombariRPelosiG. 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

  • 5

    NagtegaalIDOdzeRDKimstraDParadisVRuggeMSchirmacherPet al. The 2019 WHO Classification of Tumours of the Digestive System. Histopathology (2020) 76:182–8. doi: 10.1111/his.13975

  • 6

    LugliAKirschRAjiokaYBosmanFCathomasGDawsonHet 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

  • 7

    WatanabeTMuroKAjiokaYHashiguchiYItoYSaitoYet 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:134. doi: 10.1007/s10147-017-1101-6

  • 8

    MitrovicBSchaefferDFRiddellRHKirschR. Tumor Budding in Colorectal Carcinoma: Time to Take Notice. Mod Pathol (2012) 25:1315–25. doi: 10.1038/modpathol.2012.94

  • 9

    UenoHIshiguroMNakataniEIshikawaTUetakeHMurotaniKet 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

  • 10

    ZlobecILugliA. Epithelial Mesenchymal Transition and Tumor Budding in Aggressive Colorectal Cancer: Tumor Budding. Oncotarget (2010) 1:651–61. doi: 10.18632/oncotarget.199

  • 11

    TommeleinJVersetLBoterbergTDemetterPBrackeMDe WeverO. Cancer-Associated Fibroblasts Connect Metastasis-Promoting Communication in Colorectal Cancer. Front Oncol (2015) 5:63. doi: 10.3389/fonc.2015.00063

  • 12

    SugaiTYamadaNEizukaMSugimotoRUesugiNOsakabeMet 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

  • 13

    SugaiTUesugiNKitadaYYamadaNOsakabeMEizukaMet 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

  • 14

    HashiguchiYMuroKSaitoYYamadaNOsakabeMEizukaMet 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:142. doi: 10.1007/s10147-019-01485-z

  • 15

    SasakiKSugaiTIshidaKOsakabeMAmanoHKimuraHet 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:18. doi: 10.1016/j.humpath.2018.03.006

  • 16

    ZhangWShiXPengYWuMZhangPXieRet 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

  • 17

    MansourRNEnderamiSEArdeshirylajimiAFooladsazKFathiMGanjiSM. 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

  • 18

    van WykHCRoseweirAAlexanderPParkJHHorganPGMcMillanDCet 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

  • 19

    SunZSchwenzerARuppTMurdamoothooDVeglianteRLefebvreOet 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

  • 20

    YangZZhangCQiWCuiCCuiYXuanY. 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

  • 21

    MurakamiTKikuchiHIshimatsuHIinoIHirotsuAMatsumotoTet 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

  • 22

    YangZZhangCFengYQuanMCuiYXuanY. 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

  • 23

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

  • 24

    GulinoAFerrettiEDe SmaeleE. Hedgehog Signalling in Colon Cancer and Stem Cells. EMBO Mol Med (2009) 1:300–2. doi: 10.1002/emmm.200900042

  • 25

    SteccaBMasCClementVZbindenMCorreaRPiguetVet 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

  • 26

    OrendGChiquet-EhrismannR. Tenascin-C Induced Signaling in Cancer. Cancer Lett (2006) 244:143–63. doi: 10.1016/j.canlet.2006.02.017

  • 27

    SariINPhiLTHJunNWijayaYTLeeSKwonHY. Hedgehog Signaling in Cancer: A Prospective Therapeutic Target for Eradicating Cancer Stem Cells. Cells (2018) 7:208. doi: 10.3390/cells7110208

  • 28

    WuJYXuXFXuLNiuPQWangFHuGYet 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

  • 29

    TakahashiYSawadaGKurashigeJMatsumuraTUchiRUeoHet al. Tumor-Derived Tenascin-C Promotes the Epithelial-Mesenchymal Transition in Colorectal Cancer Cells. Anticancer Res (2013) 33:192.

  • 30

    SagaertXVanstapelAVerbeekS. Tumor Heterogeneity in Colorectal Cancer: What Do We Know So Far? Pathobiology (2018) 85:7284. doi: 10.1159/000486721

  • 31

    JangBGKimHSChangWYBaeJMOhHJWenXet al. Prognostic Significance of Stromal GREM1 Expression in Colorectal Cancer. Hum Pathol (2017) 62:5665. doi: 10.1016/j.humpath.2016.12.018

  • 32

    OhHJBaeJMWenXYChoNYKimJHKangGH. 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

Summary

Keywords

cancer-associated fibroblast, colorectal cancer, cluster analysis, prognostic marker, tenascin-C

Citation

Hashimoto M, Uesugi N, Osakabe M, Yanagawa N, Otsuka K, Kajiwara Y, Ueno H, Sasaki A and Sugai T (2021) Expression Patterns of Microenvironmental Factors and Tenascin-C at the Invasive Front of Stage II and III Colorectal Cancer: Novel Tumor Prognostic Markers. Front. Oncol. 11:690816. doi: 10.3389/fonc.2021.690816

Received

04 April 2021

Accepted

02 August 2021

Published

19 August 2021

Volume

11 - 2021

Edited by

Alex Nicolas Gordon-Weeks, University of Oxford, United Kingdom

Reviewed by

Juan Carlos De Vicente Rodríguez, University of Oviedo, Spain; Eva Andreuzzi, Aviano Oncology Reference Center (IRCCS), Italy

Updates

Copyright

*Correspondence: Tamotsu Sugai,

This article was submitted to Surgical Oncology, a section of the journal Frontiers in Oncology

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics