Expression Pattern and Prognostic Value of EPHA/EFNA in Breast Cancer by Bioinformatics Analysis: Revealing Its Importance in Chemotherapy

Liang, Zheng; Wang, Xu; Dong, Kaiti; Li, Xinhua; Qin, Chenge; Zhou, Huifang

doi:https://doi.org/10.1155/2021/5575704

BioMed Research International

On this page

Abstract Introduction Methods and Materials Results Discussion Abbreviations Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Chemo-radiotherapy for Multimodality Treatment of Tumours

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 5575704 | https://doi.org/10.1155/2021/5575704

Expression Pattern and Prognostic Value of EPHA/EFNA in Breast Cancer by Bioinformatics Analysis: Revealing Its Importance in Chemotherapy

Zheng Liang,¹Xu Wang,²Kaiti Dong,¹Xinhua Li,¹Chenge Qin,¹and Huifang Zhou¹

Academic Editor: Qiang Liu

Received17 Feb 2021

Revised26 Mar 2021

Accepted10 Apr 2021

Published23 Apr 2021

Abstract

The activities of the ephrin family in breast cancer (BrCa) are complex. Family A receptors (EPHA) and ligands (EFNA) can act as oncogenes or tumor suppressors and are implicated in chemoresistance. Here, we examined the expression pattern and prognostic value of the EPHA/EFNA family in patients with breast cancer, including patients with different subtypes or different chemotherapy cohorts. In the UALCAN database, the mRNA expression of EPHA1, EPHA10, EFNA1, EFNA3, and EFNA4 was significantly higher, whereas that of EPHA2, EPHA4, EPHA5, and EFNA5 was significantly lower in breast cancer tissues than in paracancerous tissues. The transcriptional levels of EPHA/EFNA family members were correlated with intrinsic subclasses of breast cancer. The relationship between EPHA/EFNA and the clinicopathological parameters of BrCa was analyzed using bc-GenExMiner V4.5. EPHA1, EPHA2, EPHA4, EPHA7, EFNA3, EFNA4, and EFNA5 were upregulated in estrogen receptor- (ER-) and progesterone receptor- (PR-) negative tumors, whereas EPHA3, EPHA6, and EFNA1 were upregulated in ER- and PR-positive tumors. EPHA1, EPHA2, EFNA3, and EFNA4 mRNA expression was significantly higher in human epidermal growth factor receptor 2- (HER2-) positive tumors than in HER2-negative tumors. Triple-negative status was positively correlated with EPHA1, EPHA2, EPHA4, EPHA7, EFNA3, EFNA4, and EFNA5 and negatively correlated with EPHA3 and EPHA10 mRNA expression. Genetic alterations of EPHA/EFNA in breast cancer varied from 1.1% to 10% for individual genes, as determined by the cBioPortal database. The Kaplan–Meier plotter indicated that high EphA7 mRNA expression was associated with poor overall survival (OS) and recurrence-free survival (RFS), especially in the HER2 and luminal A subtypes. EFNA4 was predicted to have poor OS and RFS in breast cancers, especially in luminal B, basal-like subtype, and patients treated with adjuvant chemotherapy. High EPHA3 expression was significantly associated with better OS and RFS, especially in the luminal A subtype, but with poor RFS in BrCa patients receiving chemotherapy. Our findings systematically elucidate the expression pattern and prognostic value of the EPHA/EFNA family in BrCa, which might provide potential prognostic factors and novel targets in BrCa patients, including those with different subtypes or treated with chemotherapy.

1. Introduction

Breast cancer is the most commonly diagnosed female cancer and the cause of 685,000 cancer mortality in 2020 worldwide [1]. Tumor recurrence and metastasis contribute to the high death rate [2]. Despite extensive research into the treatment of BrCa, chemotherapy resistance is an important issue limiting the efficacy of treatment. Novel biomarkers to predict prognosis or the sensitivity to chemotherapy are urgently needed.

Receptor tyrosine kinases (RTKs) play an important role in a variety of cellular processes in cancer [3]. Ephrins, also known as ephrin ligands, and Eph receptors (Ephs), which are RTKs, are key regulators of physiological and pathological processes involved in development and disease, such as cellular motility, cell repulsion, and cell adhesion [4]. The Ephrin family consists of multiple Ephs and ephrins. Both receptors and ligands are membrane-bound proteins that require direct cell-cell interaction for activation. Eph/ephrin signal transduction occurs not only in the receptor-expressing cell but also in the ligand-expressing cell via bidirectional signaling [5]. The ligands can have a glycosylphosphatidylinositol anchor (A type) or a membrane-spanning protein domain (B type). The receptors are also categorized as A or B according to the type of ligand they bind to. Ephrin family A includes ten receptors named EPHA (1–10) and five ligands designated as EFNA (1–5) [6]. The interaction between ligands and receptors via bidirectional signaling and its involvement in cancer biology are mediated by complex processes [7, 8]. Several Ephrin A (EPHA/EFNA) family members are overexpressed or downregulated in a variety of tumors, suggesting that they act as oncogenes or as tumor suppressors according to the cellular context [9]. Ephrin A family members that are overexpressed in cancer including EFNA1 in melanoma [10]; EPHA2 in prostate cancer [11], nasopharyngeal carcinoma [12], and squamous-cell carcinoma of the head and neck [13]; and EPHA3 in non-small-cell lung cancer [14]. EPHA7 acts as a tumor suppressor in follicular lymphoma and is a potential therapeutic target [15]. EphA2 expression levels are associated with the invasiveness and aggressive behavior of BrCa [16]. The role of EPHA/EFNA as tumor suppressors in breast carcinogenesis was also demonstrated by targeting EPHA2 [17]. EPHA2 regulates the sensitivity to paclitaxel in nasopharyngeal carcinoma via the phosphoinositide 3-kinase/Akt signaling pathway [18]. Cisplatin chemotherapy-induced ERK1/2-RSK1/2-EphA2-GPRC5A signaling is related to acquired chemoresistance in ovarian cancer [19]. Previous findings indicate that EPHA/EFNA family members may serve as biomarkers for predicting prognosis or the response to treatment, which prompted us to analyze the expression pattern and prognostic role of the EPHA/EFNA family in BrCa.

Online databases provide access to a wealth of information, such as microarray RNA chips from the Gene Expression Omnibus (GEO) [20] and RNA sequences from The Cancer Genome Atlas (TCGA) [21]. In this study, we compared the transcriptional levels of EPHA/EFNA in BrCa and paracancerous tissues using the UALCAN database [22]. The relationship between members of the EPHA/EFNA family and clinicopathologic characteristics of BrCa patients was analyzed using breast cancer gene expression miner (bc-GenExMiner) v4.5 [23, 24]. Genetic alterations of EPHA/EFNA family members, including mutations and putative copy number alterations (CNAs), were analyzed in cBioPortal [25]. Moreover, the association of EPHA/EFNA gene expression with clinical outcomes was assessed in patients with BrCa, including patients with different subtypes or those undergoing chemotherapy, using the Kaplan–Meier plotter database [26].

2. Methods and Materials

2.1. Gene Expression Analysis in UALCAN

UALCAN (http://ualcan.path.uab.edu/) is an online comprehensive and interactive platform based on RNA-seq data and clinical information from TCGA database [22]. In this study, UALCAN was used to investigate the transcriptional levels of EPHA/EFNA family members in primary BrCa tissues and their association with the subtype.

2.2. Analysis of the Relationship between EPHA/EFNA and Clinicopathologic Characteristics in Bc-GenExMiner

Bc-GenExMiner v4.5 (http://bcgenex.centregauducheau.fr/BC-GEM), a statistical mining tool of published BrCa transcriptomic data, was used to evaluate the relationship between EPHA/EFNA family members and the clinicopathologic characteristics of BrCa patients, including age, nodal status, estrogen receptor (ER) status, progesterone receptor (PR) status, human epidermal growth factor receptor 2 (HER2) status, lymph node status, triple-negative tumors, and P53 status (sequence). Information on ER, PR, HER2, lymph node status, and pathological grade was not available for all 3,996 patients.

2.3. Genetic Alterations Analysis in cBioPortal

OncoPrint is a feature of cBioPortal, an open-source web application that allows researchers to explore and analyze cancer genomic datasets (http://www.cbioportal.org/) [27, 28]. In this study, genetic alterations of EPHA/EFNA family members including mutations and putative CNAs were analyzed. A total of 1084 tumor samples with RNA-seq data on cBioPortal were included in the study.

2.4. The Prognostic Value of EPHA/EFNA mRNA Expression in BrCa Patients including Those with Different Subtypes and Undergoing Chemotherapy

The Kaplan–Meier plotter is an online database that facilitates predicting the effect of gene expression on survival in cancer. Sources for the database include GEO, EGA (European Genome-Phenome Archive), and TCGA. This platform, which contains gene expression information and survival data of BrCa patients, was used to perform a meta-analysis to verify the prognostic value of EPHA/EFNA family members for predicting overall survival (OS) and relapse-free survival (RFS). Additional analyses were restricted to cohorts according to subtype and including the patients treated with chemotherapy. Kaplan–Meier survival plots were used to compare the prognosis of all cohorts. Hazard ratios (HRs), 95% confidence intervals (CIs), and log-rank values were calculated and displayed online.

2.5. Statistical Analysis

The differential mRNA expression of EPHA/EFNA in BrCa tissues was compared by Student’s -test. The log-rank test was used to compare Kaplan-Meier survival plots. A value of <0.05 was considered statistically significant.

3. Results

3.1. Transcriptional Levels of EPHA/EFNA in BrCa

To evaluate the expression pattern of EPHA/EFNA in BrCa, the transcriptional levels of EPHA/EFNA family members were compared between BrCa and paracancerous tissues using the UALCAN database. As shown in Figure 1, EPHA2 (Figure 1(b), ), EPHA4 (Figure 1(d), ), and EPHA5 (Figure 1(e), ) expression was significantly lower in BrCa tissues than in paracancerous tissues. EPHA1 (Figure 1(a), ) and EPHA10 (Figure 1(h), ) expression was significantly higher in BrCa than in paracancerous tissues (Figure 1(f), ). The expression of EPHA3 (Figure 1(c), ), EPHA6 (Figure 1(f), ), and EPHA7 (Figure 1(g), ) did not differ significantly between BrCa and paracancerous tissues. Analysis of EFNA expression showed that the transcriptional levels of EFNA1 (Figure 1(i), ), EFNA3 (Figure 1(j), ), and EFNA4 (Figure 1(k), ) were significantly higher in BrCa than in paracancerous tissues, whereas the transcriptional level of EFNA5 (Figure 1(l), ) was significantly lower in BrCa than in paracancerous tissues. EPHA8 and EFNA2 were expressed at markedly low levels according to the UALCAN database, and they were not included in the analysis.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Figure 1

Transcriptional levels of EPHA/EFNA in breast cancer and paracancerous tissues from the TCGA dataset based on data mining via UALCAN. The transcriptional levels of EPHA1, EPHA10, EFNA1, EFNA3, and EFNA4 (a, h, i, j, and k) were upregulated, whereas the transcriptional levels of EPHA2, EPHA4, EPHA5, and EFNA5 (b, d, e, and l) were downregulated in breast cancer tissues compared with paracancerous tissues. The transcriptional levels of EPHA3, EPHA6, and EPHA7 (c, f, and g) did not differ significantly between breast cancer and paracancerous tissues.

3.2. Transcriptional Levels of EPHA/EFNA in Different BrCa Subtypes

The classification of BrCa into subtypes is helpful for predicting the therapeutic response and prognosis of patients [29]. Here, we analyzed the transcriptional levels of EPHA/EFNA family members according to BrCa subtype using the UALCAN database. As shown in Figure 2, EPHA3, EPHA5, EPHA6, and EPHA10 mRNA levels were low in HER2-positive and triple-negative BrCa patients, whereas other EPHA family members did not show this trend. EPHA3, EPHA5, EPHA6, and EPHA10 were expressed at high levels in the luminal subtype (Figures 2(c), 2(e), 2(f), and 2(h)). The highest mRNA expression levels of EFNA were detected in triple-negative tissues except EFNA3, whereas the luminal subtype showed the lowest EFNA mRNA levels (Figures 2(i)–2(l)) except EFNA1. Taken together, these findings indicate that the transcriptional levels of EPHA/EFNA family members are correlated with intrinsic subclasses in BrCa patients.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Figure 2

The transcriptional level of EPHA/EFNA in breast cancer patients with different subtypes. HER2-positive and triple-negative breast cancer patients tended to express lower levels of EPHA3, EPHA5, EPHA6, and EPHA10 mRNA. The highest mRNA expression of EPHA3, EPHA5, EPHA6, and EPHA10 was found in the luminal subtype. The highest mRNA expression of EFNA was detected in triple-negative tissues except EFNA3, and the lowest mRNA expression of EFNA was detected in the luminal subtype, except for EFNA1.

3.3. Association between EPHA/EFNA mRNA Expression and the Clinicopathological Features of Patients with BrCa

We used bc-GenExMiner v4.5 to examine the relationship between EPHA/EFNA and the clinicopathologic characteristics of patients. For the age parameter, EPHA2 (), EPHA3 (), EFNA4 (), and EFNA5 () were expressed at high levels in patients aged ≤51 years (Table 1). EFNA1 () and EFNA3 () were expressed at high levels in older patients, whereas the mRNA expression levels of the other EPHA/EFNA family members were not significantly correlated with age. EPHA2 () mRNA was higher in BrCa patients with negative lymph nodes than in those with positive lymph nodes, whereas EFNA3 () showed the opposite trend. The mRNA expression of the other EPHA/EFNA family members was not significantly associated with nodal status. ER- and PR-negative patients had higher levels of EPHA1 (), EPHA2 (), EPHA4 (), EPHA7 ( and =0.0046, respectively), EFNA3 (), EFNA4 (), and EFNA5 (). On the other hand, EPHA3 (), EPHA6 ( and =0.0066, respectively), and EFNA1 () were higher in ER- and PR-positive patients. The EPHA10 mRNA level was higher in ER-positive patients () but not significantly associated with PR status (). EPHA5 mRNA expression was not significantly associated with ER () and PR () status. EPHA1, EPHA2, EFNA3, and EFNA4 (, all) mRNA levels were significantly higher in the HER2-positive group than in the HER2-negative group. Only EPHA6 () and EFNA1 () were significantly increased in the HER2-negative group. The mRNA expression of the other EPHA/EFNA family members was not associated with HER2 status. Triple-negative status was positively correlated with the mRNA expression of EPHA1 (), EPHA2 (), EPHA4 (), EPHA7 (), EFNA3 (), EFNA4 (), and EFNA5 () and negatively correlated with EPHA3 () and EPHA10 (). The mRNA expression of the other EPHA/EFNA family members was not associated with triple-negative status. P53 mutant status (sequence) was positively correlated with EPHA1, EPHA2, EFNA3, EFNA4, and EFNA5 ( for all) and P53 wild-type status was positively correlated with EPHA3 (), EPHA6 (), and EFNA1 () mRNA expression. The mRNA expression of the other EPHA/EFNA family members was not associated with P53 mutant status.

3.4. Genetic Alterations of EPHA/EFNA in BrCa

Different kinds of genetic alterations, such as missense mutations, amplification, and deep deletions, regulate cancer-related gene expression and participate in oncogenesis. We speculated that genetic alterations may regulate the transcriptional levels of EPHA/EFNA. To investigate mutations and CNAs of the EPHA/EFNA family in BrCa, the OncoPrint feature of cBioPortal (http://www.cbioportal.org) was used to investigate the proportion and percentage of specimens with genetic alterations in EPHA/EFNA. The frequency of alterations in these genes among BrCa samples varied from 1.1% to 10% for individual genes as shown in Figure 3.

3.5. The Prognostic Value of EPHA/EFNA mRNA Expression in BrCa Patients

The prognostic value of the mRNA expression of 12 EPHA/EFNA family members in BrCa patients was examined using the Kaplan–Meier plotter. High mRNA expression of EPHA7 (Figure 4(g), , 95% CI: 1.09–2.05, ) and EFNA4 (Figure 4(k), , 95% CI: 1.04–1.64, ) was associated with poor OS. EPHA6 mRNA expression (Figure 4(f), , 95% CI: 0.98–1.84, ) was moderately associated with poor OS. High mRNA expression of EPHA3 (Figure 4(c), , 95% CI: 0.57–0.89, ) and EPHA4 (Figure 4(d), , 95% CI: 0.52–0.99, ) was significantly associated with better OS. The mRNA expression levels of the other EPHA/EFNA family members were not significantly correlated with OS.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Regarding RFS, high mRNA expression of EPHA7 (Figure 5(g), , 95% CI: 1.04–1.42, ), EFNA3 (Figure 5(j), , 95% CI: 1.04–1.33, ), and EFNA4 (Figure 5(k), , 95% CI: 1.14–1.44, ) was associated with worse RFS. High mRNA expression of other EPHA/EFNA family members was significantly associated with better RFS except EFNA1. The RFS curves are shown in Figure 5 (EPHA1: , 95% CI: 0.61–0.76, ; EPHA2: , 95% CI: 0.72–0.89, ; EPHA3: , 95% CI: 0.7–0.87, ; EPHA4: , 95% CI: 0.62–0.87, ; EPHA5: , 95% CI: 0.61–0.86, ; EPHA6: , 95% CI: 0.69–0.97, ; EPHA10: , 95% CI: 0.52–0.71, ; and EFNA5: , 95% CI: 0.66–0.82, ).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

3.6. The Prognostic Value of EPHA/EFNA mRNA Expression in BrCa Patients with Different Subtypes

To further analyze the effect of EPHA/EFNA according to the BrCa subtype, the prognostic value of EPHA/EFNA family members was assessed in BrCa patients with different molecular subtypes, including basal-like, luminal A, luminal B, and HER2+ subtypes according to the 2011 St. Gallen criteria [30]. Because OS data were lacking for some patients, this analysis was limited to RFS. In the basal-like subtype, high mRNA expression of EFNA1 (, 95% CI: 1.06–1.76, ) and EFNA4 (, 95% CI: 1.16–2.01, ) predicted an unfavorable RFS, whereas high mRNA expression levels of EPHA1 (, 95% CI: 0.48–0.79, ), EPHA4 (, 95% CI: 0.49–0.96, ), EPHA5 (, 95% CI: 0.29–0.68, ), EPHA7 (, 95% CI: 0.49–0.95, ), and EFNA5 (, 95% CI: 0.54–0.89, ) were correlated with better RFS. The remaining EPHA/EFNA members were not associated with prognosis in basal-like BrCa. (Table 2).

In the luminal A subtype, high mRNA expression levels of EPHA1 (, 95% CI: 0.53–0.75, ), EPHA2 (, 95% CI:0.51–0.72, ), EPHA3 (, 95% CI: 0.61–0.87, ), EPHA4 (, 95% CI: 0.54–0.88, ), EPHA10 (, 95% CI: 0.42–0.7, ), EFNA3 (, 95% CI: 0.68–0.96, ), EFNA4 (, 95% CI: 0.6–0.86, ), and EFNA5 (, 95% CI: 0.51–0.73, ) were associated with better RFS, whereas the high mRNA expression level of EPHA7 (49, 95% CI: 1.16–1.91, ) was associated with unfavorable RFS. The remaining EPHA/EFNA members were not associated with prognosis in luminal A BrCa.

In the luminal B subtype, high mRNA expression levels of EPHA1 (, 95% CI: 0.52–0.76, ), EPHA2 (, 95% CI: 0.63–0.93, ), EPHA4 (, 95% CI: 0.42–0.79, ), EPHA5 (, 95% CI: 0.48–0.89, ), EPHA6 (, 95% CI: 0.45–0.85, ), EPHA10 (, 95% CI: 0.4–0.75, ), and EFNA5 (, 95% CI: 0.55–0.83, ) were associated with better RFS, whereas high mRNA expression levels of EFNA1 (, 95% CI: 1.09–1.75, ), EFNA3 (, 95% CI: 1.1–1.76, ), and EFNA4 (, 95% CI: 1.05–1.6, ) were associated with worse RFS. The remaining EPHA/EFNA members were not associated with prognosis in luminal B BrCa.

In HER2+ BrCa patients, high mRNA expression levels of EPHA1 (, 95% CI: 0.37–0.8, ), EPHA6 (, 95% CI: 0.29–0.94, ), and EFNA5 (, 95% CI: 0.45–0.97, ) were correlated with better RFS. Only high mRNA expression of EPHA7 (, 95% CI: 1.45–3.72, ) was associated with worse RFS. The remaining EPHA/EFNA members were not associated with prognosis in HER2-overexpressing BrCa. These results indicate that EPHA/EFNA may serve as potential prognostic predictors in BrCa patients with different subtype.

3.7. The Prognostic Value of EPHA/EFNA mRNA Expression in BrCa Patients Treated with Chemotherapy

The prognostic value of EPHA/EFNA expression was analyzed in BrCa patients receiving different chemotherapy regimens, including adjuvant chemotherapy, neoadjuvant chemotherapy, and no chemotherapy. As shown in Table 3, high expression of EPHA2 (, 95% CI: 1.08–2.06, ), EPHA3 (, 95% CI: 1.12–2.15, ), EPHA4 (, 95% CI: 1.2–3.18, ), EFNA3 (, 95% CI: 1.44–2.65, ), and EFNA4 (, 95% CI: 1.02–1.93, ) and low expression of EPHA5 (, 95% CI: 0.25–0.93, ), EPHA6 (, 95% CI: 0.27–0.95, ), and EPHA10 (, 95% CI: 0.29–0.8, ) were significantly correlated with poor RFS in BrCa patients treated with adjuvant chemotherapy. High expression levels of EPHA3 (, 95% CI: 1–3.8, ) and EPHA10 (, 95% CI: 1.16–5.9, ) and low expression of EPHA1 (, 95% CI: 0.22–0.68, ), EPHA2 (, 95% CI: 0.31–0.97, ), and EPHA5 (, 95% CI: 0.16–0.77, ) were significantly correlated with poor RFS in BrCa patients treated with neoadjuvant chemotherapy. These results indicate that EPHA/EFNA may serve as potential prognostic factors in BrCa patients treated with chemotherapy, suggesting that these genes are potential targets for the treatment of BrCa.

4. Discussion

The activities of the EPHA/EFNA family in BrCa are complex and paradoxical. The expression and prognostic value of EPHA/EFNA in BrCa have not been extensively investigated. In the present study, two cancer databases were used to analyze the transcriptional levels of EPHA/EFNA family members in BrCa and paracancerous tissues, as well as their association with the BrCa subtype and clinicopathological features. The genetic alterations of EPHA/EFNA family members, including mutations and putative CNAs, were analyzed using cBioPortal. The Kaplan–Meier plotter was used to analyze the association between the expression levels of EPHA/EFNA and OS or RFS in BrCa patients, as well as RFS in different BrCa subtypes, including patients undergoing chemotherapy. The results may be valuable in identifying new BrCa biomarkers for predicting prognosis or sensitivity to chemotherapy and suggest that EPHA and EFNA play both oncogenic and tumor suppressor roles in BrCa.

Expression analysis showed that the transcriptional levels of EPHA2, EPHA4, and EPHA5 were significantly lower in BrCa tissues than in nontumor tissues and EPHA1 and EPHA10 were significantly upregulated in BrCa tissues. Brantley et al. investigated EPHA/EFNA protein expression in BrCa and found that EPHA2, EPHA4, and EPHA7 were significantly upregulated in BrCa samples relative to normal controls. The discrepancy between the protein and mRNA levels of EPHA2 may be due to the fact that high levels of EPHA2 in tumor cells are the result of increased protein stability [31]. In certain malignant breast cell models, EPHA2 protein levels are 50- to 500-fold higher despite comparable levels of EphA2 mRNA [32, 33]. Discrepancies in EPHA2 expression may also result from the inclusion of both invasive and noninvasive breast tumors in the TCGA database. The mRNA level of EPHA2 is higher in invasive tumors than in normal breast cells; however, EPHA2 expression is lower in noninvasive breast tumors than in normal breast cells [34]. In this study, EFNA was upregulated in BrCa tissues except for EFNA5. EPHA2 is the dominant and the best characterized EPHA receptor in the BrCa. The role of EPHA2 in breast tumor progression is controversial, and conflicting data on the clinical significance of EPHA2 have been reported in different studies [35, 36]. For example, data show that EPHA2 is overexpressed in BrCa clinical samples; however, there is also evidence that EPHA2 acts as a tumor suppressor in breast carcinogenesis [17]. The malignant behavior of EPHA2 is mediated by ligand-independent signaling, and its antioncogenic properties are attributed to ligand-dependent signaling [37]. The crosstalk between EphA2 and BrCa oncogenic pathways promotes tumor cell malignancy in ligand-independent signaling [38, 39]. We found an inverse correlation between EPHA2 and the mRNA expression of A-type ligands in the database, as shown by the downregulation of EPHA2 and EFNA overexpression in BrCa cells. Ligand upregulation in the tumor indicates that the ligand-dependent pathway is dominant in the database. Downregulation of the EPHA2 receptor by the ephrin ligand may involve ligand-mediated receptor internalization. In general, EPHA2 negatively regulates tumor growth and migration after canonical ligand-induced EPHA2 signaling, which inhibits the AKT-mTORC1 and MAPK pathways [17]. The interaction between EPHA2 and its ligand activates a negative feedback pathway mediated by growth factor-activated RAS signaling [40]. EFNA1 is upregulated in noninvasive breast cells, thereby inhibiting invasiveness, whereas EFNA1 is downregulated in invasive tumors, allowing EPHA2 to participate in invasion [34].

We also compared the differential transcriptional levels of EPHA/EFNA family members according to the different intrinsic subtypes of BrCa. The results showed differences in the expression patterns between the BrCa subtypes. EPHA3, EPHA5, EPHA6, and EPHA10 were expressed at the highest levels in luminal tissues, whereas HER2-positive and triple-negative tissues tended to express lower levels of EPHA3, EPHA5, EPHA6, and EPHA10. The highest mRNA expression levels of EFNA were found in triple-negative tissues except EFNA3, whereas the lowest mRNA expression levels of EFNA were found in the luminal subtype except EFNA1. Analysis of the bc-GenExMiner database showed that EPHA1, EPHA2, EPHA4, EPHA7, EFNA3, EFNA4, and EFNA5 were upregulated in ER- and PR-negative patients and positively correlated with triple-negative status. Consistent with this, some studies have shown that EPHA2 overexpression in BrCa is negatively correlated with ER and PR status [32, 41, 42]. EPHA2 overexpression decreases estrogen dependence and tamoxifen sensitivity [43], and EPHA2 is preferentially expressed in the basal-like phenotype [44]. In contrast, EPHA2 expression is also negatively regulated by ERα [41] and wild-type p53 [45]. The correlation of EPHA2 with HER2-positive status in the present study is consistent with the results of previous studies [39, 46], and EPHA2 is associated with resistance to trastuzumab therapy [47].

Analysis of the other EPHA/EFNA family members showed that EPHA1 was significantly upregulated in BrCa tissues compared to paracancerous tissues. A previous report showed that EPHA1 downregulation is associated with the invasiveness of breast carcinoma cells [34]. Consistent with the present profiling analyses, EPHA4 expression is associated with basal-like BrCa [48] and EPHA5 has been reported to act as a tumor suppressor, which may be related to abrrant promoter methylation [49]. Liu et al. [50] showed that EPHA7 mRNA was downregulated in BrCa specimens and loss of EPHA7 expression is more common in high-grade, early TNM-stage patients, without lymph node metastasis and correlation with negative HER2 status. Correspondingly, the present data also found downregulation of EPHA7 expression in BrCa specimens. EPHA10 is the only kinase-deficient Eph receptor [51]. The present data indicate that EPHA10 is upregulated in BrCa tissues, which is similar with a previous study showing that EPHA10 expression is significantly lower in invasive than in noninvasive breast tumors, and is absent in normal cells [34]. Similar with the previous study, the present data showed EFNA4 was upregulated in TNBC [52]. EFNA4 is required for proper differentiation and polarization of mammary epithelial cells, signifying a biologic basis for the overexpression of EFNA4 in BrCa [53].

The most common cancer-related genetic alteration is the DNA CNAs. The OncoPrint feature of cBioPortal was used to determine the frequency of genetic alterations in the EPHA/EFNA family. The results showed that EPHA/EFNA was not frequently amplified. This finding suggests that the EPHA/EFNA family does not affect BrCa survival through DNA alterations, whereas it may affect BrCa through alterations of the interaction network. The copy-gain frequency of EFNA1, EFNA3, and EFNA4 accounts for a large proportion of BrCa samples, which may explain the significant upregulation of EFNA1, EFNA3, and EFNA4 in BrCa tissues from UALCAN.

The Kaplan–Meier plotter is a survey of public microarray data repositories of survival from 3955 patients with BrCa. Besides, in general patients, we analyzed the prognostic value of EPHA/EFNA mRNA in different BrCa subtypes and in BrCa patients treated with chemotherapy. Large-scale expression profiling studies revealed a negative association between the overexpression of EPHA2, EPHA4, and EPHA7 and overall and disease-free survival in BrCa [54]. Consistent with this report and its positive correlation with the triple-negative (TNBC) status in the bc-GenExMiner, the Kaplan–Meier plotter indicated that high EPHA7 mRNA expression predicted a poor OS and RFS in BrCa, especially in the HER2+ and luminal A subtypes. EFNA4 was associated with poor OS and RFS in BrCa, which is consistent with its positive correlation with the TNBC status. In particular, EFNA4 was associated with poor RFS in luminal B and basal-like subtypes and in BrCa patients treated with adjuvant chemotherapy. High mRNA expression of EPHA3 was significantly associated with better OS and RFS, especially in the luminal A subtype, which is consistent with its negative correlation with the TNBC status. However, in patients receiving adjuvant chemotherapy and neoadjuvant chemotherapy, EPHA3 is a risk factor. These results suggest that EPHA3, EPHA7 and EFNA4 are involved in BrCa development and could predict the prognosis of patients with BrCa in various subtypes and in those receiving different chemotherapy regimens. BrCa is a heterogeneous disease with subtype-dependent histopathological and clinical significance. In the subgroup analysis, the triple-negative (basal-like) subtype is a unique subtype of BrCa with a poor prognosis and more likely to develop chemoresistance [55]. In this study, we showed that low expression of EPHA1, EPHA4, EPHA5, and EPHA7 and high expression of EFNA1, EFNA4, and EFNA5 predicted an unfavorable prognosis in basal-like patients. Although we found that high mRNA expression of EPHA2 was associated with a longer RFS in patients with BrCa, it was a risk factor in the cohort receiving adjuvant chemotherapy. In fact, the RSK1/2-EphA2-GPRC5A oncogenic signaling association with platinum chemotherapy resistance has been reported in recent study [19] . Furthermore, in clinical practice, these biomarkers coupled with the specific role of EPHA/EFNA signaling in cancer could promote targeted therapeutics, as reported before [52, 56].

The present study has several limitations that should be addressed in future studies. First, mRNA levels are not always indicative of a functional protein and thus may not fully represent the protein expression of EPHA/EFNA [25]. Future studies should include protein detection techniques to accurately assess the protein levels. Second, the function of EPHA3 and EFNA7 in BrCa has not been studied in BrCa and needs to be investigated in the future. Third, multivariate analysis could not be used in the database to correct the associations between different clinicopathological features.

In conclusion, we showed that the EPHA/EFNA family is widely expressed in BrCa tumor cells. EPHA/EFNA were identified as prognostic factors and potential targets for BrCa, which may improve our understanding of the complexity and heterogeneity of BrCa at the molecular level. The present results may help develop tools to accurately predict prognosis and design customized therapies.

Abbreviations

EPHA:	Ephrin A
HR:	Hazard ratio
BrCa:	Breast cancer
KM plotter:	Kaplan–Meier plotter
RFS:	Relapse-free survival
RTK:	Tyrosine kinase receptor
CI:	Confidence interval
ER:	Estrogen receptor
PR:	Progesterone receptor
HER2:	Human epidermal growth factor receptor 2
SPSS:	Statistical product and service solutions
OS:	Overall survival.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no potential conflicts of interest.

Authors’ Contributions

ZL conceived and designed the study, analyzed the data, and wrote the manuscript. XW conceived and designed the study. ZL and KD performed the statistical analysis and analyzed the data. XL and CQ participated in data preparation, analysis, and figure preparation. ZHF revised the final manuscript. All authors have read and approved the manuscript for publication.

Acknowledgments

This study was financially supported by the Natural Science Foundation of Tianjin City, China (Grant no. 17JCYBJC25600), and the project of the Tianjin Health Bureau (grant no. 2014KZ124).

References

H. Sung, J. Ferlay, R. Siegel, I. Soerjomataram, A. Jemal, and F. Bray, “Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: a Cancer Journal for Clinicians, pp. 1–41, 2021.
View at: Google Scholar
D. A. Berry, K. A. Cronin, S. K. Plevritis et al., “Effect of screening and adjuvant therapy on mortality from breast cancer,” The New England Journal of Medicine, vol. 353, no. 17, pp. 1784–1792, 2005.
View at: Publisher Site | Google Scholar
Z. Du and C. M. Lovly, “Mechanisms of receptor tyrosine kinase activation in cancer,” Molecular Cancer, vol. 17, no. 1, p. 58, 2018.
View at: Publisher Site | Google Scholar
E. M. Lisabeth, G. Falivelli, and E. B. Pasquale, “Eph receptor signaling and ephrins,” Cold Spring Harbor Perspectives in Biology, vol. 5, no. 9, 2013.
View at: Publisher Site | Google Scholar
E. B. Pasquale, “Eph-ephrin bidirectional signaling in physiology and disease,” Cell, vol. 133, pp. 38–52, 2008.
View at: Publisher Site | Google Scholar
Eph Nomenclature Committee, “Unified nomenclature for Eph family receptors and their ligands, the ephrins,” Cell, vol. 90, no. 3, pp. 403-404, 1997.
View at: Publisher Site | Google Scholar
E. B. Pasquale, “Eph receptors and ephrins in cancer: bidirectional signalling and beyond,” Nature Reviews. Cancer, vol. 10, no. 3, pp. 165–180, 2010.
View at: Publisher Site | Google Scholar
C.-T. J. Kou and R. P. Kandpal, “Differential expression patterns of Eph receptors and Ephrin ligands in human cancers,” BioMed Research International, vol. 2018, Article ID 7390104, 23 pages, 2018.
View at: Publisher Site | Google Scholar
B. Mosch, B. Reissenweber, C. Neuber, and J. Pietzsch, “Eph receptors and ephrin ligands: important players in angiogenesis and tumor angiogenesis,” Journal of Oncology, vol. 2010, Article ID 135285, 12 pages, 2010.
View at: Publisher Site | Google Scholar
D. J. Easty, S. P. Hill, M. Y. Hsu et al., “Up-regulation of ephrin-A1 during melanoma progression,” International Journal of Cancer, vol. 84, no. 5, pp. 494–501, 1999.
View at: Publisher Site | Google Scholar
J. Walker-Daniels, K. Coffman, M. Azimi et al., “Overexpression of the EphA2 tyrosine kinase in prostate cancer,” Prostate, vol. 41, no. 4, pp. 275–280, 1999.
View at: Publisher Site | Google Scholar
Y.-P. Xiang, T. Xiao, Q. G. Li et al., “Y772 phosphorylation of EphA2 is responsible for EphA2-dependent NPC nasopharyngeal carcinoma growth by Shp2/Erk-1/2 signaling pathway,” Cell Death & Disease, vol. 11, no. 8, p. 709, 2020.
View at: Publisher Site | Google Scholar
Y. Liu, X. Zhang, Y. Qiu et al., “Clinical significance of EphA2 expression in squamous-cell carcinoma of the head and neck,” Journal of Cancer Research and Clinical Oncology, vol. 137, no. 5, pp. 761–769, 2011.
View at: Publisher Site | Google Scholar
R. Govindan, L. Ding, M. Griffith et al., “Genomic landscape of non-small cell lung cancer in smokers and never-smokers,” Cell, vol. 150, no. 6, pp. 1121–1134, 2012.
View at: Publisher Site | Google Scholar
E. Oricchio, G. Nanjangud, A. L. Wolfe et al., “The Eph-receptor A7 is a soluble tumor suppressor for follicular lymphoma,” Cell, vol. 147, no. 3, pp. 554–564, 2011.
View at: Publisher Site | Google Scholar
J. Walker-Daniels, A. R. Hess, M. J. C. Hendrix, and M. S. Kinch, “Differential regulation of EphA2 in normal and malignant cells,” The American Journal of Pathology, vol. 162, no. 4, pp. 1037–1042, 2003.
View at: Publisher Site | Google Scholar
L. W. Noblitt, D. S. Bangari, S. Shukla et al., “Decreased tumorigenic potential of EphA2-overexpressing breast cancer cells following treatment with adenoviral vectors that express EphrinA1,” Cancer Gene Therapy, vol. 11, no. 11, pp. 757–766, 2004.
View at: Publisher Site | Google Scholar
Y. Wang, Y. Liu, G. Li et al., “Ephrin type-A receptor 2 regulates sensitivity to paclitaxel in nasopharyngeal carcinoma via the phosphoinositide 3-kinase/Akt signalling pathway,” Molecular Medicine Reports, vol. 11, no. 2, pp. 924–930, 2015.
View at: Publisher Site | Google Scholar
L. Moyano-Galceran, E. A. Pietilä, S. P. Turunen et al., “Adaptive RSK-EphA2-GPRC5A signaling switch triggers chemotherapy resistance in ovarian cancer,” EMBO Molecular Medicine, vol. 12, no. 4, article e11177, 2020.
View at: Publisher Site | Google Scholar
T. Barrett, S. E. Wilhite, P. Ledoux et al., “NCBI GEO: archive for functional genomics data sets—update,” Nucleic Acids Research, vol. 41, no. D1, pp. D991–D995, 2012.
View at: Publisher Site | Google Scholar
The Cancer Genome Atlas Research Network, J. N. Weinstein, E. A. Collisson et al., “The Cancer Genome Atlas Pan-Cancer analysis project,” Nature Genetics, vol. 45, no. 10, pp. 1113–1120, 2013.
View at: Publisher Site | Google Scholar
D. S. Chandrashekar, B. Bashel, S. A. H. Balasubramanya et al., “UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses,” Neoplasia, vol. 19, no. 8, pp. 649–658, 2017.
View at: Publisher Site | Google Scholar
P. Jézéquel, M. Campone, W. Gouraud et al., “bc-GenExMiner: an easy-to-use online platform for gene prognostic analyses in breast cancer,” Breast Cancer Research and Treatment, vol. 131, no. 3, pp. 765–775, 2012.
View at: Publisher Site | Google Scholar
P. Jezequel, J. S. Frenel, L. Campion et al., “bc-GenExMiner 3.0: new mining module computes breast cancer gene expression correlation analyses,” Database, vol. 2013, article bas060, 2013.
View at: Publisher Site | Google Scholar
J. Mei, L. Hao, X. Liu et al., “Comprehensive analysis of peroxiredoxins expression profiles and prognostic values in breast cancer,” Biomarker Research, vol. 7, no. 1, p. 16, 2019.
View at: Publisher Site | Google Scholar
B. Györffy, A. Lanczky, A. C. Eklund et al., “An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients,” Breast Cancer Research and Treatment, vol. 123, no. 3, pp. 725–731, 2010.
View at: Publisher Site | Google Scholar
E. Cerami, J. Gao, U. Dogrusoz et al., “The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data,” AACR., vol. 2, no. 5, pp. 401–404, 2012.
View at: Publisher Site | Google Scholar
J. Gao, B. A. Aksoy, U. Dogrusoz et al., “Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal,” Science Signaling, vol. 6, no. 269, p. pl1, 2013.
View at: Publisher Site | Google Scholar
W. Goto, S. Kashiwagi, K. Takada et al., “Significance of intrinsic breast cancer subtypes on the long-term prognosis after neoadjuvant chemotherapy,” Journal of Translational Medicine, vol. 16, no. 1, p. 307, 2018.
View at: Publisher Site | Google Scholar
A. Goldhirsch, W. C. Wood, A. S. Coates et al., “Strategies for subtypes--dealing with the diversity of breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011,” Annals of Oncology, vol. 22, no. 8, pp. 1736–1747, 2011.
View at: Publisher Site | Google Scholar
J. Walker-Daniels, D. J. Riese, and M. S. Kinch, “c-Cbl-dependent EphA2 protein degradation is induced by ligand binding,” Molecular Cancer Research, vol. 1, no. 1, pp. 79–87, 2002.
View at: Google Scholar
D. P. Zelinski, N. D. Zantek, J. C. Stewart, A. R. Irizarry, and M. S. Kinch, “EphA2 overexpression causes tumorigenesis of mammary epithelial cells,” Cancer Research, vol. 61, no. 5, pp. 2301–2306, 2001.
View at: Google Scholar
N. D. Zantek, J. Walker-Daniels, J. Stewart et al., “MCF-10A-NeoST: a new cell system for studying cell-ECM and cell-cell interactions in breast cancer,” Clinical Cancer Research, vol. 7, no. 11, pp. 3640–3648, 2001.
View at: Google Scholar
B. P. Fox and R. P. Kandpal, “Invasiveness of breast carcinoma cells and transcript profile: Eph receptors and ephrin ligands as molecular markers of potential diagnostic and prognostic application,” Biochemical and Biophysical Research Communications, vol. 318, no. 4, pp. 882–892, 2004.
View at: Publisher Site | Google Scholar
M. Pan, “Overexpression ofEphA2gene in invasive human breast cancer and its association with hormone receptor status,” Journal of Clinical Oncology, vol. 23, 16_suppl, p. 9583, 2005.
View at: Publisher Site | Google Scholar
K. Ogawa, R. Pasqualini, R. A. Lindberg, R. Kain, A. L. Freeman, and E. B. Pasquale, “The ephrin-A1 ligand and its receptor, EphA2, are expressed during tumor neovascularization,” Oncogene, vol. 19, no. 52, pp. 6043–6052, 2000.
View at: Publisher Site | Google Scholar
H. Miao, D. Q. Li, A. Mukherjee et al., “EphA2 mediates ligand-dependent inhibition and ligand-independent promotion of cell migration and invasion via a reciprocal regulatory loop with Akt,” Cancer Cell, vol. 16, no. 1, pp. 9–20, 2009.
View at: Publisher Site | Google Scholar
D. M. Brantley-Sieders, G. Zhuang, D. Hicks et al., “The receptor tyrosine kinase EphA2 promotes mammary adenocarcinoma tumorigenesis and metastatic progression in mice by amplifying ErbB2 signaling,” The Journal of Clinical Investigation, vol. 118, no. 1, pp. 64–78, 2008.
View at: Publisher Site | Google Scholar
G. Zhuang, D. Brantley-Sieders, D. Hicks, W. Fang, Y. Hwang, and J. Chen, “Interactions of EphA2 and HER2 promote tumor progression and anti-HER2 resistance,” AACR., vol. 68, 2008.
View at: Google Scholar
M. Macrae, R. M. Neve, P. Rodriguez-Viciana et al., “A conditional feedback loop regulates Ras activity through EphA2,” Cancer Cell, vol. 8, no. 2, pp. 111–118, 2005.
View at: Publisher Site | Google Scholar
D. P. Zelinski, N. D. Zantek, J. Walker-Daniels, M. A. Peters, E. J. Taparowsky, and M. S. Kinch, “Estrogen and Myc negatively regulate expression of the EphA2 tyrosine kinase,” Journal of Cellular Biochemistry, vol. 85, no. 4, pp. 714–720, 2002.
View at: Publisher Site | Google Scholar
A. A. Kamat, D. Coffey, W. M. Merritt et al., “EphA2 overexpression is associated with lack of hormone receptor expression and poor outcome in endometrial cancer,” Cancer, vol. 115, no. 12, pp. 2684–2692, 2009.
View at: Publisher Site | Google Scholar
M. Lu, K. D. Miller, Y. Gokmen-Polar, M. H. Jeng, and M. S. Kinch, “EphA2 overexpression decreases estrogen dependence and tamoxifen sensitivity,” Cancer Research, vol. 63, no. 12, pp. 3425–3429, 2003.
View at: Google Scholar
F. Huang, K. Reeves, X. Han et al., “Identification of candidate molecular markers predicting sensitivity in solid tumors to dasatinib: rationale for patient selection,” Cancer Research, vol. 67, no. 5, pp. 2226–2238, 2007.
View at: Publisher Site | Google Scholar
M. Dohn, J. Jiang, and X. Chen, “Receptor tyrosine kinase EphA2 is regulated by p53-family proteins and induces apoptosis,” Oncogene, vol. 20, no. 45, pp. 6503–6515, 2001.
View at: Publisher Site | Google Scholar
V. M. Youngblood, L. C. Kim, D. N. Edwards et al., “The Ephrin-A1/EPHA2 signaling axis regulates glutamine metabolism in HER2-positive breast cancer,” Cancer Research, vol. 76, no. 7, pp. 1825–1836, 2016.
View at: Publisher Site | Google Scholar
G. Zhuang, D. M. Brantley-Sieders, D. Vaught et al., “Elevation of receptor tyrosine kinase EphA2 mediates resistance to trastuzumab therapy,” Cancer Research, vol. 70, no. 1, pp. 299–308, 2010.
View at: Publisher Site | Google Scholar
I. Y. Hachim, M. Villatoro, L. Canaff et al., “Transforming growth factor-beta regulation of ephrin type-A receptor 4 signaling in breast cancer cellular migration,” Scientific Reports, vol. 7, no. 1, p. 14976, 2017.
View at: Publisher Site | Google Scholar
D. Y. Fu, Z. M. Wang, B. L. Wang et al., “Frequent epigenetic inactivation of the receptor tyrosine kinase EphA5 by promoter methylation in human breast cancer,” Human Pathology, vol. 41, no. 1, pp. 48–58, 2010.
View at: Publisher Site | Google Scholar
Z. Liu, Q. Zhang, X. Li, and Z. Tao, “Aberrant expression of receptor tyrosine kinase EphA7 in breast cancers,” International Journal of Clinical and Experimental Pathology, vol. 9, no. 7, pp. 7352–7358, 2016.
View at: Google Scholar
K. K. Murai and E. B. Pasquale, “‘Eph’ective signaling: forward, reverse and crosstalk,” Journal of Cell Science, vol. 116, no. 14, pp. 2823–2832, 2003.
View at: Publisher Site | Google Scholar
M. Damelin, A. Bankovich, A. Park et al., “Anti-EFNA4 calicheamicin conjugates effectively target triple-negative breast and ovarian tumor-initiating cells to result in sustained tumor regressions,” Clinical Cancer Research, vol. 21, no. 18, pp. 4165–4173, 2015.
View at: Publisher Site | Google Scholar
A. Burleigh, S. McKinney, J. Brimhall et al., “A co-culture genome-wide RNAi screen with mammary epithelial cells reveals transmembrane signals required for growth and differentiation,” Breast Cancer Research, vol. 17, no. 1, p. 4, 2015.
View at: Publisher Site | Google Scholar
D. M. Brantley-Sieders, A. Jiang, K. Sarma et al., “Eph/ephrin profiling in human breast cancer reveals significant associations between expression level and clinical outcome,” PLoS One, vol. 6, no. 9, article e24426, 2011.
View at: Publisher Site | Google Scholar
S. Boichuk, A. Galembikova, A. Sitenkov et al., “Establishment and characterization of a triple negative basal-like breast cancer cell line with multi-drug resistance,” Oncology Letters, vol. 14, no. 4, pp. 5039–5045, 2017.
View at: Publisher Site | Google Scholar
J. W. Lee, H. D. Han, M. M. K. Shahzad et al., “EphA2 immunoconjugate as molecularly targeted chemotherapy for ovarian carcinoma,” Journal of the National Cancer Institute, vol. 101, no. 17, pp. 1193–1205, 2009.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Zheng Liang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

737

Downloads

1104

Citations