Ligand-Independent Epidermal Growth Factor Receptor Overexpression Correlates with Poor Prognosis in Colorectal Cancer

Article information

Cancer Res Treat. 2018;50(4):1351-1361
Publication date (electronic) : 2018 January 17
doi : https://doi.org/10.4143/crt.2017.487
1Department of Diagnostic Pathology, Samkwang Medical Laboratories, Seoul, Korea
2Department of Pathology, Seoul National University College of Medicine, Seoul, Korea
3Department of Pathology, Seoul National University Bundang Hospital, Seongnam, Korea
4Department of Pathology, Kyungpook National University Hospital, School of Medicine, Kyungpook National University, Daegu, Korea
5Department of Surgery, Seoul National University Bundang Hospital, Seongnam, Korea
Correspondence: Hye Seung Lee, MD, PhD Department of Pathology, Seoul National University Bundang Hospital, 82 Gumi-ro 173beon-gil, Bundang-gu, Seongnam 13620, Korea Tel: 82-31-787-7714 Fax: 82-31-787-4012 E-mail: hye2@snu.ac.kr
*Sumi Yun and Yoonjin Kwak contributed equally to this work.
Received 2017 October 13; Accepted 2018 January 15.

Abstract

Purpose

Molecular treatments targeting epidermal growth factor receptors (EGFRs) are important strategies for advanced colorectal cancer (CRC). However, clinicopathologic implications of EGFRs and EGFR ligand signaling have not been fully evaluated. We evaluated the expression of EGFR ligands and correlation with their receptors, clinicopathologic factors, and patients’ survival with CRC.

Materials and Methods

The expression of EGFR ligands, including heparin binding epidermal growth factor-like growth factor (HBEGF), transforming growth factor (TGF), betacellulin, and epidermal growth factor (EGF), were evaluated in 331 consecutive CRC samples using mRNA in situ hybridization (ISH). We also evaluated the expression status of EGFR, human epidermal growth factor receptor 2 (HER2), HER3, and HER4 using immunohistochemistry and/or silver ISH.

Results

Unlike low incidences of TGF (38.1%), betacellulin (7.9%), and EGF (2.1%), HBEGF expression was noted in 62.2% of CRC samples. However, the expression of each EGFR ligand did not reveal significant correlations with survival. The combined analyses of EGFR ligands and EGFR expression indicated that the ligands‒/EGFR+ group showed a significant association with the worst disease-free survival (DFS; p=0.018) and overall survival (OS; p=0.005). It was also an independent, unfavorable prognostic factor for DFS (p=0.026) and OS (p=0.007). Additionally, HER4 nuclear expression, regardless of ligand expression, was an independent, favorable prognostic factor for DFS (p=0.034) and OS (p=0.049), by multivariate analysis.

Conclusion

Ligand-independent EGFR overexpression was suggested to have a significant prognostic impact; thus, the expression status of EGFR ligands, in addition to EGFR, might be necessary for predicting patients' outcome in CRC.

Introduction

Recently, considerable progress has been made in molecular studies of carcinogenesis and cancer progression. This deeper understanding has led to advances in new treatment strategies, including the use of molecularly targeted agents. One of the therapeutic targets in clinical use is the epidermal growth factor receptor (EGFR) [1-3]. EGFR, a member of the EGFR family, is known to be overexpressed in various cancers including colorectal cancer (CRC) [1,2]. Several EGFR signaling inhibitors, such as cetuximab and panitumumab, against the EGFR’s extracellular domain have an important clinical role in treatment and provide significant survival benefits in advanced CRC patients [4]. Unfortunately, a significant number of CRC patients with activated EGFR signaling suffered resistance to these treatments [4]. Several studies suggested that other genetic alterations influenced the therapeutic resistance in these cases, and more studies focused on the underlying mechanisms of resistance to EGFR inhibitors and predictive biomarkers to discriminate the patients who may receive greater benefits from treatment with EGFR inhibitors [4].

Various ligands are known to activate EGFR signaling, including heparin-binding epidermal growth factor–like growth factor (HBEGF), transforming growth factor (TGF), betacellulin, and epidermal growth factor (EGF) [5]. These ligands bind to the extracellular domain of EGFR and activate downstream signaling pathways [5]. Several studies have investigated the protein and mRNA expressions of some EGFR ligands, such as amphiregulin and epiregulin [6-8]. These ligands are commonly detected in CRC, suggesting their potential prognostic and predictive impact for anti-EGFR inhibition therapy in CRC [6,7,9-11]. Khelwatty et al. [8] also demonstrated that the evaluation of not only EGFR expression, but also the expression of its ligands may be used to predict the effect of anti-EGFR treatments. Furthermore, some EGFR ligands including HBEGF and betacellulin interact with human epidermal growth factor receptor 4 (HER4)— another member of the EGFR family—that may have close associations with favorable features and prognosis in cancer [12-14]. However, the clinicopathologic impact of various EGFR ligands and underlying mechanisms in ligand-dependent EGFR signaling remains unclear.

Therefore, wqe evaluated the mRNA expression of EGFR ligands including HBEGF, TGF, betacellulin, and EGF in a consecutive CRC cohort, and explored the clinicopathologic implications of EGFR ligands in relation to EGFR families.

Materials and Methods

1. Patients characteristics and tissue microarray construction

