Prediction and Prognostic Significance of BCAR3 Expression in Patients
Total Page:16
File Type:pdf, Size:1020Kb
Zhang et al. J Transl Med (2018) 16:363 https://doi.org/10.1186/s12967-018-1728-8 Journal of Translational Medicine RESEARCH Open Access Prediction and prognostic signifcance of BCAR3 expression in patients with multiple myeloma Weilong Zhang1†, Yuansheng Lin2†, Xiaoni Liu3†, Xue He4, Ye Zhang5, Wei Fu6, Zuozhen Yang1, Ping Yang1, Jing Wang1, Kai Hu1, Xiuru Zhang4, Weiyou Liu3*, Xiaoliang Yuan3* and Hongmei Jing1* Abstract Background: Multiple myeloma (MM) is the plasma cell tumor, which is characterized by clonal proliferation of tumor cells, with high risk of progression to renal impairment, bone damage and amyloidosis. Although the survival rate of patients with MM has improved in the past decade, most people inevitably relapse. The treatment and progno- sis of MM are still urgent problems. Breast Cancer Antiestrogen Resistance 3 (BCAR3) is a protein-coding gene that is associated with many tumors. However, there have been few studies on the relationship of BCAR3 and MM. Methods: We analyzed 1878 MM patients (1930 samples) from 7 independent datasets. First, we compared the BCAR3 expression level of MM patients in diferent stages and MM patients with diferent amplifcation of 1q21. Sec- ond, we analyzed BCAR3 expression levels in MM patients with diferent molecular subtypes. Finally, we explored the event-free survival rate (EFS) and overall survival rate (OS) of MM patients with high or low BCAR3 expression, includ- ing patients before and after relapse, and their therapeutic responses to bortezomib and dexamethasone. Results: The expression of BCAR3 showed a decreasing trend in stages I, II and III (P 0.00068). With the increase of 1q21 amplifcation level, the expression of BCAR3 decreased (P 0.022). Patients with= high BCAR3 expression had higher EFS and OS (EFS: P < 0.0001, OS: P < 0.0001). The expression= of BCAR3 gene before relapse was higher than that after relapse (P 0.0045). BCAR3 is an independent factor afecting prognosis (EFS: P 5.17E 03; OS: P 3.33E 04). = = − = − Conclusion: We found that high expression level of BCAR3 predicted better prognosis of MM patients. Low expres- sion of BCAR3 at diagnosis can predict early relapse. BCAR3 is an independent prognostic factor for MM. BCAR3 can be used as a potential biomarker. Keywords: BCAR3, Multiple myeloma, Prognosis, Gene expression profle Background heterogeneous disease with diferent clinical characteris- MM is a B cell diferentiated tumor characterized tics [4]. By recognizing genetic mechanism and mutation, by clonal proliferation of tumor cells [1–3]. MM is a a normal plasma cell transited into the following disease stages: monoclonal gammopathy of undetermined sig- nifcance, smouldering myeloma, myeloma and plasma *Correspondence: [email protected]; [email protected]; [email protected] cell leukaemia [5]. Te International Staging System (ISS) †Weilong Zhang, Yuansheng Lin and Xiaoni Liu contributed equally to uses the staging criteria to divide MM into three phases, this work combining serum albumin levels with β2-microglobulin 1 Department of Hematology, Lymphoma Research Center, Peking University Third Hospital, No. 49 North Garden Road, Haidian District, to determine the prognosis of MM patients [6]. Revised Beijing 100191, People’s Republic of China International Staging System (R-ISS) is a simple and 3 Department of Respiratory Medicine, The First Afliated Hospital efective prognostic staging system that combines of Gannan Medical University, No. 23 Qingnian Road, Zhanggong District, Ganzhou 341000, People’s Republic of China three prognostic tools (ISS: International Staging Sys- Full list of author information is available at the end of the article tem, CA: chromosomal abnormalities, LDH: lactate © The Author(s) 2018. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/ publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Zhang et al. J Transl Med (2018) 16:363 Page 2 of 9 dehydrogenase) to better assess patient prognosis [7]. criteria for patient selection in our research were stated. Gene expression profling (GEP) is important for reveal- 1) All MM patients with the published high throughput ing MM molecular heterogeneity of diferent patients. gene expression data. 2) All the patients should have Te molecular classifcation established by UAMS (Te some information such as clinical features, biochemical University of Arkansas for Medical Sciences) based on examination, karyotype, therapy or therapy response. GEP data [8]. Te molecular basis of MM is defned by Te study was approved by the Human Research Eth- unsupervised clustering analysis of mRNA expression ics Committee of Peking University third hospital. Te profles, which are divided into seven molecular subtypes research was conducted in accordance International [9]. Recurrence is a major problem in most MM patients. Conference on and the Declaration of Helsinki. Among all MM patients, early recurrence and poor prog- nosis are about 10% to 15% [10]. 1q21 amplifcation is the most common chromosomal Microarray analysis aberration in MM and is considered a high-risk genetic All microarray data were analyzed, and the signif- feature [11]. Te function of the BCAR3 (Breast Cancer cantly abnormal expression genes were systemati- Antiestrogen Resistance 3) gene is varied. BCAR3 can cally screened as predictive biomarkers. Te diferent promote cell migration, proliferation, and BCAR3 was expression of genes between BCAR3-low to BCAR3- identifed as a molecular linkage between PTP and Cas high group were also analyzed and ranked by fold- [12–14]. BCAR3 protein participates in the signaling change values (log2, FC > 0.8 or < − 0.8, P < 0.05). pathway of EGF through its SH2 domain, leading to cell We retrieved GSE24080 (559 samples) from the NCBI cycle progression, and BCAR3 itself belong to a mitogenic GEO database and analyzed the expression of BCAR3 signaling pathway [15, 16]. BCAR3 and p130Cas were among diferent ISS stages, diferent molecular sub- associated with anti-estrogen in breast cancer, Rac activa- types, diferent 1q21 amplifcation levels, and difer- tion. BCAR3 can regulate the Src signal transmission of ent survival (EFS and OS). We retrieved GSE82307 (66 BCAR3-p130 (cas) complex dependence mode [17–19]. samples) from the NCBI GEO database. We analyzed During the initial stage of gonadal development, BCAR3 the expression of BCAR3 before and after recurrence in in gonad development is very important and is expressed the same patient. We retrieved GSE19554 (38 samples) in XY gonads [20]. Studies have shown that the BCAR3 in 19 MM patients from the NCBI GEO database, and gene is the frst spontaneous mutation associated with we analyzed BCAR3 expression between the baseline cataracts caused by lens compression. Tis new cataract (before chemotherapy) and pre-1st (after induction of model could provide further knowledge about the func- chemotherapy, before bone marrow transplantation) in tion of BCAR3 protein [21]. In the study of tamoxifen in the same sample. We retrieved GSE19784 (308 samples) the treatment of metastatic breast cancer, high expres- from the NCBI GEO database and analyzed the expres- sion of BCAR3 is related to good progression-free sur- sion of BCAR3 in diferent molecular subtypes. We vival, and the expression level of BCAR3 in primary breast retrieved GSE83503 (585 samples) from the NCBI GEO tumors is relatively low, which is related to the survival database and analyzed the expression of BCAR3 in the rate of distant metastasis [14]. Ovarian cancer is a dis- recurrence and non-recurrence groups. We retrieved ease characterized by tumor heterogeneity, which is dif- GSE9782 (238 samples) from the NCBI GEO database cult to be diagnosed and treated. Studies have shown that and analyzed the expression of BCAR3 between the inhibiting the expression of BCAR3 gene can inhibit the diferent treatment response of bortezomib and dexa- cell proliferation of ovarian cancer [22]. However, there is methasone. All patients from GSE39754 (136 samples) no study reporting the relationship between BCAR3 and were treated with vincristine, adriamycin, and dexa- MM so far. By integrating data of 1878 MM patients, we methasone (VAD), followed by autologous stem cell found that BCAR3 gene is closely related to MM. transplantation (ASCT). Ten measuring “therapeutic response” after ASCT. Analysis of BCAR3 expression in Methods each treatment response compared to the average of all Data source treatment responses. Terapeutic response: complete In our study, gene expression microarrays of 1878 response (CR), very good partial response (VGPR), par- MM patients (1930 samples) were derived from Gene tial response (PR), no response, stable disease (NR), no Expression Omnibus database, including datasets response, progressive disease (Prog). GSE24080 (559 samples) [23], GSE82307 (66 samples) Gene ontology (GO) analysis [24], GSE19784 (308 samples) [25], GSE83503 (585 samples) [26], GSE9782 (238 samples) [27], GSE39754 Use the DAVID to analyze the 559 samples (dataset (136 samples) and GSE19554 (38 samples) [28, 29]. Te GSE24080), and fnd out the enrichment pathways for Zhang et al. J Transl Med (2018) 16:363 Page 3 of 9 diferent expressed genes between BCAR3-low and Result BCAR3-high group [30]. Te results were ranked by the BCAR3 is expressed in multiple myeloma of diferent molecular types P value (− log10, P < 0.05). We compared BCAR3 gene expression level of diferent 1q21 amplifcation. With the increase of the amplifcation Statistical analyses level of 1q21, the expression level of BCAR3 showed an Statistical analyses were performed by R software v3.1.3 overall downward trend (Fig. 1a, P = 0.022, Kruskal–Wal- (ggplot2 and survminer package). Te log-rank test and lis test).