Next Article in Journal
Exosomal miRNAs in the Tumor Microenvironment of Multiple Myeloma
Next Article in Special Issue
Pharmacogenomic Analyses Implicate B Cell Developmental Status and MKL1 as Determinants of Sensitivity toward Anti-CD20 Monoclonal Antibody Therapy
Previous Article in Journal
Bta-miR-206 and a Novel lncRNA-lncA2B1 Promote Myogenesis of Skeletal Muscle Satellite Cells via Common Binding Protein HNRNPA2B1
Previous Article in Special Issue
Molecular Drivers of Myelodysplastic Neoplasms (MDS)—Classification and Prognostic Relevance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The TT Genotype of the KIAA1524 rs2278911 Polymorphism Is Associated with Poor Prognosis in Multiple Myeloma

1
Department of Haematooncology and Bone Marrow Transplantation, Medical University of Lublin, 20-059 Lublin, Poland
2
Body Composition Research Laboratory, Department of Preclinical Sciences, Medical University of Lublin, 20-059 Lublin, Poland
3
Department of Human Physiology, Medical University of Lublin, 20-059 Lublin, Poland
4
Division of Mathematics, Department of Production Computerisation and Robotisation, Lublin University of Technology, 20-618 Lublin, Poland
5
Collegium Medicum, University of Warmia and Mazury in Olsztyn, 10-719 Olsztyn, Poland
6
Department of Pneumonology, Oncology and Allergology, Medical University of Lublin, 20-059 Lublin, Poland
7
Institute of Computer Science, Maria Curie-Sklodowska University, 20-033 Lublin, Poland
*
Author to whom correspondence should be addressed.
Cells 2023, 12(7), 1029; https://doi.org/10.3390/cells12071029
Submission received: 11 January 2023 / Revised: 13 March 2023 / Accepted: 22 March 2023 / Published: 28 March 2023
(This article belongs to the Special Issue Molecular Mechanisms of Hematological Malignancies)

Abstract

:
Background: The KIAA1524 gene encodes an oncoprotein, CIP2A, which inhibits the phosphorylation of the Akt kinase B, stabilizes the c-Myc protein, and, through that, promotes cancerogenesis. An increase in CIP2A expression has been observed in numerous solid tumors and hematologic malignancies, including multiple myeloma (MM). The aim of our study was to evaluate the clinical impact of the functional single nucleotide polymorphisms (SNP) of the KIAA1524 gene (rs2278911, 686C > T) in MM patients. Methods: The study group consisted of 128 patients with de novo MM. EDTA venous blood samples were collected prior to the treatment. The SNPs were analyzed by Real-Time PCR with the use of specific Taqman probes. Results: Multivariable analysis revealed that variables independently associated with shorter progression-free survival (PFS) included thrombocytopenia, delTP53 and IGH/CCND1 translocation and the TT genotype of the KIAA1524 gene (686C > T) (median PFS: 6 vs. 25 months; HR = 7.18). On the other hand, autologous haematopoietic stem cell transplantation (AHSCT) was related to a lower risk of early disease progression. Moreover, light chain disease, International Staging System (ISS) 3, poor performance status, hypoalbuminemia, IGH/FGFR3 translocation and the TT genotype of the KIAA1524 gene (686C > T) were independent prognostic factors associated with shorter overall survival (OS) (median OS: 8 vs. 45 months; HR = 7.08). Conclusion: The evaluation of the SNP 686C > T of the KIAA1524 gene could be used as a diagnostic tool in MM patients at risk of early disease progression and death.

1. Introduction

Multiple myeloma (MM) is a malignant neoplasm characterized by the clonal proliferation of atypical plasma cells within the bone marrow, and accounts for around 10% of all hematologic malignancies [1]. The risk factors for MM include age, race, sex, family history, chronic inflammation, immunodeficiency and exposure to radiation and organic solvents [2,3,4,5]. Patients may experience a wide range of unspecific and usually mild symptoms, including fatigue, nausea, vomiting, weight loss and recurring infections [3]. At the time of presentation, 73% of the patients have anemia, almost 80% have osteolytic bone disease and 19% have acute kidney injury [6]. The imbalance between the osteoblasts’ and the osteoclasts’ activity results in the development of bone lesions [7]. The progression of the disease leads to bone fractures, osteoporosis, multi-organ failure and death [3].
The current mainstay of MM treatment includes proteasome inhibitors (PI) and immunomodulatory drugs (IMiDs) known as monoclonal antibodies. Autologous hematopoietic stem cell transplantation (AHSCT) remains the standard of care for young patients with newly diagnosed MM [7,8]. Although the introduction of new drugs has significantly improved the treatment’s results, MM remains a disease with a very unpredictable course, with the estimated survival ranging from less than six months to more than 10 years [9]. The clinical outcome of MM is known to be affected by several factors, including the patient’s performance status, age, comorbidities, therapy choice, microenvironment, tumor burden and characteristics of the tumor—including the chromosomal changes and gene expression. However, even within the same risk groups, we observe a significant heterogeneity of the outcome [4,10,11].
The MM cells derive from B-line cells, which develop in the germinal centers located within the lymph nodes and the spleen. Malignant development begins with the acquisition of hyperdiploidy or translocation with the immunoglobulin heavy chain (IGH) gene locus [4]. Based on the chromosome ploidy and other parameters, genetic abnormalities are divided into the hyperdiploidy karyotype, which is usually associated with translocations of IGH locus, and the hypodiploid karyotype, which is mostly linked to trisomies [12]. The first stage of the malignant transformation is the B cell’s immortalization. One of the most common initiating mutations, observed in nearly 50% of MM patients, is a chromosomal translocation with the IGH locus at the 14q32 region, which may involve many other loci, such as 11q13, 4p16.3, 6p21, 16q23 or 20q11 [12,13]. Hyperdiploidy is present in 55% of MM patients and is associated with trisomies of chromosomes such as 3, 5, 7, 9, 11, 15, 19 and 21 [12,14]. While trisomies 3 and 5 are beneficial for patients’ survival, trisomy 21 is linked with poor prognosis [14]. The genomic instability promotes further molecular changes. The most common somatic mutations, which are present in around 50% of patients, occur in genes of the RAS/MAPK pathway, such as KRAS, NRAS and BRAF [4]. Other cytogenic changes have been described as secondary genetic aberrations, among which the most frequent are MYC rearrangements, del (13q), del (17p), (del)1p and amp(1q) [10,12]. These have been connected to MM progression.
The signal pathways JAK2/STAT3, Ras/Raf, PI3K/AKT and NFκB have been identified as those in control of MM cells’ apoptosis, migration and proliferation, and may be linked to MM cells’ drug resistance [15,16]. The interactions between MM cells and non-malignant cells from the bone marrow microenvironment—such as stromal and endothelial cells—contribute to elevated levels of interleukin 6 (Il-6), vascular endothelial growth factor (VEGF) and insulin-like growth factor 1 (IGF-1), which promote the PI3K/Akt/mTOR pathway’s activation [17]. As that signaling pathway regulates a wide range of MM’s cellular processes, it is under intense research for targeting treatment [16,17]. So far, its role in the development of MMs has not been fully explained; however, in most other neoplasms, its activity is controlled by phosphatases such as PP2A (protein phosphatase 2).
The KIAA1524 gene encodes the cancerous inhibitor of PP2A (CIP2A)—also known as a p90 tumor-associated antigen—which is a human oncoprotein whose overexpression has been observed in patients with solid tumors of the head and neck, breast, stomach and large intestine [18]. The c-Myc stabilization in the tumor cells is promoted by the CIP2A thanks to its ability to inhibit the c-Myc protein’s dephosphorylation, which is mediated by the activity of PP2A towards c-Myc serine 62 [18,19]. Additionally, CIP2A is also responsible for the downregulation of protein kinase B (Akt 9s) phosphorylation [20]. Myc and CIP2A appear to promote each other’s expression in a positive regulatory loop [21]. As the inhibition of CIP2A may allow for the PP2A-mediated proteolytic degradation of c-Myc in cancer cells and thereby limit the cancer cells’ growth and transformation, the CIP2A may be a promising candidate for targeted treatment [18,19].

2. Materials and Methods

2.1. The Study Population