A total of 331 CRC samples were collected from 2005 to 2006 at Seoul National University Bundang Hospital. All patents had undergone curative surgical resection without preoperative chemoradiation therapy. Patients’ characteristics and clinicopathologic details including age, sex, tumor location, tumor differentiation, pathologic stage (American Joint Committee on Cancer TNM stage, seventh edition), time of death, tumor recurrence, and adjuvant chemoradiation treatment status were obtained from medical and pathologic records. Disease-free survival (DFS) was defined as the time from the date of curative resection to tumor recurrence and overall survival (OS) was calculated to the date of death. Of these 331 CRC cases, patients were treated with chemotherapy alone in 234 cases and both chemotherapy and radiotherapy in 15 cases. Among 188 patients with stage III/IV CRC, 155 and nine patients received treatments with additional chemotherapy alone and chemotherapy plus radiotherapy, respectively. After a microscopic review of each case, 2-mm-sized tissue cores from the most representative tumor areas were extracted to construct tissue microarray (TMA) blocks, as previously described [15].

2. mRNA in situ hybridization

Next, the expression levels of EGFR ligands were measured by mRNA in situ hybridization (ISH) using the RNA scope 2.0 HD Reagent kit (Advanced Cell Diagnostics Inc., Hayward, CA) on TMA blocks, according to the manufacturer’s instructions. After pretreatment with heat and protease, 4-μm tissue sections were incubated and hybridized with the target probes for EGFR ligands. Each hybridized slide was then hybridized with preamplifier, amplifier, and horseradish peroxidase‒labeled probes in order. The colorimetric reaction with 3,3′-diaminobenzidine followed. As a positive control, ubiquitin C staining was applied under otherwise identical conditions, and cases with negative or weak ubiquitin C signals were excluded from further analysis. The slides were scored according to the manufacturer’s scoring guideline: 0, no staining or less than 1 dot/cell; 1, 1-3 dots/cell; 2, 4-10 dots/cell; 3, > 10 dots/cell and dot clusters visible in less than 10% positive cells under ×20 objective; and 4, > 10 dots/cell and dot clusters visible in more than 10% positive cells under ×20 objective.

3. Immunohistochemistry and interpretation

Immunohistochemical staining was performed on the 4-μm-thick sections of TMA blocks with EGFR (pharmDX kit, prediluted, DAKO, Glostrup, Denmark), HER3 (1:20, DAKO), and HER4 (1:8,000, Thermo Scientific, Fremont, CA) antibodies. Staining was carried out according to the manufacturer’s protocols. In brief, the slides from the TMA blocks with 4-μm thickness were deparaffinized and rehydrated. Then, these slides were placed for antigen retrieval and incubated with antibodies [13,16,17]. EGFR-positive expression was defined as membranous staining in > 0% of the tumor cells, according to the DAKO EGFR pharmDX assay interpretation guideline. For evaluation of HER3 and HER4 expression, a previously described 4-tiered scoring system based on the staining intensity was applied [16]. The cytoplasmic and/or membranous expression of HER3 was categorized into four grades, as follows: 0, negative; 1, weak staining; 2, moderate staining; and 3, strong staining. In addition, cytoplasmic and nuclear expression of HER4 was evaluated separately. Positivity for HER3 and HER4 expression was considered as ≥ 10% of tumor cells with moderate to strong immunostaining.

4. Assessment of HER2 status

The evaluation of HER2 status was performed by immunohistochemistry and silver in situ hybridization (SISH), according to the DAKO HercepTest interpretation manual for gastric cancer [16]. The slides were stained with an antibody to HER2/neu (predilution, Ventana Medical Systems, Tucson, AZ) using an automated BenchMark XT stainer. HER2 SISH was also performed with INFORM HER2 DNA and INFORM chromosome 17 (CEP17) probes (Ventana Medical Systems). The entire core was scanned and the signals for 20 non-overlapping tumor cell nuclei were calculated in the hot spot area. HER2 status was considered positive if ≥ 10% of tumor cells revealed 2+ and 3+ immunostaining and a ratio between HER2 and CEP17 of ≥ 2.0 upon SISH.

5. Microsatellite instability determination

The microsatellite instability (MSI) status was determined using five markers of NCI recommended microsatellite loci (BAT-25, BAT-26, D5S346, D17S250, and D2S123) [13]. Briefly, tissue sections obtained from the tumor area and matching normal colonic mucosa were analyzed for MSI status using an automated DNA sequencer (ABI3731 Genetic Analyzer, Applied Biosystems, Foster City, CA). The MSI status was divided into three groups as followed: microsatellite stable (MSS; no unstable markers), microsatellite instabilitylow (MSI-L; one unstable marker), and microsatellite instability-high (MSI-H; two or more unstable markers).

6. Statistical analysis

We implemented the chi-square test and Fisher exact test to compare the differences between the expression of EGFR ligands and EGFR families, and various clinicopathologic characteristics. The Kaplan-Meier method with the log-rank test was applied for survival analysis. To analyze the prognostic significance of the expression of EGFR ligands and EGFR families, we performed univariate and multivariate Cox regression analyses. The SPSS ver. 21 (IBM Corp., Armonk, NY) was used for all statistical analyses, and we considered p < 0.05 as statistically significant.

7. Ethical statement

The study was approved by the Institutional Review Board of The Seoul National University Bundang Hospital (B-1210/174-301) and performed in accordance with the principles of the Declaration of Helsinki. The informed consent was waived.

Results

1. Patient characteristics

A total of 331 CRC cases consisting of 184 men and 147 women were included in this study. The median age of these patients was 66 years (range, 20 to 95 years). Regarding tumor histology, well to moderate differentiation was observed in 300 cases (90.6%), and poorly differentiated and mucinous-type CRCs were observed in 31 cases (9.4%). Of the 320 cases with available results for MSI status, 27 cases (8.4%) were categorized to MSI-H and 293 cases (91.6%) were determined as MSS/MSI-L. S1 Table summarizes the patients’ characteristics.

