1. Introduction
Multiple Myeloma (MM) is a plasma cell malignancy characterized by monoclonal plasma cell proliferation in the bone marrow, leading to significant morbidity owing to bone lesions, anemia, renal failure, and hypercalcemia. Recent advances in MM treatment, characterized by the introduction of novel therapeutics, have markedly improved patient survival rates. These developments, along with innovations in cytogenetics and risk-stratification models, have progressively refined the prognostic assessment of MM. Despite these advances, prognosis remains suboptimal for a subset of patients not adequately classified by current risk models, highlighting the need for novel prognostic indicators [
1,
2].
MicroRNAs (miRNAs) are a class of non-coding RNAs that regulate gene expression, predominantly post-transcriptionally influencing various cellular processes including proliferation, apoptosis, and cancer pathogenesis. Post- transcriptionally, miRNAs bind to untranslated regions of mRNA and regulate specific gene expression and act as proto-oncogenesor suppressor miRNAS [
3,
4]. Epigenetic processes such as DNA methylation, histone modification, and RNA modification influence the expression levels of miRNAs, while as a feedback function, miRNAs influence epigenetic enzymes such as DNA methyltransferases , histone deacetylases and histone methyltransferases and operate as epigenetic regulators [
5]. Although, several studies indicate that miRNA expression patterns are related to the biology and clinical behavior of solid and hematologic malignancies, the origin of miRNAs, the mechanisms of transit (secreted as free miRNAs actively or during cell apoptosis) and the most prominent role (intracellularly or extracellularly), are not well known [
6]. Due to their stability as molecules, extracellular expression levels have been investigated for their potential prognostic impact in many tumors. However, limited research is available regarding the prognostic significance of circulating miRNAs in MM; miRNAs can be found in the serum in two forms, free miRNA and exosomal miRNA.
MiR-16 was initially associated with hematological neoplasms, such as chronic lymphocytic leukemia [
7] and was among the first miRNAs examined in MM as part of Al Masri et al.'s study on miRNA expression analysis in multiple myeloma [
8,
9]. Few studies in MM bone marrow have examined miR-16, showing downregulation and suggesting its potential as a prognostic marker [
10].
Recent studies have shown that miR-21 is a proto-oncogene that plays a significant role in bone marrow-mediated growth, survival, and resistance of MM cells to therapy. Although, miR-21 has been studied in many cancers [
11], data of serum levels in MM are extremely limited.
The role of miRNAs in MM represents an appealing approach for advancing our understanding of the molecular underpinnings of the disease and identifying novel prognostic tools.
Over more than a decade, the therapeutic landscape of MM has been revolutionized by the inclusion of lenalidomide in treatment regimens, initially in relapsed/refractory MM (RRMM), and then rapidly in both the first-line and maintenance therapy strategies, which were significantly improved. However, resistance to lenalidomide emerged as a challenge as these patients often have poor outcomes [
12].
This backdrop sets the stage for our investigation of the prognostic relevance of miR-16 and miR-21 in newly diagnosed and relapsed MM, and their potential as biomarkers for lenalidomide treatment response. This study aimed to enlighten the prognostic landscape of MM and explore novel avenues for personalized therapy, potentially guiding more effective treatment strategies and improving patient outcomes.
2. Patients and Methods
2.1. Patients and Samples
In this study, we retrospectively enrolled a total of 48 MM patients at diagnosis and 35 at relapse prior to Lenalidomide-Dexamethasone (LD) treatment. All patients were treated in our department with LD, 13 at first line and 35 at relapse, of whom 21 were tested twice at diagnosis and before LD initiation. Systematic collection of serum samples was performed from both, the patient and control group of healthy individuals, ensuring that the samples were stored in a frozen state for subsequent analysis. Serum samples were collected at the time of diagnosis and before the first administration of lenalidomide and dexamethasone. Our study’s inclusion criteria were based on the International Myeloma Working Group (IMWG) criteria for MM diagnosis, and all patients received a therapeutic protocol consisting of Lenalidomide with dexamethasone at either diagnosis or relapse. None of our patients received Lenalidomide as maintenance therapy. The therapeutic protocol comprised lenalidomide (d1-d21) in combination with dexamethasone 20 or 40 md (d1, d8, d15, and d22) in a 28-day cycle (LD). The Glomerular Filtration Rate (GFR) was determined using the Cockcroft-Gault equation.
2.2. MicroRNA Quantification
2.2.1. RNA Isolation
Total RNA enriched for miRNAs was isolated from 100 μl of blood serum using the MagMAX™ mirVana™ Total RNA Isolation Kit (Thermo Fisher Scientific Inc.) according to the manufacturer’s protocol. Magnetic-bead technology was used to uniformly recover high-quality RNA. The Implen P330 nanophotometer (Implen GmbH) was used to evaluate the concentration and purity of the extracted total RNA
2.2.2. cDNA Synthesis
MiRNA quantification was performed for miR-16-5p and miR21-5p. A TaqMan® Advanced miRNA cDNA Synthesis Kit (Thermo Fisher Scientific Inc.) was used to generate cDNA, which uses universal primers to consistently amplify all targets, including low-expression miRNA targets, thereby increasing assay sensitivity.
2.2.3. Serum microRNA Expression Quantification by Quantitative Realtime PCR (qRT-PCR)
TaqMan® Advanced miRNA Assays (for hsa-miR-16-5p (assay ID #477860_mir) for hsa-miR-21-5p (assay ID #477975_mir) and for cel-miR-39-3p (assay ID #478293_mir)) were used after cDNA synthesis to perform qRT-PCR with TaqMan® Fast Advanced Master Mix providing high specificity and the StepOnePlus™ Real-Time PCR System (Thermo Fisher Scientific Inc.). Each sample was analyzed in duplicate by qPCR. ExpressionSuite™ Software was used for the analysis and calculation of relative gene expression using the comparative Ct (ΔΔCt) method and normalization of sample-to-sample variation to an exogenous control. For accurate and precise quantification of circulating miRNA levels, the Caenorhabditis elegans miRNA Cel-miR-39-3p was used as an exogenous control. Finally, the relative levels of miRNA in patient samples were compared to the lowest value of healthy individuals’ miRNA (due to exceptionally elevated values of miRNA expression in HI compared to MM patients), which was used as a reference sample, and the final results were presented as fold change in expression using the 2ΔΔCt formula.
2.3. Statistical Analysis
Mann–Whitney test was used to analyze differential miRNAs expression between MM patients and HIs. We used the median of miRNA expression levels as a cut off point in the correlation and survival analysis.
To assess the impact of varying serum miRNA 16 and 21 levels at diagnosis (Dx) and LD initiation, the ratio of these two variables was calculated (miR-LD/DX). For the prognostic utility of this analysis, a ratio value greater than two was selected as the cutoff standard.
Clinical characteristics of the different groups and their associations with miRNAs were evaluated using Spearman correlation analysis. The Kaplan-Meier method was used to calculate the overall survival curves according to miRNA-16 and miRNA-21 expression. Survival differences according to gene expression were analyzed using the log-rank test. All statistical analyses were performed using SPSS software (version 28.0, IBM, USA). Statistical significance was set at P < 0.05.
2.4. Ethical Approval
The study was approved by the institutional review board (IRB) of our hospital (protocol number 191 and date of approval 21 February 2019). Informed consent was obtained from each patient prior to participation in the study, followed by a comprehensive review of their medical records. Patient confidentiality was maintained throughout the study.
4. Discussion
Recent research has focused on miRNAs function and their involvement in epigenetic regulation of tumors. Studies have shown that miR-16 and miR-21 regulate cell cycle and are linked to the pathogenesis of MM.
Downregulation of miR-16 enhances cell proliferation through AKT3 kinase, MAP kinases, ribosomal protein S6, and the NF-kB activator MAP3KIP3 inhibition, and influences angiogenesis and apoptosis by increasing the expression of VEGF, BCL2, CCND1, WNT3A, and MCL1 genes [
10]. In MM, miRNAs-15a and -16 may disable disease progression by inhibiting various kinase pathways [
2]. Moreover, miR-16 regulates the cell cycle and is involved in drug resistance in MM through induction of IL-6 and cell cycle progression [
13]. Recent research showed that serum miR-16 expression levels differed amongst malignancies. MiR-16 had lower differential expression levels in esophageal and stomach cancers compared to normal controls, indicating that miR-16 may operate as a tumor suppressor [
14]. In MM, few studies in the bone marrow, have demonstrated that miR-16 expression was downregulated compared with normal control cells and their expression levels had prognostic impact [
15]. We studied serum levels of miR-16 which were significantly down regulated in newly diagnosed MM patients compared to normal controls. (
Table 4) Similar to our findings, circulating miR-16 was found to be downregulated both as free serum miR-16 [
16,
17], as well as in exosomal serum form [
4] when compared to healthy individuals miR-16 expression.
MiR-21 is involved in several signal transduction pathways, important to MM pathophysiology, such as the JAK/STAT-3 and NF-KB. Although, circulating miR-21 has been shown to operate as a proto-oncogene in many cancers [
9,
18], the studies of serum miR-21 in MM are extremely limited with conflicting results [
16,
19]. In one study, serum expression of mir-21 studied in 30 MM patients was significantly higher than in the MGUS group (14 patients) and the control group [
16]. On the contrary, miR-21 in serum of MM, SMM and MGUS patients was shown to be downregulated only in MM patients, indicating that considerable down-regulation of this miRNA may be associated with subsequent events in disease progression. Furthermore, serum exosomes of miR-21 have been found downregulated in MM patients[
20]. In our study, serum expression miR-21 was significantly downregulated in newly diagnosed MM patients compared to normal controls. (Table 5).
The expression levels of serum miRNA-16 and 21 measured right before LD treatment in the same patients, remained significantly suppressed compared to healthy individuals, but a significant up-regulation compared to the levels at diagnosis was observed. (
Figure 1). Clinical characteristics, important in MM for the diagnosis and prognosis, were significantly correlated to serum miR-16 in our patients. This was the case for renal impairment, hypercalcemia, B2M and CRP serum levels. Serum Levels of B2M, renal impairment, CRP, and neu/lymph ratio were all associated with miR-21 at the time of diagnosis. MiR-21 serum levels at relapse were associated with ISS stage. From previous studies, plasma levels of miRNA 16 were associated with B2M [
21], ISS disease stage and PFS in myeloma [
15]. In addition, serum miR-21 was significantly different according to Durie and Salmon stages, correlates with B2M, IgG and albumin levels [
16]
Increased interest has been observed regarding the value of miRNAs in drug resistance in MM. In this context, serum levels of distinct miRNAs were shown by Robak et al. [
22] to be expressed differently in MM patients based on Bortezomib sensitivity. Additional investigation has shown that miR-16 levels increased after melphalan and bortezomib treatment and decreased following the impact of IL-6 fluctuations, indicating their potential value as a response biomarker to therapies as well [
23]. Likewise, increased expression of miR-21 in bone marrow stromal MM cells, promoted resistance to dexamethasone, doxorubicin, or bortezomib through upregulation of nuclear factor-κB (NF-κB) [
24]. In contrast, the role of miRNAs in the modulation of the activity of IMiDs has been poorly investigated. Recently, argonaute 2 (AGO2), which plays a pivotal role in miRNAs’ maturation and function, was identified as a cereblon binding partner and it was found that the steady-state levels of AGO2 are regulated by cereblon [
25]. This led to the concept that miRNAs had a synergistic impact with anti-MM factors and might be used to enhance anti-MM therapy or used as prognostic implicators. In an attempt to determine the potential association of miRNA expression to Lenalidomide, researchers showed that serum miRNA expression levels, miR-26a-5p, miR-29c-3p, miR-30b-5p, miR-30c-5p, and miR-331-3p were significantly down-regulated in poor responders to treatment [
26].
In addition, we observed that higher levels of miR-16 (x2- twice the expression of median miR-16) correlated with response to treatment greater than PR and higher miR-21 levels with a response better than VgPR. Additionally, an increased ratio between serum at diagnosis and before LD treatment (DX/LD) of both miR-16 and miR-21 showed significant association with late relapse to LD treatment (>24 months). Indeed, this is a very interesting finding taking into consideration that poor survival outcome has been linked to early relapse, irrespectively of the depth of response to LD achieved [
27].
Finally, our findings showed that miR-16 and miR-21 levels have prognostic value for MM patients survival. Higher serum levels of miR-16 were associated with improved OS. Furthermore, a difference in serum levels between diagnosis and relapse MM was indicated by a lower than two miR-16 Dx/Rel ratio linked with shorter TTP and a miR-21 Dx/Rel ratio related with shorter TTNT, when evaluated prior to LD treatment. Lower expression of other miRNAs (miR-30b5p, miR-30c-5p, miR-193a-5p, miR-331-3p) has been related to poorer overall survival or significantly shorter TTP after lenalidomide treatment [
26].
One report by Rocci et al 2014 [
6] found that lower miR-16 levels are significantly related to shorter OS. Interestingly, when they compared, intracellular miRNAs (from MM cells) and serum miRNAs from the same patients, no correlation was observed, thus suggesting a different origin for circulating miRNAs, such as normal immune/inflammatory or stromal cells. This may be the case also for miR-21 difference between low serum levels and intracellular increased expression levels of miR-21 in MM but this remains to be elucidated with further research.
Author Contributions
M- C K and GVZ D; designed the study, supervised and corrected the manuscript, P.P. S; supervised the study, A-I Gkioka; performed the study and the measurements, wrote the paper and followed patients, M T; performed the experiments and wrote the paper, AK; wrote the paper, followed patients, C- AM2, AP, T –MT; performed experiments, AG, AA, MP, VB; followed patients.