The study group included 128 (64 male and 64 female) newly diagnosed MM patients who received first-line chemotherapy with thalidomide (CTD—cyclophosphamide + thalidomide + dexamethasone; n = 27; 21.1%), bortezomib (V(C)D—bortezomib + cyclophosphamide + dexamethasone; n = 58; 45.3%) or bortezomib and thalidomide (VTD—bortezomib + thalidomide + dexamethasone; n = 43; 33.6%). Patients were recruited in the Department of Hemato-oncology and Bone Marrow Transplantation in Lublin from 2015 to 2020.
MM diagnosis was based on SLiM CRAB criteria, according to the International Myeloma Working Group (IMWG) recommendations. Disease staging was determined according to Durie–Salmon and the current ISS staging system. Treatment response was assessed according to IMWG 2016 consensus criteria [22,23,24,25]. The severity of adverse effects was evaluated using CTCAE (version 5.0) [26].
Detailed patients’ characteristics are presented in Table 1.

2.2. DNA Isolation

The DNA isolation from 200 μL of whole peripheral blood was performed using the column method with a dedicated kit, according to the manufacturer’s recommendations (DNA Blood Mini Kit, Qiagen, Toronto, ON, Canada). Subsequently, the spectrophotometric evaluation of the concentration and quality of the obtained DNA was carried out using a NanoDrop Lite Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).

2.3. Single Nucleotide Polymorphism Analysis

A StepOnePlus Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) was used for the SNP genotyping. The real-time PCR was performed according to the manufacturer-established protocol with the use of the Genotyping Master Mix and TaqMan probes specific for the KIAA1524 SNP (rs2278911; 686C > T) (C__11237075_30, Thermo Fisher Scientific, USA). A selected probe was functionally tested by the manufacturer. The qPCR was performed in 10 µL reaction volume on 96-well plates. A rection mix containing 5 µL of Taqman Genotyping Master Mix, 0.5 µL of Taqman SNP assay (20×) and 4.5 µL of 0.2 ng/µL diluted DNA template was added to each well. The thermal cycling protocol was as follows: enzyme activation at 95 °C, 40 cycles of denaturation at 95 °C for 15 s and annealing at 60 °C for 1 min. All sample tests were performed in triplicate. After the amplification, the obtained genotypes were analyzed with StepOne Software v2.3 (Applied Biosystems, Waltham, MA, USA). Additionally, 10% of the samples were randomly selected and re-analyzed in a sequencing device (3500 Genetic Analyzer, Life Technologies). A 100% consistency of the results was obtained.

2.4. Plasma Cells Isolation

The bone marrow mononuclear cells (BMMC) were isolated from 3 mL of bone marrow (diluted 10×) using the density gradient centrifugation method (Biocoll, AG Biochrom, Berlin Germany), then subjected to washing and a viability assessment followed by the cell separation. The cell separation was performed after the incubation with CD138 antibodies and conducted according to the manufacturer’s instructions with the use of the magnetic-activated cell sorting technique (MACS® Cell Separation, Bergisch Gladback, Germany). Subsequently, the assessment of viability and quantity of plasma cells was performed.

2.5. Gene Expression in the Bone Marrow Plasma Cells

Total RNA was isolated from the plasma cells (the maximum of 107 cells were used) with the use of an RNeasy Mini Kit (Qiagen, Toronto, ON, Canada). The expression of the KIAA and MYC genes was measured with the use of a StepOnePlus device and a dedicated High-Capacity cDNA Reverse Transcription Kit, TaqMan™ Fast Advanced Master Mix and TaqMan probes (Thermo Fisher Scientific, USA). For the KIAA1524 (assay ID: Hs00405413) and MYC (assay ID: hs00153408) genes, the ACTB (assay ID: hs99999903) gene was used as the reference gene. The relative expression of the investigated molecular targets was calculated using the 2−ΔCt and 2−ΔΔCt formulas.

2.6. Cytogenetic Assessment

The bone marrow samples were subjected to cytogenetic analysis. Abnormalities, such as the TP53 gene deletion, IGH gene rearrangements and CKS1B gene amplification, were detected by simultaneous staining of cytoplasmic immunoglobulin with the fluorescence in situ hybridization (cIg-FISH) technique according to the recommendations of Ross et al., with some modifications [27]. The following probes, all from Abbott Molecular (Abbott Park, IL, USA), were used: Vysis TP53/CEP 17 FISH Probe Kit for detection of del(17p13.1), Vysis IGH/FGFR3 Dual Colour, Dual Fusion Translocation Probe for detection of t(4;14)(p16;q32), Vysis IGH/MYC/CEP 8 Tri-colur, Dual Fusion Translocation Probe for detection of t(8;14)(q24;q32), Vysis IGH/CCND1 Dual Colour, Dual Fusion Translocation Probe for detection of t(11;14)(q13;q32), Vysis IGH/MAF Dual Colour, Dual Fusion Translocation Probe for detection of t(14;16)(q32;q23) and Vysis 1q21 CKS1B SpectrumOrange/1p32 CDKN2C SpectrumGreen FISH Probe Kit for the detection of amp(1q32). Fluorescent microscopic analysis was performed by scoring 100 AMCA-positive (anti-mitochondrial antibody-positive) plasma cells to determine the frequency of each aberration. The cut-off values were set at 20% for deletion/amplification probes and 10% for dual fusion probes, according to the recommendations of the European Myeloma Network [15,16].

2.7. Statistical Analysis

MedCalc software (v.15.8) was used for statistical analysis of the acquired data. A chi-square test was performed to assess the KIAA1524 genotypes’ distribution depending on the selected demographic, as well as clinical and molecular variables. Due to the abnormal distribution of the studied continuous data in the D’Agostino–Pearson test, either the non-parametric Mann–Whitney U Test (for comparing two independent groups) or the ANOVA Kruskal–Wallis test (for comparing more than two independent groups) was performed in order to compare the KIAA1524′s expression levels depending on the distribution of SNP located in this gene. The Spearman’s rank correlation coefficient test was applied to assess the correlation between the expression of the KIAA1524 gene and the MYC gene.
The influence of selected demographic, clinical and molecular variables on the risk of no response to chemotherapy was estimated by the odds ratio (OR) and corresponding 95% confidence intervals (CI). In terms of the survival evaluation, the log-rank test was applied for the univariable analysis, whereas Cox regression models were used in the multivariable analysis. The influence of selected demographic, clinical and molecular variables on the risk of death was described by the hazard ratio (HR) and corresponding 95% CI. The Kaplan–Meier estimation method was utilized to generate the survival curves, whereas forest plot graphs were used to illustrate the results of the multivariable analysis. In cases in which both the composite (e.g., disease stage according to ISS) and its constituent variables (e.g., albumin or B2M levels) were statistically significant according to the univariable analysis, only the first one of them was included in the multivariate models. The backward elimination method was used for the selection of variables to be included in the results adjustment. Finally, in the multivariable PFS analysis, all results were adjusted by the following variables: gender, AHSCT, PLT, del17/TP53, t(11;14) IGH/CCND1 and KIAA gene genotype. On the other hand, in the multivariable OS analysis, all results were adjusted by diagnosis, performance status, anemia before treatment, t(4;14) IGH/FGFR3 and KIAA gene genotype. The progression-free survival (PFS) was defined as the period between the start of the treatment and the occurrence of disease progression (complete data), or of the last documented follow-up (censored data). The overall survival (OS) was calculated from the start of treatment until the patient’s death (complete data) or the last documented follow-up (censored).
The data from the MMRF CoMMpass project were retrieved from “https://portal.gdc.cancer.gov/ (accessed on 21 March 2023)” through the R language wrapper library GenomicDataCommons (installed by BiocManager). Based on the obtained data, Pearson correlations between genes of interest were computed, and further selected results were illustrated with the use of the corrplot library (with heat maps presenting correlation matrix).
In all analyses, results with a p-value < 0.05 was considered statistically significant. Additionally, p-values in the range of 0.05–0.06 were considered to show a trend toward statistical significance.

3. Results

3.1. The KIAA1524 Genotypes’ Distribution in Relation to the Demographic, Clinical and Molecular Factors

The only significant difference observed between the patients with various KIAA1524 genotypes was the prevalence of the t(4;14) IGH/FGFR3, whose absence was more common in patients with the CC genotype than in the carriers of the T allele (80% vs. 20%, p = 0.0095). The detailed data is presented in the Supplementary Table S1.

3.2. Response to Treatment

