Abstract
Background/Aim: Early-stage colorectal cancer (CRC) carries a wide range of survival probabilities. Novel biomarkers in this setting are eagerly awaited. Cancer stem cells (CSCs) are considered one of the reasons for treatment failure. This study sought to determine whether activation of pathways governing the function of CSC's could correlate with treatment outcomes. Materials and Methods: Tumor specimens from 325 patients were analyzed with immunohistochemistry (IHC) for Hedgehog and Notch pathway activation and results were correlated with prognosis. Results: Positive Notch3 protein expression was an unfavorable prognostic factor for disease-free survival (DFS) and overall survival (OS) (HR=2.43, p=0.024 and HR=2.56, p=0.028, respectively). Activation of the Shh pathway showed univariately longer DFS (HR=0.49, p=0.032). Possible crosstalk between the two pathways was indicated. No further associations between pathway activation and outcome were evident. Conclusion: Apart from Notch 3, activation of the pathways, as indicated by IHC expression of their components, did not result in differences in terms of DFS or OS.
Colorectal cancer (CRC) is the second leading cause of cancer-related deaths worldwide, with an estimate of 881,000 deaths occurring in 2018 (1). Despite major advances in early diagnosis, surgical treatment and chemotherapy agents, still, a significant proportion of early-stage cases will develop metastases and ultimately succumb from their disease. One of the reasons for this unfortunate evolution is the presence of cancer stem cells (CSCs), which, among other properties, seem to be resistant to current treatments. Fundamental to the discovery of new treatments targeting CSCs is our full and in depth understanding of the biology behind them (2).
Among the molecular pathways governing the function of CSCs are Hedgehog and Notch. The Hedgehog pathway, first described in 1980 (3), plays an important role during embryogenesis and the development of several organs and tissues, among which the gastrointestinal tract. Deregulation of this pathway is associated with developmental deficits in embryos and cancer development in adult patients. Its role is well established in cases of holoprosencephaly, Gorlin syndrome, basal cell carcinomas and medulloblastomas (4). Mutations or overexpression of the pathway components can result in cancer development (5-8). The Hedgehog family consists of three proteins, sonic (Shh), indian (Ihh) and desert (Dhh), which are the ligands of 12-pass transmembrane protein receptors known as Patched. In the absence of ligand interaction, patched receptors keep the 7-pass transmembrane protein smoothened (Smo) inactive. Upon ligand binding, Smo activates the downstream glioma-associated oncogene homologue isoforms (Gli-1,-2,-3) in the cytoplasm, which then translocate to the nucleus, inducing the expression or downregulation of their target genes (mainly regulatory of the cell cycle or proapoptotic genes) (9). Hedgehog, besides other functions, has been found to be involved in tumorigenic development and cancer metastases both in vitro and in vivo (10). It also seems to play an important role in stem cell survival and self-renewal. Unfortunately, the results are not consistent and findings from several studies are contradictory (11-15).
The Notch Pathway genes were first described in 1917 upon an observation of a genetic anomaly causing the formation of notches in the wings of Drosophila melanogaster (16). Subsequent studies revealed that this pathway plays a crucial role in pattern formation during embryonic development, through specific cell to cell interactions that determine the fate of each cell and inducing the formation of different cells (17). The pathway, in brief, consists of five ligands [Jagged-1, Jagged-2, Delta-like-1 (Dll-1), Delta-like-3 (Dll-3) and Delta-like-4 (Dll-4)], and four receptors (Notch1, Notch2, Notch3 and Notch4) (18). Deregulation of the pathway results in developmental deficits in humans, like the Allagile syndrome (19), cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) (20), Hajdu-Cheney syndrome (21) and spondylocostal dysostosis (22), but also in neoplasia, like T-ALL (23) and diffuse large B-cell lymphoma (24).
Selected patient and tumor characteristics.
Activation of the pathway starts upon the interaction of the ligands with the receptors. This leads to a proteolytic cleavage of the extracellular part of the receptor followed by a second cleavage in the intracellular part that releases the Notch Intracellular Domain (NICD), which then enters the nucleus and forms a transcriptional activation complex that promotes the expression of a number of genes, like Myc, CyclinD, hes1 etc. 1 (25-28). Notch can act both as a major oncogene, associated with tumor progression and metastasis, as well as a tumor suppressor (10). It has been also shown that Notch pathway activation contributes to the development and survival of CSCs (29, 30). Especially in colon cancer, an interplay between Notch and Wnt has been documented and it seems that Notch's activation is imperative for the adenoma-carcinoma progression of APC mutant patients (31).
REMARK diagram.
The role of Hedgehog and Notch pathways in solid tumors, among which in CRC, has been studied and reviewed quite extensively (32, 33). These are master developmental pathways, commonly activated in many tumors in which they not only play a crucial role in tumor initiation, but also in primary and metastatic tumor development and cancer stem cell regulation (10).
This study sought to correlate the expression of Notch2, Notch3, Jag1, Shh, Smo, Patched-1 and Gli1 with various clinico-pathological parameters and patients' outcome with the intention to identify meaningful associations.
Materials and Methods
Clinical study. A total of 325 colorectal cancer patients treated with adjuvant chemotherapy within a randomized trial setting (34) were included in the present study (Figure 1). The clinical protocol was approved by the Institutional Review Boards (IRBs) in participating institutions and by the National Organization for Medicines (53386/14-10-05). The trial was registered in the Australian New Zealand Clinical Trials Registry and allocated the following Registration Number: ANZCTR 12610000509066. The present translational project was approved by the Bioethics Committee of the Aristotle University of Thessaloniki School of Health Sciences, Faculty of Medicine (1/8-11-2012). Written informed consent for participation in the trial was obtained from all the patients and optionally a separate informed consent was obtained for providing biological material for research purposes.
Selected patient and tumor characteristics at baseline are presented in Table I. More than half of the patients were men (55.4%) and the median age at diagnosis was 64.5 years. The majority of the tumors were of G2 histological grade (N=235, 72.3%) and of stage III (N=202, 62.2%). Regarding the primary site, 137 tumors were located in the left colon (42.2%), 98 in the right colon (30.2%) and 87 in the rectum (26.8%). Formalin-fixed paraffin-embedded (FFPE) tumor tissues were prospectively collected and retrospectively analyzed for the expression of important biomarkers and the identification of mutations in the KRAS and BRAF genes.
KRAS/BRAF tumor genotyping. Screening of tumor DNA for KRAS mutations (coding exons 1, 2 and 3, CCDS 8702.1) and BRAF V600E mutation (coding exon 15, CCDS 5863.1) was performed with dd-sequencing of PCR products, amplified with appropriate M13-coupled nested primers, as previously described in detail (35, 36).
Immunohistochemical (IHC) evaluation and scoring system. For Notch2, Notch3 and Jag1 staining intensity and percentage of positive cells were evaluated by an experienced pathologist (TK). Regarding the cutoffs, staining intensity was graded as 1+, 2+ and 3+, while the percentage of positive cells was scored into four categories (0 for 0%, 1 for 1-33%, 2 for >33-66% and 3 for >66-100%). The product of the intensity and percentage categories (range=0-9) was used as the final score and classified as negative (0-4) and positive (5-9). Similarly, for Gli, Patch, Shh and Smo intensity and positivity rates were evaluated. Staining intensity was scored as 1+, 2+ and 3+, while the percentage of positive cells was graded as 0 for <5%, 1 for 5-25%, 2 for >26-50%, 3 for >51-75% and 4 for >75%. The product of the intensity and percentage categories (range=0-12) was used as the final score and classified as negative (0-4) and positive (5-12). Representative staining images of Notch and Hedgehog pathway markers are presented in Figure 2.
Representative staining images of Notch and Hedgehog pathway markers in colorectal carcinomas. All magnifications ×400. Smo: Smoothened; SHH: Sonic Hedgehog. For each marker: left microphotograph, positive; right microphotograph, negative.
Tumors were classified, according to IHC positivity for MLH1, PMS2, MSH2 and MSH6, in MMR proficient (pMMR) if all proteins were expressed and MMR deficient (dMMR) in case of null expression of at least one protein.
Statistical analysis. Continuous variables are presented as medians with the corresponding range and categorical variables as frequencies with percentages. The Chi-square or Fisher's exact (where appropriate) test was used for group comparisons of categorical data.
Overall survival (OS) was defined as the time (in months) from the date of colorectal cancer diagnosis to the date of death or last contact, while disease-free survival (DFS) was defined as the time (in months) from the date of diagnosis to the first documented progression, death without prior documented progression or last contact, whichever occurred first. Surviving patients (for OS and DFS) and patients without progression (for DFS) were censored at the date of last contact. Survival distributions were estimated using the Kaplan–Meier method and compared across groups with the log-rank test. The associations between the factors of interest (see below) and progression/mortality rates were evaluated with hazard ratios estimated with univariate and multivariate Cox proportional hazard regression models.
The following parameters were studied in relation to DFS/OS: 1) clinicopathological: age, body mass index (BMI), gender (man, woman), ECOG performance status (0, 1-2), TNM stage (I-II, III), histological grade (G1, G2, G3), obstruction (yes, no), perforation (yes, no), primary site (left colon, right colon, rectum); 2) IHC markers considered as 2-level categorical variables (negative vs. positive) using the cut-offs previously described in the “IHC evaluation and scoring system” section: Jag1, Notch2, Notch3, Gli, Patched, Smo, Shh; MMR status (deficient, proficient); 3) mutational markers considered as 2-level categorical variables (mutant, wild-type): BRAF and KRAS. In addition, combined variables were examined with the following categories: Jag1-Notch pathway: (i) Jag1-positive or Noth3-positive and Notch2-positive, indicating the activated state of the pathway; (ii) Notch2-positive only; (iii) other. Shh-pathway: (i) (Gli-positive and Shh-positive) or Smo-positive or Patched-positive, indicating the activated state of the pathway; (ii) other.
Bar plots showing the association of Jag1 protein expression with histological grade (A); Patched protein expression with primary site (B).
In the multivariate analyses, each of the markers of interest that was found to be significant or revealed a trend towards significance in the univariate analyses (p<0.10) was adjusted for the significant clinicopathological parameters identified by univariate analyses with respect to DFS and OS.
All analyses were performed in the entire cohort of patients. No adjustment for multiple comparisons was performed based on the exploratory nature of the study. Analyses were performed using the SAS (SAS for Windows, version 9.3, SAS Institute Inc., Cary, NC) and R statistical software (37). Statistical significance was set at 2-sided p=0.050.
Results
Among 310 tumors informative for Jag1, 321 informative for Notch2 and 316 informative for Notch3, 63 (20.3%), 241 (75.1%) and 11 (3.5%) were positive, respectively. Similarly, among 303 tumors informative for Gli, 318 informative for Patched, 310 informative for Shh and 294 informative of Smo, 139 tumors (45.9%), 25 (7.9%), 292 (94.2%) and 15 (5.1%) were positive, respectively. In total, 59 patients had the Jag-Notch pathway activated, while more than half of the patients with informative tumors (149 patients) had the Shh pathway activated. Regarding MMR status, among 288 tumors informative for all MMR proteins, 254 (88.2%) were pMMR. In addition, among 312 NGS-informative tumors, 107 had KRAS mutations and 15 carried BRAF mutations.
Jag1 was significantly associated with Gli (chi-square, p=0.018), Notch2 (p=0.024) and Notch3 (p=0.017). Jag1-positive tumors, compared to Jag1-negative, were more frequently Gli-positive (59.3% vs. 42.1%), Notch2-positive (87.1% vs. 73.5%) and Notch3-positive (8.1% vs. 2%). Additionally, Shh-positive tumors were more frequently Gli-positive compared to Shh-negative tumors (47% vs. 17.6%, p=0.018), while Smo-positive tumors were more frequently Patched-positive compared to Smo-negative tumors (26.7% vs. 6.1%, p=0.003). No significant association was observed between any of the IHC markers and BRAF/KRAS mutational status.
Jag1 IHC expression was significantly associated with histological grade (p=0.002). Jag1-positive tumors were more frequently of G1 histological grade compared to Jag1-negative tumors (14.8% vs. 3.3%, Figure 3A). Additionally, Patched IHC expression was associated with primary site (p=0.018). Patched-positive tumors, compared to Patched-negative, were more frequently located in the right colon (44.0% vs. 29.7%) and less frequently located in the left colon (16.0% vs. 45.2%) (Figure 3B). No other significant associations were observed. No significant associations were observed between selected clinicopathological parameters and Jag1-Notch or Shh pathway.
Kaplan–Meier curves with respect to DFS (A) and OS (B) based on Notch3 protein expression.
At a median follow-up of 122.6 months (95%CI=120.9-124.4), 87 patients had died (26.8%) and 106 had progressed (32.6%). The median OS and DFS had not been reached yet at the time of the analysis. Regarding the effect of the study markers on patients' outcome, only positive Notch3 protein expression had a significant, unfavorable effect on both DFS and OS (Figure 4). Notch3-positive patients presented with shorter DFS and OS compared to those with Notch3- negative tumors (DFS HR=2.53, 95%CI=1.17-5.45, Wald's p=0.018; OS HR=2.67, 95%CI=1.16-6.13, p=0.021). Upon adjustment for age and stage (parameters that showed statistical significance in the univariate analyses for DFS and OS), positive Notch3 expression retained its unfavorable prognostic significance for DFS and OS (Table II). However, the power of these results is limited by the small number of patients with Notch3-positive tumors and should, therefore, be interpreted with care until further validated in larger cohorts. Similarly, when examining the effect of Jag1-Notch pathway and Shh-pathway activation on DFS and OS no significant associations were observed (DFS log-rank p=0.91 and p=0.14 and OS log-rank p=0.86 and p=0.21, respectively). Of note, patients with activated Shh pathway seemed to fair better in terms of DFS after the first two years. The hazard rate functions of each category of the Shh pathway parameter (activated, other) for DFS were assessed. Patients in both groups seemed to have the same hazard ratio for the first 20 months of follow-up (i.e. activated vs. other, HR is approximately 1). However, between 20 and 70 months, patients with activated Shh pathway seemed to have a lower hazard ratio compared to patients without an activated pathway (i.e., activated vs. other, HR <1). In comparison, in the period between 70 and 100 months the hazard ratio of the activated group increased while, simultaneously, the hazard ratio in the non-activated group decreased and the relationship was reversed (i.e., activated vs. other, HR >1). Therefore, the effect of Shh pathway activation was further examined, excluding patients with a DFS lower than 20 and greater than 70 months. A significant association was observed between Shh pathway activation and reduced risk for progression (HR=0.49, 95%CI=0.26-0.94, p=0.032), which however was not retained in the multivariate analysis adjusting for age and stage (HR=0.54, 95%CI=0.26-1.10, p=0.088, Table II).
Hazard ratios and 95%CIs estimated by univariate and multivariate Cox regression models with respect to DFS and OS.
Discussion
In this observational study, 325 CRC patients and their respective CRC tumor specimens were analyzed retrospectively for important and relevant to CRC biomarkers (KRAS, BRAF, MMR status), as well as, for the evaluation of two major developmental pathways, Notch and Hedgehog, and the results were correlated with clinicopathological parameters and outcome.
The strongest point of our study regards Notch3. Patients whose tumors expressed high levels of Notch3 were shown to have worse outcomes, meaning, approximately 2.5 times reduced DFS and OS as this is evident from the respective HRs. These findings agree with current evidence in the international literature, where an association between high levels of Notch3 and a more aggressive malignant cell phenotype has been described (38). Nevertheless, due to the small number of patients in the Notch3-positive group (N=11), our results should be interpreted with caution. Furthermore, our analysis showed a statistically significant association between Notch3 positivity and the existence of BRAF mutation in CRC tumors, but again, due to the very limited numbers, this observation should be interpreted cautiously.
In our study 75.1% of the patients were IHC Notch2 positive. Indeed, high Notch2 expression has been observed in the majority of GI malignant cells in various publications (39-41). Unfortunately, this observation was not correlated with outcome in our series, as it was elsewhere (42, 43).
The role of Jagged-1 in CRC still remains elusive. Overexpression has been correlated with bad prognosis, especially when combined with high levels of Notch2, and with increased probability for relapse and lymph node metastasis (44). Furthermore, Jagged-1 expression intensity has been correlated with higher grade, higher TNM stage, depth of infiltration and lymph node metastasis (45), something that was not replicated in our analysis, where more patients were included and Jagged-1 expression was statistically significantly associated with a lower grade at diagnosis.
In several studies, Gli1 expression level in CRC tissues varies from 30% to 79%, with staining being either cytoplasmic or nuclear (46-51). In our study, 45,87% of the tumors stained positive for Gli-1. Expression of Gli-1 protein in colon cancer tissues seems to play a key role in the occurrence and development of colon cancer, while high Gli-1 expression may promote postoperative liver metastasis development (52). In our study, no significant associations between Gli-1 and prognosis were found.
Crosstalk of Notch pathway with other embryonic pathways, such as WNT, has been implied in the literature (39, 53), but little is known about any associations with the Hedgehog pathway. In our study, statistically significant associations were observed between Jagged-1 expression and Gli-1 expression. Since Gli-1 is the primary effector molecule of the Hedgehog pathway activation, this association may indicate a possible crosstalk between the two pathways, a notion that needs further exploration.
Patched-1 overexpression has been correlated with moderately differentiated tumors (47) and nodal status (50), as well as, with worse DFS and OS in patients with early-stage cancer that underwent curative treatment (47, 50). Contradictory results were reported in a small study of 19 patients, where low levels of Patched-1 mRNA expression were associated with a higher risk of metastasis (54). In our study, an association between Patched-1 expression and location of the primary tumor was observed, with the right sided tumors showing high expression, whereas the left sided tumors showed low expression. No association of Patched-1 expression with stage, DFS or OS was detected.
Smo, like Patched-1, is one of the two transmembrane proteins of the Hedgehog pathway that along with the downstream glioma-associated transcription factors are responsible for the activity of the pathway. High levels of this molecule have been described in colorectal cancer cells. Higher Smo expression has been correlated with an increased probability of lymph node metastasis and higher T stage (48). Apart from the expected correlation of Smo-positive, with Gli- and Patched-1-positive tumors, no other significant associations of Smo were observed in our series.
It has been reported that epithelial malignancies of the adult gut are associated with abnormal SHH/SHH expression (11, 55, 56). In concordance with this observation, in our study, 94.2% of the CRC tumors overexpressed Shh. Furthermore, in a retrospective analysis it was shown that Shh overexpression is associated with poorer DFS and OS in patients with non-metastatic colorectal cancer that underwent curative surgery (50). In our series, patients with an activated Shh pathway seemed to fair better in terms of DFS after the first two years, during which the majority of relapses have occurred. In contrast, in the follow-up period of 70 to 100 months, the HR of the Shh activated group was increased. When excluding from the analysis patients with a DFS lower than 20 months (poor prognosis patients) and greater than 70 months (good prognosis patients) Shh activation characterized the patients that would experience a reduced risk for relapse; this is an observation that could justify further exploration.
Conclusion
Although the Notch and Hedgehog pathways have been shown in several studies to play a role in CRC development and progression, data are still contradictory and unclear (57). In our study, activation of the two examined pathways, as this is indicated by IHC overexpression of their components, did not result in significant associations in terms of outcome in early-stage CRC patients, with the exception of Notch3. Due to the retrospective nature of our analysis and the limited number of Notch3-positive tumors, our data are not enough to promote Notch3 to a putative biomarker indicative of early progression. Nevertheless, it is possible that potential biomarkers may be “hiding” alongside these two pathways, and further research is warranted in order to depict their exact role in the metastatic evolution in CRC patients.
Acknowledgements
The Authors are indebted to all patients and their families for the provision of biological material for research purposes. The Authors would also like to thank Emily Daskalaki for excellent technical assistance, Helen Sgouramalli, MSc for block sectioning, Maria Moschoni for coordinating the data management and Stella Dallidou for secretarial assistance.
This study was supported by an internal Hellenic Cooperative Oncology Group (HeCOG) translational research grant (HE TRANS_CRC) and by a Hellenic Society for Medical Oncology (HeSMO) grant. The funders played no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Footnotes
Authors' Contributions
Study Conception: GR, TK, ND, GF; Study Design: GR, TK, KM, GAK, IK, VK, ND, GF; Acquisition of data: GR, TM, HPK, GR, ES, KP, EP, PP, GK, DM, APB, CC, KP, DP, GF; Analysis and interpretation of data: GR, TK, ZS, KM, GAK, IK, VK, SC, ND, GF; Drafted and revised the manuscript: GR, TK, ZS, KM, GAK, IK, VK, SC, ND, GF; All Authors have revised and approved the final version.
This article is freely accessible online.
- Received February 14, 2019.
- Revision received March 6, 2019.
- Accepted March 14, 2019.
- Copyright© 2019, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved