m+Medelevator.RU
← Поиск публикаций
Journal of Cancer · 2020

Clinical prognostic implications of EPB41L4A expression in multiple myeloma.

Название и аннотация: Перевод в очереди
Обработка статьи · этапы, ошибки и управление
Загрузка состояния…

История · было → стало

Состояние обновляется каждые 5 секунд, пока блок открыт. Процент показывается только там, где известен объём работы. Пауза останавливает продолжение этапа; текущий запрос к источнику может завершиться.

Zhang W, Lai R, He X, Liu X, Zhang Y, Yang Z, Yang P, Wang J, Hu K, Yuan X, Zhang X, Liu W, Jing H

Аннотация

Background: Multiple myeloma (MM) is one of the most common incurable malignancies in malignant plasma cell disease. EPB41L4A is a target gene for the Wnt/β-catenin pathway, which is closely related to the survival of multiple myeloma cells. However, there is currently no research report on the prognostic significance of the EPB41L4A gene in MM. Methods: We studied the biological significance and prognostic significance of EPB41L4A expression in MM by integrating 1956 MM samples from 7 datasets, and explored the relationship between EPB41L4A expression and MM ISS stage, molecular type, therapeutic response and survival. Results: We found that the expression level of EPB41L4A is inversely proportional to the copy number of 1q21 (P = 3.4e-13). EPB41L4A was low expressed in MAF, MMSET and proliferating molecular typing patients (P Conclusions: We found that the high expression of EPB41L4A predicts good survival level in MM.

Перевод полного текста: Перевод в очереди

Читать найденный полный текст

https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/4.0/ ). See http://ivyspring.com/terms for full terms and conditions.

Background

Multiple myeloma (MM) is one of the most common incurable malignancies in malignant plasma cell disease. 1 It is characterized by malignant proliferation of monoclonal plasma cells in the bone marrow and secretion of a large number of monoclonal immunoglobulin. 2 The main clinical manifestations were infection, anemia, immunosuppression, bone destruction and renal failure. 3 In recent years, the incidence of MM has gradually increased; accounting for about 1% of all cancers, and about 10% of hematopoietic tumors. 4 The median overall survival can be several months or more than ten years. 5 At present, the specific mechanism of plasma cell disorders and progression to symptomatic MM disease has not been fully elucidated. However, there has been some progress in the study of the staging of MM disease, which can predict the prognosis of patients. In the early years, staging of tumors was performed by using the Durie-Salmon staging system. 6 In recent years, the International Myeloma Working Group (IMWG) has proposed an international staging system (ISS) for staging multiple myeloma with β2-MG and serum albumin levels. 7 It helps to predict the prognosis of the disease. With the development of molecular diagnostic techniques, the IMWG proposed the revised International Classification System (R-ISS) in the original ISS staging system in 2015. 8 , 9 The treatment options of MM vary, but at present it is chemotherapy-based and autologous stem cell transplantation (ASCT). 10 , 11 At present, some studies have shown that some genes can predict the clinical outcome of MM patients. For example, 1q21 amplification is the most common chromosomal mutation in MM patients and indicates a poor prognosis. 12 And two GEP scores for MM patient classification and prognosis (based on the Arkansas GEP70 score and the HOVON SYK92 score) have been commercialized. We can still find some "new genes" that can also predict the clinical outcome of MM patients. The "new gene" can be integrated with genes that have been found to predict the prognosis of MM patients to better predict the prognosis of MM patients, and can also provide a reference for opening up a new therapeutic pathway.

EPB41L4A ( e rythrocyte p rotein b and 4.1-l ike 4a, also called Nbl4) which is a member of the band FERM (Four-point-one, Ezrin, Radixin, Moesin) protein superfamily were discovered in almost all cellular organism in molecular ecology, especially high expressions in brain, liver, thymus, and peripheral blood leukocytes. 13 The protein superfamily members are important regulators of embryonic development. They are a group of membrane-associated proteins whose major biological functions are involved in the regulation of cytoskeletal rearrangements, intracellular trafficking and signal transduction. Studies have suggested that mutations in several genes in the FERM protein family are associated with human cancers and blood cell diseases. 14 Studies have shown that EPB41L4A is a target gene for the Wnt/β-catenin pathway and is associated with cell polarity or proliferation. (13)Beta-catenin is a key regulator of cell proliferation and is frequently mutated in various human cancer types. Recent advances in cancer research have revealed that β-catenin plays a role in intercellular adhesion and Wnt signaling pathways and is also a major player in carcinogenesis in various tissues (colon, liver, ovary, skin, and blood); 15 - 18 the Wnt signaling pathway plays an important role in regulating cell proliferation, migration, differentiation, and stem cell self-renewal, which is closely related to the survival of multiple myeloma cells. 19

Only a few studies have shown that EPB41L4A has genetic susceptibility to colon, nerve/brain muscle and congenital keratosis. 20 However, there is currently no clinical prognostic study on the EPB41L4A gene in any types of cancer. The study integrates multiple data from MM to achieve three goals. First, determine the clinical relationship between EPB41L4A and MM. Second, to understand whether EPB41L4A acts on MM is related to the cell cycle. Finally, understand the relationship between this gene and the recurrence before and after treatment.

Methods

Data source

In our study, gene expression microarrays of 1956 MM samples were derived from Gene Expression Omnibus (GEO) database. We selected patients in our research with the criteria as follows. 1) All MM patients should have the whole transcriptome gene expression data. 2) All the patients should include some clinical information such as clinical biochemical examination, therapy regimen or therapy response. GSE24080 (559 samples from 559 patients), 21 the relationship between EPB41L4A expression and ISS stage, serotype, survival (EFS and OS), 1q21 amplification, molecular subtype, and pathway was analyzed. GSE9782 (264 samples from 264 patients), 22 the relationship between EPB41L4A expression and survival (OS) was analyzed. GSE82307 (66 samples from 33 patients), 23 we analyzed the relationship between EPB41L4A expression before and after relapse or treatment. GSE19554 (38 samples from 19 patients), 24 we analyzed the relationship between EPB41L4A expression before and after treatment. GSE19784 (308 samples from 308 patients), 25 the relationship between the expression of EPB41L4A and the 9 molecular types was analyzed. GSE83503 (585 samples from 585 patients), 26 the relationship between expression of EPB41L4A and relapse was analyzed. GSE9782 (238 samples from 238 patients), 22 the relationship between the expression of EPB41L4A and the drug treatment response (dexamethasone and bortezomib) was analyzed using the U133A and U133B array test methods, respectively. GSE39754 (136 samples from 136 patients), 27 the relationship between the expression of EPB41L4A and ASCT after drug treatment induction therapeutic response was analyzed. Our research is in line with the Helsinki Declaration. And our research was approved by the Ethics Committee of Peking University Third Hospital and the Ethics Committee of Gannan Medical University. These patients have signed informed consent and obtained approval from the Ethics Committee. 21 - 27

Microarray analysis

Gene expression of each dataset was calculated by RMA method (robust multiarray averaging). The relative RNA expression value of each gene was log-transformed using log2. We analyzed expression of each gene and the survival (EFS and OS) using survcomp package (hazard.ratio function) with the GSE24080 datasest. Hazard ratio and P-value was estimated through Cox regression model. EPB41L4A gene is in the high rank among all genes by ranking with the P-value.

The cut point of EPB41L4A gene expression was conducted using survminer package (surv_cutpoint function) with maximally selected rank statistics method. The MM patients with higher expression than this cut point constituted the EPB41L4A high group, and the MM patients with expression lower than this cut point constituted the EPB41L4A low group. The gene expression profiles of EPB41L4A high group and EPB41L4A low group were analyzed. P < 0.05 in unpaired t test and Foldchange (FC, log2) > 0.8 or < -0.8 was used to indicate expression levels of different genes.

Gene Ontology (GO) analysis

By using the default parameters of the DAVID tool, 28 different expressed genes between EPB41L4A high group and EPB41L4A low group was analyzed of the pathway enrichment. The results were ranked by the P value (-log10).

Statistics

This study performed statistical analysis by R software v3.1.3 (ggplot2 and survminer package). More than two sets of samples were compared using the Kruskal-Wallis test. The log-rank test, Cox regression analysis and Fisher's exact test are used for survival analysis. The average values of two samples were compared using the unpaired t test or Wilcoxom test. The average values of more than two samples were compared using the Anova.

Result

Low expression of EPB41L4A in MAF, MMSET and proliferating molecular typing

We compared the expression levels of EPB41L4A at different amplification levels of 1q21. This result shows that the higher the amplification levels of 1q21, the lower the expression level of EPB41L4A gene. There was a significant difference in the expression level of the gene between each amplification level (dataset GSE24080 , total 559 samples, P = 3.4e-13, Kruskal-Wallis Test, Figure 1 A). The comparison of the expression level of EPB41L4A in all 7 molecular subtypes with the average value of all molecular subtypes showed that the gene was highly expressed in the hyper diploid type with statistical significance (dataset GSE24080 , total 559 samples, P <= 0.0001, unpaired t test, two sided, Figure 1 B). The expression level of EPB41L4A was significantly lower and statistically significant in patients with MAF, MMSET and proliferating molecular types (dataset GSE24080 , total 559 samples, P <= 0.001, unpaired t test, two sided, Figure 1 B). However, comparing the expression levels of EPB41L4A in CD1, CD2 and LB molecular types and the average level of all 7 molecular types was not statistically significant (dataset GSE24080 , total 559 samples, P > 0.05, unpaired t test, two sided, Figure 1 B). We also compared the expression levels of the 9 molecular types (the other method of molecular typing) with the average of all molecular types. The data was derived from another dataset GSE19784 , a total of 308 MM specimens. The expression level of EPB41L4A in 9 molecular typing patients is approximately the same as that in 7 molecular typing (P = 2.4e-07, Anova, Figure S1 and Table S1 ). Similarly, the expression level of EPB41L4A was decreased in proliferative molecular type, although there was not statistical significant (P > 0.05, unpaired t test, two sided, Figure S1 ).

To compare the expression level of EPB41L4A gene in different stages of MM, we analyze 559 MM expression profiles from the GSE24080 dataset. There was a significant difference in the expression of EPB41L4A gene between ISS I and III phase of MM (P = 0.019, Wilcoxon Test, Figure S2 A). However, comparing this gene expression level between ISS I and ISS II Phase was not statistically significant (P = 0.46, Wilcoxon Test, Figure S2 A), and the comparison of gene expression levels between ISS II and ISS III was also not statistically significant (P = 0.15, Wilcoxon Test, Figure S2 A). We also compared the expression levels of EPB41L4A at different stages of different serotypes. In the serotype IgA group, we found that the expression level of EPB41L4A was lower in the ISS III than in the first two stages, and the difference in the expression level of EPB41L4A in ISS II and ISS III was statistically significant. (P = 0.045, Wilcoxon Test, Figure S2 B). However, there was no significant difference in the expression levels of this gene between ISS I and ISS II (P = 0.46, Wilcoxon Test, Figure S2 B), and similarly no significant difference between ISS I and ISS III (P = 0.084, Wilcoxon Test, Figure S2 B). The expression of EPB41L4A gene in serum-free light chain (FLC) and serotype IgG groups was not statistically significant at all Phases. (FLC: P = 0.8, IgG: P = 0.15, Kruskal-Wallis Test, Figure S2 B)

EPB41L4A expression predicts the survival level in MM

Since EPB41L4A is low expressed in proliferative patients, we further analyzed the relationship between the expression level of EPB41L4A and survival. We compared event-free survival (EFS) and overall survival (OS) in 559 MM patients (dataset GSE24080 ) with the EPB41L4A-high group and the EPB41L4A-low group. We found that the high expression of EPB41L4A predicts significantly good EFS (P < 0.0001, log-rank test, Figure 2 A) and OS (P < 0.0001, log-rank test, Figure 2 A). We also analyzed the relationship between EPB41L4A expression level and survival through an independent dataset (dataset GSE9782 , total 264 samples). Similarly, we found that EPB41L4A-high group had better OS (P < 0.0001, log-rank test, Figure S4 ) than EPB41L4A-low group. In addition, we also further analyzed the expression of this gene in the early (ISS I) and mid-late (ISS II and III) of MM patients with OS and EFS. The results showed that the expression level of EPB41L4A is closely related to OS and EFS in various stages of the disease. Especially in the early stage of the disease, the expression level of the gene has more significant difference in the survival period (ISS I: EFS P = 1e-04, OS P = 0.00036; ISS II and III EFS P = 0.00048, OS P = 0.0064; log-rank test; Figure 2 B and Figure S5 ).

EPB41L4A expression as an independent prognostic factor in MM patient

We found that EPB41L4A was an independent prognostic factor in 559 MM patients (EPS: P = 8.80e-06, OS P = 7.54e-05, Cox regression analysis, Table S2 ). The EFS hazard ratio for EFS41L4A was 0.50 (95% CI: 0.37 -0.68, Cox regression analysis, Table S2 ) and the OS risk ratio for EFS41L4A was 0.49 (95% CI: 0.34-0.70, Cox regression analysis, Table S2 ). The Table S2 shows that EPB41L4A is a prognostic factor independent of β-2 microglobulin (B2M), LDH, albumin 21 , MRI and bone marrow biopsy plasma cells (BMPC). In addition, to ensure that the comparison between the EPB41L4A-high group and the EPB41L4A-low group is significant and reduces bias, we excluded confounding factors. The basic clinical features of SEX, RACE and ISOTYPE of MM patients enrolled in these two groups were essentially the same (P > 0.05, Fisher's exact test, Table S3 ). In addition, baseline characteristics of B2M, CRP, CREAT, ALB, ASPC, BMPC, and MRI were also consistent (P > 0.05; unpaired t-test, bilateral; Table S3 ). However, the clinical features of LDH and hemoglobin (HGB) were inconsistent between the two groups of MM patients (LDH: P < 0.001; HGB: P = 0.025; unpaired t test, two sided; Table S3 ). We also analyzed the clinical features of MM patients based on EPB41L4A expression levels in other datasets, except for ISOTYPE, ALB (dataset GSE9782 , P <0.001, Table S4 ), other clinical features (SEX, RACE, Age, B2M, treatment and treatment response) were essentially identical (P > 0.05, Table S4 ).

Both passive regulation of cell division and DNA replication were meaningfully enriched pathways in EPB41L4A-high group MM

We compared the gene expression profiles of the EPB41L4A-high group and the EPB41L4A-low group MM (Figure 3 A). In conclusion, 227 up-regulated genes and 233 down-regulated genes were observed between EPB41L4A-high group and EPB41L4A-low group in MM. Figure 4 shows top 12 up-regulated genes and top 12 down-regulated genes. We found that ISL2 and CCND1 gene expression was up-regulated in EPB41L4A-high group (dataset GSE24080 , total 559 samples, ISL2: P < 0.001; CCND1 P < 0.001; Figure 3 A). However, FGFR3 and CCND2 are down-regulated genes, that is the expression of FGFR3 and CCND2 is elevated in patients with low expression of EPB41L4A (dataset GSE24080 , total 559 samples, FGFR3: P < 0.001; CCND2 P < 0.001; Figure 3 A). It is further shown that high expression of EPB41L4A is beneficial to prognosis. Through Gene Ontology (GO) analysis, we found that the top 15 pathway is related to the cell cycle, and the most significant enrichment pathway through all the expressed genes is cell division, mitotic nuclear division and cell adhesion (dataset GSE24080 , total 559 samples, P < 0.0001, Figure 3 B). We were surprised to find that in the DNA replication pathway all of the different gene expressions are down-regulated (dataset GSE24080 , total 559 samples, P < 0.0001, Figure 4 ).

The relationship between EPB41L4A and MM recurrence

Although current high-dose therapy, autologous stem cell transplantation, chemotherapy and supportive care can significantly improve MM survival, almost all MM patients will eventually relapse, so we explore whether this gene is associated with recurrence. We analyzed 585 MM expression profiles from the GSE83503 dataset. We found that patients with relapsed MM had lower expression levels of EPB41L4A than those without recurrence (P = 0.0039, unpaired t test, two sided, Figure S3 and Table S5 ). We compared the expression of EPB41L4A before and after recurrence in 33 pairs of MM patient samples from the GSE82307 dataset. We found no difference in the expression of EPB41L4A before and after relapse (P = 0.13, Wilcoxon Test, Figure 5 A and Table S6 ). We also compared EPB41L4A expression before and after treatment in 19 pairs of MM patient samples from the GSE19554 dataset. We found no difference in the expression of EPB41L4A before treatment and pre-1st (P = 0.7, Wilcoxon Test, Figure 5 B).

The relationship between EPB41L4A expression and therapeutic response

To further understand the relationship between EPB41L4A expression and therapeutic response, we compared the relationship between EPB41L4A expression and dexamethasone (Dex) and bortezomib treatment responses, respectively. We found that EPB41L4A expression was elevated and statistically significant in dexamethasone-treated partial response (PR) patients (dataset GSE9782 , total 238 samples, P <= 0.05, Anova, Figure S6 ). However, the expression of EPB41L4A was not statistically significant in the other treatment responses of Dex and bortezomib (dataset GSE9782 , 238 samples, P > 0.05, Anova, Figure S6 ). We also compared the relationship between the expression of EPB41L4A and the therapeutic response of autologous stem cell transplantation (ASCT) after three drug-induced treatments (vincristine, doxorubicin, and dexamethasone (VAD)). We found that EPB41L4A expression was not statistically significant in all of its therapeutic responses (dataset GSE39754 , 136 samples, P > 0.05, Anova, Figure S7 ).

Discussion

As the understanding of the pathogenesis of MM patients continues to improve, 29 the heterogeneity of the disease becomes more and more obvious. These "high-risk" patients often have one or more characteristics of cytogenetic abnormalities, but these abnormalities do not always lead to poor prognosis. 30 EPB41L4A is a key target of Wnt signaling pathway, 13 and Wnt signaling pathway can regulate cell proliferation and differentiation, 19 and is related to renal cell carcinoma complicated with complex of colon, nerve/brain muscle, congenital keratoses and tuberous sclerosis. 31 It is closely related to the survival of MM patients. 19 At present, there is no research report on prognostic significance and biological implication of EPB41L4A gene in MM. Therefore, by integrating several datasets, we found that the expression level of EPB41L4A is inversely proportional to the MM stage, and the prognosis of MM patients with high expression of EPB41L4A is better.

A number of previous studies have shown that the long non-coding RNA EPB41L4A-AS1 and the long non-coding RNA EPB41L4A-AS2 are associated with a variety of cancers. For example: 1) EPB41L4A-AS1 is a p53-regulated gene. And the consumption of EPB41L4A-AS1 in tumor treatment can increase the sensitivity of tumors to glutaminase inhibitors. 32 2) EPB41L4A-AS1 plays an important role in the diagnosis and prognosis evaluation of various cancers including cervical cancer, liver cancer, breast cancer, bladder cancer, glioblastoma, and colorectal cancer. 33 - 35 That is, the low expression of EPB41L4A-AS1 is associated with a low survival rate of these cancer types. 3) EPB41L4A-AS2 is an important target for anti-metastasis treatment in patients with head and neck squamous cell carcinoma (HNSCC). 36 It can also effectively evaluate the prognosis of anti-metastasis treatment in patients with HNSCC. 4) Long non-coding RNA EPB41L4A-AS2 can inhibit cell proliferation, invasion and promote apoptosis of hepatocellular carcinoma, renal cancer, non-small cell lung cancer and breast cancer. 37 - 40 However, there is currently no clinical survival study on the EPB41L4A gene in any types of cancer. There are only a few studies related to genetics of the EPB41L4A. For example: 1) Studies have revealed genetic alterations in two patients with tuberous sclerosis complex-associated renal cell carcinoma (TSC-RCC) by whole exome sequencing (WES). However, only one patient with TSC-RCC was found to have a genetic alteration in EPB41L4A. 31 2) It was found by WES that EPB41L4A has genetic susceptibility to colon, nerve/brain muscle and congenital keratosis. 20 To understand the role of EPB41L4A in cancer, studies have explored the role of EPB41L4A in cell lines. This study found that when β-catenin was depleted in SW480 cells, the RNA expression of EPB41L4A was significantly reduced (SW480 cells are colorectal cancer cell lines). 41 To the best of our knowledge, we have not found any study about the clinical prognostic relationship between EPB41L4A and tumor. In our study, we obtained 1956 samples by integrating 7 datasets and found a relationship between the expression level of EPB41L4A and the clinical prognosis of MM.

By analyzing the dataset of 559 MM patients, we found that EPB41L4A is a "good" gene, and high expression of EPB41L4A suggests a better prognosis. There are 4 points of support for EPB41L4A low expression suggest poor prognosis. 1) 1q21 gene amplification suggests a low survival rate and a poor prognosis, and the higher the amplification, the worse the prognosis. 42 , 43 We found that the expression level of EPB41L4A is opposite to the 1q21 amplification copy number. 2) The molecular typing by MAF, MMSET, and Proliferation suggested that MM had a poor prognosis. 44 We found that the expression of EPB41L4A in these molecular typing was significantly lower than the average of all 7 molecular typing. 3) EFS and OS in EPB41L4A high expression group were significantly higher than EPB41L4A low expression group. 4) The expression level of EPB41L4A in patients with relapsed MM is lower than that in patients without recurrence.

Our study found that some "good" genes are also highly expressed in patients with high expression of EPB41L4A, while some "bad" genes are also highly expressed in patients with low expression of EPB41L4A. For example, these genes are ISL2, KIT, CCND1 (highly expressed in patients with high expression of EPB41L4A), CCND2 and FGFR3 (highly expressed in patients with low expression of EPB41L4A) and the like. KIT (CD117) is a transmembrane glycoprotein of the class III receptor tyrosine kinase family. 45 Several studies have shown that c-Kit expression in MM cells is functional and associated with survival pathways (Akt pathway). High expression of KIT suggests a good prognosis. 46 - 48 High expression of CCND1 can make MM patients have a better prognosis by highly specific inhibition of translation of myeloma cells. 49 , 50 However, CCND1 is also involved in the most common translocation in myeloma 11 ; 14 . The translocation 11 ; 14 placed CCND1 under the transcriptional control of the immunoglobulin heavy chain (IgH) enhancer, resulting in a dysregulation of CCND1, thereby accelerating the G1 to S phase transition in the plasma cell. 51 , 52 Arjun Lakshman et al. found that patients with high expression of CCND1 with 11 ; 14 translocation had a poor prognosis. 53 Overexpression of CCND2 reduces survival of MM patients by affecting RNA transcription. 54 , 55 Fibroblast growth factor receptor 3 (FGFR3) is a transmembrane tyrosine kinase receptor. 56 , 57 It is involved in intracellular signaling pathways and promotes cell proliferation, migration and differentiation. Overexpression of FGFR3 can lead to proliferation of myeloma cells leading to poor prognosis in MM patients, but it is often used as a therapeutic target to improve MM prognosis. 58 , 59 Not only is the expression of EPB41L4A highly consistent with “good” or “bad” genes in MM, but EPB41L4A expression is an independent prognostic factor for MM patients (see results).

EPB41L4A is a key target of the wnt signaling pathway, 13 which can affect cell proliferation, migration and differentiation. 19 EPB41L4A improves the prognosis of patients with MM and is related to the cell cycle. We found that the top 15 pathways of EPB41L4A are related to the cell cycle, and the high expression of EPB41L4A will reduce the expression of most genes in DNA replication and cell proliferation pathway. We compared the expression of EPB41L4A before and after relapse in the same MM patient and compared the expression of EPB41L4A in the same patient after untreated and pre-1st chemotherapy, and found that there was no difference in the expression of EPB41L4A before and after, but we also found that MM patients with relapsed had lower expression levels of EPB41L4A than MM patients with non-relapsed. Therefore, we hypothesize that the prognosis and recurrence levels of MM patients may be determined by the initial state of EPB41L4A expression level, but this requires further experiments to verify.

However, our research has some shortcomings. The detailed mechanism of EPB41L4A in MM is not deep enough conducted. Multiple markers were not combined to observe clinical predictions of MM patients. There is also a lack of prospective studies to demonstrate whether the expression level of EPB41L4A can alter clinical decisions for early treatment.

In summary, the high expression of EPB41L4A predicts good survival level in MM. The high expression of EPB41L4A was show as a good classifier in MM. The expression level of the EPB41L4A can predict the prognosis of MM patients, but is determined by the initial expression state of EPB41L4A, and the expression level of EPB41L4A does not change due to treatment and relapse. The EPB41L4A may be related to cell division and DNA replication pathway which can cause the better survival level and lower recurrence level of MM.

Supplementary Material

Supplementary figures and tables.

Click here for additional data file.

We thank 7 datasets ( GSE24080 , GSE9782 , GSE19784 , GSE83503 , GSE39754 , GSE82307 and GSE19554 ) which have previously been deposited by others in GEO database.

Funding

This work was funded by National Natural Science Foundation of China (81800195 and 81460315), interdisciplinary medicine Seed Fund of Peking University (BMU2018MB004), Beijing Natural Science Foundation (7132183 and 7182178), China Health Promotion Foundation (CHPF-zlkysx-001), Scientific Research Foundation (20141114) from Health Commission of Jiangxi Province, and Science and Technology Research Foundation (GJJ14676) from Educational Commission of Jiangxi Province, China.

Ethics approval and consent to participate

This study was approved by the Helsinki declaration of 1964 and its subsequent amendments

Author contributions

HMJ, WYL and XRZ conceived the project. WLZ, RL and XH analyzed the data. WLZ, RL, XH, XNL, YZ, ZZY, PY, JW, KH, XLY, XRZ, WYL and HMJ contributed towards the interpretation of the data. All authors wrote and approved the final manuscript.

Abbreviations

Multiple myeloma

International Myeloma Working Group

International Staging System

Revised International Classification System

autologous stem cell transplantation

erythrocyte protein band 4.1-like 4a, also called Nbl4

Four-point-one, Ezrin, Radixin, Moesin

Gene Expression Omnibus

Event-free survival time, defined from date of registration to the occurrence of death from any cause, disease progression or relapse, or censored at the date of last contact

Overall survival time, defined from date of registration to the date of death from any cause or censored at the date of last contact

Beta-2 microglobulin

Albumin

Haemoglobin

Number of magnetic resonance imaging (MRI)-defined focal lesions (skull, spine, pelvis)

Bone marrow biopsy plasma cells (%)

Hazard ratio

Confidence interval

C-reactive protein

Creatinine

Aspirate plasma cells (%)

head and neck squamous cell carcinoma

tuberous sclerosis complex-associated renal cell carcinoma

whole exome sequencing

1 Palumbo A Anderson K Multiple myeloma The New England journal of medicine 2011 364 11 1046 60 21410373 10.1056/NEJMra1011442

2 Kumar SK Rajkumar SV The multiple myelomas - current concepts in cytogenetic classification and therapy Nature reviews Clinical oncology 2018 15 7 409 21 10.1038/s41571-018-0018-y 29686421

3 Kyle RA Rajkumar SV An overview of the progress in the treatment of multiple myeloma Expert review of hematology 2014 7 1 5 7 24433105 10.1586/17474086.2014.870030

4 Siegel RL Miller KD Jemal A Cancer statistics, 2018 CA Cancer J Clin 2018 68 1 7 30 29313949 10.3322/caac.21442

5 Cossu G Terrier LM Benboubker L Destrieux C Velut S Francois P Spinal metastases in multiple myeloma: A high-risk subgroup for ISS III patients Surgical oncology 2018 27 2 321 6 29937188 10.1016/j.suronc.2018.05.005

6 Durie BG Salmon SE A clinical staging system for multiple myeloma. Correlation of measured myeloma cell mass with presenting clinical features, response to treatment, and survival Cancer 1975 36 3 842 54 1182674 10.1002/1097-0142(197509)36:3<842::aid-cncr2820360303>3.0.co;2-u

7 Greipp PR San Miguel J Durie BG Crowley JJ Barlogie B Blade J International staging system for multiple myeloma Journal of clinical oncology: official journal of the American Society of Clinical Oncology 2005 23 15 3412 20 15809451 10.1200/JCO.2005.04.242

8 Palumbo A Avet-Loiseau H Oliva S Lokhorst HM Goldschmidt H Rosinol L Revised International Staging System for Multiple Myeloma: A Report From International Myeloma Working Group Journal of clinical oncology: official journal of the American Society of Clinical Oncology 2015 33 26 2863 9 26240224 10.1200/JCO.2015.61.2267 PMC4846284

9 Fonseca R Bergsagel PL Drach J Shaughnessy J Gutierrez N Stewart AK International Myeloma Working Group molecular classification of multiple myeloma: spotlight review Leukemia 2009 23 12 2210 21 19798094 10.1038/leu.2009.174 PMC2964268

10 Scalzulli E Grammatico S Vozella F Petrucci MT Proteasome inhibitors for the treatment of multiple myeloma Expert opinion on pharmacotherapy 2018 19 4 375 86 29478351 10.1080/14656566.2018.1441287

11 Kohler M Greil C Hudecek M Lonial S Raje N Wasch R Current developments in immunotherapy in the treatment of multiple myeloma Cancer 2018 124 10 2075 85 29409124 10.1002/cncr.31243

12 Marchesini M Ogoti Y Fiorini E Aktas Samur A Nezi L D'Anca M ILF2 Is a Regulator of RNA Splicing and DNA Damage Response in 1q21-Amplified Multiple Myeloma Cancer cell 2017 32 1 88 100.e6 28669490 10.1016/j.ccell.2017.05.011 PMC5593798

13 Guo Y Christine KS Conlon F Gessert S Kuhl M Expression analysis of epb41l4a during Xenopus laevis embryogenesis Development genes and evolution 2011 221 2 113 9 21556855 10.1007/s00427-011-0362-8 PMC3612984

14 Tepass U FERM proteins in animal morphogenesis Current opinion in genetics & development 2009 19 4 357 67 19596566 10.1016/j.gde.2009.05.006

15 Rosenbluh J Wang X Hahn WC Genomic insights into WNT/beta-catenin signaling Trends in pharmacological sciences 2014 35 2 103 9 24365576 10.1016/j.tips.2013.11.007 PMC3917662

16 Anastas JN Moon RT WNT signalling pathways as therapeutic targets in cancer Nature reviews Cancer 2013 13 1 11 26 23258168 10.1038/nrc3419

17 Wang R Sun Q Wang P Liu M Xiong S Luo J Notch and Wnt/beta-catenin signaling pathway play important roles in activating liver cancer stem cells Oncotarget 2016 7 5 5754 68 26735577 10.18632/oncotarget.6805 PMC4868719

18 Ghahhari NM Babashah S Interplay between microRNAs and WNT/beta-catenin signalling pathway regulates epithelial-mesenchymal transition in cancer European journal of cancer (Oxford, England: 1990) 2015 51 12 1638 49 10.1016/j.ejca.2015.04.021 26025765

19 Spaan I Raymakers RA van de Stolpe A Peperzak V Wnt signaling in multiple myeloma: a central player in disease with therapeutic potential J Hematol Oncol 2018 11 1 67 85 29776381 10.1186/s13045-018-0615-3 PMC5960217

20 Mathur P Medicherla KM Chaudhary S Patel M Bagali P Suravajhala P Whole exome sequencing reveals rare variants linked to congenital pouch colon Sci Rep 2018 8 1 6646 54 29703930 10.1038/s41598-018-24967-y PMC5923232

21 Shi L Campbell G Jones WD Campagne F Wen Z Walker SJ The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models Nature biotechnology 2010 28 8 827 38 10.1038/nbt.1665 PMC3315840 20676074

22 Mulligan G Mitsiades C Bryant B Zhan F Chng WJ Roels S Gene expression profiling and correlation with outcome in clinical trials of the proteasome inhibitor bortezomib Blood 2007 109 8 3177 88 17185464 10.1182/blood-2006-09-044974