A significantly higher risk of non-response after the 2nd cycle of chemotherapy was observed in patients with poor performance status (PS: 2–4: OR = 3.95), those who were treated with a CTD regimen (OR = 5.35), those with elevated CRP (OR = 3.32) and those with the TT genotype of the KIAA1524 gene (OR = 8.33). On the other hand, the risk of non-response after the 2nd cycle of chemotherapy was significantly lower in patients treated with the VTD regimen (OR = 0.05) and those with the CC genotype of the KIAA1524 gene (OR = 0.23). In turn, a significantly higher risk of non-response after the 4th cycle of chemotherapy was observed in patients with elevated CRP (OR = 6.70), those with the presence of t(11.14)/CCND1 translocation (a trend toward significance; OR = 6.42) and those who were carriers of the TT genotype of the KIAA1524 gene (OR = 20.20). On the other hand, a significantly higher risk of non-response after the 6th cycle of therapy was observed in patients with thrombocytopenia (OR = 9.22), those with the presence of 17p/TP53 deletion (OR = 12.46) and those with the presence of t(11.14)/CCND1 translocation (OR = 55.20). Moreover, a significantly higher risk of non-response after the 8th cycle of therapy was observed in patients with the presence of t(11.14)/CCND1 translocation (OR = 16.00). Detailed data on the influence of the selected demographic, clinical and molecular variables on chemotherapy response are presented in Supplementary Table S2.

3.3. Progression-Free Survival

3.3.1. Univariable Analysis

According to the univariable analysis, the factors significantly correlated with an increased risk of PFS reduction were: male sex (HR = 1.65), stage 3 disease according to the ISS classification (HR = 1.74), higher stage of chronic kidney disease (>G2: HR = 1.85), poor performance status (PS: 2–4; HR = 1.60), anemia before treatment (HR = 2.40) low platelet count (HR = 2.08), low albumin level (HR = 2.29), high LDH level (HR = 1.95), high β2-microglobulin level (B2M) (HR = 2.58), high level of creatinine (HR = 1.64), the presence of del17p/TP53 (HR = 1.84) and the presence of t(11;14) IGH/CCND1 (HR = 2.49). Additionally, a higher risk of PFS reduction was observed in patients with the TT variant of the KIAA1524 gene (rs2278911) as compared to patients with other variants (median PFS: 6 vs. 25 months; HR = 4.67) (Figure 1A). A decreased risk of PFS reduction was observed in patients treated with VTD as compared to those treated with CTD or V(C)D schemes (HR = 0.56), those in whom AHSCT was applied (HR = 0.37) and those with the CC genotype (compared to carriers of the T allele, HR = 0.40).

3.3.2. Multivariable Analysis

Multivariable analysis confirmed that thrombocytopenia (HR = 3.07), high LDH level (HR = 2.59), presence of del17p/TP53 (HR = 2.25), presence of t(11;14) IGH/CCND1 (a trend toward significance; HR = 2.34) and the TT genotype of the KIAA1524 gene (rs2278911) (HR = 7.57) were independent prognostic factors for a higher risk of progression. However, AHSCT was the only independent prognostic factor for a lower risk of progression (HR = 0.48) [Figure 2A].

3.4. Overall Survival

3.4.1. Univariable Analysis

The univariable analysis identified the following factors as significantly related to higher risk of OS death: light chain disease (HR = 2.37), stage 3 disease according to the ISS classification (stage 3; HR = 2.03), higher stage of chronic kidney disease (>G2: HR = 1.78), poor performance status (PS: 2–4; HR = 2.01), the presence of anemia before treatment (HR = 2.43), low albumin levels (HR = 2.28), high B2M levels (HR = 2.73), high creatinine levels (HR = 1.74), the presence of del17p/TP53 (HR = 1.84) and the presence of t(11;14) IGH/CCND1 (HR = 2.62). Moreover, we observed a higher risk of death in patients with the TT variant of the KIAA1524 gene (rs2278911) as compared to those with other genotypes (median OS: 8 months vs. 45 months; HR = 4.39) (Figure 1B). A significantly lower risk of OS reduction was observed in the subjects in whom AHSCT was applied (HR = 0.53) and those with the CC variant of the KIAA1524 gene (rs2278911) as compared to carriers of the T allele (HR = 0.41).

3.4.2. Multivariable Analysis

The multivariable analysis confirmed the independent prognostic value of the following factors: light chain disease (HR = 3.52), the stage 3 disease according to the ISS classification (HR = 2.11), performance status (a trend toward significance; HR = 2.19), low albumin levels (HR = 4.08), the presence of t(4;14), IGH/FGFR3 (HR = 3.22) and the presence of the TT genotype of the KIAA1524 gene (rs2278911) (HR = 6.89) [Figure 2B].
A graphical representation of the results of the multivariable analysis of PFS and OS is shown in Figure 2. Detailed data on the influence of the selected demographic, clinical and molecular variables on PFS and OS are presented in Supplementary Table S3.
Figure 1. Kaplan–Meier curves presenting the influence of KIAA gene genotypes on progression-free survival (A) and overall survival (B).
Figure 1. Kaplan–Meier curves presenting the influence of KIAA gene genotypes on progression-free survival (A) and overall survival (B).
Cells 12 01029 g001
Figure 2. Forest plot showing the results of multivariate analysis for progression-free survival (A) and overall survival (B).
Figure 2. Forest plot showing the results of multivariate analysis for progression-free survival (A) and overall survival (B).
Cells 12 01029 g002

3.5. Relationship between the Expression of the KIAA1524 Gene, Its SNP (rs2278811) and the Expression of the MYC Gene

We noted that the expression of the KIAA1524 gene was significantly lower in the carriers of the CC genotype of this gene compared to the carriers of the CT and TT genotypes (in both cases p < 0.05) (Figure 3A). Moreover, we found an intermediate, positive correlation between expression of the KIAA1524 gene and the MYC gene (rho = 0.515; p = 0.0345) (Figure 3B).

3.6. Relationship between the Expression of the KIAA1524 Gene and the Expression of Other Genes Based on Data from the CoMMpass Project

Based on data from the CoMMpass project, we assessed the correlation between the expression of the KIAA1524 gene (known in the above database as CIP2A) and other genes known to be involved in the biology and pathomechanism of MM (Figure 4). A statistically significant and positive correlation between the expression of the CIP2A gene and the TP53, IL6ST, MCL1, BCL-2, AKT1 and E2F1 genes was found. In addition, a positive correlation between the expression of CIP2A and the MYC gene, which was demonstrated in our study, was further validated in the data retrieved from this project. Moreover, we found that the studied gene was strongly correlated with genes involved, among others, in the regulation of the cell cycle and differentiation, DNA replication, DNA damage detection and repair, mRNA processing and transport, spindle dynamics, chromosome processing and segregation (Figure 5).

4. Discussion

MYC affects a number of metabolic pathways and plays a crucial key role in the development of MM. It plays a role in cycle regulation, proliferation, differentiation and metabolism of the cells, and its dysregulation has been shown to promote the transition from monoclonal gammopathy of undetermined significance (MGUS) to MM. Its activity can be regulated by its mutations, as well as by changes in pathways involved in its regulation (e.g., it may be controlled by mTORC1 and RAS on a post-transcriptional level) [28,29,30].
The two most important phosphorylation sites on MYC, responsible for its stabilization and degradation are Serine62 (S62), which stabilizes MYC and participates in the target gene selection when phosphorylated, and Threonine 58 (T58), which runs the MYC degradation when phosphorylated [30,31]. S62 is phosphorylated by the extracellular-regulated kinases (ERKs) 1 and 2, while the glycogen synthase kinase (GSK-3β) is responsible for the T58′s phosphorylation, which can only be carried out after the S62′s phosphorylation [30]. Before MYC can be delivered to and degraded by the ubiquitin–proteasome system, it must be dephosphorylated in position S62 by PP2A [32,33]. CIP2A inhibits the PP2A, and, through that, stabilizes MYC [30,32]. Moreover, PP2A’s inhibition increases ERK’s activation and promotes cell proliferation and survival, as well as the expression of invasive proteases [30].
Since CIP2A also influences the RAS pathway, it may indirectly regulate the MYC pathway through its elements [34,35,36] The Ras proto-oncogene has been shown to enhance the ERK’s activity, restrict the GSK-3β’s activity and affect the PI3K pathway-mediated T58′s phosphorylation, making it a MYC stabilizer [30]. In addition, the mutations in MYC ubiquitin ligase (FBW7) may prevent the MYC’s degradation, and the ubiquitin-specific protease USP28 has been shown to contribute to the MYC stabilization in neoplastic cells. RalA is a Ras effector, and its phosphorylation is also an important element of PP2A-mediated tumor suppression [30]. Altogether, MYC’s activity is regulated not only by CIP2A, but also by a number of other factors. The complexity of these regulatory systems would justify the intermediate, positive correlation we observed between the expression of the KIAA1524 gene and the MYC gene (rho = 0.515; p = 0.0345) (Figure 2B).
Nevertheless, CIP2A is not the only factor capable of changing the PP2A’s activity. GFI-1 has been shown (growth factor independence-1) to inhibit the SGPP1 (sphingosine-1-phosphate phosphatase 1) gene transcription, leading to increased S1P (sphingosine-1-phosphate) synthesis. This resulted in a S1P/Ceramides ratio imbalance, triggering the PP2A inactivation and, as a consequence, c-Myc stabilization in the MM cells [33]. Such interactions must be considered when judging the prognostic roles of KIAA1524 and CIP2A.
In the case of the 686 SNP of the KIAA1524 (rs2278911) gene, when cystein becomes replaced by tymin, the 229th amino acid coded by the gene changes from arginine (Arg) to glutamine (Gln) [37]. Such an alteration has been shown to lower the charge and change the surface properties of the CIP2A’s N-terminal domain, affecting the protein-to-protein interactions [34]. In addition, Preuss et al. found that certain modifications within the B subunit affected PP2A’s activity towards its targets [38]. A research paper by Wang et al. [35] revealed that, with its N-terminal (amino acids 1–560), CIP2A could bind directly to the B56α and B56γ subunits of the PP2A protein—both of which are tumor suppressor units—and that the region necessary for this was located between the amino acids 159 and 245 and positively charged. Moreover, the researchers found that changing the 230th (asparagine—neutrally charged) nucleotide to a negatively charged glutamic acid led to decreased binding of CIP2A to the B56α, and concluded that the B56 subunit was most likely to bind to a positively charged surface of the CIP2A [35].
In our study, the presence of the TT genotype of the KIAA1524 gene (rs2278911) had an independent prognostic value, with the afflicted patients being over seven times more likely to experience PFS and OS reduction compared to the carriers of other variants. Since the 686C > T alteration resulted in an exchange of a positively charged amino acid for a neutral one, it seems that in order to explain the TT genotype’s correlation with worse prognosis, we must search beyond the currently known direct CIP2A-PP2A interactions. In a paper published in 2014, three years before Wang’s publication, Puustinen et al. wrote that the mechanism in which CIP2A modulated the mTORC1-related PP2A activity was most likely that of an allosteric inhibitor, and stipulated that it was more likely to affect the substrate availability rather than the enzyme’s activity [39]. In 2021, Laine et al. discovered that poor clinical outcomes in basal-like breast cancer (BLBC) patients was related to high CIP2A expression and not to the alterations in its gene’s sequence. They also established that CIP2A expression was upregulated by the DNA-PF/CHK1 activity, TP53 inactivation and EGFR pathway activation [40]. As we have also observed a significantly lower expression of CIP2A among patients with the CC genotype, we would stipulate that that the consequences of indirect CIP2A-mediated PP2A modulation may outweigh the effects of the direct PP2A-binding, and that the most likely cause for the reported outcome–SNP correlation is, in fact, altered CIP2A expression. Our results would also align with a finding by Li et al., who speculated that Arg229Gln altered the CIP2A expression [37].
A number of drugs [41,42,43,44,45] have shown the ability to influence KIAA1524, CIP2A or their affected pathways. Bedewy and Elmaghraby observed that MM patients had higher CIP2A expression, which decreased during bortezomib treatment. Furthermore, higher CIP2A levels were correlated with a worse response to bortezomib, as well as with shorter PFS [45]. The research on KIAA1524 and CIP2A as potential treatment targets is ongoing, but in order for it to progress, many obstacles must be overcome.
One of the biggest issues we are currently facing is a very convoluted net of CIP2A–PP2A interactions, which may be affected by a number of factors and has yet to be fully explored [34,46,47,48], as well as the partially contradicting data on the matter. A high CIP2A concentration has been found to serve as a negative prognostic factor in a number of human cancers [35], and CIP2A knockdown has been shown to inhibit the activity of MYC, E2F1 and Akt [30]; restrict tumor growth in oral cancer mice models [47] and result in decreased proliferation and increased apoptosis of human MM cells [48]. However, such results could not be achieved in all tumor types [46]. Additionally, although there seems to be a general consensus that CIP2A-mediated PP2A inhibition supports MM development, we did come across some research which pointed towards an oncogenic role of PP2A. In 2001 Kang and Chou found that PP2A regulated the ROI (reactive oxygen intermediates) levels within MM cells, which promoted the expression of IL-6 and BCL-2, thus supporting the MM’s development [49]. Three years later, Mitsuthashi et al. demonstrated that PP2A prevented the degradation of the glycoprotein 130 (also known as IL-6 signal transducer). The researchers observed that bortezomib-mediated gp130 degradation provoked the MM cells’ apoptosis by decreasing their sensitivity to the IL-6 [50]. In 2021, Slomp et al. found that PP2A reversed the dephosphorylated and stabilized MCL1 (induced myeloid leukemia cell differentiation protein) in MM [45].
Based on data from the CoMMpass, we identified 30 candidate genes, the regulation of expression of which may be related to KIAA: AC069499.1, MYH15, HNRNPR, KIF11, CENPI, USP1, CEP78, AMMECR1, RRM1, PLK4, TMPO, SMC4, SGO1, CEP55, BUB1, OSBPL11, KIF4A, CAPG, SPAG5, INCENP, POLA1, CDK1, ZWINT, SMC2, PCNA, SPATA5, NCAPG2, GPSM2, NCAPH and TOP2A. Disturbances in their expression have been the subject of many studies on the pathogenesis of solid tumors, e.g., lung, breast, endometrial, stomach and colon cancer [51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68].
Interestingly, these include genes encoding centromeric proteins (CENPI) [51]; centrosomal proteins (CEP55, CEP78) [52,62]; motor proteins (MYH15, KIF11) [53,54] involved in the construction of the karyokinetic spindle or proteins regulating cell division processes (PLK4, SGO1, BUB1, KIF4A, NCAPG, SPAG5, INCENP, CDK1, ZWINT, NCAPG2, GPSM2, NCAPH and TOP2A) [55,60,63,66,67,68,69,70,71,72,73]; and genes encoding proteins involved in DNA repair processes (USP1, RRM1, KIF4A, POLA1 and PCNA) [58,59,65,74,75], which are defined as oncogenes (SPAG5, SMC2) [67,76].
In the context of multiple myeloma, the expression of the HNRNPA2B1, USP1, RRM1, SPAG5, PCNA and TOP2A genes has been studied. It has been proven that HNRNPA2B1 increases the stability of mRNA transcripts for interleukin enhancer-binding factor 3 (ILF3), and its overexpression induces cell proliferation through the AKT3 kinase. Increased HNRNPA2B1 gene expression is associated with an unfavorable prognosis [57]. The USP1 protein leads to regulated DNA repair through deubiquitination and homologous recombination pathways, and inhibits cell differentiation by stabilizing tumor-promoting inhibitors of DNA-binding proteins. Inhibition of USP1 has been shown to reduce myeloma cell viability, inhibit the growth of multiple myeloma cells and overcome bortezomib resistance [58]. Higher expression of the RRM1 gene is associated with shorter overall survival. A blockade of the RRM1 gene results in the upregulation of DNA damage response genes and p53-regulated genes, thereby inhibiting the growth and promoting apoptosis of myeloma cells and cells of the bone marrow microenvironment [59]. SPAG5 gene expression correlates with myeloma cell malignancy, and SPAG5 overexpression is associated with poorer clinical outcomes. In vitro, knockdown of the SPAG5 gene leads to anti-MM effects, cell growth arrest and apoptosis via the PI3K/AKT/mTOR pathway [67]. The PCNA gene encodes a multifunctional protein which is essential for DNA replication and repair and is overexpressed in myeloma cells. PCNA inhibition induces apoptosis and sensitizes cells to melphalan, a common anti-myeloma drug [75]. TOP2A has been demonstrated to be another potential target for anti-myeloma therapy. It has been observed that overexpression of TOP2A topoisomerase in myeloma cells is associated with resistance to proteasome inhibitors [77].

5. Conclusions

The varied clinical course of MM, differences in the efficacy of specific therapeutic strategies, and different lengths of time needed for chemoresistance development imply a need to identify risk stratification factors that would enable the personalization of therapy and the improvement of treatment outcomes. For this purpose, we can use prognostic markers, which assess the risk of an unfavorable course of a disease and its recurrence, as well as predictive factors, which determine the probability of obtaining remission with a specific treatment method. Although our results are preliminary, they indicate that the assessment of the SNP 686C > T of the KIAA1524 gene may be a useful biomarker for determining the prognosis of MM patients.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells12071029/s1, Table S1: Distribution of KIAA1524 genotypes depending on selected demographic, clinical and molecular variables. Table S2: Chemotherapy response depending on selected demographic, clinical and molecular variables. Table S3: Survival of MM patients depending on selected demographic, clinical and molecular variables.

Author Contributions

Conceptualization, A.S.-S. and R.M.; methodology, R.M., M.M., S.C., I.H.-M., T.K. and A.K.; formal analysis, R.-M. and A.K.; investigation, A.S.-S.; data curation, A.S.-S. and S.C.; writing—original draft preparation, A.S.-S., R.M., P.B. and M.J.; writing—review and editing, A.S.-S. and M.S.; visualization, A.S.-S.; supervision, T.M.-M. and M.H.; project administration, A.S.-S.; funding acquisition, A.S.-S. and M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Statutory Funds of the Medical University of Lublin No. DS177, provided by the Polish Ministry of Science and Higher Education.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki. Local ethical approval was obtained from the Bioethical Committee at the Medical University of Lublin (protocol code KE-0254/77/2019).

Informed Consent Statement

Informed consent was obtained from all patients involved in the study.

Data Availability Statement

