- Open Access
Evaluation of the prognostic significance of HER family mRNA expression in high-risk early breast cancer: a Hellenic Cooperative Oncology Group (HeCOG) validation study
Journal of Translational Medicine volume 13, Article number: 171 (2015)
The aim of the study was to evaluate the prognostic ability of the transcriptional profiling of the HER family genes in early breast cancer, as a validation analysis of another previously published HeCOG study.
RNA was extracted from 663 formalin-fixed paraffin-embedded (FFPE) tumor tissue samples of high-risk early breast cancer patients enrolled in the randomized HE10/00 trial. Relative mRNA expression of all four HER family members was assessed by quantitative reverse transcription-polymerase chain reaction (qRT-PCR).
In compliance with our previous study, the overall agreement between qRT-PCR and IHC/FISH for HER2 status determination was good (69%). Likewise, the overall concordance between qRT-PCR and IHC for EGFR status was high (81%). In line with our previously reported data, we demonstrated a positive association between HER2 and HER3 mRNA expression. Similarly, mRNA expression of HER3 and HER4 was positively associated with each other and negatively associated with EGFR. Regarding relationships with clinico-pathological parameters, our findings are also in agreement with our previous results. Generally, increased EGFR and HER2 mRNA expression was related to unfavorable, whereas high HER3 and HER4 mRNA expression was associated with favorable clinico-pathological parameters. In univariate analysis, no significant association between EGFR, HER2 and HER3 mRNA expression and overall survival (OS) or disease-free survival (DFS) was demonstrated. However, high EGFR protein expression was associated with significantly shorter OS (log-rank, p = 0.015). In compliance with our previously published data, increased HER4 mRNA expression had a significantly favorable prognostic value in terms of OS (p = 0.044) and DFS (p = 0.047). In multivariate analysis, among all HER receptors, only EGFR protein expression was found to affect OS (Wald’s p = 0.028) and DFS (p = 0.015) independently. Concerning the combined expression of all four HER family receptors, the combination of high EGFR, high HER2, low HER3 and low HER4 mRNA expression was associated with a trend for shorter OS (log-rank, p = 0.065) and significantly worse DFS (p = 0.033), compared with all other co-expression profiles.
These data indicate that qRT-PCR may represent a valid alternative method for evaluating the expression of HER family members in FFPE breast carcinoma tissue samples.
Australian New Zealand Clinical Trials Registry ACTRN12609001036202
The human epidermal growth factor receptor (HER) family consists of four homologous members: ErbB-1 (epidermal growth factor receptor [EGFR] or HER1), ErbB-2 (HER2), ErbB-3 (HER3) and ErbB-4 (HER4) . The HER family constitutes a promising area for the development of targeted treatments in patients with breast cancer and considerable therapeutic progress has already been achieved. Apart from HER2, the significance of lateral signaling partners is also increasingly recognized, given the role of dimerization among HER receptors. Therefore, the evaluation of all HER family members as a whole is considered important.
Overexpression and/or amplification of the HER2 receptor occurs in 15–30% of breast cancer cases and is associated with aggressive course of the disease and unfavorable clinical outcome . HER2 status can be evaluated at the DNA, the mRNA or the protein level. Although various assays are available, immunohistochemistry (IHC) is the technique of choice in the routine practice. Fluorescence in situ hybridization (FISH) analysis is required in cases with (2+) HER2 immunostaining . However, despite efforts to standardize these techniques, there is a considerable intra- and inter-laboratory variability of the results [4, 5]. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) represents an alternative test for the determination of HER2 status. qRT-PCR does not necessitate experienced pathologists to interpret the results and is associated with reproducible and quantitative findings. Previous reports have demonstrated that the qRT-PCR can be applied in archival formalin-fixed paraffin-embedded (FFPE) tissues .
In a previous study , we investigated the potential prognostic value of the transcriptional profiling of all four HER family genes in patients with high-risk early breast cancer, receiving dose-dense anthracycline-based sequential adjuvant chemotherapy with or without paclitaxel. The results of this analysis, using qRT-PCR, suggested that EGFR and HER2 are unfavorable prognostic markers, while HER3 and HER4 mRNA expression is related to better clinical outcome. In order to validate the above findings in the current study, we evaluated the prognostic ability of HER family mRNA expression using qRT-PCR, in a larger series of high-risk patients with early breast cancer. These patients were treated with dose-dense sequential or concurrent epirubicin and paclitaxel, followed by ‘intensified’ CMF, within the context of the Hellenic Cooperative Oncology Group (HeCOG) 10/00 randomized phase III trial .
The HE10/00 trial  was a randomized phase III study (ACTRN12609001036202). Patients were treated with three cycles of epirubicin (E), followed by three cycles of paclitaxel (T, Taxol, Bristol Myers-Squibb, Princeton, NJ), followed by three cycles of intensified CMF (cyclophosphamide, methotrexate and 5-fluorouracil) (E–T–CMF, all cycles given every 2 weeks) or with four cycles of epirubicin/paclitaxel (ET) combination (given on the same day every 3 weeks) followed by three cycles of intensified CMF every 2 weeks (ET–CMF). By study design, the cumulative doses and the chemotherapy duration were identical in the two arms but dose intensity of epirubicin and paclitaxel was double in the E–T–CMF arm. A total of 1,086 eligible patients with node-positive operable breast cancer were accrued in a period of 5 years (2000–2005). HER2-positive patients received trastuzumab upon relapse, as previously described . Treatment schedules, baseline characteristics and clinical outcomes have already been described in detail . Primary tumor diameter, axillary nodal status and tumor grade were obtained from the pathology report. The clinical protocol was approved by local regulatory authorities, while the present translational research studies were approved by the Bioethics Committee of the Aristotle University of Thessaloniki School of Medicine. All patients signed a study-specific written informed consent before randomization, which in addition to giving consent for the trial allowed the use of biological material for future research purposes.
Tissue microarray (TMA) construction
Formalin-fixed paraffin-embedded (FFPE) tumor tissue samples were prospectively collected from 663 patients (61% of 1,086 randomized patients) who were part of the HE10/00 study population. For 39% of the randomized patients, we were either not successful in obtaining an FFPE tumor tissue sample or the sample obtained was not of adequate quality or quantity for biomarker evaluation. The REMARK diagram  for the study is shown in Figure 1. Hematoxylin–eosin stained sections from the tissue blocks were reviewed by two experienced breast cancer pathologists and the most representative tumor areas were marked for the construction of the ΤΜΑ blocks with the use of a manual arrayer (Model I, Beecher Instruments, San Prairie, WI), as previously described [11, 12]. Each case was represented by two tissue cores, 1.5 mm in diameter, obtained from the most representative areas of primary invasive tumors or in some cases (7.2%) from synchronous axillary lymph node metastases and re-embedded in 35 microarray blocks. Each TMA block contained 38–66 tissue cores from the original tumor tissue blocks, while cores from various neoplastic, non-neoplastic and reactive tissues were also included, serving as controls for slide-based assays. Cases not represented, damaged or inadequate on the TMA sections were re-cut from the original blocks and these sections were used for protein and gene analysis. Histological grade was evaluated according to the Scarff, Bloom and Richardson system.
Immunohistochemical labeling was performed according to standard protocols on serial 2.5 μm thick sections from the TMA blocks or the original blocks. All cases were also stained for vimentin (clone V9, Dako, Glostrup, Denmark) and cytokeratin 8/18 (clone 5D3, Novocastra™, Leica Biosystems, Newcastle, UK), which were used as control stains for tissue immunoreactivity and fixation, as well as identification of tumor cells. Tissue samples negative for the above antibodies were excluded from the study. To assure optimal reactivity, immunostaining was done 7–10 days after sectioning at the Laboratory of Molecular Oncology of the Hellenic Foundation for Cancer Research, Aristotle University of Thessaloniki School of Medicine. The staining procedures for EGFR (clone 31G7, Invitrogen, Carlsbad, CA), HER2 (A0485 polyclonal antibody, Dako), estrogen receptor (ER, clone 6F11, Novocastra™, Leica Biosystems), progesterone receptor (PgR, clone 1A6, Novocastra™, Leica Biosystems) and Ki67 (clone MIB-1, Dako) were performed using a Bond Max™ autostainer (Leica Microsystems, Wetzlar, Germany), as previously described .
Interpretation of the IHC results
The evaluation of all IHC sections was done by two experienced breast cancer pathologists, blinded as to the patients’ clinical characteristics and survival data, according to existing established criteria, as previously described . Briefly, EGFR protein expression was scored in a scale from 0 to 3+ and was considered to be positive if membraneous staining of 1+, 2+ or 3+ was present in ≥1% of tumor cells ; HER2 protein expression was scored in a scale from 0 to 3+, the latter corresponding to uniform, intense membrane staining in >30% of invasive tumor cells ; ER and PgR were evaluated using the Histoscore method (max score 400) and were considered positive if staining was present in ≥1% of tumor cell nuclei  and, for Ki67, the expression was defined as low (<14%) or high (≥14%) based on the percentage of stained tumor cell nuclei . The mean percentage of stained cells from the two cores was calculated, while in cases with different intensities, the higher intensity score obtained from the two cores was used. If one of the tissue cores was lost or damaged the overall score was determined from the remaining one. When whole tissue sections were used, the entire tumor area was evaluated.
Fluorescence in situ hybridization (FISH)
TMA sections or whole sections (5 μm thick) were cut for FISH analysis, using the ZytoLight® SPEC HER2/TOP2A/CEP17 triple-color probe kit (ZytoVision, Bremerhaven, Germany). The FISH was performed according to the manufacturer’s protocol with minor modifications. Four carcinoma cell lines (MDA-MB-231, MDA-MB-175, MDA-MB-453 and SK-BR-3) from the Oracle HER2 Control Slide (Leica Biosystems), with a known HER2 gene status, were also used as a control for the FISH assays and analyzed for HER2 genomic status.
For all probes, sequential (5 planes at 1.0 μm) digital images were captured using the Plan Apo VC 100×/1.40 oil objective (Nikon, Kanagawa, Japan), using specific filters for each probe. The resulting images were reconstructed using specifically developed software for cytogenetics (XCyto-Gen, ALPHELYS, Plaisir, France). Processed sections were considered eligible for FISH evaluation according to the ASCO/CAP criteria . For the evaluation of HER2/TOP2A/CEP17 status, non-overlapping nuclei from the invasive part of the tumor were randomly selected, according to morphological criteria using DAPI staining, and scored. The virtual slides of HER2, ER or PgR stains, created as previously described  were used for selecting the invasive part of the tumor in each TMA. Twenty tumor nuclei were counted according to Press et al. . The HER2 gene was considered to be amplified when the HER2/CEP17 ratio was ≥2.2 , or the mean HER2 copy number was >6 .
In cases with ratios at or near the cut-off (1.8–2.2), additional 20 or 40 nuclei were counted and the ratio was recalculated. In cases with a borderline ratio at 60 nuclei, additional FISH assays were performed in whole sections . All primary image data of the TMA and whole tumor sections have been digitally scanned and made publicly available at: http://www.hecog-images.gr/HER2/TOP2A/CEN17/FISH_HE10/97_HE10/00.
RNA isolation from FFPE tissue and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) assessment
Prior to RNA isolation, macrodissection of tumor areas was performed in most of the FFPE sections with <50% tumor cell content. More than one FFPE section (10 μm thick) was used for RNA extraction when the tumor surface of a given sample was less than 0.25 cm2. From each FFPE section or macrodissected tissue fragments, RNA was extracted using a standardized fully automated isolation method for total RNA from FFPE tissue, based on germanium-coated magnetic beads (XTRAKT kit, STRATIFYER Molecular Pathology GmbH, Cologne, Germany) in combination with a liquid handling robot (XTRAKT XL, STRATIFYER Molecular Pathology GmbH), as previously described . The method involves extraction-integrated deparaffinization and DNase I digestion steps. The quality and quantity of RNA was checked by measuring CALM2 expression as a surrogate for amplifiable mRNA by qRT-PCR. CALM2 was used as endogenous reference, since it had previously been identified as being highly expressed among breast cancer tissue samples.
qRT-PCR primers and labeled hydrolysis probes were selected using Primer Express® Software, Version 2.2 and 3 (Applied Biosystems/Life Technologies, Karlsruhe, Germany), according to the manufacturer’s instructions, and were controlled for single nucleotide polymorphisms. All primers, probes and amplicons were checked for their specificity against nucleotide databases at NCBI using basic local alignment search tool (BLAST). Primers and probes were purchased from Eurogentec S.A. (Seraing, Belgium). For each primer/probe set, the amplification efficiency was tested, aiming to reach comparable efficiency of >90% (efficiency range from 91 to 100%). Primers and hydrolysis probes were diluted to 100 µM, using a stock solution with nuclease-free water (Life Technologies GmbH, Darmstadt, Germany). qRT-PCR was applied for the relative quantification (RQ) of the HER family genes EGFR, HER2, HER3 and HER4. The Primer/probe (FAM/TAMRA-labeled) sets used for amplification of the target and reference genes are shown in Table 1.
For PCR, 0.5 µM of each primer and 0.25 µM of each probe were used. All quantitative reverse-transcription PCRs were performed in triplicates using the SuperScript® III Platinum® One-Step qRT-PCR kit (Invitrogen/Life Technologies, Darmstadt, Germany) according to the manufacturer’s instructions. Experiments were performed on a Stratagene Mx3005p (Agilent Technologies, Waldbronn, Germany) with 30 min at 50°C and 2 min at 95°C followed by 40 cycles of 15 s at 95°C and 30 s at 60°C. The lengths of the amplicons detected by the EGFR, HER2, HER3, HER4 and CALM2 assays were 93, 61, 81, 75 and 72 bp, respectively, with PCR efficiencies [E = 1(10 − slope)] of 89.1, 97.2, 96.0, 92.0 and 99.7%, respectively. Samples were considered eligible for further investigation when the cycle threshold (CT) values of the housekeeping gene were <32 (triplicate mean values). Relative expression levels (relative quantification, RQ) of the target transcripts were calculated as 40—DCT values (DCT = mean CT target gene − mean CT housekeeping gene) to yield positively correlated numbers and to facilitate comparisons. A commercially available human reference RNA (Stratagene qPCR Human Reference Total RNA, Agilent Technologies, Waldbronn, Germany) was used as positive control. No-template controls were assessed in parallel to exclude contamination.
Continuous data are presented as medians and corresponding ranges, while categorical data are presented as counts and percentages. Separation of continuous markers into high and low expression was performed using predetermined cut-offs based on a previously published study . These cut-offs were the 75th percentile for EGFR and the 50th percentile, or the median as it is usually referred to, for HER2, HER3 and HER4. The Spearman’s correlation coefficient method was used to assess correlations among continuous variables.
Comparisons of categorical with continuous variables were made using the Mann–Whitney or the Kruskal–Wallis tests, while the Chi square test was used for testing associations between categorical variables. Cohen’s Kappa and calculation of sensitivity and specificity were used as assessing tools in order to determine the concordance between the qRT-PCR and IHC/FISH methods.
Disease-free survival (DFS) was measured from the date of diagnosis until verified disease progression, death or last contact (whichever occurred first), while overall survival (OS) was measured from the date of diagnosis until death from any cause or date of last contact. Time-to-event distributions were estimated using Kaplan–Meier curves, while log-rank tests and univariate Cox analyses were used for assessing differences statistically. Possible prognostic significance of the markers among breast tumor subtypes was assessed using univariate Cox analyses with interactions. Univariate tests were determined at the level of 0.1%, controlling for multiple comparisons and keeping the overall type I error rate at the level of 5%. Due to the exploratory nature of this study, p values less than 0.05 were presented in the results section.
In the multivariate Cox regression analysis, significance was determined at the level of 15%, while variable selection was performed based on the likelihood ratio test, among the following factors: number of positive nodes, tumor size, type of surgery, adjuvant radiotherapy, adjuvant hormonal therapy and subtype classification. All tests were two-sided. The statistical analysis complied with the reporting recommendations for tumor marker prognostic studies (REMARK)  and was performed using the SAS software (SAS for Windows, version 9.3, SAS Institute Inc., Cary, NC).
Basic patient and tumor characteristics according to treatment arm and in the entire study population are shown in Table 2. The differences among these characteristics between the cohort of 663 patients included in the analysis and the remaining 423 patients of the HE10/00 study that were not included are shown in Table 3. The cohort of patients used in the current analysis included cases with more aggressive characteristics, such as higher histological grade, higher number of positive nodes and tumor size, with higher histological grade being found to be associated with HER2 status (66% in HER2-positive vs. 46% in HER2-negative patients, p < 0.001).
Relative mRNA expression values of HER family receptors
The distribution of relative expression values (2 to the power of DCT) of mRNA encoding for HER family members is shown in Figure 2. The median value for EGFR was 421.7 (range 1–38,431), for HER2 1.3 (range 0.02–2,256), for HER3 1.2 (range 0.01–73.1) and for HER4 10.7 (range 0.12–53,231).
Concordance between qRT-PCR and IHC
The total number of samples with data available from both IHC/FISH and qRT-PCR was 625 and 632 for EGFR and HER2, respectively. For EGFR, 93 of the 625 tumors (14.9%) were IHC positive, whereas 153 tumors (24.5%) had EGFR mRNA expression at or above the 75th percentile, as assessed by qRT-PCR. For HER2, 143 of the 632 tumors (22.6%) were IHC/FISH positive, whereas 316 tumors (50.0%) had HER2 mRNA expression at or above the median, as assessed by qRT-PCR. For these tumors, we found a statistically significant association between the evaluations obtained by the two methods, for the EGFR (Chi square test, p < 0.001) and the HER2 (Chi square test, p < 0.001) receptors. The observed overall concordance between the determination of HER2 by qRT-PCR and IHC/FISH was 69.1%. The levels for sensitivity and specificity were 92.3 and 62.4%, respectively, while Cohen’s kappa was 0.38 (95% CI 0.32–0.44). The overall agreement between qRT-PCR and IHC for EGFR was 80.8%. Sensitivity and specificity were 67.7 and 83.1%, respectively, while Cohen’s kappa was 0.40 (95% CI 0.32–0.49) (Table 4).
Relationships among HER family receptor mRNA expression
A positive correlation was demonstrated between HER2 and HER3 mRNA expression levels (r = 0.42, p < 0.0001), as well as between HER2 and HER4 mRNA values (r = 0.11, p = 0.007). Furthermore, HER3 and HER4 mRNA expression levels were positively correlated to each other (r = 0.56, p < 0.0001) and negatively so with EGFR (r = −0.16, p < 0.0001 and r = −0.15, p < 0.0001, respectively). Scattered plots of the mRNA expression values of different HER family members are presented in Figure 3.
Association of HER family receptor mRNA expression with clinicopathological parameters
EGFR mRNA expression was inversely associated with the presence of ER (Mann–Whitney, p < 0.001) and PgR (p = 0.005). Moreover, EGFR mRNA expression was inversely related to age (p = 0.003), post-menopausal status (p = 0.003) and tumor size (p = 0.001). A negative correlation between HER2 mRNA expression and PgR was found (p = 0.003). HER3 mRNA expression was associated with ER (p < 0.001) and PgR positivity (p < 0.001), as well as with histological grade I/II (p = 0.001). In addition, HER3 mRNA expression was inversely related to the number of positive lymph nodes (p = 0.017). HER4 mRNA expression was correlated with ER and PgR positivity (p < 0.001 and p < 0.001, respectively), as well as with histological grade I/II (Kruskal–Wallis, p < 0.001). Furthermore, HER4 was associated with lobular and mixed histology (p = 0.002).
Regarding the associations of HER family receptor mRNA expression with breast cancer subtypes, high EGFR mRNA expression was associated with triple-negative breast cancer (TNBC) (p < 0.001), while high HER2 mRNA expression with luminal-HER2 and HER2-enriched subtypes (both, p < 0.001). High HER3 mRNA expression was associated with luminal A and luminal B subtypes (both, p < 0.001), whereas high HER4 mRNA expression with luminal A tumors (p < 0.001).
Prognostic value of HER family receptor mRNA expression
Survival status of the patients was updated in March 2012. The median follow-up time was 98.9 months (range 0.1–132.5 months). During this time, 218 patients (32.9%) had developed a relapse and 151 patients (22.8%) had died. The 8-year OS was 78.2%, whereas the 8-year DFS was 68.5%. The median OS and DFS have not been reached yet.
Concerning patients with high EGFR mRNA expression, a trend for significantly reduced OS (HR = 1.38, 95% CI 0.98–1.96, Wald’s p = 0.068, log-rank p = 0.074) was observed. However, with respect to protein expression, patients whose tumors were positive for EGFR had significantly reduced OS (HR = 1.65, 95% CI 1.10–2.48, Wald’s p = 0.017, log-rank p = 0.015) (Figure 4). No significant associations between HER2 (either mRNA expression or HER2 status by IHC/FISH) and OS were found. Similarly, HER3 mRNA expression was not associated with OS. In contrast, high HER4 mRNA expression (using the median value as a cut-off point) had a favorable prognostic value in terms of OS (HR = 0.72, 95% CI 0.52–0.99, Wald’s p = 0.045, log-rank p = 0.044) (Figure 5). In the multivariate analysis that included 612 patients, EGFR protein expression (HR = 1.66, 95% CI 1.06–2.60, Wald’s p = 0.028), tumor size (HR = 1.60, 95% CI 1.07–2.40, p = 0.023) and the number of positive nodes (HR = 2.74, 95% CI 1.86–4.04, p < 0.001) retained their prognostic significance for OS (Table 5).
With respect to DFS, no significant associations were demonstrated for EGFR, HER2 or HER3 mRNA expression. On the other hand, high HER4 mRNA expression was associated with lower risk for relapse (HR = 0.76, 95% CI 0.59–1.00, Wald’s p = 0.048, log-rank p = 0.047) (Figure 5). However, in multivariate analysis only EGFR protein expression (HR = 1.56, 95% CI 1.09–2.24, Wald’s p = 0.015), type of surgery (HR = 1.49, 95% CI 1.07–2.06, p = 0.017) and the number of positive nodes (HR = 2.13, 95% CI 1.58–2.87, p < 0.001) independently affected DFS (Table 5).
Prognostic value of HER family receptor mRNA expression among subtypes
The prognostic significance of the HER family receptor mRNA expression among breast cancer subtypes was assessed by the use of univariate interaction tests. High EGFR mRNA expression was found to have unfavorable prognostic value in HER2-enriched cases in terms of OS (HR = 5.30, 95% CI 1.71–16.45, interaction p = 0.031). Adjusting for basic clinical and treatment characteristics, the prognostic ability of EGFR mRNA expression remained significant for OS, not only for the HER2-enriched tumors (HR = 9.12, 95% CI 2.87–29.03, Wald’s p = 0.0002), but also for the triple-negative patients (HR = 2.32, 95% CI 1.02–5.27, p = 0.045), with an overall interaction p value of 0.003 (data not shown).
Prognostic value of HER family member co-expression
Regarding the prognostic significance of specific co-expression patterns of all four HER family receptors, we found that the combination of high EGFR, high HER2, low HER3, and low HER4 mRNA expression was associated with significantly worse DFS (log-rank, p = 0.033) and a trend for decreased OS (p = 0.065), compared with all other possible co-expression profiles.
In the current study, we used qRT-PCR to evaluate the transcriptional profiling of all four HER family receptor genes, in a large series of high-risk (predominantly T2-3, node-positive) early breast cancer patients, with a substantially long follow-up of 8 years. The aim of this analysis was to validate the findings of a previously published similar study, conducted from our group . Most of the available studies have assessed protein expression and/or gene amplification of specific HER family members. Therefore, the prognosis of patients with breast cancer relating to the transcriptional profiling of all HER family receptors has not been extensively evaluated.
In accord with our previous results, the overall agreement between qRT-PCR (using the 50th percentile cut-off) and IHC/FISH for the determination of HER2 status was good (69%), with a Cohen’s kappa value of only 0.38, due to inherent differences in the cut-offs used in the two methods. These findings are in line with other studies [6, 22, 23]. Similarly, in agreement with our previous study the overall concordance between qRT-PCR (using the 75th percentile cut-off) and IHC for the evaluation of EGFR status was high (81%). Moreover, in complete agreement with our previous results, the qRT-PCR technique for HER2 (using the median value as a threshold) was associated with a high level of sensitivity (92%) and acceptable specificity (62%). The limitations regarding qRT-PCR as an alternative method for assessing HER family members in FFPE samples has previously been described .
Concerning the associations with clinicopathological parameters, our results are in compliance with our previous findings, as well. In general, HER3 and HER4 mRNA expression was associated with favorable, while EGFR and HER2 mRNA expression was associated with unfavorable parameters.
In the prognostic analyses, no significant associations between EGFR mRNA expression and OS were found in the total population, irrespectively of the cut-off points that we used. In our previous study , we demonstrated a reduced OS in patients with increased EGFR mRNA expression, using the 75th percentile as a cut-off. In the present analysis using the same cut-off, a trend (p = 0.074) towards a decreased OS was shown. Evaluating treatment arms separately, the difference in OS using the 75th percentile was significant in patients receiving the E–T–CMF regimen (p = 0.044). Moreover, EGFR mRNA expression was found to be a significant negative prognostic factor for OS in HER2-enriched and TNBC cases.
With regard to protein expression, patients with tumors positive for EGFR had significantly shorter OS in the total study population. Moreover, the negative prognostic significance of EGFR protein expression for OS and DFS was demonstrated in the multivariate analysis of the study. EGFR is considered to be a negative prognostic factor in breast cancer and such an association has been shown in ours, as well as in other studies [24–27].
In contrast to our previously reported data  and other studies [23, 28] no significant associations between HER2 status (either by mRNA expression or IHC/FISH) and OS or DFS were demonstrated. A definitive explanation for this unexpected finding cannot be given. Although adjuvant treatment with trastuzumab could be proposed to be a possible reason, none of the patients enrolled in the HE10/00 phase III trial received adjuvant therapy with trastuzumab. On the other hand, all patients had been treated with paclitaxel. So far, data concerning the interaction between HER2 receptor status and sensitivity to taxanes are not conclusive. However, the CALGB 9344 adjuvant study  has shown that patients with HER2-positive tumors had significant benefit from the addition of paclitaxel to doxorubicin/cyclophosphamide irrespectively of ER status, while there was no benefit in HER2-negative, ER-positive patients . Moreover, a meta-analysis  of three trials [30, 32, 33] reported a significant interaction between HER2 status and adjuvant taxane therapy in terms of DFS. Such an interaction could possibly explain (at least in part) the absence of negative prognostic significance of HER2 status in the current analysis.
In our previous study  a positive association between HER3 mRNA expression and OS was reported. In contrast, no significant association between HER3 mRNA expression and clinical outcome was found in the present analysis. However, HER3 mRNA expression was associated with several favorable clinicopathological parameters such as ER and PgR positivity, histological grade I/II and fewer positive nodes. The available data in the literature regarding the prognostic role of HER3 in patients with breast cancer are contradictory . Although a number of studies support a negative prognostic value of HER3 [25, 35, 36], other reports suggest a positive prognostic value [24, 37, 38]. So far, a conclusive explanation for the discrepancies among studies concerning the prognostic value of HER3 expression in breast cancer cannot be given. Recent data indicate that the sub-cellular distribution of HER receptors considerably affects their biological activities . Therefore, the effect of HER3 on clinical outcome may be better evaluated taking into consideration, not only the expression of HER3, but also the sub-cellular distribution of the receptor and the expression levels of HER3 ligands .
In complete agreement with our previously reported results , an association of high HER4 mRNA expression (using the median value as a cut-off point) with increased OS and DFS was found. However, the positive prognostic value of high HER4 mRNA expression was not maintained in the multivariate analysis. Similarly, other studies have also demonstrated the positive prognostic ability of HER4 in patients with breast cancer, both at the mRNA and the protein level [24, 25, 41]. Existing evidence suggests that HER4 signaling promotes the differentiation and inhibition of growth in breast cancer cells . In cell line experiments, when HER2-positive cancer cells were transfected to overexpress HER4, a reduction in proliferation and an increase in apoptosis was noted .
With regard to the prognostic potential of the combined expression patterns of the different HER family members, our findings indicate that it may be the co-expression patterns, rather than the expression of individual HER family receptors, that should be taken into account when assessing the prognosis of patients with breast cancer. This observation may prove to have important therapeutic implications.
The present study suggests that EGFR protein overexpression, as assessed by IHC, is a negative prognostic factor for OS and DFS in patients with high-risk operable breast cancer. In addition, high HER4 mRNA expression was associated with a better clinical outcome in terms of OS and DFS in the univariate analysis. Since the clinical study was conducted in the pre-trastuzumab era, the above findings appear to have purely prognostic and not predictive significance. Furthermore, qRT-PCR may be a valid alternative technique for the determination of co-expression patterns of HER family receptors in FFPE breast tumor samples with possible important therapeutic implications.
epidermal growth factor receptor
fluorescence in situ hybridization
Hellenic Cooperative Oncology Group
human epidermal growth factor receptor
polymerase chain reaction
quantitative reverse transcription-polymerase chain reaction
Mosesson Y, Yarden Y (2004) Oncogenic growth factor receptors: implications for signal transduction therapy. Semin Cancer Biol 14:262–270
Slamon DJ, Clark GM, Wong SG, Levin WJ, Ullrich A, McGuire WL (1987) Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science 235:177–182
Bartlett J, Mallon E, Cooke T (2003) The clinical evaluation of HER-2 status: which test to use? J Pathol 199:411–417
Gancberg D, Jarvinen T, di Leo A, Rouas G, Cardoso F, Paesmans M et al (2002) Evaluation of HER-2/NEU protein expression in breast cancer by immunohistochemistry: an interlaboratory study assessing the reproducibility of HER-2/NEU testing. Breast Cancer Res Treat 74:113–120
Thomson TA, Hayes MM, Spinelli JJ, Hilland E, Sawrenko C, Phillips D et al (2001) HER-2/neu in breast cancer: interobserver variability and performance of immunohistochemistry with four antibodies compared with fluorescent in situ hybridization. Mod Pathol 14:1079–1086
Gjerdrum LM, Sorensen BS, Kjeldsen E, Sorensen FB, Nexo E, Hamilton-Dutoit S (2004) Real-time quantitative PCR of microdissected paraffin-embedded breast carcinoma: an alternative method for HER-2/neu analysis. J Mol Diagn 6:42–51
Koutras AK, Kalogeras KT, Dimopoulos MA, Wirtz RM, Dafni U, Briasoulis E et al (2008) Evaluation of the prognostic and predictive value of HER family mRNA expression in high-risk early breast cancer: a Hellenic Cooperative Oncology Group (HeCOG) study. Br J Cancer 99:1775–1785
Gogas H, Dafni U, Karina M, Papadimitriou C, Batistatou A, Bobos M et al (2012) Postoperative dose-dense sequential versus concomitant administration of epirubicin and paclitaxel in patients with node-positive breast cancer: 5-year results of the Hellenic Cooperative Oncology Group HE 10/00 phase III Trial. Breast Cancer Res Treat 132:609–619
Razis E, Bobos M, Kotoula V, Eleftheraki AG, Kalofonos HP, Pavlakis K et al (2011) Evaluation of the association of PIK3CA mutations and PTEN loss with efficacy of trastuzumab therapy in metastatic breast cancer. Breast Cancer Res Treat 128:447–456
McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, Clark GM (2006) REporting recommendations for tumor MARKer prognostic studies (REMARK). Breast Cancer Res Treat 100:229–235
Kononen J, Bubendorf L, Kallioniemi A, Barlund M, Schraml P, Leighton S et al (1998) Tissue microarrays for high-throughput molecular profiling of tumor specimens. Nat Med 4:844–847
Skacel M, Skilton B, Pettay JD, Tubbs RR (2002) Tissue microarrays: a powerful tool for high-throughput analysis of clinical specimens: a review of the method with validation data. Appl Immunohistochem Mol Morphol 10:1–6
Fountzilas G, Kourea HP, Bobos M, Televantou D, Kotoula V, Papadimitriou C et al (2011) Paclitaxel and bevacizumab as first line combined treatment in patients with metastatic breast cancer: the Hellenic Cooperative Oncology Group experience with biological marker evaluation. Anticancer Res 31:3007–3018
Schippinger W, Dandachi N, Regitnig P, Hofmann G, Balic M, Neumann R et al (2007) The predictive value of EGFR and HER-2/neu in tumor tissue and serum for response to anthracycline-based neoadjuvant chemotherapy of breast cancer. Am J Clin Pathol 128:630–637
Wolff AC, Hammond ME, Schwartz JN, Hagerty KL, Allred DC, Cote RJ et al (2007) American Society of Clinical Oncology/College of American Pathologists guideline recommendations for human epidermal growth factor receptor 2 testing in breast cancer. J Clin Oncol 25:118–145
Hammond ME, Hayes DF, Dowsett M, Allred DC, Hagerty KL, Badve S et al (2010) American Society of Clinical Oncology/College Of American Pathologists guideline recommendations for immunohistochemical testing of estrogen and progesterone receptors in breast cancer. J Clin Oncol 28:2784–2795
Cheang MC, Chia SK, Voduc D, Gao D, Leung S, Snider J et al (2009) Ki67 index, HER2 status, and prognosis of patients with luminal B breast cancer. J Natl Cancer Inst 101:736–750
Press MF, Sauter G, Buyse M, Bernstein L, Guzman R, Santiago A et al (2011) Alteration of topoisomerase II-alpha gene in human breast cancer: association with responsiveness to anthracycline-based chemotherapy. J Clin Oncol 29:859–867
Vanden Bempt I, Van Loo P, Drijkoningen M, Neven P, Smeets A, Christiaens MR (2008) Polysomy 17 in breast cancer: clinicopathologic significance and impact on HER-2 testing. J Clin Oncol 26:4869–4874
Sauter G, Lee J, Bartlett JM, Slamon DJ, Press MF (2009) Guidelines for human epidermal growth factor receptor 2 testing: biologic and methodologic considerations. J Clin Oncol 27:1323–1333
Pentheroudakis G, Kotoula V, Eleftheraki AG, Tsolaki E, Wirtz RM, Kalogeras KT et al (2013) Prognostic significance of ESR1 gene amplification, mRNA/protein expression and functional profiles in high-risk early breast cancer: a translational study of the Hellenic Cooperative Oncology Group (HeCOG). PLoS One 8:e70634
Ginestier C, Charafe-Jauffret E, Penault-Llorca F, Geneix J, Adelaide J, Chaffanet M et al (2004) Comparative multi-methodological measurement of ERBB2 status in breast cancer. J Pathol 202:286–298
Vinatzer U, Dampier B, Streubel B, Pacher M, Seewald MJ, Stratowa C et al (2005) Expression of HER2 and the coamplified genes GRB7 and MLN64 in human breast cancer: quantitative real-time reverse transcription-PCR as a diagnostic alternative to immunohistochemistry and fluorescence in situ hybridization. Clin Cancer Res 11:8348–8357
Pawlowski V, Revillion F, Hebbar M, Hornez L, Peyrat JP (2000) Prognostic value of the type I growth factor receptors in a large series of human primary breast cancers quantified with a real-time reverse transcription-polymerase chain reaction assay. Clin Cancer Res 6:4217–4225
Witton CJ, Reeves JR, Going JJ, Cooke TG, Bartlett JM (2003) Expression of the HER1-4 family of receptor tyrosine kinases in breast cancer. J Pathol 200:290–297
Tsutsui S, Ohno S, Murakami S, Hachitanda Y, Oda S (2002) Prognostic value of epidermal growth factor receptor (EGFR) and its relationship to the estrogen receptor status in 1,029 patients with breast cancer. Breast Cancer Res Treat 71:67–75
Tzaida O, Gogas H, Dafni U, Kyroudi A, Papaspyrou I, Kyriakou V et al (2007) Evaluation of the prognostic and predictive value of HER-1/EGFR in breast cancer patients participating in a randomized study with dose-dense sequential adjuvant chemotherapy. Oncology 72:388–396
Bergqvist J, Ohd JF, Smeds J, Klaar S, Isola J, Nordgren H et al (2007) Quantitative real-time PCR analysis and microarray-based RNA expression of HER2 in relation to outcome. Ann Oncol 18:845–850
Henderson IC, Berry DA, Demetri GD, Cirrincione CT, Goldstein LJ, Martino S et al (2003) Improved outcomes from adding sequential paclitaxel but not from escalating doxorubicin dose in an adjuvant chemotherapy regimen for patients with node-positive primary breast cancer. J Clin Oncol 21:976–983
Hayes DF, Thor AD, Dressler LG, Weaver D, Edgerton S, Cowan D et al (2007) HER2 and response to paclitaxel in node-positive breast cancer. N Engl J Med 357:1496–1506
Dhesy-Thind B, Pritchard KI, Messersmith H, O’Malley F, Elavathil L, Trudeau M (2008) HER2/neu in systemic therapy for women with breast cancer: a systematic review. Breast Cancer Res Treat 109:209–229
Martin M, Pienkowski T, Mackey J, Pawlicki M, Guastalla JP, Weaver C et al (2005) Adjuvant docetaxel for node-positive breast cancer. N Engl J Med 352:2302–2313
Kostopoulos I, Arapantoni-Dadioti P, Gogas H, Papadopoulos S, Malamou-Mitsi V, Scopa CD et al (2006) Evaluation of the prognostic value of HER-2 and VEGF in breast cancer patients participating in a randomized study with dose-dense sequential adjuvant chemotherapy. Breast Cancer Res Treat 96:251–261
Amin DN, Campbell MR, Moasser MM (2010) The role of HER3, the unpretentious member of the HER family, in cancer biology and cancer therapeutics. Semin Cell Dev Biol 21:944–950
Bieche I, Onody P, Tozlu S, Driouch K, Vidaud M, Lidereau R (2003) Prognostic value of ERBB family mRNA expression in breast carcinomas. Int J Cancer 106:758–765
Sassen A, Rochon J, Wild P, Hartmann A, Hofstaedter F, Schwarz S et al (2008) Cytogenetic analysis of HER1/EGFR, HER2, HER3 and HER4 in 278 breast cancer patients. Breast Cancer Res 10:R2
Quinn CM, Ostrowski JL, Lane SA, Loney DP, Teasdale J, Benson FA (1994) c-erbB-3 protein expression in human breast cancer: comparison with other tumour variables and survival. Histopathology 25:247–252
Lee Y, Cho S, Seo JH, Shin BK, Kim HK, Kim I et al (2007) Correlated expression of erbB-3 with hormone receptor expression and favorable clinical outcome in invasive ductal carcinomas of the breast. Am J Clin Pathol 128:1041–1049
Offterdinger M, Schofer C, Weipoltshammer K, Grunt TW (2002) c-erbB-3: a nuclear protein in mammary epithelial cells. J Cell Biol 157:929–939
Koutras AK, Fountzilas G, Kalogeras KT, Starakis I, Iconomou G, Kalofonos HP (2010) The upgraded role of HER3 and HER4 receptors in breast cancer. Crit Rev Oncol Hematol 74:73–78
Suo Z, Risberg B, Kalsson MG, Willman K, Tierens A, Skovlund E et al (2002) EGFR family expression in breast carcinomas. c-erbB-2 and c-erbB-4 receptors have different effects on survival. J Pathol 196:17–25
Naresh A, Long W, Vidal GA, Wimley WC, Marrero L, Sartor CI et al (2006) The ERBB4/HER4 intracellular domain 4ICD is a BH3-only protein promoting apoptosis of breast cancer cells. Cancer Res 66:6412–6420
Sartor CI, Zhou H, Kozlowska E, Guttridge K, Kawata E, Caskey L et al (2001) Her4 mediates ligand-dependent antiproliferative and differentiation responses in human breast cancer cells. Mol Cell Biol 21:4265–4275
AK conceived of the study, participated in its design as well as in the acquisition, analysis and interpretation of the data and drafted the manuscript. KTK conceived of the study, participated in its design as well as in the analysis and interpretation of data and drafted the manuscript. RMW conceived of the study, participated in the acquisition, analysis and interpretation of the data and revised critically the manuscript. ZA performed the statistical analysis. MB participated in the interpretation of the data, drafted the manuscript and revised it critically. FZ participated in the acquisition of the data and revised critically the manuscript. EV participated in the interpretation of the data, drafted the manuscript and revised it critically. ET, HG, GP, NP, CM, CC, DB, PP, GA, DP, HPK participated in the acquisition of data and revised critically the manuscript. GF conceived of the study, participated in its design as well as in the analysis and interpretation of the data, drafted the manuscript and revised it critically. All authors read and approved the final manuscript.
The authors are indebted to all patients and their families for their trust and participation in the HE10/00 trial and for the provision of biological material for research purposes. The authors also wish to thank all HeCOG personnel (data managers, research assistants and monitors) for their dedication, M. Moschoni for data coordination and T. Spinari for collection of FFPE tissue blocks. Translational research was supported by a HeCOG research grant: HE TRANS_BR.
Compliance with ethical guidelines
Competing interests The authors declare they have no competing interests.
Angelos Koutras, Konstantine T Kalogeras and Ralph M Wirtz have contributed equally to this work
About this article
Cite this article
Koutras, A., Kalogeras, K.T., Wirtz, R.M. et al. Evaluation of the prognostic significance of HER family mRNA expression in high-risk early breast cancer: a Hellenic Cooperative Oncology Group (HeCOG) validation study. J Transl Med 13, 171 (2015). https://doi.org/10.1186/s12967-015-0530-0
- HER family
- Prognostic value
- Breast cancer