2. Expression of EGFR ligands and EGFR families in CRC

The expression of EGFR ligands was not detected in the non-neoplastic colorectal mucosa. In CRC, EGFR ligands were expressed in the cytoplasm of tumor cells. The CRC cases were divided into two groups according to the expression of EGFR ligands: positive (1-4) or negative (0). Although betacellulin (7.9%, 26/331) and EGF (2.1%, 7/331) were faintly expressed in a small number of CRC cases, expression of HBEGF and TGF were observed in 62.2% (206/331) and 38.1% (126/331) with various staining intensities, respectively (Table 1, Fig. 1).

Distribution of EGFR ligands expression in CRC

Fig. 1.

Representative images of mRNA expression of epidermal growth factor receptor ligands in colorectal cancer. (A) Heparin binding epidermal growth factor–like growth factor (×400). (B) Transforming growth factor (×400). (C) Betacellulin (×400). (D) Epidermal growth factor (×400).

Immunohistochemical expression for EGFR families were also performed in 331 CRC tissue samples, using TMA blocks. EGFR and HER3 were overexpressed in 129 (39.0%) and 228 (68.9%) out of 331 cases, respectively. Twelve cases (3.6%) showed HER2 positivity. Regarding to the sublocalization, HER4 cytoplasmic and nuclear expression were found in 207 (62.5%) and 36 (10.9%) cases, respectively.

3. Association between the expression of EGFR ligands and EGFR families

We explored the correlation between the expression of EGFR ligands and EGFR families. Our results revealed that HBEGF and TGF expressions were significantly associated with HER3 expression (p=0.003 and p=0.024, respectively). However, we did not observe other significant associations between the expression of EGFR ligands and EGFR families (Table 2).

Association between expression of EGFR ligands and EGFR families

4. Association between EGFR ligands and clinicopathologic variables

HBEGF expression revealed a positive correlation with low grade differentiation (p < 0.001), but not with other various clinicopathologic characteristics (all p > 0.05). TGF expression revealed a tendency to be detected in low-grade CRCs with borderline statistical significance (p=0.062). Although the frequencies were low, EGF positivity was associated with vascular invasion (p=0.025) and tended to be detected in CRC with expanding tumor borders (p=0.066). However, no significant difference between betacellulin expression and various clinicopathologic variables was noted. Regarding the association with MSI status, TGF expression was significantly associated with MSS/MSI-L status (p=0.009) (S1 Table). Upon Kaplan-Meier survival analysis with the log-rank test, no significant relationships with either DFS or OS were revealed according to the expression status of each EGFR ligand in CRC (S2 Fig.).

5. Combined analysis of EGFR ligands and EGFR in CRC

We assessed the correlation with clinicopathologic parameters depending on the combined status of EGFR ligands. A total of 241 cases (72.8%) revealed positivity for at least one EGFR ligand in CRC. EGFR ligand positivity was significantly correlated with low grade differentiation (p < 0.001) and the presence of lymphatic invasion (p=0.030), but not with other clinicopathologic factors. The expression of EGFR ligands was more frequently detected in MSS/MSI-L CRC (p=0.009) (Table 3).

Expression of EGFR ligands, EGFR, HER4, and their correlation with clinicopathologic characteristics

Next, we analyzed the clinicopathologic significance depending on the combined status of EGFR ligands and EGFR expression. Interestingly, the EGFR-positive CRC group without ligand expression revealed a significant association with poorly differentiated histology (p=0.023) (Table 4). We explored survival differences using Kaplan-Meier analysis. Patients with ligand–/EGFR+ revealed the worst DFS compared to other groups, but these differences were not statistically significant (p=0.123). The ligand–/EGFR+ group also revealed the poorest OS (p=0.039). After dichotomization, we compared the prognostic significance in the CRC with ligand–/EGFR+ group compared to the CRC with any other expression group; patients with ligand–/EGFR+ revealed significantly shorter DFS and OS (p=0.018 and p=0.005, respectively) (Fig. 2). Moreover, we performed multivariate analysis to determine the prognostic impact. In the Cox model, dichotomization into EGFR ligands–/EGFR+ and others predicted unfavorable DFS and OS (p=0.026 and p=0.007, respectively). Tumor differentiation, invading borders, neural invasion, and pathologic stage also indicated their prognostic significance for patients’ survival (Table 5).

Combined expression of EGFR ligands and EGFR expression

Fig. 2.

Kaplan-Meier survival curves according to the combined expression status of ligands and epidermal growth factor receptor (EGFR). (A) Disease-free survival. (B) Overall survival.

Univariate and multivariate survival analysis

6. Association of combined status of HBEGF, betacellulin, and HER4 expression with clinicopathologic paracmeters

Since HBEGF and betacellulin have also been reported as ligands for HER4, we explored the clinicopathologic significance of HBEGF, betacellulin, and HER4 expression. HER4 cytoplasmic expression was found in 207 cases (62.5%), while its nuclear expression was observed in 36 cases (10.9%) (S3 Fig.). HER4 cytoplasmic expression was correlated with lower T category (p=0.004). HER4 nuclear immunoreactivity was significantly associated with lower T category (p=0.001) and the absence of lymphatic invasion (p=0.043). Other clinicopathologic parameters were not correlated with HER4 cytoplasmic or nuclear expression (Table 3). Patients with HER4 nuclear expression, regardless of its ligands, HBEGF and betacellulin expression status, revealed superior DFS (p=0.035) and OS (p=0.030), as determined in the survival analysis. However, we failed to determine a significant correlation between HER4 cytoplasmic expression or its co-expression with its ligands with regard to DFS and OS (S4 Fig.). HER4 nuclear expression predicted favorable DFS and OS upon multivariate analysis (p=0.034 and p=0.049, respectively) (Table 5).

Discussion

In this study, we found that EGFR ligands including HBEGF, TGF, betacellulin, and EGF were expressed with various intensities in CRC tissue. While we failed to determine the significant prognostic values of each EGFR ligand, CRC patients with EGFR overexpression and no other ligand expression presented with significantly unfavorable outcomes, compared to others. This result indicated its independent prognostic value for poor survival. Furthermore, HER4 nuclear expression, regardless of its ligands, HBEGF, and betacellulin expressions, is also a better prognostic indicator for CRC. Our results may provide the basis for understanding EGFR signaling pathways and the development of new treatments targeting EGFR signaling.

EGFR, a well-known transmembrane receptor kinase, is composed of three structural and functional domains: the extracellular ligand binding domain, transmembrane region, and cytoplasmic domain containing tyrosine kinase [18]. Previous studies identified two different EGFR downward signaling pathways. Various ligands activate EGFR signaling by binding to the extracellular domain, which results in conformational changes of EGFR and activation of several canonical signals, such as ERK and Akt [5]. Contrary to ligand-dependent signaling, ligand-independent constitutive activation of EGFR induces a distinct signaling pathway through the transcription factor interferon regulatory factor 3 [18,19]. Interestingly, these two signaling pathways are exclusively activated [18,19].

Previous studies have mainly focused on the prognostic and predictive impacts of amphiregulin and epiregulin, which are commonly expressed in CRC and might be associated with favorable prognostic effects in CRC [6,7,9-11]. However, the clinicopathologic implications of other ligands, including HBEGF, TGF, betacellulin, and EGF, are insufficiently investigated. For HBEGF, a pivotal role in cancer progression has been proposed in various cancers, including cervical, breast, and gastric cancers [20-23]. In CRC, Kawasaki et al. [24] suggested that HBEGF has roles in CRC progression through the recruitment of intestinal myofibroblasts. In our cohort, although HBEGF expression is commonly detected in CRC, we were unable to determine a significant correlation with various clinicopathologic factors except for tumor differentiation. Interestingly, stromal HBEGF expression was observed in 21.5% (71/331) of our cohort. Among these, 63 cases revealed HBEGF expression in both stromal and tumor cells (S5 Fig.). Recent studies by Schrevel et al. [20] and Murata et al. [25] described the possible role of HBEGF from fibroblasts and macrophages, which are associated with cancer cell proliferation in cervical cancer. Based on these results of HBEGF expression in the tumor microenvironment, we analyzed the association between stromal HBEGF expression and clinicopathologic parameters, but no significant differences including patients’ survival depending on stromal HBEGF expression was observed in CRC (data not shown).

A recent study by Khelwatty et al. [8] conducted in metastatic CRC using immunohistochemistry, demonstrated that betacellulin was the most commonly expressed ligand and was associated with unfavorable prognosis in metastatic CRC. However, HBEGF expression was not detected in that study. Inconsistent with their results, our study found different expression frequencies, wherein HBEGF and TGFwere relatively commonly expressed, while betacellulin expression presented in only a few subsets of the consecutive CRC cohort, using mRNA ISH. Immunohistochemistry is considered as a reproducible and efficient tool for the characterization of proteins of interest on formalin fixed paraffin embedded slides. However, there is a paucity of reliable commercially available immunohistochemical antibodies against various EGFR ligands. In addition, the cutoff value of immunohistochemical expression for EGFR ligands remains variable and controversial. Yoshida et al. [26] observed the predictive impact of ≥ 30% expression of four EGFR ligands including amphiregulin, HBEGF, TGF-α, and epiregulin. Khelwatty et al. [8] adopted various cutoff points (> 5-50% of tumor cells) for each EGFR ligand to avoid the effect of low specific staining. Considering the above matters, we evaluated the expression of four EGFR ligands, using mRNA ISH, but mRNA levels were not completely overlapping with the levels of their translated protein expression. It was expected that several factors including different cohorts and detection methods may contribute to explain these discordant results. In addition, previous studies have evaluated the mRNA expression of several EGFR ligands including EGF, TGF in colon cancer cell lines and human colon cancer tissue [27,28]. As these studies, EGFmRNA was significantly expressed in colon cancer cell lines and fresh-frozen human colon cancer specimen. Rajagopal et al. [28] also detected the EGFmRNA expression in non-neoplastic colon specimen, by reverse transcription–polymerase chain reaction analysis, in contrast with our results that no EGFmRNA expression was found in non-neoplastic colonic mucosa. In present study, we examined EGR ligand expression in formalin fixed paraffin embedded specimens by mRNA ISH. The difference of samples used for mRNA measurements and detection methods might explain these conflicting results. Further studies regarding the expression of EGFR ligands are warranted to validate these results.