The research data are available from the corresponding author upon reasonable request. The data are not publicly available due to the fact that the data sheet contains information that exceeds the scope of this study, and which may be used for other research papers in the future.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rajkumar, S.V. Multiple myeloma: Every year a new standard? Hematol. Oncol. 2019, 37 (Suppl. 1), 62–65. [Google Scholar] [CrossRef] [Green Version]
  2. Fairfield, H.; Falank, C.; Avery, L.; Reagan, M.R. Multiple myeloma in the marrow: Pathogenesis and treatments. Ann. N. Y. Acad. Sci. 2016, 1364, 32–51. [Google Scholar] [CrossRef] [Green Version]
  3. Michels, T.C.; E Petersen, K. Multiple Myeloma: Diagnosis and Treatment. Am. Fam. Physician 2017, 95, 373–383. [Google Scholar] [PubMed]
  4. Van de Donk, N.W.C.J.; Pawlyn, C.; Yong, K.L. Multiple myeloma. Lancet 2021, 397, 410–427. [Google Scholar] [CrossRef] [PubMed]
  5. Sergentanis, T.N.; Zagouri, F.; Tsilimidos, G.; Tsagianni, A.; Tseliou, M.; Dimopoulos, M.A.; Psaltopoulou, T. Risk Factors for Multiple Myeloma: A Systematic Review of Meta-Analyses. Clin. Lymphoma Myeloma Leuk. 2015, 15, 563–577.e3. [Google Scholar] [CrossRef]
  6. Cowan, A.J.; Green, D.J.; Kwok, M.; Lee, S.; Coffey, D.G.; Holmberg, L.A.; Tuazon, S.; Gopal, A.K.; Libby, E.N. Diagnosis and Management of Multiple Myeloma. JAMA 2022, 327, 464–477. [Google Scholar] [CrossRef] [PubMed]
  7. Palumbo, A.; Anderson, K. Multiple Myeloma. N. Engl. J. Med. 2011, 364, 1046–1060. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Al Hamed, R.; Bazarbachi, A.H.; Malard, F.; Harousseau, J.-L.; Mohty, M. Current status of autologous stem cell transplantation for multiple myeloma. Blood Cancer J. 2019, 9, 44. [Google Scholar] [CrossRef] [Green Version]
  9. Greipp, P.R.; Miguel, J.S.; Durie, B.G.; Crowley, J.J.; Barlogie, B.; Bladé, J.; Boccadoro, M.; Child, J.A.; Avet-Loiseau, H.; Kyle, R.A.; et al. International Staging System for Multiple Myeloma. J. Clin. Oncol. 2005, 23, 3412–3420. [Google Scholar] [CrossRef] [Green Version]
  10. Schavgoulidze, A.; Cazaubiel, T.; Perrot, A.; Avet-Loiseau, H.; Corre, J. Multiple Myeloma: Heterogeneous in Every Way. Cancers 2021, 13, 1285. [Google Scholar] [CrossRef]
  11. Avet-Loiseau, H.; Attal, M.; Campion, L.; Caillot, D.; Hulin, C.; Marit, G.; Stoppa, A.-M.; Voillat, L.; Wetterwald, M.; Pegourie, B.; et al. Long-Term Analysis of the IFM 99 Trials for Myeloma: Cytogenetic Abnormalities [t(4;14), del(17p), 1q gains] Play a Major Role in Defining Long-Term Survival. J. Clin. Oncol. 2012, 30, 1949–1952. [Google Scholar] [CrossRef] [PubMed]
  12. Hu, L.; Ru, K.; Zhang, L.; Huang, Y.; Zhu, X.; Liu, H.; Zetterberg, A.; Cheng, T.; Miao, W. Fluorescence in situ hybridization (FISH): An increasingly demanded tool for biomarker research and personalized medicine. Biomark. Res. 2014, 2, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Takimoto, M.; Ogawa, K.; Kato, Y.; Saito, T.; Suzuki, T.; Irei, M.; Shibuya, Y.; Suzuki, Y.; Kato, M.; Inoue, Y.; et al. Close relation between 14q32/IGH translocations and chromosome 13 abnormalities in multiple myeloma: A high incidence of 11q13/CCND1 and 16q23/MAF. Int. J. Hematol. 2008, 87, 260–265. [Google Scholar] [CrossRef] [PubMed]
  14. Corre, J.; Munshi, N.C.; Avet-Loiseau, H. Risk factors in multiple myeloma: Is it time for a revision? Blood 2021, 137, 16–19. [Google Scholar] [CrossRef]
  15. Podar, K.; Chauhan, D.; Anderson, K.C. Bone marrow microenvironment and the identification of new targets for myeloma therapy. Leukemia 2008, 23, 10–24. [Google Scholar] [CrossRef] [Green Version]
  16. Ramakrishnan, V.; Kimlinger, T.; Haug, J.; Painuly, U.; Wellik, L.; Halling, T.; Rajkumar, S.V.; Kumar, S. Anti-Myeloma Activity of Akt Inhibition Is Linked to the Activation Status of PI3K/Akt and MEK/ERK Pathway. PLoS ONE 2012, 7, e50005. [Google Scholar] [CrossRef] [Green Version]
  17. Ramakrishnan, V.; Kumar, S. PI3K/AKT/mTOR pathway in multiple myeloma: From basic biology to clinical promise. Leuk. Lymphoma 2018, 59, 2524–2534. [Google Scholar] [CrossRef]
  18. Côme, C.; Laine, A.; Chanrion, M.; Edgren, H.; Mattila, E.; Liu, X.; Jonkers, J.; Ivaska, J.; Isola, J.; Darbon, J.-M.; et al. CIP2A Is Associated with Human Breast Cancer Aggressivity. Clin. Cancer Res. 2009, 15, 5092–5100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Junttila, M.R.; Puustinen, P.; Niemelä, M.; Ahola, R.; Arnold, H.; Böttzauw, T.; Ala-Aho, R.; Nielsen, C.; Ivaska, J.; Taya, Y.; et al. CIP2A Inhibits PP2A in Human Malignancies. Cell 2007, 130, 51–62. [Google Scholar] [CrossRef] [Green Version]
  20. Soofiyani, S.R.; Hejazi, M.S.; Baradaran, B. The role of CIP2A in cancer: A review and update. Biomed. Pharm. 2017, 96, 626–633. [Google Scholar] [CrossRef]
  21. Khanna, A.; Böckelman, C.; Hemmes, A.; Junttila, M.R.; Wiksten, J.-P.; Lundin, M.; Junnila, S.; Murphy, D.; Evan, G.I.; Haglund, C.; et al. MYC-Dependent Regulation and Prognostic Role of CIP2A in Gastric Cancer. J. Natl. Cancer Inst. 2009, 101, 793–805. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Rajkumar, S.V.; Dimopoulos, M.A.; Palumbo, A.; Blade, J.; Merlini, G.; Mateos, M.-V.; Kumar, S.; Hillengass, J.; Kastritis, E.; Richardson, P.; et al. International Myeloma Working Group updated criteria for the diagnosis of multiple myeloma. Lancet Oncol. 2014, 15, e538–e548. [Google Scholar] [CrossRef] [PubMed]
  23. Durie, B.G.; Salmon, S.E. 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, 842–854. [Google Scholar] [CrossRef] [PubMed]
  24. D’Agostino, M.; Cairns, D.A.; Lahuerta, J.J.; Wester, R.; Bertsch, U.; Waage, A.; Zamagni, E.; Mateos, M.-V.; Dall’Olio, D.; van de Donk, N.W.C.J.; et al. Second Revision of the International Staging System (R2-ISS) for Overall Survival in Multiple Myeloma: A European Myeloma Network (EMN) Report Within the HARMONY Project. J. Clin. Oncol. 2022, 40, 3406–3418. [Google Scholar] [CrossRef]
  25. International Myeloma Working Group. Criteria for the classification of monoclonal gammopathies, multiple myeloma and related disorders: A report of the International Myeloma Working Group. Br. J. Haematol. 2003, 121, 749–757. [Google Scholar] [CrossRef] [Green Version]
  26. U.S. Department of Health and Human Services. Common Terminology Criteria for Adverse Events (CTCAE); Version 5; U.S. Department of Health and Human Services: Washington, DC, USA, 2017. [Google Scholar]
  27. Abdallah, N.; Rajkumar, S.V.; Greipp, P.; Kapoor, P.; Gertz, M.A.; Dispenzieri, A.; Baughn, L.B.; Lacy, M.Q.; Hayman, S.R.; Buadi, F.K.; et al. Cytogenetic abnormalities in multiple myeloma: Association with disease characteristics and treatment response. Blood Cancer J. 2020, 10, 82. [Google Scholar] [CrossRef]
  28. Saadoune, C.; Nouadi, B.; Hamdaoui, H.; Chegdani, F.; Bennis, F. Multiple Myeloma: Bioinformatic Analysis for Identification of Key Genes and Pathways. Bioinform. Biol. Insights 2022, 16. [Google Scholar] [CrossRef]
  29. Misund, K.; Bruinink, D.H.O.; Coward, E.; Hoogenboezem, R.M.; Rustad, E.H.; Sanders, M.A.; Rye, M.; Sponaas, A.-M.; van der Holt, B.; Zweegman, S.; et al. Clonal evolution after treatment pressure in multiple myeloma: Heterogenous genomic aberrations and transcriptomic convergence. Leukemia 2022, 36, 1887–1897. [Google Scholar] [CrossRef]
  30. Junttila, M.R.; Westermarck, J. Mechanisms of MYC stabilization in human malignancies. Cell Cycle 2008, 7, 592–596. [Google Scholar] [CrossRef] [Green Version]
  31. Xing, M.-L.; Lu, Y.-F.; Wang, D.-F.; Zou, X.-Y.; Zhang, S.-X.; Yun, Z. Clinical significance of sCIP2A levels in breast cancer. Eur. Rev. Med. Pharmacol. Sci. 2016, 20, 82–91. [Google Scholar]
  32. Gomes, L.R.; Menck, C.F.M.; Cuervo, A.M. Chaperone-mediated autophagy prevents cellular transformation by regulating MYC proteasomal degradation. Autophagy 2017, 13, 928–940. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Petrusca, D.N.; Mulcrone, P.L.; Macar, D.A.; Bishop, R.T.; Berdyshev, E.; Suvannasankha, A.; Anderson, J.L.; Sun, Q.; Auron, P.E.; Galson, D.L.; et al. GFI1-Dependent Repression of SGPP1 Increases Multiple Myeloma Cell Survival. Cancers 2022, 14, 772. [Google Scholar] [CrossRef] [PubMed]
  34. Dahlström, K.M.; Salminen, T.A. 3D model for Cancerous Inhibitor of Protein Phosphatase 2A armadillo domain unveils highly conserved protein–protein interaction characteristics. J. Theor. Biol. 2015, 386, 78–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wang, J.; Okkeri, J.; Pavic, K.; Wang, Z.; Kauko, O.; Halonen, T.; Sarek, G.; Ojala, P.M.; Rao, Z.; Xu, W.; et al. Oncoprotein CIP 2A is stabilized via interaction with tumor suppressor PP 2A/B56. EMBO Rep. 2017, 18, 437–450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Myant, K.; Qiao, X.; Halonen, T.; Come, C.; Laine, A.; Janghorban, M.; Partanen, J.I.; Cassidy, J.; Ogg, E.-L.; Cammareri, P.; et al. Serine 62-Phosphorylated MYC Associates with Nuclear Lamins and Its Regulation by CIP2A Is Essential for Regenerative Proliferation. Cell Rep. 2015, 12, 1019–1031. [Google Scholar] [CrossRef] [Green Version]
  37. Li, Y.; Wang, K.; Dai, L.; Wang, P.; Song, C.; Shi, J.; Ren, P.; Ye, H.; Zhang, J. HapMap-based study of CIP2A gene polymorphisms and HCC susceptibility. Oncol. Lett. 2012, 4, 358–364. [Google Scholar] [CrossRef] [Green Version]
  38. Preuss, K.-D.; Fadle, N.; Regitz, E.; Held, G.; Pfreundschuh, M. Inactivation of protein-phosphatase 2A causing hyperphosphorylation of autoantigenic paraprotein targets in MGUS/MM is due to an exchange of its regulatory subunits. Int. J. Cancer 2014, 135, 2046–2053. [Google Scholar] [CrossRef]
  39. Puustinen, P.; Rytter, A.; Mortensen, M.; Kohonen, P.; Moreira, J.; Jäättelä, M. CIP2A oncoprotein controls cell growth and autophagy through mTORC1 activation. J. Cell Biol. 2014, 204, 713–727. [Google Scholar] [CrossRef]
  40. Laine, A.; Nagelli, S.G.; Farrington, C.; Butt, U.; Cvrljevic, A.N.; Vainonen, J.P.; Feringa, F.M.; Grönroos, T.J.; Gautam, P.; Khan, S.; et al. CIP2A Interacts with TopBP1 and Drives Basal-Like Breast Cancer Tumorigenesis. Cancer Res. 2021, 81, 4319–4331. [Google Scholar] [CrossRef]
  41. Liu, C.-Y.; Hu, M.-H.; Hsu, C.-J.; Huang, C.-T.; Wang, D.-S.; Tsai, W.-C.; Chen, Y.-T.; Lee, C.-H.; Chu, P.-Y.; Chen, M.-H.; et al. Lapatinib inhibits CIP2A/PP2A/p-Akt signaling and induces apoptosis in triple negative breast cancer cells. Oncotarget 2016, 7, 9135–9149. [Google Scholar] [CrossRef] [Green Version]
  42. Lin, Y.-C.; Chen, K.-C.; Chen, C.-C.; Cheng, A.-L.; Chen, K.-F. CIP2A-mediated Akt activation plays a role in bortezomib-induced apoptosis in head and neck squamous cell carcinoma cells. Oral Oncol. 2012, 48, 585–593. [Google Scholar] [CrossRef] [PubMed]
  43. Bedewy, A.M.L.; Elmaghraby, S.M. CIP2A expression in Bortezomib-treated multiple myeloma. JBUON 2020, 25, 395–400. [Google Scholar]
  44. Nath, S.; Ohlmeyer, M.; Salathe, M.A.; Poon, J.; Baumlin, N.; Foronjy, R.F.; Geraghty, P. Chronic Cigarette Smoke Exposure Subdues PP2A Activity by Enhancing Expression of the Oncogene CIP2A. Am. J. Respir. Cell Mol. Biol. 2018, 59, 695–705. [Google Scholar] [CrossRef] [PubMed]
  45. Slomp, A.; Moesbergen, L.M.; Eldering, E.; Kersten, M.J.; Minnema, M.C.; Peperzak, V. Phosphatase PP2A enhances MCL-1 protein half-life in multiple myeloma cells. Cell Death Dis. 2021, 12, 229. [Google Scholar] [CrossRef] [PubMed]
  46. McDonald, W.J.; Sangster, S.M.; Moffat, L.D.; Henderson, M.J.; Too, C.K. α4 phosphoprotein interacts with EDD E3 ubiquitin ligase and poly(A)-binding protein. J. Cell. Biochem. 2010, 110, 1123–1129. [Google Scholar] [CrossRef]
  47. Alexander-Bryant, A.A.; Dumitriu, A.; Attaway, C.C.; Yu, H.; Jakymiw, A. Fusogenic-Oligoarginine peptide-mediated silencing of the CIP2A oncogene suppresses oral cancer tumor growth in vivo. J. Control. Release 2015, 218, 72–81. [Google Scholar] [CrossRef] [Green Version]
  48. Yang, X.; Zhang, Y.; Liu, H.; Lin, Z. Cancerous Inhibitor of PP2A Silencing Inhibits Proliferation and Promotes Apoptosis in Human Multiple Myeloma Cells. BioMed Res. Int. 2016, 2016, 6864135. [Google Scholar] [CrossRef] [Green Version]
  49. Kang, H.S.; Choi, I. Protein Phosphatase 2A Modulates the Proliferation of Human Multiple Myeloma Cells via Regulation of the Production of Reactive Oxygen Intermediates and Anti-Apoptotic Factors. Cell. Immunol. 2001, 213, 34–44. [Google Scholar] [CrossRef]
  50. Mitsuhashi, S.; Shima, H.; Tanuma, N.; Sasa, S.; Onoe, K.; Ubukata, M.; Kikuchi, K. Protein phosphatase type 2A, PP2A, is involved in degradation of gp130. Mol. Cell. Biochem. 2005, 269, 183–187. [Google Scholar] [CrossRef] [PubMed]
  51. Thangavelu, P.U.; Lin, C.-Y.; Vaidyanathan, S.; Nguyen, T.H.; Dray, E.; Duijf, P.H. Overexpression of the E2F target gene CENPI promotes chromosome instability and predicts poor prognosis in estrogen receptor-positive breast cancer. Oncotarget 2017, 8, 62167–62182. [Google Scholar] [CrossRef] [Green Version]
  52. Brunk, K.; Zhu, M.; Bärenz, F.; Kratz, A.-S.; Haselmann-Weiss, U.; Antony, C.; Hoffmann, I. Cep78 is a novel centriolar protein involved in Plk4-induced centriole overduplication. J. Cell Sci. 2016, 129, 2713–2718. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Gladilin, E.; Ohse, S.; Boerries, M.; Busch, H.; Xu, C.; Schneider, M.; Meister, M.; Eils, R. TGFβ-induced cytoskeletal remodeling mediates elevation of cell stiffness and invasiveness in NSCLC. Sci. Rep. 2019, 9, 7667. [Google Scholar] [CrossRef] [Green Version]
  54. Garcia-Saez, I.; Skoufias, D.A. Eg5 targeting agents: From new anti-mitotic based inhibitor discovery to cancer therapy and resistance. Biochem. Pharmacol. 2020, 184, 114364. [Google Scholar] [CrossRef] [PubMed]
  55. Maniswami, R.R.; Prashanth, S.; Karanth, A.V.; Koushik, S.; Govindaraj, H.; Mullangi, R.; Rajagopal, S.; Jegatheesan, S.K. PLK4: A link between centriole biogenesis and cancer. Expert Opin. Ther. Targets 2017, 22, 59–73. [Google Scholar] [CrossRef]
  56. Wang, J.; Liu, X.; Chu, H.-J.; Li, N.; Huang, L.-Y.; Chen, J. Centromere Protein I (CENP-I) Is Upregulated in Gastric Cancer, Predicts Poor Prognosis, and Promotes Tumor Cell Proliferation and Migration. Technol. Cancer Res. Treat. 2021, 20. [Google Scholar] [CrossRef] [PubMed]
  57. Jiang, F.; Tang, X.; Tang, C.; Hua, Z.; Ke, M.; Wang, C.; Zhao, J.; Gao, S.; Jurczyszyn, A.; Janz, S.; et al. HNRNPA2B1 promotes multiple myeloma progression by increasing AKT3 expression via m6A-dependent stabilization of ILF3 mRNA. J. Hematol. Oncol. 2021, 14, 54. [Google Scholar] [CrossRef]
  58. Das, D.S.; Das, A.; Ray, A.; Song, Y.; Samur, M.K.; Munshi, N.C.; Chauhan, D.; Anderson, K.C. Blockade of Deubiquitylating Enzyme USP1 Inhibits DNA Repair and Triggers Apoptosis in Multiple Myeloma Cells. Clin. Cancer Res. 2017, 23, 4280–4289. [Google Scholar] [CrossRef] [Green Version]
  59. Sagawa, M.; Ohguchi, H.; Harada, T.; Samur, M.K.; Tai, Y.-T.; Munshi, N.C.; Kizaki, M.; Hideshima, T.; Anderson, K.C. Ribonucleotide Reductase Catalytic Subunit M1 (RRM1) as a Novel Therapeutic Target in Multiple Myeloma. Clin. Cancer Res. 2017, 23, 5225–5237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. He, H.; Zheng, C.; Tang, Y. Overexpression of SMC4 predicts a poor prognosis in cervical cancer and facilitates cancer cell malignancy phenotype by activating NF-κB pathway. Hum. Cell 2021, 34, 1888–1898. [Google Scholar] [CrossRef]
  61. Chen, Q.; Wan, X.; Chen, Y.; Liu, C.; Gu, M.; Wang, Z. SGO1 induces proliferation and metastasis of prostate cancer through AKT-mediated signaling pathway. Am. J. Cancer Res. 2019, 9, 2693–2705. [Google Scholar]
  62. Zhang, X.; Xu, Q.; Li, E.; Shi, T.; Chen, H. CEP55 predicts the poor prognosis and promotes tumorigenesis in endometrial cancer by regulating the Foxo1 signaling. Mol. Cell. Biochem. 2022. [Google Scholar] [CrossRef] [PubMed]
  63. Jiang, N.; Liao, Y.; Wang, M.; Wang, Y.; Wang, K.; Guo, J.; Wu, P.; Zhong, B.; Guo, T.; Wu, C. BUB1 drives the occurrence and development of bladder cancer by mediating the STAT3 signaling pathway. J. Exp. Clin. Cancer Res. 2021, 40, 378. [Google Scholar] [CrossRef]
  64. Chou, C.-W.; Hsieh, Y.-H.; Ku, S.-C.; Shen, W.-J.; Anuraga, G.; Ta, H.D.K.; Lee, K.-H.; Lee, Y.-C.; Lin, C.-H.; Wang, C.-Y.; et al. Potential Prognostic Biomarkers of OSBPL Family Genes in Patients with Pancreatic Ductal Adenocarcinoma. Biomedicines 2021, 9, 1601. [Google Scholar] [CrossRef] [PubMed]
  65. Pan, J.; Lei, X.; Mao, X. Identification of KIF4A as a pan-cancer diagnostic and prognostic biomarker via bioinformatics analysis and validation in osteosarcoma cell lines. PeerJ 2021, 9, e11455. [Google Scholar] [CrossRef] [PubMed]
  66. Shi, Y.; Ge, C.; Fang, D.; Wei, W.; Li, L.; Wei, Q.; Yu, H. NCAPG facilitates colorectal cancer cell proliferation, migration, invasion and epithelial–mesenchymal transition by activating the Wnt/β-catenin signaling pathway. Cancer Cell Int. 2022, 22, 119. [Google Scholar] [CrossRef]
  67. Zeng, X.; Xu, W.; Tong, J.; Liu, J.; Zhang, Z.; Liu, M.; Wu, C.; Yu, Q.; Ye, C.; Wu, C.; et al. SPAG5 as a novel biomarker and potential therapeutic target via regulating AKT pathway in multiple myeloma. Leuk. Lymphoma 2022, 63, 2565–2572. [Google Scholar] [CrossRef]
  68. Li, B.; Xiao, Q.; Shan, L.; Song, Y. NCAPH promotes cell proliferation and inhibits cell apoptosis of bladder cancer cells through MEK/ERK signaling pathway. Cell Cycle 2022, 21, 427–438. [Google Scholar] [CrossRef]
  69. Xie, B.; Wang, S.; Jiang, N.; Li, J.J. Cyclin B1/CDK1-regulated mitochondrial bioenergetics in cell cycle progression and tumor resistance. Cancer Lett. 2018, 443, 56–66. [Google Scholar] [CrossRef]
  70. Kim, J.H.; Youn, Y.; Lee, J.-C.; Kim, J.; Hwang, J.-H. Involvement of the NF-κB signaling pathway in proliferation and invasion inhibited by Zwint-1 deficiency in Pancreatic Cancer Cells. J. Cancer 2020, 11, 5601–5611. [Google Scholar] [CrossRef]
  71. Kim, J.H.; Shim, J.; Ji, M.-J.; Jung, Y.; Bong, S.M.; Jang, Y.-J.; Yoon, E.-K.; Lee, S.-J.; Kim, K.G.; Kim, Y.H.; et al. The condensin component NCAPG2 regulates microtubule–kinetochore attachment through recruitment of polo-like kinase 1 to kinetochores. Nat. Commun. 2014, 5, 4588. [Google Scholar] [CrossRef] [Green Version]
  72. Yang, D.; Ji, F.; Li, Y.; Jiao, Y.; Fang, X. GPSM2 Serves as an Independent Prognostic Biomarker for Liver Cancer Survival. Technol. Cancer Res. Treat. 2020, 19. [Google Scholar] [CrossRef]
  73. Adams, R.; Wheatleya, S.; Gouldsworthy, A.; Kandels-Lewis, S.; Carmena, M.; Smythe, C.; Gerloff, D.; Earnshaw, W. INCENP binds the Aurora-related kinase AIRK2 and is required to target it to chromosomes, the central spindle and cleavage furrow. Curr. Biol. 2000, 10, 1075–1078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Wang, T.; Liu, D.; Wang, L.; Liu, M.; Duan, W.; Yi, J.; Yi, Y. DNA Repair Genes Are Associated with Subtype Classification, Prognosis, and Immune Infiltration in Uveal Melanoma. J. Oncol. 2022, 2022, 1965451. [Google Scholar] [CrossRef]
  75. Müller, R.; Misund, K.; Holien, T.; Bachke, S.; Gilljam, K.M.; Våtsveen, T.K.; Rø, T.B.; Bellacchio, E.; Sundan, A.; Otterlei, M. Targeting Proliferating Cell Nuclear Antigen and Its Protein Interactions Induces Apoptosis in Multiple Myeloma Cells. PLoS ONE 2013, 8, e70430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Han, Y.H.; Wan, Y.; Xiong, H.; Sun, G.L. Structural maintenance of chromosomes 2 is identified as an oncogene in bladder cancer in vitro and in vivo. Neoplasma 2020, 67, 364–370. [Google Scholar] [CrossRef]
  77. Reale, A.; Khong, T.; Mithraprabhu, S.; Savvidou, I.; Hocking, J.; Bergin, K.; Ramachandran, M.; Chen, M.; Dammacco, F.; Ria, R.; et al. TOP2A expression predicts responsiveness to carfilzomib in myeloma and informs novel combinatorial strategies for enhanced proteasome inhibitor cell killing. Leuk. Lymphoma 2020, 62, 337–347. [Google Scholar] [CrossRef] [PubMed]