23 Weinhold N Ashby C Rasche L Chavan SS Stein C Stephens OW Clonal selection and double-hit events involving tumor suppressor genes underlie relapse in myeloma Blood 2016 128 13 1735 44 27516441 10.1182/blood-2016-06-723007 PMC5043128

24 Zhou W Yang Y Xia J Wang H Salama ME Xiong W NEK2 induces drug resistance mainly through activation of efflux drug pumps and is associated with poor prognosis in myeloma and other cancers Cancer cell 2013 23 1 48 62 23328480 10.1016/j.ccr.2012.12.001 PMC3954609

25 Broyl A Hose D Lokhorst H de Knegt Y Peeters J Jauch A Gene expression profiling for molecular classification of multiple myeloma in newly diagnosed patients Blood 2010 116 14 2543 53 20574050 10.1182/blood-2009-12-261032

26 Miannay B Minvielle S Roux O Drouin P Avet-Loiseau H Guerin-Charbonnel C Logic programming reveals alteration of key transcription factors in multiple myeloma Sci Rep 2017 7 1 9257 69 28835615 10.1038/s41598-017-09378-9 PMC5569101

27 Chauhan D Tian Z Nicholson B Kumar KG Zhou B Carrasco R A small molecule inhibitor of ubiquitin-specific protease-7 induces apoptosis in multiple myeloma cells and overcomes bortezomib resistance Cancer cell 2012 22 3 345 58 22975377 10.1016/j.ccr.2012.08.007 PMC3478134

28 Huang da W Sherman BT Lempicki RA Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources Nature protocols 2009 4 1 44 57 19131956 10.1038/nprot.2008.211

29 Kuiper R van Duin M van Vliet MH Broijl A van der Holt B El Jarari L Prediction of high- and low-risk multiple myeloma based on gene expression and the International Staging System Blood 2015 126 17 1996 2004 26330243 10.1182/blood-2015-05-644039 PMC4616233

30 Dispenzieri A Myeloma: management of the newly diagnosed high-risk patient Hematology American Society of Hematology Education Program 2016 2016 1 485 94 27913520 10.1182/asheducation-2016.1.485 PMC6142458

31 Park JH Lee C Chang MS Kim K Choi S Lee H Molecular Characterization and Putative Pathogenic Pathways of Tuberous Sclerosis Complex-Associated Renal Cell Carcinoma Translational oncology 2018 11 4 962 70 29925043 10.1016/j.tranon.2018.05.010 PMC6008483

32 Liao M Liao W Xu N Li B Liu F Zhang S LncRNA EPB41L4A-AS1 regulates glycolysis and glutaminolysis by mediating nucleolar translocation of HDAC2 EBioMedicine 2019 41 200 13 30796006 10.1016/j.ebiom.2019.01.035 PMC6444057

33 Korivi Jyothiraj S Meenakumari B Sundersingh S Sridevi V Rajkumar T Pandey A Construction of lncRNA-mediated ceRNA network to reveal clinically relevant lncRNA biomarkers in glioblastomas Molecular oncology 2019 17 5 4369 74 10.3892/ol.2019.10114 PMC6444437 30944630

34 He M Lin Y Xu Y Identification of prognostic biomarkers in colorectal cancer using a long non-coding RNA-mediated competitive endogenous RNA network Oncology letters 2019 17 3 2687 94 30854042 10.3892/ol.2019.9936 PMC6365949

35 Deva Magendhra Rao AK Patel K Identification of lncRNAs associated with early-stage breast cancer and their prognostic implications Molecular oncology 2019 13 6 1342 55 30959550 10.1002/1878-0261.12489 PMC6547626

36 Huang T Huang W Lu H Zhang BY Ma J Zhao D Identification and validation a TGF-β-associated long non-coding RNA of head and neck squamous cell carcinoma by bioinformatics method Journal of Translational Medicine 2018 16 1 46 60 29490660 10.1186/s12967-018-1418-6 PMC5831574

37 Xu S Wang P You Z Meng H Mu G Bai X The long non-coding RNA EPB41L4A-AS2 inhibits tumor proliferation and is associated with favorable prognoses in breast cancer and other solid tumors Oncotarget 2016 7 15 20704 17 26980733 10.18632/oncotarget.8007 PMC4991486

38 Wang YG Wang T Shi M Zhai B Long noncoding RNA EPB41L4A-AS2 inhibits hepatocellular carcinoma development by sponging miR-301a-5p and targeting FOXL1 Journal of experimental & clinical cancer research: CR 2019 38 1 153 66 30971290 10.1186/s13046-019-1128-9 PMC6458726

39 Shu J Li S Chen YB Zhu QF Yu XH Long non-coding RNA EPB41L4A-AS2 inhibited non-small cell lung cancer proliferation, invasion and promoted cell apoptosis Neoplasma 2018 65 5 664 72 30249102 10.4149/neo_2018_170713N480

40 Pang B Wang Q Ning S Wu J Zhang X Chen Y Landscape of tumor suppressor long noncoding RNAs in breast cancer Journal of experimental & clinical cancer research: CR 2019 38 1 79 97 30764831 10.1186/s13046-019-1096-0 PMC6376750

41 Ishiguro H Furukawa Y Daigo Y Miyoshi Y Nagasawa Y Nishiwaki T Isolation and characterization of human NBL4, a gene involved in the beta-catenin/tcf signaling pathway Japanese journal of cancer research: Gann 2000 91 6 597 603 10874211 10.1111/j.1349-7006.2000.tb00987.x PMC5926395