Because the expression of each EGFR ligand was not a prognostic marker in CRC according to our results, we adopted the combined expression status of four ligands to evaluate the clinicopathologic significance of ligand associated EGFR overexpression in CRC. Our results revealed that ligand-independent EGFR overexpression was related to poor prognosis in CRC patients. In addition, different from low-affinity ligands such as amphiregulin and epiregulin, the ligands that we tested, including HBEGF, TGF, betacellulin, and EGF, are known as high-affinity ligands and these ligands have relatively high levels of signaling potential and activating effects for canonical pathways [6]. From these findings, we infer that ligand-dependent EGFR signaling is not a decisive factor for poor prognosis. Several recent studies have described that there were two distinct signaling pathway, ligand dependent and ligand independent constitutive pathway and these two signaling showed the mutually exclusive nature [18,19]. However, the downstream signaling of EGFR activation cooperatively or independently with EGFR ligands was not clear in CRC. In present study, we evaluated the expression status of EGFR and its ligands and showed the prognostic impact depending on their expression phenotypes in CRC, but we could not show the association with EGFR downstream signaling deeply. Furthermore, immunohistochemical staining using EGFR pharmDX antibodies, which were used to determine the expression status of EGFR in this study, have also been identified to be reactive in cases of EGFRvIII mutants, which trigger ligand-independent EGFR signaling [19]. Given these points, our study has certain limitations and further studies will be necessary to elucidate EGFR downstream signaling, considering the association of its ligands in CRC.

For HER4 expression, immunohistochemically, cytoplasmic localization was the predominant pattern and some cases revealed nuclear immunostaining in our study. HER4 activation generates two structurally and functionally different isoforms by alternative splicing: CYT1 and CYT2. While the CYT2 isoform can enter the nucleus easily and regulates transcription related to favorable features including cell differentiation, the CYT1 isoform is insufficient to translocate into the nucleus and contains a phosphatidylinositol-3 kinase binding motif to activate its signaling pathway to induce tumor cell proliferation [13,14]. These findings suggested that, as for staining patterns, nuclear accumulation of HER4 might indicate better prognosis, in contrast to cytoplasmic localization that is associated with tumor progression [12]. We previously reported that HER4 nuclear expression might be regarded as a favorable prognostic marker in gastric cancer [13]. However, clinicopathologic significance depending on the subcellular localization of the HER4 protein was not fully evaluated in CRC. Based on the results of previous studies, cytoplasmic HER4 indicated a positive association with poor survival in CRC [12,29]. In our study, the prognostic impact of HER4 cytoplasmic expression was not significant, which was inconsistent with previous results. Moreover, HER4 nuclear expression was significantly associated with the absence of lymphatic invasion and lower pathologic stage, and was determined to be an independent factor for improved prognosis in CRC. These findings are generally consistent with the results demonstrated in gastric cancer [13]. Interestingly, these prognostic impacts were not influenced by the status of its ligands, or HBEGF and betacellulin expressions. Therefore, our results suggest that HER4 nuclear expression, not cytoplasmic localization, and regardless of ligand expression, might be a biomarker for predicting better prognosis in CRC patients. Because there were more known ligands for HER4 that we did not test, larger and more comprehensive studies for the subcellular localization of HER4 expression are needed to explain the underlying mechanism and the different prognostic impacts in CRC.

In conclusion, we identified that EGFR overexpression without ligand expression could be a biomarker for unfavorable prognosis. Moreover, the evaluation of EGFR ligand expression, in addition to EGFR expression, could provide meaningful information for predicting patients’ survival in CRC. Together, HER4 nuclear expression, regardless of its ligands expression, is suggested as an independent favorable prognostic factor in CRC.

Electronic Supplementary Material

Notes

Conflict of interest relevant to this article was not reported.

Acknowledgements

This study was funded through a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number: HI14C1813).

References