Figure 3. Box–whisker graph presenting a comparison of expression of the KIAA1524 gene depending on genotypes of the SNP of this gene (A). Distribution plot presenting correlation between expression of KIAA1524 gene expression of the MYC gene (B).
Figure 3. Box–whisker graph presenting a comparison of expression of the KIAA1524 gene depending on genotypes of the SNP of this gene (A). Distribution plot presenting correlation between expression of KIAA1524 gene expression of the MYC gene (B).
Cells 12 01029 g003
Figure 4. Heat map presenting correlation matrix, including CIP2A/KIAA1524 and selected genes involved in the biology or pathomechanism of MM, created based on data retrieved from the CoMMpass project.
Figure 4. Heat map presenting correlation matrix, including CIP2A/KIAA1524 and selected genes involved in the biology or pathomechanism of MM, created based on data retrieved from the CoMMpass project.
Cells 12 01029 g004
Figure 5. Heat map presenting correlation matrix including 30 genes demonstrating the strongest correlation with CIP2A/KIAA1524, created based on data retrieved from the CoMMpass project.
Figure 5. Heat map presenting correlation matrix including 30 genes demonstrating the strongest correlation with CIP2A/KIAA1524, created based on data retrieved from the CoMMpass project.
Cells 12 01029 g005
Table 1. Baseline characteristics of the study group.
Table 1. Baseline characteristics of the study group.
Variablen = 128 (100%)
Sex
Men64 (50%)
Women64 (50%)
Age
<6558 (45.3%)
≥6570 (54.7%)
Myeloma type
IgG72 (56.2.%)
IgA29 (22.7%)
Light chains25 (19.5%)
Non-secretory1 (0.8%)
Plasmocytoma1 (0.8%)
Light chain type
Kappa 80 (63%)
Lambda 47 (37%)
No data: n = 1
Durie–Salmon stage
I 12 (9.4%)
II 16 (12.5%)
III 100 (78.1%)
ISS stage
133 (25.8%)
236 (28.1%)
359 (46.1%)
Renal function
A—creatinine < 2 mg/dL102 (79.7%)
B—creatinine ≥ 2 mg/dL26 (20.3%)
Performance status
011 (8.6%)
149 (38.3%)
251 (39.8%)
315 (11.7%)
42 (1.6%)
Body weight loss
No 47 (50.5%)
Yes 46 (49.5%)
No data: n = 35
5%14 (30.4%)
10%32 (69.6%)
Anemia grade before treatment (WHO)
Absent or I°65 (50.8%)
II°, III° or IV°63 (49.2%)
Treatment protocol
CTD27 (21.1%)
V(C)D58 (45.3%)
VTD43 (33.6%)
AHSCT82 (64.1%)
No46 (35.9%)
Yes
del 17p/TP53
Absent65 (74.7%)
Present22 (25.3%)
No data: n = 41
t(4;14) IGH/FGFR3
Absent75 (85.2%)
Present13 (14.8%)
No data: n = 40
t(11;14) IGH/CCND1
Absent87 (88.8%)
Present11 (11.2%)
No data: n = 30
t(14;16) IGH/MAF
Absent82 (83.2%)
Present6 (6.8%)
No data: n = 40
Other IGH rearrangement
Absent84 (85.7%)
Present14 (14.3%)
No data: n = 30
amp 1q21/CKS1B
Absent11 (47.8%)
Present12 (52.2%)
No data: n = 105
CTD—cyclophosphamide, thalidomide, dexamethsone; ISS—Multiple Myeloma International Staging System; WHO—World Health Organization; V(C)D—bortezomib, (cyclophosphamide), dexamethasone; VTD—bortezomib, thalidomide, dexamethasone.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Szudy-Szczyrek, A.; Mlak, R.; Mazurek, M.; Krajka, T.; Chocholska, S.; Bitkowska, P.; Jutrzenka, M.; Szczyrek, M.; Homa-Mlak, I.; Krajka, A.; et al. The TT Genotype of the KIAA1524 rs2278911 Polymorphism Is Associated with Poor Prognosis in Multiple Myeloma. Cells 2023, 12, 1029. https://doi.org/10.3390/cells12071029