42 Grzasko N Hus M Pluta A Jurczyszyn A Walter-Croneck A Morawska M Additional genetic abnormalities significantly worsen poor prognosis associated with 1q21 amplification in multiple myeloma patients Hematological oncology 2013 31 1 41 8 22674819 10.1002/hon.2018

43 Shaughnessy JD Jr Zhan F Burington BE Huang Y Colla S Hanamura I A validated gene expression model of high-risk multiple myeloma is defined by deregulated expression of genes mapping to chromosome 1 Blood 2007 109 6 2276 84 17105813 10.1182/blood-2006-07-038430

44 Decaux O Lode L Magrangeas F Charbonnel C Gouraud W Jezequel P Prediction of survival in multiple myeloma based on gene expression profiles reveals cell cycle and chromosomal instability signatures in high-risk patients and hyperdiploid signatures in low-risk patients: a study of the Intergroupe Francophone du Myelome Journal of clinical oncology: official journal of the American Society of Clinical Oncology 2008 26 29 4798 805 18591550 10.1200/JCO.2007.13.8545

45 Ullrich A Schlessinger J Signal transduction by receptors with tyrosine kinase activity Cell 1990 61 2 203 12 2158859 10.1016/0092-8674(90)90801-k

46 Montero JC Lopez-Perez R San Miguel JF Pandiella A Expression of c-Kit isoforms in multiple myeloma: differences in signaling and drug sensitivity Haematologica 2008 93 6 851 9 18443272 10.3324/haematol.12171

47 Bataille R Pellat-Deceunynck C Robillard N Avet-Loiseau H Harousseau JL Moreau P CD117 (c-kit) is aberrantly expressed in a subset of MGUS and multiple myeloma with unexpectedly good prognosis Leukemia research 2008 32 3 379 82 17767956 10.1016/j.leukres.2007.07.016

48 Miettinen M Lasota J KIT (CD117): a review on expression in normal and neoplastic tissues, and mutations and their clinicopathologic correlation Applied immunohistochemistry & molecular morphology: AIMM 2005 13 3 205 20 16082245 10.1097/01.pai.0000173054.83414.22

49 Manier S Huynh D Inhibiting the oncogenic translation program is an effective therapeutic strategy in multiple myeloma Science translational medicine 2017 9 389 eaal2668 (017) 28490664 10.1126/scitranslmed.aal2668 PMC5718051

50 Walker BA Mavrommatis K Wardell CP Identification of novel mutational drivers reveals oncogene dependencies in multiple myeloma Blood 2018 132 6 587 97 29884741 10.1182/blood-2018-03-840132 PMC6097138

51 Chesi M Bergsagel PL Brents LA Smith CM Gerhard DS Kuehl WM Dysregulation of cyclin D1 by translocation into an IgH gamma switch region in two multiple myeloma cell lines Blood 1996 88 2 674 81 8695815

52 Specht K Haralambieva E Bink K Kremer M Mandl-Weber S Koch I Different mechanisms of cyclin D1 overexpression in multiple myeloma revealed by fluorescence in situ hybridization and quantitative analysis of mRNA levels Blood 2004 104 4 1120 6 15090460 10.1182/blood-2003-11-3837

53 Lakshman A Alhaj Moustafa M Rajkumar SV Dispenzieri A Gertz MA Buadi FK Natural history of t(11;14) multiple myeloma Leukemia 2018 32 1 131 8 28655925 10.1038/leu.2017.204

54 Misiewicz-Krzeminska I Sarasquete ME Vicente-Duenas C Krzeminski P Wiktorska K Corchete LA Post-transcriptional Modifications Contribute to the Upregulation of Cyclin D2 in Multiple Myeloma Clinical cancer research: an official journal of the American Association for Cancer Research 2016 22 1 207 17 26341922 10.1158/1078-0432.CCR-14-2796

55 Gutierrez NC Sarasquete ME Misiewicz-Krzeminska I Delgado M De Las Rivas J Ticona FV Deregulation of microRNA expression in the different genetic subtypes of multiple myeloma and correlation with gene expression profiling Leukemia 2010 24 3 629 37 20054351 10.1038/leu.2009.274

56 Turner N Grose R Fibroblast growth factor signalling: from development to cancer Nature reviews Cancer 2010 10 2 116 29 20094046 10.1038/nrc2780

57 Zhang X Ibrahimi OA Olsen SK Umemori H Mohammadi M Ornitz DM Receptor specificity of the fibroblast growth factor family. The complete mammalian FGF family The Journal of biological chemistry 2006 281 23 15694 700 16597617 10.1074/jbc.M601252200 PMC2080618

58 Kalff A Spencer A The t(4;14) translocation and FGFR3 overexpression in multiple myeloma: prognostic implications and current clinical strategies Blood Cancer Journal 2012 2 e89 (2012) 22961061 10.1038/bcj.2012.37 PMC3461707

59 St-Germain JR Taylor P Tong J Jin LL Nikolic A Stewart II Multiple myeloma phosphotyrosine proteomic profile associated with FGFR3 expression, ligand activation, and drug inhibition Proceedings of the National Academy of Sciences of the United States of America 2009 106 47 20127 32 19901323 10.1073/pnas.0910957106 PMC2775037

ИП Урсолов Владлен Сергеевич

ИНН 500118482352 · ОГРНИП 314500117700049

Реквизиты и контакты