1. Spano JP, Lagorce C, Atlan D, Milano G, Domont J, Benamouzig R, et al. Impact of EGFR expression on colorectal cancer patient prognosis and survival. Ann Oncol 2005;16:102–8.
2. Saijo N. Present status and problems on molecular targeted therapy of cancer. Cancer Res Treat 2012;44:1–10.
3. Kim D, Hong YS, Kim JE, Kim KP, Lee JL, Chun SM, et al. Use of a High-throughput genotyping platform (OncoMap) for RAS mutational analysis to predict cetuximab efficacy in patients with metastatic colorectal cancer. Cancer Res Treat 2017;49:37–43.
4. Karapetis CS, Khambata-Ford S, Jonker DJ, O'Callaghan CJ, Tu D, Tebbutt NC, et al. K-ras mutations and benefit from cetuximab in advanced colorectal cancer. N Engl J Med 2008;359:1757–65.
5. Schneider MR, Wolf E. The epidermal growth factor receptor ligands at a glance. J Cell Physiol 2009;218:460–6.
6. Khambata-Ford S, Garrett CR, Meropol NJ, Basik M, Harbison CT, Wu S, et al. Expression of epiregulin and amphiregulin and K-ras mutation status predict disease control in metastatic colorectal cancer patients treated with cetuximab. J Clin Oncol 2007;25:3230–7.
7. Hobor S, Van Emburgh BO, Crowley E, Misale S, Di Nicolantonio F, Bardelli A. TGFα and amphiregulin paracrine network promotes resistance to EGFR blockade in colorectal cancer cells. Clin Cancer Res 2014;20:6429–38.
8. Khelwatty S, Essapen S, Bagwan I, Green M, Seddon A, Modjtahedi H. The impact of co-expression of wild-type EGFR and its ligands determined by immunohistochemistry for response to treatment with cetuximab in patients with metastatic colorectal cancer. Oncotarget 2017;8:7666–77.
9. Ohchi T, Akagi Y, Kinugasa T, Kakuma T, Kawahara A, Sasatomi T, et al. Amphiregulin is a prognostic factor in colorectal cancer. Anticancer Res 2012;32:2315–21.
10. Jacobs B, De Roock W, Piessevaux H, Van Oirbeek R, Biesmans B, De Schutter J, et al. Amphiregulin and epiregulin mRNA expression in primary tumors predicts outcome in metastatic colorectal cancer treated with cetuximab. J Clin Oncol 2009;27:5068–74.
11. Jonker DJ, Karapetis CS, Harbison C, O'Callaghan CJ, Tu D, Simes RJ, et al. Epiregulin gene expression as a biomarker of benefit from cetuximab in the treatment of advanced colorectal cancer. Br J Cancer 2014;110:648–55.
12. Fujiwara S, Hung M, Yamamoto-Ibusuk CM, Yamamoto Y, Yamamoto S, Tomiguchi M, et al. The localization of HER4 intracellular domain and expression of its alternately-spliced isoforms have prognostic significance in ER+ HER2- breast cancer. Oncotarget 2014;5:3919–30.
13. Yun S, Koh J, Nam SK, Park JO, Lee SM, Lee K, et al. Clinical significance of overexpression of NRG1 and its receptors, HER3 and HER4, in gastric cancer patients. Gastric Cancer 2018;21:225–36.
14. Wang J, Yin J, Yang Q, Ding F, Chen X, Li B, et al. Human epidermal growth factor receptor 4 (HER4) is a favorable prognostic marker of breast cancer: a systematic review and meta-analysis. Oncotarget 2016;7:76693–703.
15. Lee KS, Nam SK, Koh J, Kim DW, Kang SB, Choe G, et al. Stromal expression of microRNA-21 in advanced colorectal cancer patients with distant metastases. J Pathol Transl Med 2016;50:270–7.
16. Seo AN, Kwak Y, Kim WH, Kim DW, Kang SB, Choe G, et al. HER3 protein expression in relation to HER2 positivity in patients with primary colorectal cancer: clinical relevance and prognostic value. Virchows Arch 2015;466:645–54.
17. Atkins D, Reiffen KA, Tegtmeier CL, Winther H, Bonato MS, Storkel S. Immunohistochemical detection of EGFR in paraffin-embedded tumor tissues: variation in staining intensity due to choice of fixative and storage time of tissue sections. J Histochem Cytochem 2004;52:893–901.
18. Chakraborty S, Li L, Puliyappadamba VT, Guo G, Hatanpaa KJ, Mickey B, et al. Constitutive and ligand-induced EGFR signalling triggers distinct and mutually exclusive downstream signalling networks. Nat Commun 2014;5:5811.
19. Guo G, Gong K, Wohlfeld B, Hatanpaa KJ, Zhao D, Habib AA. Ligand-independent EGFR signaling. Cancer Res 2015;75:3436–41.
20. Schrevel M, Osse EM, Prins FA, Trimbos JB, Fleuren GJ, Gorter A, et al. Autocrine expression of the epidermal growth factor receptor ligand heparin-binding EGF-like growth factor in cervical cancer. Int J Oncol 2017;50:1947–54.
21. Iwamoto R, Takagi M, Akatsuka J, Ono K, Kishi Y, Mekada E. Characterization of a novel anti-human HB-EGF monoclonal antibody applicable for paraffin-embedded tissues and diagnosis of HB-EGF-related cancers. Monoclon Antib Immunodiagn Immunother 2016;35:73–82.
22. Zhou ZN, Sharma VP, Beaty BT, Roh-Johnson M, Peterson EA, Van Rooijen N, et al. Autocrine HBEGF expression promotes breast cancer intravasation, metastasis and macrophage-independent invasion in vivo. Oncogene 2014;33:3784–93.
23. Byeon SJ, Lee HS, Kim MA, Lee BL, Kim WH. Expression of the ERBB family of ligands and receptors in gastric cancer. Pathobiology 2017;84:210–7.
24. Kawasaki H, Saotome T, Usui T, Ohama T, Sato K. Regulation of intestinal myofibroblasts by KRas-mutated colorectal cancer cells through heparin-binding epidermal growth factor-like growth factor. Oncol Rep 2017;37:3128–36.
25. Murata T, Mizushima H, Chinen I, Moribe H, Yagi S, Hoffman RM, et al. HB-EGF and PDGF mediate reciprocal interactions of carcinoma cells with cancer-associated fibroblasts to support progression of uterine cervical cancers. Cancer Res 2011;71:6633–42.
26. Yoshida M, Shimura T, Sato M, Ebi M, Nakazawa T, Takeyama H, et al. A novel predictive strategy by immunohistochemical analysis of four EGFR ligands in metastatic colorectal cancer treated with anti-EGFR antibodies. J Cancer Res Clin Oncol 2013;139:367–78.
27. Ito M, Yoshida K, Kyo E, Ayhan A, Nakayama H, Yasui W, et al. Expression of several growth factors and their receptor genes in human colon carcinomas. Virchows Arch B Cell Pathol Incl Mol Pathol 1990;59:173–8.
28. Rajagopal S, Huang S, Moskal TL, Lee BN, el-Naggar AK, Chakrabarty S. Epidermal growth factor expression in human colon and colon carcinomas: anti-sense epidermal growth factor receptor RNA down-regulates the proliferation of human colon cancer cells. Int J Cancer 1995;62:661–7.
29. Mota JM, Collier KA, Barros Costa RL, Taxter T, Kalyan A, Leite CA, et al. A comprehensive review of heregulins, HER3, and HER4 as potential therapeutic targets in cancer. Oncotarget 2017;8:89284–306.

Article information Continued

Fig. 1.

