Abstract
Purpose
Triple-negative breast cancer (TNBC) is an aggressive subtype of breast cancer with a poor prognosis and a lack of targeted therapy. Overexpression of FRAT1 is thought to be associated with this aggressive subtype of cancer. Here, we performed a comprehensive analysis and assessed the association between overexpression of FRAT1 and TNBC.
Methods
First, using different web-based bioinformatics platforms (TIMER 2.0, UALCAN, and GEPIA 2), the expression of FRAT1 was assessed. Then, the expression of the FRAT1 protein and hormone receptors and HER2 status were assessed by immunohistochemical analysis. For samples of tumors with equivocal immunoreactivity, we performed silver in situ hybridization of the HER2 gene to determine an accurate HER2 status. Next, we used the R package and bc-GenExMiner 4.8 to analyze the relationship between FRAT1 expression and clinicopathological parameters in breast cancer patients. Finally, we determined the relationship between FRAT1 overexpression and prognosis in patients.
Results
The expression of FRAT1 in breast cancer tissues is significantly higher than in normal tissue. FRAT1 expression was significantly related to worse overall survival (P < 0.05) and was correlated with these clinicopathological features: T stage, N stage, age, high histologic grade, estrogen receptor status, progesterone receptor status, Her-2 status, TNBC status, basal-like status, CK5/6 status, and Ki67 status.
Breast cancer is the most common cancer and is the leading cause of cancer death among women [1]. According to Korean National Cancer Registration data, the incidence rate of breast cancer is increasing annually [2]. Breast cancer is a heterogeneous disease in terms of histology, therapeutic response, and patient outcomes [3]. Triple-negative breast cancer (TNBC) is a highly heterogeneous and clinically aggressive breast cancer that accounts for 15%–20% of all breast cancers and has the highest death rate of all breast cancer subtypes [456]. Because of its heterogeneity, various combined therapeutic modalities are implemented. Clinicopathological factors such as TNM stage, histological grade, estrogen receptor (ER), progesterone receptor (PR) expression, and human epidermal growth factor receptor 2 (HER-2) gene amplification are considered prognostic factors which provide the key for therapeutic strategies. Therefore, it is essential to define specific tumor markers to reflect the prognosis.
FRAT1 (frequently rearranged in advanced T-cell lymphomas 1), a member of the Frat family, is located on human chromosome 10q24.1, encoding a 29-kDa protein comprising 279 amino acids [7]. FRAT1 has been considered to be a positive activator of the Wnt/beta-catenin (β-catenin) signaling pathway, which is a family of proteins that are involved in many vital cellular functions such as regeneration and organogenesis [8910]. FRAT1 has been reported to be overexpressed in several tumors, such as non-small-cell lung cancer (NSCLC) [8], gliomas [11], ovarian cancers [12], esophageal squamous cell carcinomas [13], and colon cancer [14]. This data suggests that FRAT1 is associated with an aggressive phenotype. The relevance of the expression status of FRAT1 in breast cancer has not been determined, and its role in biological behavior is also unknown.
In the present study, we investigated the prognostic validity of FRAT1 expression and analyzed correlations with clinicopathological parameters and overall survival (OS) in breast cancer patients. These results may have a great contribution to the prediction of prognosis.
This study includes not only the breast cancer dataset from The Cancer Genome Atlas (TCGA) database (http://dna00.bio.kyutech.ac.jp/PrognoScan/index.html) [15], but we also collected 226 cases of invasive breast cancer at Konkuk University Medical Center between January 2009 and December 2011. Clinicopathological data such as age, tumor size, nodal stage, histologic grade, ER status, PR status, HER-2 amplification status, and CK5/6 expressions were obtained from the patients’ medical records and pathological reports. All patients were of the female gender, and their treatment was either modified radical mastectomy or breast-conserving surgery, and tumors were completely removed by surgery with pathological evaluation. After the completion of the surgery, adjuvant treatment with radiotherapy, chemotherapy, and hormone therapy was administered based on international guidelines. All the patients were staged according to the TNM stage, which was based on the criteria of American Joint Committee on Cancer, 8th edition. The patients were classified into 4 molecular subtypes, as follows: luminal A: ER-positive (ER+) and/or PR-positive (PR+), HER-2-negative (HER-2-) and low Ki67 level (<14.0%); luminal B: ER+ and/or PR+, Ki67 ≥14.0%, HER-2- or HER-2-positive (HER-2+); HER-2: ER-negative (ER-) and PR-negative (PR-), HER-2+; and TNBC: ER-, PR-, and HER-2-. We also analyzed the survival times of patients. The median follow-up time was 24 months.
The study was approved by the Institutional Review Board of Konkuk University Medical Center (No. KUH1210054).
TIMER (Tumor Immune Estimation Resource) 2.0 (http://timer.cistrome.org) is a comprehensive resource for systematical analysis of immune infiltrates across diverse cancer types [16]. In this study, TIMER 2.0 was used to analyze the expression of FRAT1 in breast cancer tissues and normal tissues across all TCGA datasets.
UALCAN (http://ualcan.path.uab.edu) is a website for the online analysis and mining of cancer data based mainly on the analysis of relevant cancer data in the TCGA database [17]. In this study, we used it to assess FRAT1 expression in breast cancer tissues and normal breast tissues.
GEPIA (Gene Expression Profiling Interactive Analysis) 2 (http://gepia2.cancer-pku.cn/#index) is a web server for analyzing the RNA sequencing expression data of tumors and normal samples from the TCGA and GTEx (Genotype-Tissue Expression) projects, using a standard processing pipeline [18]. The GEPIA box-plot tool was used for tumor/normal differential expression analysis for FRAT1 in breast cancers.
Four-micrometer sections were made from the formalin-fixed, paraffin-embedded blocks, which were then deparaffinized in xylene and rehydrated in a graded series of alcohol solutions. For antigen retrieval, paraffin tissue sections were cooked with 10-mM sodium citrate buffer, pH 6.0, at a sub-boiling temperature for 15 minutes and cooled for 20 minutes at room temperature. The sections were washed twice with Tris-buffered saline for 10 minutes. The antibody incubations were carried out at room temperature for 1 hour. FRAT1 was detected with standard avidin-biotin-peroxidase techniques using diaminobenzidine as the chromogen. Afterward, the slides were briefly counterstained with hematoxylin, dehydrated, and mounted. The immunostained samples were evaluated by 2 pathologists. The expression of FRAT1 was graded in terms of the percentage of cytoplasmic membrane staining in each block as well as the intensity of staining. However, there is no common cutoff value for the patterns of staining. In this study, the grade was designated according to the following criteria: negative, staining in less than 1% and positive, staining in more than 1% of tumor cells (Table 1).
ER and PR statuses were determined by immunohistochemistry (IHC) using the American Society of Clinical Oncology/College of American Pathologists (ASCO/CAP) guidelines. Tumors with >1% positively stained tumor cells were classified as positive for ER and PR. HER-2 status was also determined by IHC. HER-2 positive tumors were defined as 3+ and HER-2 negative tumors were defined as 0 or 1+, using the ASCO/CAP guidelines as the criteria for immunohistochemical staining. We performed silver in situ hybridization (SISH) to determine an accurate HER-2 status in tumors with equivocal immunoreactivity (2+).
We used PrognoScan (http://dna00.bio.kyutech.ac.jp/PrognoScan/index.html) to analyze the overexpressed FRAT1 and OS in breast cancer. The Kaplan-Meier plotter, and R2: Kaplan-Meier scanner with the log-rank test was used to examine the data. The cohorts from GSE (GSE9893) and TCGA were divided into 2 groups according to FRAT1 expression (high, >67%; low, ≤33%). All statistical tests were 2-sided.
The Breast Cancer Gene-Expression Miner ver. 4.8 (bc-GenExMiner 4.8, http://bcgenex.ico.unicancer.fr/BC-GEM/GEM-Accueil.php?js=1) is developed as a web-based tool including a MySQL relational database [19]. It was used to analyze the relationship between FRAT1 expression and clinicopathological parameters in breast cancer patients.
R software ver. 3.6.3 (R Foundation for Statistical Computing, Vienna, Austria) was used to perform the statistical analysis. Differences between groups were compared using the Wilcoxon rank-sum test or Student t-test, as appropriate. The chi-square test was used to analyze the possible correlations between FRAT1 expression and clinicopathologic factors. OS was defined as the time from surgery until the time of death. OS was estimated using the PrognoScan, Kaplan-Meier plotter, and Kaplan-Meier scanner. A P-value of <0.05 was considered statistically significant. IBM SPSS Statistics ver. 20.0 (IBM Corp, Armonk, NY, USA) was used for statistical analyses.
TIMER 2.0 was used to analyze the messenger RNA (mRNA) levels of FRAT1 in tumor tissues and normal tissues. The results showed that FRAT1 was overexpressed in breast cancer as compared to its expression levels in normal tissues, P < 0.05 (Fig. 1). To further evaluate the differential expression of FRAT1, we compared its expression level in the TCGA dataset using UALCAN. As shown in Fig. 2A, the expression of FRAT1 in breast cancer tissues is significantly higher than in normal tissue, P < 0.001 (Fig. 2A). In addition, it was discovered that the expression of FRAT1 mRNA in breast cancer tissues was significantly higher than that in normal tissues using GEPIA2 web analysis (P < 0.01) (Fig. 2B).
The median age was 56 years (range, 27–79 years). Among them, 66 patients underwent modified radical mastectomy, and 130 underwent breast-conserving surgery. The distribution of tumor and node stage was as follows: T1 in 104, T2 in 77, T3 in 12, N0 in 124, N1 in 39, N2 in 24, and N3 in 4 cases.
The expression of FRAT1 in breast cancer tissue was correlated with the patients’ clinicopathologic features. The clinicopathological characteristics of the patients according to FRAT1 positivity and negativity are shown in Table 2. In this study, FRAT1 was positively and markedly correlated with an increase in the TNM stages. FRAT1 positivity was correlated with tumor size (P = 0.004) and nodal status (P = 0.008). The FRAT1 positivity group had an association with a high histologic grade (P < 0.001) compared to the normal group. The expression of FRAT1 was shown to be increased with histological grade. Interestingly, CK5/6 is the most useful and important marker for the identification of the basal subgroup of TNBC [20]. Eighteen out of 21 cases with positive staining of CK5/6 showed FRAT1 expression, which is statistically significant (P < 0.001). Furthermore, the FRAT1 positivity group also had an association with ER negativity (P < 0.001) and PR negativity (P < 0.001). Although there was no statistically significant result, FRAT1 positivity was observed more in HER-2 negativity (P = 0.363). There was a significant difference in the subtype and the FRAT1 positivity group had more TNBC (11 out of 21). There was no significant correlation between FRAT1 expression and age and HER-2 status. Kaplan-Meier survival analysis demonstrated that the group with FRAT1 expression showed significantly worse OS than the one without FRAT1 expression (Fig. 4).
In the TCGA database, differential expression of FRAT1 in breast cancer was correlated with T stage, lymph node stage, and age (Table 3). In Fig. 5, we used bc-GenExMiner 4.8 to compare the transcription levels of FRAT1 among groups of patients according to different clinicopathological characteristics. The results showed that the overexpression of FRAT1 was significantly related to ER (IHC) status, PR (IHC) status, Her-2 (IHC) status, TNBC (IHC) status, basal-like (PAM 50) status, and Ki67 (IHC) status, P < 0.001 (Fig. 5).
Next, the prognostic value of FRAT1 mRNA expression in human breast cancer was further evaluated, using these online databases, including PrognoScan, and R2: Kaplan-Meier scanner. As shown in Fig. 6 (PrognoScan), we found that overexpression of FRAT1 was associated with poorer OS (P = 0.010) (Fig. 6). The result of R2: Kaplan-Meier scanner showed that overexpression of FRAT1 was associated with a poor prognosis for patients with breast cancer (P = 0.025) (Fig. 7). The results revealed that overexpression of FRAT1 was significantly associated with poor prognosis. All these findings indicate that FRAT1 may serve as a poor predictor of the prognosis of patients with breast cancer.
TNBC is a heterogeneous and deadly disease in women with limited therapy choices, necessitating the development of novel targeted therapies [21]. The FRAT1 mapped to human chromosome 10q24.1, which is a human homolog of mouse proto-oncogene FRAT1 [7]. Proto-oncogene FRAT1 was first identified as a T-cell lymphoma proto-oncogene that encodes a 29KD protein. FRAT1 regulates Wnt/β-catenin signaling through inhibiting the acidification of GSK3 on β-catenin and translocation into the nucleus, where it binds to the TCF/LEF signal transduction pathway [1322]. Dysregulation of Wnt/β-catenin signaling has been thought to play an important role in different kinds of cancers. Extensive studies have shown that the Wnt signal transduction pathway is associated with different kinds of human cancers, such as colon cancer, acute leukemia, lung cancer, and glioma [1223]. FRAT1 has been identified as a component of the Wnt signaling pathway in previous research [2425]. FRAT1 is a positive regulator of the Wnt/β-catenin pathway [12]. Overexpression of FRAT1 leads to β-catenin stabilization through dissociation of GSK-3β from Axin and inhibition of β-catenin phosphorylation [2526]. High FRAT1 expression is associated with increased malignancy and a higher clinical stage [27]. FRAT1 was associated with invasion, metastasis, and malignant phenotypes. In this study, we observed FRAT1 overexpression in breast cancer. FRAT1 positivity was correlated with tumor size, nodal status, histologic grade, and CK5/6. Furthermore, the FRAT1 positivity group also had an association with ER negativity and PR negativity. There was a significant difference in the subtype and the FRAT1 positivity group had more TNBC. TNBC is one of the subtypes of malignant breast cancer with the worst prognosis [28]. Most patients have been diagnosed at an advanced stage or experience a recurrence within a short period of time. Thus, a special approach to treatment and strategy are needed. Importantly, FRAT1 positivity is correlated with poor OS in breast cancer patients. Furthermore, we assessed the expression pattern of FRAT1 in human breast cancer by bioinformatics validation. Then, the prognostic significance of FRAT1 expression in breast cancer was also analyzed. Higher expression of FRAT1 indicated poorer OS in breast cancer in different kinds of online databases. Furthermore, previous studies on different types of malignancies have reported the role of high expression of FRAT1 protein but not in breast cancer. Interestingly, previous studies demonstrated that Wnt7a protein expression was significantly correlated with ER expression, and loss of Wnt7a expression was associated with poor disease-free survival in breast cancer patients [29]. Wnt7a is a member of the Wnt family, which is regulated by FRAT1. The exact function that FRAT1 plays in the pathway in controlling carcinogenesis is unknown. Therefore, further research will be needed. In conclusion, FRAT1 protein is abnormally expressed in breast cancer. FRAT1 might be involved in the malignancy progression of breast cancer. High expression of FRAT1 is correlated with clinical and pathological stages and poorer OS. Therefore, positive FRAT1 expression is a biomarker for cancer severity and a predictive marker of poor prognosis in breast cancer.
Notes
References
1. Althuis MD, Dozier JM, Anderson WF, Devesa SS, Brinton LA. Global trends in breast cancer incidence and mortality 1973-1997. Int J Epidemiol. 2005; 34:405–412. PMID: 15737977.
2. Schnitt SJ. Classification and prognosis of invasive breast cancer: from morphology to molecular taxonomy. Mod Pathol. 2010; 23 Suppl 2:S60–S64. PMID: 20436504.
4. Liao C, Zhang Y, Fan C, Herring LE, Liu J, Locasale JW, et al. Identification of BBOX1 as a therapeutic target in triple-negative breast cancer. Cancer Discov. 2020; 10:1706–1721. PMID: 32690540.
5. Anders CK, Carey LA. Biology, metastatic patterns, and treatment of patients with triple-negative breast cancer. Clin Breast Cancer. 2009; 9(Suppl 2):S73–S81. PMID: 19596646.
6. Metzger-Filho O, Tutt A, de Azambuja E, Saini KS, Viale G, Loi S, et al. Dissecting the heterogeneity of triple-negative breast cancer. J Clin Oncol. 2012; 30:1879–1887. PMID: 22454417.
7. Saitoh T, Katoh M. FRAT1 and FRAT2, clustered in human chromosome 10q24.1 region, are up-regulated in gastric cancer. Int J Oncol. 2001; 19:311–315. PMID: 11445844.
8. Zhang Y, Yu JH, Lin XY, Miao Y, Han Y, Fan CF, et al. Overexpression of Frat1 correlates with malignant phenotype and advanced stage in human non-small cell lung cancer. Virchows Arch. 2011; 459:255–263. PMID: 21818639.
9. Hagen T, Cross DA, Culbert AA, West A, Frame S, Morrice N, et al. FRAT1, a substrate-specific regulator of glycogen synthase kinase-3 activity, is a cellular substrate of protein kinase A. J Biol Chem. 2006; 281:35021–35029. PMID: 16982607.
10. Krishnamurthy N, Kurzrock R. Targeting the Wnt/beta-catenin pathway in cancer: update on effectors and inhibitors. Cancer Treat Rev. 2018; 62:50–60. PMID: 29169144.
11. Guo G, Mao X, Wang P, Liu B, Zhang X, Jiang X, et al. The expression profile of FRAT1 in human gliomas. Brain Res. 2010; 1320:152–158. PMID: 20096670.
12. Wang Y, Hewitt SM, Liu S, Zhou X, Zhu H, Zhou C, et al. Tissue microarray analysis of human FRAT1 expression and its correlation with the subcellular localisation of beta-catenin in ovarian tumours. Br J Cancer. 2006; 94:686–691. PMID: 16479254.
13. Wang Y, Liu S, Zhu H, Zhang W, Zhang G, Zhou X, et al. FRAT1 overexpression leads to aberrant activation of beta-catenin/TCF pathway in esophageal squamous cell carcinoma. Int J Cancer. 2008; 123:561–568. PMID: 18498136.
14. Zhu K, Guo J, Wang H, Yu W. FRAT1 expression regulates proliferation in colon cancer cells. Oncol Lett. 2016; 12:4761–4766. PMID: 28101222.
15. Tomczak K, Czerwińska P, Wiznerowicz M. The Cancer Genome Atlas (TCGA): an immeasurable source of knowledge. Contemp Oncol (Pozn). 2015; 19(1A):A68–A77. PMID: 25691825.
16. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020; 48(W1):W509–W514. PMID: 32442275.
17. Chandrashekar DS, Bashel B, Balasubramanya SA, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BV, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017; 19:649–658. PMID: 28732212.
18. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017; 45(W1):W98–W102. PMID: 28407145.
19. Jézéquel P, Campone M, Gouraud W, Guérin-Charbonnel C, Leux C, Ricolleau G, et al. bc-GenExMiner: an easy-to-use online platform for gene prognostic analyses in breast cancer. Breast Cancer Res Treat. 2012; 131:765–775. PMID: 21452023.
20. Hashmi AA, Naz S, Hashmi SK, Hussain ZF, Irfan M, Bakar SM, et al. Cytokeratin 5/6 and cytokeratin 8/18 expression in triple negative breast cancers: clinicopathologic significance in South-Asian population. BMC Res Notes. 2018; 11:372. PMID: 29884220.
21. Bianchini G, Balko JM, Mayer IA, Sanders ME, Gianni L. Triple-negative breast cancer: challenges and opportunities of a heterogeneous disease. Nat Rev Clin Oncol. 2016; 13:674–690. PMID: 27184417.
22. Saitoh T, Mine T, Katoh M. Molecular cloning and expression of proto-oncogene FRAT1 in human cancer. Int J Oncol. 2002; 20:785–789. PMID: 11894125.
23. Zhang Y, Han Y, Zheng R, Yu JH, Miao Y, Wang L, et al. Expression of Frat1 correlates with expression of β-catenin and is associated with a poor clinical outcome in human SCC and AC. Tumour Biol. 2012; 33:1437–1444. PMID: 22528942.
24. Yost C, Farr GH 3rd, Pierce SB, Ferkey DM, Chen MM, Kimelman D. GBP, an inhibitor of GSK-3, is implicated in Xenopus development and oncogenesis. Cell. 1998; 93:1031–1041. PMID: 9635432.
25. Thomas GM, Frame S, Goedert M, Nathke I, Polakis P, Cohen P. A GSK3-binding peptide from FRAT1 selectively inhibits the GSK3-catalysed phosphorylation of axin and beta-catenin. FEBS Lett. 1999; 458:247–251. PMID: 10481074.
26. Li L, Yuan H, Weaver CD, Mao J, Farr GH 3rd, Sussman DJ, et al. Axin and Frat1 interact with dvl and GSK, bridging Dvl to GSK in Wnt-mediated regulation of LEF-1. EMBO J. 1999; 18:4233–4240. PMID: 10428961.
27. Walf-Vorderwülbecke V, de Boer J, Horton SJ, van Amerongen R, Proost N, Berns A, et al. Frat2 mediates the oncogenic activation of Rac by MLL fusions. Blood. 2012; 120:4819–4828. PMID: 23074275.
28. Jiang Y, Liu Y, Tan X, Yu S, Luo J. TPX2 as a novel prognostic indicator and promising therapeutic target in triple-negative breast cancer. Clin Breast Cancer. 2019; 19:450–455. PMID: 31494045.
29. Yi K, Min KW, Wi YC, Kim Y, Shin SJ, Chung MS, et al. Wnt7a deficiency could predict worse disease-free and overall survival in estrogen receptor-positive breast cancer. J Breast Cancer. 2017; 20:361–367. PMID: 29285041.