BioMed Research International

BioMed Research International / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 7629456 | https://doi.org/10.1155/2020/7629456

S. Zmorzyński, S. Popek-Marciniec, W. Styk, M. Wojcierowska-Litwin, I. Korszeń-Pilecka, A. Szudy-Szczyrek, S. Chocholska, M. Hus, A. A. Filip, "The Impact of the NOD2/CARD15 Variant (3020insC) and PSMA6 Polymorphism (-8C>G) on the Development and Outcome of Multiple Myeloma", BioMed Research International, vol. 2020, Article ID 7629456, 15 pages, 2020. https://doi.org/10.1155/2020/7629456

The Impact of the NOD2/CARD15 Variant (3020insC) and PSMA6 Polymorphism (-8C>G) on the Development and Outcome of Multiple Myeloma

Academic Editor: Francesco Di Raimondo
Received27 Mar 2020
Accepted25 May 2020
Published08 Jun 2020

Abstract

Introduction. Multiple myeloma (MM) is a hematological malignancy characterized by genetic variety. The 3020insC variant of the NOD2/CARD15 gene results in the upregulation of proinflammatory cytokines. Chronic inflammation and abnormal function of the proteasome system may lead to MM development. The polymorphism (-8C>G) in the PSMA6 gene affects proteasome activity. The aim of our study was to analyze the possible relationship of NOD/CARD15 and PSMA6 genes with the risk of development and outcome of MM, as well as the sensitivity to bortezomib (proteasome inhibitor) in cell cultures derived from MM patients. Objects and Methods. Genomic DNA from 100 newly diagnosed MM patients and 100 healthy blood donors was analyzed by methods such as PCR-RFLP (for PSMA6 genotyping) and automated DNA sequencing (for NOD2/CARD15 genotyping). In a subgroup of 50 MM patients, nucleated bone marrow cells were treated with bortezomib in vitro. Results. Patients with PSMA6 CG+GG genotypes had higher chances for progressive disease (, 95% CI 1.07-23.16, ), shorter overall survival taking into account the type of treatment (), and increased risk of death due to MM at the level of tendency (, 95% CI 1.02-21.97, ). The presence of NOD2/CARD15 3020insC decreased the risk of renal dysfunction in MM (, 95% CI 0.07-0.74, ). The analyzed changes in NOD2/CARD15 and PSMA6 genes did not impact the MM risk. In an in vitro study, bortezomib increased the number of apoptotic cells at 8 nM and 12 nM between wild-type and 3030insC variants of NOD2/CARD15 ( and , respectively). Conclusion. The presented results suggest a possible impact of PSMA6 CG+GG genotypes on the MM outcome and the association of the NOD2/CARD15 variant with bortezomib in vitro sensitivity.

1. Introduction

Multiple myeloma (MM) is a hematological, heterogeneous malignancy characterized by infiltration of bone marrow by clonal plasma cells [1]. Inflammatory mediators, such as nuclear factor κB (NF-κB), are critical in the malignant transformation [2, 3]. Chronic inflammation with elevated NF-κB can support the development of hematological malignancies, including MM [4]. Activation of this pathway is enhanced by nucleotide oligomerization domain 2 (NOD2 also known as CARD15) [5]. This protein is encoded by the NOD2/CARD15 gene (locus 16q12.1), is a component of the host defence against intracellular pathogens, and plays an important role in immune system function. The NOD2/CARD15 gene consists of 18 exons, and its expression is observed in monocytes, macrophages, granulocytes, and some epithelial cells [5, 6]. The most common variant in the NOD2/CARD15 gene is present in the form of 3020insC, which leads to a stop codon in exon 11. The result of a nonsense change is the production of truncated NOD2 protein with the last 33 amino acids missing (Leu1007fsinsC) [5]. The variant leads to loss-of-gene function, and truncated protein has a reduced ability to recognize bacterial muramyl dipeptide and consequently to stimulate NF-κB response [5, 7]. On the other hand, the absence of wild-type NOD2 protein results in the upregulation of proinflammatory cytokines [8]. The 3020insC variant of NOD2/CARD15 also predisposes cells to many types of common cancers, mainly solid tumors—for example, ovarian cancer, breast cancer, and colorectal cancer [7, 9, 10].

The processes like apoptosis, proliferation, and inflammation are controlled in a nondirective way by the proteasome system [11]. It is a multicatalytic proteinase complex. The 20S proteasome is composed of seven α-subunits and ten β-subunits. The human proteasome subunit α6 (or p27K) is encoded by the PSMA6 gene (locus 14q13), which contains single-nucleotide polymorphism (SNP) in 5-UTR (rs10489900). This SNP is caused by substitution (-8C>G) where nucleotide with cytosine (C allele) is replaced by guanine (G allele). The G allele is associated with a higher transcription rate of the PSMA6 gene in human B cells in comparison to the C allele [12]. The proteasome-dependent pathway and higher expression of the PSMA6 gene may activate NF-κB, which leads to inflammation and transformation to MM [13].

The proteasome system consists of the 26S proteasome (including two subunits of 20S and 19S) and immunoproteasome [14]. The 26S proteasome is involved in the degradation of proteins regulating the cell cycle, cell migration, proliferation, angiogenesis, and apoptosis [15]. These proteins play a role in malignant transformation. Therefore, proteasome inhibitors preventing angiogenesis and cell migration and decreasing inflammatory response have been developed and tested for antitumor response [14, 15]. In MM patients, significant improvements in survival have been obtained, among others, by the introduction of bortezomib, a 26S proteasome inhibitor, as well as an anti-inflammatory agent—thalidomide [16]. Inhibition of the proteasome pathway in cancer cells leads to apoptosis via the activation of unfolded or abnormal protein response [17].

Taking the above into account, we hypothesized that variants in the NOD2/CARD15 gene and PMSA6 polymorphism may affect the development and outcome of MM. In this paper, we present an analysis for the studied genetic variants individually and in different combinations to evaluate the role of these genes in the MM susceptibility, survival, and response to treatment with thalidomide and/or bortezomib. Furthermore, the aim of our research was to analyze the studied DNA changes in the sensitivity prediction in cell cultures derived from MM patients. To our knowledge, the presented results were not previously published by other authors.

2. Material and Methods

2.1. Patients and Samples

For the study, bone marrow aspirates and peripheral blood samples were collected from 100 newly diagnosed patients with MM (with mean age 65.32 years), who were hospitalized at the Chair and Department of Hematooncology and Bone Marrow Transplantation, Medical University of Lublin, in years 2013-2019. The study obtained positive opinions from the Bioethics Committee of the Medical University of Lublin (nos. KE-0254/165/2013 and KE-0254/337/2016), according to the ethical standards established by the Helsinki Declaration. The research material was collected upon all patients () and healthy blood donors () who provided written informed consent.

Control samples were made of peripheral blood obtained from 100 healthy blood donors (50 males/50 females, with mean age 34.4 years) attending the Regional Blood Donation and Blood Treatment Center in Kielce, Poland.

The general characteristics of MM patients are shown in Table 1. The inclusion and exclusion criteria are described in Table 2.


VariablesMM patients ()

Sex
 Male53
 Female47
Type of MM
 IgG56
 IgA25
 Light chain19
Stage according to the international staging system
 I28
 II30
 III42
Renal failure
 No82
 Yes18
The stage of chronic kidney disease (grade)
 G130
 G228
 G3A16
 G3B12
 G47
 G57
Cytogenetic changes
 del(17p13.1)14
 t(4;14)15
 t(14;16)2
Anemia grade before treatment (WHO)
 Absent28
 I—mild32
 II—moderate30
 III—severe10
Chemotherapy
 Cyclophosphamide, thalidomide, dexamethasone (CTD)49
 Velcade, cyclophosphamide, dexamethasone (VCD)29
 Velcade, thalidomide, dexamethasone (VTD)22

At diagnosis. M: mean; SD: standard deviation.

Inclusion criteria
MM patients(i) Newly diagnosed MM patients who are required to start active treatment according to the IMWG published in 2014
(ii) Signed informed consent
(iii) Years after 18
(iv) Patients must have measurable disease, defined as follows:
(a) For secretory MM, measurable disease is defined as the presence of quantifiable monoclonal component (≥0.5 g/dl)
(b) For poor secretory or nonsecretory MM, the level of the affected serum free light chain must be ≥10 mg/dl (≥100 mg/l, with an abnormal free light chain ratio)
(v) Eastern Cooperative Oncology Group (ECOG) performance
(vi) Life expectancy more than 3 months
Control group(i) Age: 18 and older
(ii) Signed informed consent

Exclusion criteria
MM patients(i) Active smoldering MM
(ii) Active plasma cell leukemia
(iii) Documented systemic amyloid light chain amyloidosis
(iv) Active central nervous system involvement with MM
(v) Active other hematologic malignancies or solid tumors
Control group(i) Known to be infected with HIV, syphilis, tuberculosis, hepatitis B, or hepatitis C
(ii) A condition in which repeated blood draws or injections pose more than minimal risk for the subject such as hemophilia, other severe coagulation disorders, or significantly impaired venous access
(iii) A condition that requires active medical intervention or monitoring to avert serious danger to the participant’s health or well-being

Therapeutic induction regimens consisted of thalidomide and/or bortezomib combined with steroids and/or cyclophosphamide. 37 MM patients underwent autologous hematopoietic stem cell transplantation (auto-HSCT). Response to treatment was evaluated according to the International Myeloma Working Group guidelines and classified as stringent complete remission (sCR), complete response (CR), very good partial response (VGPR), partial response (PR), minimal response (MR), stable disease (SD), or progressive disease (PD) as described elsewhere [18, 19]. Overall survival (OS) encompassed time from diagnosis until relapse, progression, death due to a tumor effect, or the last follow-up and time from diagnosis until death by any cause or the last follow-up. The median follow-up time of MM patients enrolled in the study was 18 months. Progression-free survival (PFS) was estimated as the time elapsed between treatment initiation and tumor progression or death from any cause [20].

Peripheral blood was used for DNA isolation and studied genetic variant determination.

Cell cultures were established from bone marrow aspirates to carry out the research associated with cIg-FISH () and bortezomib treatment (). In an in vitro bortezomib study, MM patients without chromosomal aberrations were included.

2.2. DNA Isolation

DNA isolation from peripheral blood was performed using a commercial kit (Qiagen, Germany) according to the manufacturer’s procedure. The concentration and quality of DNA were checked using a NanoDrop device (Thermo Fisher Scientific, USA).

2.3. PSMA6 Genotyping

For the analysis of the PSMA6 polymorphism, the PCR restriction fragment length polymorphism (PCR-RFLP) method was applied according to the validated protocol of Bachmann et al. [21]. A PSMA6 gene fragment length of 100 bp was amplified by PCR using the following primers: (i)Forward 5-CTG GTG CGG GAG CTA CGG G-3(ii)Reverse 5-AAT GGT AAT GTG GCG GTC AAA AC-3

Each PCR mixture (25 μl) contained 100 ng genomic DNA and PCR buffer (Clontech), dNTP mixture (0.25 mM), HD polymerase (Clontech), and primers (10 μM of each). The touchdown PCR method was used. The mixture was heated at 95°C for 5 min and underwent 14 cycles of amplification: denaturation at 95°C for 30 s, annealing at 64.5°C for 20 s (-0.5°C/per cycle), and elongation at 72°C for 20 s. After 14 cycles, the mixture underwent 20 cycles with a constant temperature of 57.5°C. The denaturation and elongation temperatures and times were the same as the above. The final elongation takes 5 min at 72°C. The PCR was performed in the Applied Biosystems 9700 Thermal Cycler.

The PCR product was digested with RsaI (Thermo Fisher Scientific) for 16 hours at 37°C producing two fragments of 50 bp or one fragment of 100 bp for the presence of the G or C allele, respectively. RFLP products were analyzed on 3% agarose gel and stained with SimplySafe (EURx, Poland) and visualized in G:Box (Syngene, Great Britain) (Figure 1). An independent PCR analysis was carried out for each sample.

2.4. NOD2/CARD15 Genotyping

The NOD2/CARD15 frameshift change (3020insC) was analyzed by the use of automated DNA sequencing. The following primers were used to amplify 533 bp: (i)Forward 5-CTG AGC CTT TGT TGA TGA GC -3(ii)Reverse 5-TCT TCA ACC ACA TCC CCA TT-3

Each PCR mixture (25 μl) contained 50 ng genomic DNA, PCR buffer (Clontech), dNTP mixture (0.25 mM), HD polymerase (0.31 U) (Clontech), and primers (10 μM of each). The mixture was heated at 94°C for 5 min and underwent 35 cycles of amplification: denaturation at 98°C for 15 s, annealing at 60°C for 10 s, and elongation at 72°C for 20 s. The final elongation takes 5 min at 72°C. The PCR was performed in the Applied Biosystems 9700 Thermal Cycler using the HD Advantage polymerase (Clontech). Sequencing PCR was performed with the use of the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) in a thermal cycler (as described previously). The sequencing PCR product was purified by the use of an exterminator kit (A&A Biotechnology). The sequencing run module was StdSeq50_POP7 in genetic analyzer 3130 (Applied Biosystems). The results were analyzed by the use of Applied Biosystems software (Figure 2).

2.5. Cytogenetic Analyses

Abnormalities essential for MM, such as del(17p13.1) and IgHV gene rearrangements—t(4;14) and t(14;16)—were tested by simultaneous staining of cytoplasmic immunoglobulin and FISH (cIg-FISH) according to Ross et al.’s recommendations and were identified according to the previously described protocol with modifications [2024].

2.6. Bortezomib In Vitro Treatment

Bone marrow aspirates () (mean number of plasma cells—) were used to establish cell cultures as described by Zmorzynski et al. [20]. The number of apoptotic, necrotic, and viable cells was determined by means of the Annexin V-Cy3 Apoptosis Detection Kit according to the manufacturer’s protocol (Sigma-Aldrich, USA) and following the methods of Zmorzynski et al. [24].

2.7. Statistical Analysis

Laboratory values of MM patients with polymorphisms were compared using an independent -test for continuous variables and the chi-square test for categorical variables. The association of the studied polymorphisms with prognostic factors was evaluated using the chi-square test or Fisher’s exact test (when one expected value was <5). The quantitative data was shown as frequency or percentage. Deviation of genotype frequencies in controls and cases from the Hardy-Weinberg equilibrium (HWE) was assessed by the chi-squared test with Yates’s correction for the groups with less than five patients [25]. For 95% confidence interval (CI), we assumed and ; therefore, if and the corresponding , then the population is in HWE. The logistic regression was used to evaluate the fold risk of MM. The Cox proportional hazard model was used for univariate and multivariate analyses of OS and PFS. Variables of the international staging system (ISS), auto-HSCT, and genotypes (without combinations) were included in multivariate analysis. The Kaplan-Meier method and the log-rank test were used for survival analysis. The statistical power of the study was calculated according to Bacchetti and Leung [26]. We assumed a 5% error of inference and the related level of significance pointing to the existence of statistically significant differences. Statistical analyses were performed using the Statistica ver. 12.5 (StatSoft) software.

3. Results

The presented study included 100 MM patients, 53 males and 47 females, with a median age of 65.32 years. Genotyping was successful in all the individuals investigated within the study. The HWE test confirmed that the genotypic frequencies of PSMA6 and NOD2/CARD15 genes for healthy individuals (controls) and MM patients were balanced (Table 3). The allelic frequencies of the PSMA6 or NOD2/CARD15 gene between study and control groups were statistically insignificant (Table 4). The studied PSMA6 and NOD2/CARD15 variants did not impact on the risk of MM (Table 5). The CG and GG genotypes of the PSMA6 gene were analyzed together, because the number of GG cases was very low (Table 3). Moreover, we analyzed the effect of both changes—3020insC of the NOD2/CARD15 gene and CG+GG genotypes of the PSMA6 polymorphism—on the risk of MM, which was statistically insignificant (Table 5). In the case of the 3020insC variant of the NOD2/CARD15 gene, only heterozygotes were observed.


GroupsGenotypesTotalHWE value and

PSMA6 gene
CCCGGG
Control
 E84.6414.720.64100,
 O84160100
Case
 E81181100,
 O81181100
NOD2/CARD15 gene
Wild-type homozygotesHeterozygotes3020insC homozygotes
Control
 E96.043.920.04100,
 O9640100
Case
 E92.167.680.16100,
 O9280100

If and the corresponding , then the population is in HWE.

AllelesMM, of alleles (%)Controls, of alleles (%) values

PMSA6
 C180 (90)184 (92)0.58
 G20 (10)16 (8)
 Total200 (100)200 (100)
NOD2/CARD15
 Wild type192 (96)196 (98)0.24
 Mutated8 (4)4 (2)
 Total200 (100)200 (100)


GenotypesMM patientsControlsOR95% CI value

NOD2/CARD15
 Wild type9296Reference
 Mutated840.480.14-1.640.23
PMSA6
 CC8184Reference
 CG+GG19160.810.39-1.680.57
NOD2/CARD15PSMA6
 Wild typeCC6978Reference
 MutationCG+GG410.220.02-2.020.31

A univariate Cox analysis revealed that patients at stage III according to ISS had a 1.7-fold increased risk of death () (Table 6). Similar findings were observed in the case of disease relapse or progression in MM patients at stage III (, ) and without auto-HSCT (, ) (Table 6). Moreover, the univariate Cox analysis did not show an increased risk of death or disease relapse/progression in MM patients with the PSMA6 polymorphism and/or NOD2/CARD15 variant. The multivariate Cox regression analysis confirmed that patients without auto-HSCT or at stage III according to ISS had an increased risk of death or disease relapse/progression, respectively (Table 7).


VariableUnivariate Cox analysis for OSUnivariate Cox analysis for PFS
valueHR95% CI valueHR95% CI

ISS
 I+IIRR
 III0.0031.691.18-2.41<0.0011.681.27-2.22
Auto-HSCT
 YesRR
 No<0.0015.302.11-13.340.0022.541.37-4.70
PSMA6
 CCRR
 CG+GG0.130.430.14-1.310.331.150.15-3.96
NOD2/CARD15
 Wild typeRR
 Mutated0.490.730.29-1.780.930.970.45-2.06
PSMA6 and NOD2/CARD15
 CC+wild typeRR
 CG+GG+mutated0.601.680.22-12.580.253.140.43-22.97

R: reference.

VariableMultivariate Cox analysis for OSMultivariate Cox analysis for PFS
valueHR95% CI valueHR95% CI

ISS
 I+IIReferenceReference
 III0.421.180.77-1.820.011.571.11-2.21
Auto-HSCT
 YesReferenceReference
 No0.013.821.36-10.780.141.720.82-3.61
PSMA6
 CCReferenceReference
 CG+GG0.861.420.02-82.330.415.180.095-281.82
NOD2/CARD15
 Wild typeReferenceReference
 Mutated

Too small group for analysis.

The analysis of the response rate in MM showed that patients at stage III or without auto-HSCT had an increased chance of progressive disease (PD) (Table 8). Moreover, patients with CG+GG genotypes of PSMA6 had higher chances for PD (Table 8). The statistical power of this association was 0.63. The CG+GG genotypes had an increased risk of death due to MM at the level of tendency (, 95% CI 1.02-21.97, ).


VariableResponse rate
CR+VGPR+PR+SDPD
valueOR (95% CI)

ISS
 I+IIReference
 III<0.0015.44 (2.14-13.84)
Auto-HSCT
 YesReference
 No0.0153.36 (1.22-9.25)
PSMA6
 CCReference
 CG+GG0.055.0 (1.07-23.16)
NOD2/CARD15
 Wild typeReference
 Mutated0.090.4 (0.13-1.18)
PSMA6 and NOD2/CARD15
 CC+wild typeReference
 CG+GG+mutated0.42

Too small group for analysis. CR: complete response; VGPR: very good partial response; PR: partial response; SD: stable disease; PD: progressive disease patients—according to the response criteria for MM [14, 15].

Furthermore, we analyzed the potential relationships between clinical and laboratory results and selected genotypes. We found that the NOD2/CARD15 variant decreased the risk of renal dysfunction in MM patients (, 95% CI 0.07-0.74, ). We did not confirm any significant correlations among the studied variants of NOD2/CARD15 and PSMA6 genes regarding clinical parameters (Table 9). The odds ratio was used to analyze the relationship of the PSMA6 polymorphism and NOD2/CARD15 variant with chromosomal aberrations in MM patients. The results were statistically insignificant for PSMA6 CG+GG (, 95% CI 0.36-4.13, ), as well as for the NOD2/CARD15 variant (, 95% CI 0.55-12.52, ). Taking into account the type of chromosomal aberration and studied genotypes, we did not find statistically significant results - PSMA6 CC versus CG+GG - del(17p13.1) (), t(4;14) (), t(14;16) () and NOD2/CARD15 wild type versus 3020insC variant - del(17p13.1) (), t(4;14) (), t(14;16) ().


VariablesMM patientsNOD2/CARD15PSMA6
Wild typeMutated valueCCCG+GG value

Mean age (years)65.3265.2065.930.7865.3965.00.87
Free light chain ratio303241.41617.370.08324.44217.140.59
% of plasma cells in bone marrow30.7531.1828.400.6230.0533.880.47
Albumins (g/dl)3.573.583.560.913.533.750.20
β2-Microglobulin (mg/l)5.985.796.930.306.264.810.15
Calcium (mM/l)2.452.432.530.272.452.430.81
Hemoglobin (g/dl)10.3810.3910.330.9110.3510.510.74
Creatinine (mg/dl)1.571.491.980.281.691.070.16
C-reactive protein (mg/l)15.6814.023.680.3116.5712.260.63
Estimated glomerular filtration rate (ml/min/1.73 m2)70.057.990.1566.1876.190.20
PFS (months)18.3218.5916.870.7318.2518.570.94
OS (months)3126.9421.310.4326.2425.150.87

At diagnosis.

We analyzed (by ANOVA) the association between the studied genotypes and survival of MM patients. Without taking into account the treatment, we did not find statistically shorter OS and PFS in patients with the wild-type NOD2/CARD15 or PSMA6 CC genotype in comparison to those with the 3020insC variant of NOD2/CARD15 or PSMA6 CG+GG genotypes (Table 9). In the log-rank test (Figures 3 and 4) without taking into account the type of treatment, the tendency for shorter PFS in patients with CG+GG genotypes was observed (Figure 4). Furthermore, a log-rank analysis taking into account the studied changes and the types of treatment (thalidomide versus bortezomib versus both—thalidomide and bortezomib) was performed. The NOD2/CARD15 3020insC variant did not impact the type of treatment versus OS () or PFS (). Similar results were observed in the case of the CC genotype of the PSMA6 gene (OS () and PFS ()). Individuals with CG+GG genotypes treated with bortezomib (VCD) or both—bortezomib and thalidomide (VTD)—had shorter OS in comparison to those treated with thalidomide (CTD) (). Shorter PFS (at the level of tendency) was observed in patients with CG+GG genotypes treated with bortezomib (VCD) or both (bortezomib and thalidomide (VTD)) in comparison to those treated with thalidomide (CTD) ().

In in vitro studies, bortezomib increased the number of apoptotic and necrotic cells in all the studied genotypes. The only statistically significant differences were observed in the number of apoptotic cells at concentrations of bortezomib at 8 nM and 12 nM between wild-type and 3020insC variants of NOD2/CARD15 (Table 10).


GenotypesBortezomib/DMSO
Control 0.1% DMSO (0 nM) value1 nM value2 nM value4 nM value8 nM value12 nM value

Apoptotic cells (%)
 Mean5.4515.4417.4925.5537.5247.68
NOD2/CARD15 wild type5.100.3015.510.9217.290.6626.180.3739.640.01849.480.03
 NOD2/CARD15 mutated6.9615.1018.3422.7828.2839.86
PSMA6 CC5.460.9714.790.4917.650.7624.660.3035.580.0747.150.61
PSMA6 CG+GG5.4117.3917.028.1943.3349.25
 CC+wild type5.190.5514.460.7817.520.8425.090.8137.070.3448.860.45
 CG+GG+mutated6.7512.9718.2523.7530.7743.97
Necrotic cells (%)
 Mean2.434.356.6511.6912.8414.93
NOD2/CARD15 wild type2.600.474.380.876.610.9112.090.3613.150.5215.370.31
 NOD2/CARD15 mutated1.684.216.789.9711.4313.02
PSMA6 CC2.380.854.640.236.660.9711.270.4212.760.9114.990.90
PSMA6 CG+GG2.583.486.6112.9413.0414.75
 CC+wild type2.580.934.630.566.670.8811.920.6613.090.8515.460.71
 CG+GG+mutated2.423.637.0513.3712.3014.22
Viable cells (%)
 Mean89.1779.5876.0866.1859.7153.32
NOD2/CARD15 wild type89.560.4278.910.4076.160.9165.620.4958.600.1452.830.42
 NOD2/CARD15 mutated87.4982.4775.6968.6064.5155.43
PSMA6 CC90.060.1280.650.2675.390.4866.880.4760.690.2853.930.40
PSMA6 CG+GG86.4876.3578.1364.0756.7651.47
 CC+wild type90.400.2979.990.9575.340.9566.370.9559.710.7353.050.61
 CG+GG+mutated86.8579.6175.7466.8061.6750.45

4. Discussion

To our knowledge, this is the first study elucidating the correlation of the PSMA6 polymorphism and NOD2/CARD15 variant with the risk and the outcome of MM, as well as the response to bortezomib under in vitro conditions.

Little is known about the role of the PSMA6 -8C>G polymorphism in cancer development, including MM. Most of the research was focused on the association of this SNP with the risk of cardiovascular disease and diabetes mellitus [27, 28]. Our findings suggest that CG+GG genotypes of the PSMA6 gene might be associated with the response rate of MM patients and shorter OS taking into account the type of treatment. The G allele of this SNP leads to higher promoter activity and elevated PSMA6 gene expression in human B cells in comparison to the C allele [12]. Higher proteasome activity, as a result of the GG genotype, may lead to degradation of misfolded or unfit antiapoptotic proteins [29]. Bachmann and coworkers in 116 MM patients showed shorter median survival time for individuals with the GG genotype [21]. Moreover, they revealed in Cox multivariate analysis that the G allele was an independent prognostic factor for MM patients [21]. In our Cox multivariate analysis, we did not confirm the association of the G allele with the MM outcome. However, in the log-rank test, we observed higher chances of progression in patients with CG+GG genotypes. Additionally, taking into account the type of treatment, we found shorter OS in individuals with CG+GG genotypes treated with bortezomib/bortezomib+thalidomide in comparison to those (with CG+GG genotypes too) treated with thalidomide. It has a different mechanism of action in comparison to bortezomib. Thalidomide, as an immunomodulatory drug, inhibits the NF-κB pathway and provides better survival of MM patients [30]. The inflammation development through activation of the NF-κB pathway may lead to malignant transformation, as well as disease progression [13].

The α6 subunit, encoded by the PSMA6 gene, is present in the core of the 20S proteasome. It belongs to proteasome-interacting proteins and among others with the β5 subunit plays a role as an auxiliary factor. Suppression of the proteasome leads to endoplasmic reticulum stress and apoptosis [31]. Inhibition of the proteasome using β5 subunit blockers like bortezomib results in cell cycle arrest and induction of apoptosis in MM cells [32]. However, in our study, we did not observe the relationship between the PSMA6 gene polymorphism and response to bortezomib treatment under in vitro conditions. It is possible that the α6 subunit in comparison to β5 of the 20S proteasome does not modify substantially the proteasome function.

The change in the NOD2/CARD15 gene in the form of C insertion in a region encoding the ligand recognition domain (exon 11) is described as polymorphism (rs2066847) or mutation [5, 9, 10, 33, 34]. The difference is due to the fact that in cancer research, it is called mutation and in a noncancer field, the term polymorphism is used. The current guidelines of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology recommend the use of the neutral term “variant” only [35]. The 3020insC variant of the NOD2/CARD15 gene is implicated in the development of chronic inflammation and increased risk of solid tumors such as colorectal, gastric, breast, and lung cancer, as well as non-Hodgkin’s lymphoma [10, 3639]. Kurzawski et al. found a statistically significant difference in mutated NOD2/CARD15 allele distribution in 250 nonhereditary nonpolyposis colorectal cancer patients in comparison to healthy individuals [36]. In our study, we did not observe differences in NOD2/CARD15 wild-type/3020insC allele distribution between MM patients and the control group. The frequency of these alleles in our study populations was in HWE. Moreover, we did not find a statistically significant impact of the NOD2/CARD15 variant on OS and PFS of MM patients. This result is consistent with the observation of Huszno and coworkers. In their study in breast cancer patients, they found no difference in OS between 3020insC variant carriers and wild-type homozygotes of the NOD2/CARD15 gene [9].

Furthermore, we found that the NOD2/CARD15 variant decreased the risk of renal dysfunction in MM patients, which was not previously reported by other researchers. In an animal model, Stroo and coworkers observed no difference in renal injury between wild-type Nod2 double-knockout mice with chronic kidney disease [40]. They concluded that Nod2 does not play an important role in the development of renal damage [40]. On the other hand, Kruger et al. observed an unfavorable effect of the NOD2/CARD15 wild type on the outcome of graft survival in 352 patients receiving their first renal transplant [41]. Moreover, they took into account two additional SNPs (SNP8—R702W and SNP12—G908R) present in the NOD2/CARD15 gene, which are associated with increased susceptibility to Crohn’s disease [41]. It is possible that the NOD2/CARD15 variant showed only a minor effect and cooccurrence of these SNPs may affect or modify the function of NOD2 protein. In our further research, it would be worth extending the analysis by these two SNPs.

In our in vitro study, the 3020insC variant of the NOD2/CARD15 gene was associated with statistically significant lower levels of apoptotic cells at 8 nM and 12 nM of bortezomib. The loss of function of the NOD2/CARD15 gene in MM leads to apoptosis inhibition. It is known that overexpression of NOD2/CARD15 enhances apoptosis via caspase-9 [42]. It would be desirable to examine the NOD2/CARD15 gene expression in plasma cells after in vitro bortezomib treatment.

There are some limitations of our study associated with the number of MM patients recruited to the study, as well as the used methods.

The sample size was relatively small, in part due to the low incidence of the disease. We have collected detailed clinical data for 100 MM patients. Some patients were not included into this analysis, because they are still in follow-up and their clinical data are incomplete. However, in our research, the number of 100 MM patients is enough for most analyses. Some of them were not possible as a result of low frequency of genotypes. Further analysis on a larger cohort can help better understand the significance of the studied changes in the pathobiology of MM, especially in the case of disease outcome.

In our in vitro study, we used RPMI media with 10% FCS. Currently, AIM-V media are used instead of RMPI and FCS as they induce apoptosis of primary B cells [43]. In our control samples without bortezomib, spontaneous apoptosis was observed, and therefore, in our further research, apoptosis evaluation will be done using serum-free AIM-V medium. Moreover, for the evaluation of apoptosis, fluorescent microscopy instead of flow cytometry-based apoptosis detection (FACS) was used. FACS analysis is more reliable and quantitative for apoptosis evaluation. Unfortunately, during the experiment time, the flow cytometry method was not available for us and retrospective analysis of apoptosis is not possible. The set used for apoptosis and necrosis detection was dedicated and validated to fluorescent microscopy.

In this paper, we analyzed, for the first time, the significance of the PSMA6 polymorphism and NOD2/CARD15 variant in the survival and response to therapy of MM patients. In conclusion, these data underscore the possible relationship between PSMA6 CG+GG genotypes and response rate of MM patients, as well as shorter OS taking into account the type of treatment. Moreover, we observed an association of the NOD2/CARD15 variant with the lower risk of renal dysfunction in the study group and lower apoptosis levels in cell culture with bortezomib.

Data Availability

The clinical data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors are grateful to the patients who participated in this study. This project was supported by the Medical University of Lublin, Poland (Grants MNmb 232 and DS 225).

References

  1. J. J. Keats, M. Chesi, J. B. Egan et al., “Clonal competition with alternating dominance in multiple myeloma,” Blood, vol. 120, no. 5, pp. 1067–1076, 2012. View at: Publisher Site | Google Scholar
  2. Y. Fan, R. Mao, and J. Yang, “NF-κB and STAT3 signaling pathways collaboratively link inflammation to cancer,” Protein & Cell, vol. 4, no. 3, pp. 176–185, 2013. View at: Publisher Site | Google Scholar
  3. S. K. Kumar, V. Rajkumar, R. A. Kyle et al., “Multiple myeloma,” Nature Reviews Disease Primers, vol. 3, no. 1, article 17046, 2017. View at: Publisher Site | Google Scholar
  4. D. Nagel, M. Vincendeau, A. C. Eitelhuber, and D. Krappmann, “Mechanisms and consequences of constitutive NF-κB activation in B-cell lymphoid malignancies,” Oncogene, vol. 33, no. 50, pp. 5655–5665, 2014. View at: Publisher Site | Google Scholar
  5. Y. Ogura, D. K. Bonen, N. Inohara et al., “A frameshift mutation in NOD2 associated with susceptibility to Crohn's disease,” Nature, vol. 411, no. 6837, pp. 603–606, 2001. View at: Publisher Site | Google Scholar
  6. O. Gutierrez, C. Pipaon, N. Inohara et al., “Induction of Nod2 in myelomonocytic and intestinal epithelial cells via nuclear factor-κB activation,” The Journal of Biological Chemistry, vol. 277, no. 44, pp. 41701–41705, 2002. View at: Publisher Site | Google Scholar
  7. P. Magnowski, K. Medrek, M. Magnowska et al., “The 3020insC NOD2 gene mutation in patients with ovarian cancer,” Ginekologia Polska, vol. 79, no. 8, pp. 544–549, 2008. View at: Google Scholar
  8. T. Watanabe, N. Asano, P. J. Murray et al., “Muramyl dipeptide activation of nucleotide-binding oligomerization domain 2 protects mice from experimental colitis,” The Journal of Clinical Investigation, vol. 118, pp. 545–559, 2008. View at: Publisher Site | Google Scholar
  9. J. Huszno, Z. Kołosza, K. Tęcza, J. Pamuła-Piłat, M. Mazur, and E. Grzybowska, “Comparison between NOD2 gene mutation carriers (3020insC) and non-carriers in breast cancer patients: a clinicopathological and survival analysis,” Archives of Medical Science - Civilization Diseases, vol. 3, no. 1, pp. 10–e15, 2018. View at: Publisher Site | Google Scholar
  10. I. Omrane, A. Mezlini, O. Baroudi et al., “3020insC NOD2/CARD15 polymorphism associated with treatment of colorectal cancer,” Medical Oncology, vol. 31, no. 5, p. 954, 2014. View at: Publisher Site | Google Scholar
  11. S. Dimmeler, K. Breitschopf, J. Haendeler, and A. M. Zeiher, “Dephosphorylation targets Bcl-2 for ubiquitin-dependent degradation: a link between the apoptosome and the proteasome pathway,” The Journal of Experimental Medicine, vol. 189, no. 11, pp. 1815–1822, 1999. View at: Publisher Site | Google Scholar
  12. K. Ozaki, H. Sato, A. Iida et al., “A functional SNP in _PSMA6_ confers risk of myocardial infarction in the Japanese population,” Nature Genetics, vol. 38, no. 8, pp. 921–925, 2006. View at: Publisher Site | Google Scholar
  13. Z. J. Chen, “Ubiquitin signalling in the NF-κB pathway,” Nature Cell Biology, vol. 7, no. 8, pp. 758–765, 2005. View at: Publisher Site | Google Scholar
  14. Y. Matsuki-Muramoto, K. Nozawa, K. Uomori, I. Sekigawa, and Y. Takasaki, “Bortezomib treatment prevents glomerulosclerosis associated with lupus nephritis in a murine model through suppressive effects on the immune and renin-angiotensin systems,” Modern Rheumatology, vol. 27, pp. 77–86, 2016. View at: Publisher Site | Google Scholar
  15. S. Frankland-Searby and S. R. Bhaumik, “The 26S proteasome complex: an attractive target for cancer therapy,” Biochimica et Biophysica Acta (BBA) - Reviews on Cancer, vol. 1825, pp. 64–76, 2012. View at: Publisher Site | Google Scholar
  16. S. Bringhen, M. Offidani, S. Palmieri et al., “Early mortality in myeloma patients treated with first-generation novel agents thalidomide, lenalidomide, bortezomib at diagnosis: a pooled analysis,” Critical Reviews in Oncology/Hematology, vol. 130, pp. 27–35, 2018. View at: Publisher Site | Google Scholar
  17. H. Yang, J. A. Zonder, and Q. P. Dou, “Clinical development of novel proteasome inhibitors for cancer treatment,” Expert Opinion on Investigational Drugs, vol. 18, no. 7, pp. 957–971, 2009. View at: Publisher Site | Google Scholar
  18. A. Palumbo, S. V. Rajkumar, J. F. San Miguel et al., “International Myeloma Working Group consensus statement for the management, treatment, and supportive care of patients with myeloma not eligible for standard autologous stem-cell transplantation,” Journal of Clinical Oncology, vol. 32, no. 6, pp. 587–600, 2014. View at: Publisher Site | Google Scholar
  19. B. G. Durie, J. L. Harousseau, J. S. Miguel et al., “International uniform response criteria for multiple myeloma,” Leukemia, vol. 20, no. 9, pp. 1467–1473, 2006. View at: Publisher Site | Google Scholar
  20. F. Ross, H. Avet-Loiseau, G. Ameye et al., “Report from the European myeloma network on interphase FISH in multiple myeloma and related disorders,” Haematologica, vol. 97, no. 8, pp. 1272–1277, 2012. View at: Publisher Site | Google Scholar
  21. H. S. Bachmann, J. Novotny, S. Sixt et al., “The G-allele of the PSMA6-8>G polymorphism is associated with poor outcome in multiple myeloma independently of circulating proteasome serum levels,” European Journal of Haematology, vol. 85, no. 2, pp. 108–113, 2010. View at: Publisher Site | Google Scholar
  22. G. J. Ahmann, S. M. Jalal, A. L. Juneau et al., “A novel three-color, clone-specific fluorescence in situ hybridization procedure for monoclonal gammopathies,” Cancer Genetics and Cytogenetics, vol. 101, no. 1, pp. 7–11, 1998. View at: Publisher Site | Google Scholar
  23. A. Dmoszynska and S. Chocholska, “Molecular biology methods in the diagnosis of multiple myeloma,” in Molecular Aspects of Hematologic Malignancies: Diagnostic Tools and Clinical Applications, M. Witt, M. Dawidowska, and T. Szczepanski, Eds., pp. 443–449, Springer-Verlag Berlin Heidelberg, 2012. View at: Publisher Site | Google Scholar
  24. S. Zmorzynski, A. Szudy-Szczyrek, S. Popek-Marciniec et al., “ACE insertion/deletion polymorphism (rs4646994) is associated with the increased risk of multiple myeloma,” Frontiers in Oncology, vol. 9, p. 44, 2019. View at: Publisher Site | Google Scholar
  25. K. J. Preacher, Calculation for the chi-square test: An interactive calculation tool for chi-square tests of goodness of fit and independence [Computer software], 2001.
  26. P. Bacchetti and J. M. Leung, “Sample size calculations in clinical research,” Anesthesiology, vol. 97, no. 4, pp. 1028-1029, 2002. View at: Publisher Site | Google Scholar
  27. M. Barbieri, R. Marfella, M. R. Rizzo et al., “The −8 UTR C/G polymorphism of PSMA6 gene is associated with susceptibility to myocardial infarction in type 2 diabetic patients,” Atherosclerosis, vol. 201, no. 1, pp. 117–123, 2008. View at: Publisher Site | Google Scholar
  28. K. Hinohara, T. Nakajima, T. Sasaoka et al., “Replication studies for the association of _PSMA6_ polymorphism with coronary artery disease in East Asian populations,” Journal of Human Genetics, vol. 54, no. 4, pp. 248–251, 2009. View at: Publisher Site | Google Scholar
  29. M. Schmidt and D. Finley, “Regulation of proteasome activity in health and disease,” Biochimica et Biophysica Acta (BBA) - Molecular Cell Research, vol. 1843, pp. 13–25, 2014. View at: Publisher Site | Google Scholar
  30. T. Hideshima, D. Chauhan, P. Richardson et al., “NF-κB as a therapeutic target in multiple myeloma,” The Journal of Biological Chemistry, vol. 277, no. 19, pp. 16639–16647, 2002. View at: Publisher Site | Google Scholar
  31. R. Bravo, V. Parra, D. Gatica et al., “Endoplasmic reticulum and the unfolded protein response: dynamics and metabolic integration,” International Review of Cell and Molecular Biology, vol. 301, pp. 215–290, 2013. View at: Publisher Site | Google Scholar
  32. M. Ri, “Endoplasmic-reticulum stress pathway-associated mechanisms of action of proteasome inhibitors in multiple myeloma,” International Journal of Hematology, vol. 104, no. 3, pp. 273–280, 2016. View at: Publisher Site | Google Scholar
  33. N. P. Mayor, B. E. Shaw, J. A. Madrigal, and S. G. E. Marsh, “NOD2 polymorphisms and their impact on haematopoietic stem cell transplant outcome,” Bone Marrow Research, vol. 2012, Article ID 180391, 13 pages, 2012. View at: Publisher Site | Google Scholar
  34. J. Lubinski, T. Huzarski, G. Kurzawski et al., “The 3020insC allele of NOD2 predisposes to cancers of multiple organs,” Hereditary Cancer in Clinical Practice, vol. 3, no. 2, pp. 59–63, 2005. View at: Publisher Site | Google Scholar
  35. S. Richards, N. Aziz, S. Bale et al., “Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology,” Genetics in Medicine, vol. 17, no. 5, pp. 405–423, 2015. View at: Publisher Site | Google Scholar
  36. G. Kurzawski, J. Suchy, J. Kładny et al., “The NOD2 3020insC mutation and the risk of colorectal cancer,” Cancer Research, vol. 64, no. 5, pp. 1604–1606, 2004. View at: Publisher Site | Google Scholar
  37. S. Angeletti, S. Galluzzo, D. Santini et al., “NOD2/CARD15 polymorphisms impair innate immunity and increase susceptibility to gastric cancer in an Italian population,” Human Immunology, vol. 70, no. 9, pp. 729–732, 2009. View at: Publisher Site | Google Scholar
  38. M. R. Lener, D. Oszutowska, J. Castaneda et al., “Prevalence of the NOD2 3020insC mutation in aggregations of breast and lung cancer,” Breast Cancer Research and Treatment, vol. 95, no. 2, pp. 141–145, 2006. View at: Publisher Site | Google Scholar
  39. M. S. Forrest, C. F. Skibola, T. J. Lightfoot et al., “Polymorphisms in innate immunity genes and risk of non-Hodgkin lymphoma,” British Journal of Haematology, vol. 134, no. 2, pp. 180–183, 2006. View at: Publisher Site | Google Scholar
  40. I. Stroo, D. Emal, L. M. Butter et al., “No difference in renal injury and fibrosis between wild-type and NOD1/NOD2 double knockout mice with chronic kidney disease induced by ureteral obstruction,” BMC Nephrology, vol. 19, no. 1, p. 78, 2018. View at: Publisher Site | Google Scholar
  41. B. Kruger, C. A. Boger, B. Schroppel et al., “Impact of NOD2/CARD15 haplotypes on the outcome after kidney transplantation,” Transplant International, vol. 20, no. 7, pp. 600–607, 2007. View at: Publisher Site | Google Scholar
  42. M. Gazouli, “NOD2/CARD15 mediates the inflammatory responses in inflammatory bowel diseases (IBDs),” Annals of Gastroenterology, vol. 18, pp. 16–18, 2005. View at: Google Scholar
  43. S. Sinha, J. Boysen, M. Nelson et al., “Targeted Axl inhibition primes chronic lymphocytic leukemia B cells to apoptosis and shows synergistic/additive effects in combination with BTK inhibitors,” Clinical Cancer Research, vol. 21, no. 9, pp. 2115–2126, 2015. View at: Publisher Site | Google Scholar

Copyright © 2020 S. Zmorzyński et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views129
Downloads119
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.