Representative images of mRNA expression of epidermal growth factor receptor ligands in colorectal cancer. (A) Heparin binding epidermal growth factor–like growth factor (×400). (B) Transforming growth factor (×400). (C) Betacellulin (×400). (D) Epidermal growth factor (×400).

Fig. 2.

Kaplan-Meier survival curves according to the combined expression status of ligands and epidermal growth factor receptor (EGFR). (A) Disease-free survival. (B) Overall survival.

Table 1.

Distribution of EGFR ligands expression in CRC

Ligand mRNA expression
0 1 2 3 4
HBEGF 125 (37.8) 96 (29.0) 57 (17.2) 37 (11.2) 16 (4.8)
TGF 205 (61.9) 88 (26.6) 27 (8.2) 10 (3.0) 1 (0.3)
Betacellulin 305 (92.1) 23 (6.9) 3 (0.9) 0 0
EGF 324 (97.9) 1 (0.3) 6 (1.8) 0 0

Values are presented as number (%). HBEGF, heparin binding epidermal growth factor like growth factor; TGF, transforming growth factor; EGF, epidermal growth factor.

Table 2.

Association between expression of EGFR ligands and EGFR families

Ligand EGFR-positive p-value HER2-positive p-value HER3-positive p-value HER4cyto-positive p-value HER4nu-positive p-value
HBEGF
 0 44 (35.2) 0.273 6 (4.8) 0.380 74 (59.2) 0.003 71 (56.8) 0.093 12 (9.6) 0.561
 1-4 85 (41.3) 6 (2.9) 154 (74.8) 136 (66.0) 24 (11.7)
TGF
 0 78 (38.0) 0.660 9 (4.4) 0.546 132 (64.4) 0.024 131 (63.9) 0.513 22 (10.7) 0.914
 1-4 51 (40.5) 3 (2.4) 96 (76.2) 76 (60.3) 14 (11.1)
Betacellulin
 0 121 (39.7) 0.372 12 (3.9) 0.609 209 (68.5) 0.630 189 (62.0) 0.463 33 (10.8) > 0.999
 1-4 8 (30.8) 0 19 (73.1) 18 (69.2) 3 (11.5)
EGF
 0 127 (39.2) 0.710 12 (3.7) > 0.999 223 (68.8) > 0.999 203 (62.7) > 0.999 36 (11.1) > 0.999
 1-4 2 (28.6) 0 5 (71.4) 4 (57.1) 0

Values are presented as number (%). EGFR, epidermal growth factor receptor; HER2, human epidermal growth factor receptor 2; HER3, human epidermal growth factor receptor 3; HER4cyto, human epidermal growth factor receptor 4 cytoplasmic expression; HER4nu, human epidermal growth factor receptor 4 nuclear expression; HBEGF, heparin binding epidermal growth factor like growth factor; TGF, transforming growth factor; EGF, epidermal growth factor.

Table 3.

Expression of EGFR ligands, EGFR, HER4, and their correlation with clinicopathologic characteristics

Characteristic EGFR ligands
EGFR
HER4cyto
HER4nu
Positive p-value Positive p-value Positive p-value Positive p-value
Age (yr)
 ≤ 60 76 (69.7) 0.377 33 (30.3) 0.023 68 (62.4) 0.968 11 (10.1) 0.748
 > 60 165 (74.3) 96 (43.2) 139 (62.6) 25 (11.3)
Sex
 Male 138 (75.0) 0.316 62 (33.7) 0.028 118 (64.1) 0.503 18 (9.8) 0.475
 Female 103 (70.1) 67 (45.6) 89 (60.5) 18 (12.2)
Location
 Right sided 68 (70.1) 0.476 47 (48.5) 0.023 55 (56.7) 0.158 7 (7.2) 0.169
 Left sided 173 (73.9) 82 (35.0) 152 (65.0) 29 (12.4)
Depth of tumor (pT)
 T1, T2 36 (75.0) 0.712 13 (27.1) 0.068 39 (81.3) 0.004 12 (25.0) 0.001
 T3, T4 205 (72.4) 116 (41.0) 168 (59.4) 24 (8.5)
Differentiation
 Low grade 227 (75.7) < 0.001 115 (38.3) 0.458 192 (64.0) 0.087 34 (11.3) 0.553
 High grade 14 (45.2) 14 (45.2) 15 (48.4) 2 (6.5)
Tumor border
 Expanding 31 (64.6) 0.166 13 (27.1) 0.068 29 (60.4) 0.743 5 (10.4) 0.912
 Infiltrative 210 (74.2) 116 (41.0) 178 (62.9) 31 (11.0)
Lymphatic invasion
 Absent 94 (66.7) 0.030 48 (34.0) 0.113 96 (68.1) 0.072 21 (14.9) 0.043
 Present 147 (77.4) 81 (42.6) 111 (58.4) 15 (7.9)
Vascular invasion
 Absent 196 (72.9) 0.964 107 (39.8) 0.532 167 (62.1) 0.721 30 (11.2) 0.737
 Present 45 (72.6) 22 (35.5) 40 (64.5) 6 (9.7)
Neural invasion
 Absent 165 (72.1) 0.643 83 (36.2) 0.127 143 (62.4) 0.959 21 (9.2) 0.135
 Present 76 (74.5) 46 (45.1) 64 (62.7) 15 (14.7)
pTNM stage
 I, II, III 200 (73.8) 0.389 107 (39.5) 0.686 176 (64.9) 0.055 33 (12.2) 0.106
 IV 41 (68.3) 22 (36.7) 31 (51.7) 3 (5.0)
MSI status (n=320)
 MSS/MSI-L 220 (75.1) 0.009 119 (40.6) 0.135 185 (63.1) 0.247 33 (11.3) 1.000
 MSI-H 14 (51.9) 7 (25.9) 14 (51.9) 3 (11.1)

Values are presented as number (%). EGFR, epidermal growth factor receptor; HER4, human epidermal growth factor receptor 4; HER4cyto, human epidermal growth factor receptor 4 cytoplasmic expression; HER4nu, human epidermal growth factor receptor 4 nuclear expression; MSI, microsatellite instability; MSS, microsatellite stable; MSI-L, microsatellite instability–low; MSI-H, high microsatellite instability–high.

Table 4.

Combined expression of EGFR ligands and EGFR expression

Characteristic Others Ligand–/EGFR+ p-value
Age (yr)
 ≤ 60 102 (93.6) 7 (6.4) 0.131
 > 60 196 (88.3) 26 (11.7)
Sex
 Male 164 (89.1) 20 (10.9) 0.541
 Female 134 (91.2) 13 (8.8)
Location
 Right sided 86 (88.7) 11 (11.3) 0.592
 Left sided 212 (90.6) 22 (9.4)
Depth of tumor (pT)
 T1, T2 45 (93.8) 3 (6.3) 0.444
 T3, T4 253 (89.4) 30 (10.6)
Differentiation
 Low grade 274 (91.3) 26 (8.7) 0.023
 High grade 24 (77.4) 7 (22.6)
Tumor border
 Expanding 43 (89.6) 5 (10.4) > 0.999
 Infiltrative 255 (90.1) 28 (9.9)
Lymphatic invasion
 Absent 130 (92.2) 11 (7.8) 0.257
 Present 168 (88.4) 22 (11.6)
Vascular invasion
 Absent 245 (91.1) 24 (8.9) 0.185
 Present 53 (85.5) 9 (14.5)
Neural invasion
 Absent 208 (90.8) 21 (9.2) 0.467
 Present 90 (88.2) 12 (11.8)
pTNM stage
 I, II, III 245 (90.4) 26 (9.6) 0.628
 IV 53 (88.3) 7 (11.7)
MSI status (n=320)
 MSS/MSI-L 264 (90.1) 29 (9.9) 0.742
 MSI-H 24 (88.9) 3 (11.1)

Values are presented as number (%). EGFR, epidermal growth factor receptor; MSI, microsatellite instability; MSS, microsatellite stable; MSI-L, microsatellite instability–low; MSI-H, microsatellite instability–high.

Table 5.

Univariate and multivariate survival analysis

Variable Disease-free survival
Overall survival
HR 95% CI p-value HR 95% CI p-value
Univariate analysis
 Age (> 60 yr vs. ≤ 60 yr) 1.189 0.794-1.778 0.400 1.312 0.834-2.066 0.240
 Sex (male vs. female) 1.435 0.979-2.104 0.064 1.240 0.817-1.881 0.312
 Location (left sided vs. right sided) 1.198 0.789-1.821 0.397 1.168 0.738-1.849 0.507
 Differentiation (high grade vs. low grade) 3.040 1.853-4.988 < 0.001 3.583 2.135-6.014 < 0.001
 Tumor border (infiltrative vs. expanding) 4.148 1.692-10.172 0.002 4.050 1.487-11.031 0.006
 Lymphatic invasion (present vs. absent) 2.770 1.805-4.251 < 0.001 3.256 1.983-5.346 < 0.001
 Venous invasion (present vs. absent) 3.287 2.230-4.843 < 0.001 3.799 2.495-5.783 < 0.001
 Neural invasion (present vs. absent) 3.030 2.090-4.393 < 0.001 3.520 2.331-5.314 < 0.001
 pTNM stage (IV vs. I, II, III) 4.562 3.116-6.678 < 0.001 7.592 5.009-11.506 < 0.001
 MSI status (MSS/MSI-L vs. MSI-H) 1.381 0.642-2.974 0.409 1.051 0.485-2.278 0.899
 Ligand–/EGFR+ vs. others 1.843 1.099-3.088 0.020 2.146 1.251-3.682 0.006
 HER4nu (positive vs. negative) 0.423 0.186-0.964 0.041 0.347 0.128-0.946 0.039
Multivariate analysis
 Differentiation (high grade vs. low grade) 2.964 1.762-4.986 < 0.001 3.487 1.998-6.084 < 0.001
 Tumor border (infiltrative vs. expanding) 2.888 1.161-7.184 0.023 2.446 0.878-6.815 0.087
 Lymphatic invasion (present vs. absent) 1.414 0.876-2.283 0.156 1.369 0.780-2.402 0.274
 Venous invasion (present vs. absent) 1.348 0.844-2.154 0.212 1.197 0.723-1.982 0.484
 Neural invasion (present vs. absent) 2.277 1.529-3.391 < 0.001 2.593 1.658-4.057 < 0.001
 pTNM stage (IV vs. I, II, III) 3.022 2.006-4.552 < 0.001 5.578 3.543-8.780 < 0.001
 Ligand–/EGFR+ vs. others 1.824 1.074-3.098 0.026 2.158 1.230-3.787 0.007
 HER4nu (positive vs. negative) 0.408 0.178-0.935 0.034 0.363 0.132-0.994 0.049

HR, hazard ratio; CI, confidence interval; MSI, microsatellite instability; MSS, microsatellite stable; MSI-L, microsatellite instability–low; MSI-H, microsatellite instability–high; EGFR, epidermal growth factor receptor; HER4nu, human epidermal growth factor receptor 4 nuclear expression.