AMA Style

Szudy-Szczyrek A, Mlak R, Mazurek M, Krajka T, Chocholska S, Bitkowska P, Jutrzenka M, Szczyrek M, Homa-Mlak I, Krajka A, et al. The TT Genotype of the KIAA1524 rs2278911 Polymorphism Is Associated with Poor Prognosis in Multiple Myeloma. Cells. 2023; 12(7):1029. https://doi.org/10.3390/cells12071029

Chicago/Turabian Style

Szudy-Szczyrek, Aneta, Radosław Mlak, Marcin Mazurek, Tomasz Krajka, Sylwia Chocholska, Paulina Bitkowska, Marta Jutrzenka, Michał Szczyrek, Iwona Homa-Mlak, Andrzej Krajka, and et al. 2023. "The TT Genotype of the KIAA1524 rs2278911 Polymorphism Is Associated with Poor Prognosis in Multiple Myeloma" Cells 12, no. 7: 1029. https://doi.org/10.3390/cells12071029

APA Style

Szudy-Szczyrek, A., Mlak, R., Mazurek, M., Krajka, T., Chocholska, S., Bitkowska, P., Jutrzenka, M., Szczyrek, M., Homa-Mlak, I., Krajka, A., Małecka-Massalska, T., & Hus, M. (2023). The TT Genotype of the KIAA1524 rs2278911 Polymorphism Is Associated with Poor Prognosis in Multiple Myeloma. Cells, 12(7), 1029. https://doi.org/10.3390/cells12071029

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop