Skip to main content
Have a personal or library account? Click to login
Autophagy-related genes affect the survival of multiple myeloma patients depending on chromosomal abnormality Cover

Autophagy-related genes affect the survival of multiple myeloma patients depending on chromosomal abnormality

Open Access
|Jun 2023

Figures & Tables

Table 1.

Outline of the study

Study dataset (548 MM patients, GSE24080)Validation dataset 1576 MM patients, GSE4204, GSE4452, GSE4581, and GSE2658)
Demographic properties of patients (Table 2)
Survival Time According to Patients with/without Chromosomal Abnormality(Table 3)
Univariate and multivariate Cox regression analyses and focus on the significant effect of ABN on survival of MM patients (Figure 1, Tables 4 and 5)
Prognosis-related (survival-triggering or hazardous) ARG identification in patients with or without ABN (Figures 2 and 3)
Table 6 (for patients without ABN), Table 7 (for patients with ABN)
GOrilla Analysis (Figure 4)
Validation of prognosis-related ARGs in MM patients with ABN (only for 1q21 amplification) (Table 8)

[i] ARG, autophagy-related genes; ABN, chromosomal abnormality; GOrilla, Gene Ontology enRIchmentanaLysis and visuaLizAtion; MM: multiple myeloma.

Table 2.

Demographic properties in MM patients

Parameters (n = 548)n (%)
Age (years); mean ± SD; 57.04 ± 9.37
Gender
  Male217 (39.59)
  Female331 (60.40)
ABN
  Yes349 (63.68)
  No199 (36.31)
Treatment protocol
  TT2334 (60.94)
  TT3214 (39.05)
Race
  Caucasian486 (88.68)
  Other62 (11.31)

[i] TT2 involves thalidomide treatment; TT3 involves bortezomib, thalidomide, and dexamethasone treatment. ABN, chromosomal abnormality; SD, standard deviation; MM: multiple myeloma; TT2, Total Therapy 2; TT3, Total Therapy 3.

Figure 1.

The significant effect of ABN on the EFS of MM patients. The effect of independent variables on both EFS of MM patients is shown using Kaplan–Meier survival plots as indicated in (A) ABN; EFS (B) sex and EFS; (C) age and EFS; (D) race and EFS; and (E) treatment protocol. ABN, chromosomal abnormality; EFS, event-free survival; MM: multiple myeloma.

Table 3.

Survival time according to patients with or without ABN

Patients (months)OS (mean ± SD)EFS (mean ± SD)
All patients (n = 548)48.23 ± 22.3642.61 ± 22.15
Patients with ABN (n = 349)42.70 ± 24.7937.63 ± 23.96
Patients without ABN (n = 199)50.56 ± 20.5044.93 ± 20.50

[i] ABN, chromosomal abnormality; EFS, event-free survival; OS, overall survival; SD, standard deviation.

Table 4.

Univariate Cox regression results of OS and EFS analyses (n = 548)

Dependent variablesEFSOS
Independent variablesHRLower 95% CIUpper 95% CIPHRLower 95% CIUpper 95% CIP
Age1.01260.99891.0270.07251.0222*1.0051.040.0107
Sex0.99000.76521.2810.93901.00180.73341.3690.9910
Race1.45380.93722.2550.09481.02790.63671.6570.9120
Treatment protocol1.6409**1.2032.2380.00171.19820.83161.7270.3320
ABN1.8009***1.3982.3190.00005192.2054***1.6242.9950.0000408

* P < 0.05;

** P < 0.01;

*** P < 0.001.

ABN, chromosomal abnormality; CI, confidence interval; EFS, event-free survival; HR, Hazard ratio; OS, overall survival.

Table 5.

Multivariate Cox regression results of overall and event free survival analysis

Dependent variablesEFSOS
Independent variablesHazard ratioLower 95% C.IUpper 95% C.IPHazard ratioLower 95% C.IUpper 95% C.IP
Age1.01130.99711.0260.12021.02341.00571.0410.0093
Sex0.99620.76891.2910.97731.03080.75371.410.8489
Race1.40490.89692.2010.13750.90120.55151.4730.6781
Treatment protocol1.6626**1.21682.2720.00141.20650.83571.7420.3163
Chromosomal abnormality1.7627***1.36762.2720.000122.1790***1.60382.9610.0000633

* P < 0.05;

** P < 0.01;

*** P < 0.001.

ABN, chromosomal abnormality.

Table 6.

Univariate and multivariate Cox regression results of prognosis-associated ARG in MM patients without ABN

Univariate and multivariate Cox regression results of prognosis-associated ARG in MM patients without ABNGeneUnivariate Cox regressionMultivariate Cox regression
HRLower 95% CIUpper 95% CIPHRLower 95% CIUpper 95% CIP
APOL10.80962*0.68620.95530.01240.828576*0.69960.98130.02935
ARSA0.83103**0.73540.93910.0030.845735**0.74620.95860.00873
ATG16L21.2505*1.0131.5430.03741.30139*1.05071.6120.01583
ATG9A0.6562***0.51140.84190.0009230.705951**0.55060.90510.00603
ATG1010.6862**0.51560.91320.00980.703135*0.51460.96080.02703
CDKN1A0.7468*0.58490.95360.01920.74999*0.58300.96480.02519
DNAJB10.7074*0.53160.94140.01760.72321*0.54190.96530.02781
DRAM10.7561*0.57560.99320.04460.74669*0.57290.97330.030754
FKBP1B0.82008**0.72180.93180.002330.830138**0.72930.94490.00483
FOXO30.7109**0.55850.9050.005590.70415**0.55160.8990.00488
MAPK90.6577*0.45650.94740.02440.649662*0.44910.93980.022037
NCKAP10.8088***0.71980.90890.0003620.82788**0.73600.93120.00164
RAB240.7536*0.59580.9530.01820.727826**0.57230.92560.009599

* P < 0.05;

** P < 0.01;

*** P < 0.001.

ABN, chromosomal abnormality; ARG, autophagy-related genes; CI, confidence interval; HR, Hazard ratio; MM, Multiple myeloma

Table 7.

Univariate and Multivariate Cox regression results of prognosis-associated ARG in MM patients with ABN

GeneUnivariate Cox regressionMultivariate Cox regression
HRLower 95% CIUpper 95% CIPHRLower 95% CIUpper 95% CIP
ARNT1.5304*1.0642.2020.02191.498746*1.02382.1940.0374
ATG4B0.5552*0.31720.9720.03950.554451*0.31670.97080.0391
ATG4D0.5521***0.38910.78330.00080.568641**0.39190.8250.00295
ATIC1.8684**1.2362.8240.003021.857592**1.21472.8410.00427
BIRC51.30634***1.1631.4676.22e–061.343899***1.19011.5181.87e–06
CAPN101.21016*1.0081.4530.0411.26361*1.04701.5250.0147
CASP31.4334*1.0132.0290.04221.505861*1.04722.1650.0272
CDKN1A0.7103**0.56450.89380.003530.726627**0.57870.91240.00598
CDKN2A1.4842*1.0562.0850.02281.54526*1.05252.2690.0263
CXCR40.73568***0.62890.86060.0001250.765047**0.64860.90250.00148
DNAJB90.6088*0.3890.95290.02990.593613*0.37660.93560.0246
EIF2S11.6981**1.1562.4940.006921.79448**1.18722.7120.00553
EIF4EBP11.34577***1.1281.6050.0009551.312952**1.10061.5660.00249
EIF4G12.3264***1.4933.6260.0001932.206738***1.39483.4910.000721
FADD1.5694*1.0632.3170.02341.61684*1.08942.4000.01708
FKBP1A1.8452**1.1872.8680.006471.88003**1.19972.9460.00588
FKBP1B0.82577**0.72430.94140.00420.851602*0.74540.9730.0181
FOXO10.6549*0.47390.90510.01030.645133**0.46680.89170.00794
FOXO30.7363*0.56570.95830.02280.72246*0.55010.94890.01942
GABARAP0.4739**0.28960.77530.002950.464234**0.28090.76720.00276
GNAI31.5955*1.0822.3530.01841.542438*1.04152.2840.0305
HDAC61.7322*1.0592.8340.02881.645421*1.00512.6940.0477
HGS1.8098*1.142.8730.01191.696411*1.06192.7100.027
HSP90AB11.6281*1.0882.4360.01781.594093*1.05672.4050.0262
IL240.82369*0.68690.98770.03630.795624*0.65790.96220.0184
IRGM0.81103*0.67370.97630.02690.823283*0.68340.99190.0408
ITGA60.67844***0.57090.80631.06e–050.686741**0.57700.81742.35e–05
ATG130.4507***0.28660.70870.0005580.425973**0.27210.66690.00019
LAMP20.7123*0.54950.92340.01040.73084*0.56190.95050.0193
MAPK11.496*1.0272.180.03591.481801*1.00792.1790.0455
MBTPS21.445*1.0392.0090.02871.435077*1.01822.0230.0391
NCKAP10.8715*0.78230.97090.01260.8927*0.79870.99780.0456
PARP11.8204***1.2862.5770.0007291.938906***0.85402.0940.000271
PRKCD0.7481*0.57480.97350.03080.751456*0.57410.98370.03755
SIRT20.5835**0.40890.83270.002980.662074*0.45210.96960.0341
TNFSF100.88013*0.78410.9880.03040.87052*0.77590.97670.01824
TP530.7408**0.60790.90290.002960.726943**0.59470.88870.00186
TSC21.5559*1.0992.2020.01261.4891644*1.04992.1120.0255
VAMP30.571***0.41680.78230.0004850.596231**0.43170.82350.0017
WDFY31.1406*1.0171.280.0251.132911*1.01271.2670.0293
WDR45L0.5495*0.32830.91980.02270.57379*0.33960.96960.0379

* P < 0.05;

** P < 0.01;

*** P < 0.001.

ABN, chromosomal abnormality; ARG, autophagy-related genes; CI, confidence interval; HR, Hazard ratio; MM, Multiple myeloma.

Figure 2.

Survival analysis of selected ARGs in MM patients without ABN; Kaplan–Meier plots of 6 prognosis-related ARGs in MM patients without ABN. (A) ARSA, (B) ATG9A, (C) FKBP1B, (D) FOXO3, (E) NCKAP1, and (F) RAB24.

ARSA, arylsulfatase A; ATG9A, ATG9 autophagy-related 9 homolog A (Saccharomyces cerevisiae); FKBP1B, FK506 binding protein 1B, 12.6 kDa; FOXO3, forkhead box O3; NCKAP1, NCK-associated protein 1; RAB24, RAB24, member RAS oncogene family. ABN, chromosomal abnormality; ARGs, autophagy-related genes.

Figure 3.

Survival analysis of selected ARGs in MM patients with ABN; Kaplan–Meier plots of 6 prognosis-related ARGs in MM patients with ABN. (A) ATG4D (S. cerevisiae); (B) ATIC; (C) BIRC5; (D) CDKN1A (p21, Cip1); (E) CXCR4; (F) EIF2S1, 35 kDa; (G) EIF4EBP1; (H) EIF4G1; (I) FKBP1A, 12 kDa; (J) FOXO1; (K) GABARAP; (L) ITGA6; (M) ATG13, KIAA0652; (N) PARP1; (O) TP53; (P) VAMP3 (cellubrevin).

ABN, chromosomal abnormality; ARGs, autophagy-related genes; ATG4D, TG4 autophagy-related 4 homolog D; ATG13, Autophagy-Related Protein 13; ATIC, 5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase–IMP cyclohydrolase; BIRC5, baculoviral IAP repeat-containing 5; CDKN1A, cyclin-dependent kinase inhibitor 1A; CXCR4, chemokine (C-X-C motif) receptor 4; EIF2S1, eukaryotic translation initiation factor 2, subunit 1 alpha; EIF4EBP1, eukaryotic translation initiation factor 4E binding protein 1; EIF4G1, eukaryotic translation initiation factor 4 gamma, 1; FKBP1A, FK506 binding protein 1A; FOXO1, forkhead box O1; GABARAP, GABA(A) receptor-associated protein; ITGA6, integrin, alpha 6; MM: multiple myeloma; PARP1, poly (ADP-ribose) polymerase 1; TP53, tumor protein p53; VAMP3, vesicle-associated membrane protein 3.

Figure 4.

GOrilla biological process in MM patients with ABN.

ABN, chromosomal abnormality; GOrilla, Gene Ontology enRIchmentanaLysis and visuaLizAtion; MM: multiple myeloma.

Table 8.

Validation of prognosis-related ARGs in MM patients with ABN

GeneUnivariate Cox regression
HRLower 95% CIUpper 95% CIP
ARNT0.78874**0.66720.93240.00543
ATG4B2.1821***1.513.1543.29e–05
ATG4D0.2839***0.19070.42265.56e–10
ATIC3.5622***2.4855.1074.72e–12
CAPN100.4837***0.36650.63832.87e–07
CASP30.84905*0.74090.9730.0186
CDKN1A0.76380***0.66580.87620.00012
CDKN2A1.51399***1.3561.691.55e–13
CXCR40.5576***0.45680.68079.5e–09
DNAJB91.7560**1.1932.5840.00427
EIF4EBP10.81096***0.7260.90590.000207
EIF4G10.6512**0.49110.86330.00287
FADD1.6059**1.1792.1880.00268
FKBP1A0.3130***0.20930.46791.5e–08
FKBP1B0.6616***0.52280.83740.000589
FOXO10.3332***0.240.46265.25e–11
FOXO30.7879*0.62650.99080.0415
GNAI30.79854**0.67320.94720.00981
HGS0.70614***0.62540.79731.98e–08
HSP90AB10.76746***0.69080.85268.12e–07
IL242.5083***1.6753.7567.98e–06
IRGM2.3205***1.6643.2367.07e–07
ITGA64.0435***2.8475.7446.11e–15
ATG130.3571***0.22540.56581.16e–05
LAMP22.4952***1.7743.5091.48e–07
MAPK12.0933***1.3953.1420.000363
TNFSF101.2571*1.0211.5480.0314
TP531.25520*1.0381.5180.019
TSC21.4152*1.041.9250.027

* P < 0.05;

** P < 0.01;

*** P < 0.001.

ABN, chromosomal abnormality; ARGs, autophagy-related genes; CI, confidence interval; HR, Hazard ratio; MM, Multiple myeloma

DOI: https://doi.org/10.2478/abm-2022-0028 | Journal eISSN: 1875-855X | Journal ISSN: 1905-7415
Language: English
Page range: 249 - 264
Published on: Jun 16, 2023
Published by: Chulalongkorn University
In partnership with: Paradigm Publishing Services

© 2023 Gizem Ayna Duran, Yasemin Benderli Cihan, published by Chulalongkorn University
This work is licensed under the Creative Commons Attribution 4.0 License.