Next Article in Journal
Elevated Level of Nerve Growth Factor (NGF) in Serum-Derived Exosomes Predicts Poor Survival in Patients with Breast Cancer Undergoing Neoadjuvant Chemotherapy
Previous Article in Journal
Skin Lesion Classification Based on Surface Fractal Dimensions and Statistical Color Cluster Features Using an Ensemble of Machine Learning Techniques
Previous Article in Special Issue
Sampling, Logistics, and Analytics of Urine for RT-qPCR-based Diagnostics
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 21
10.3390/cancers13215258
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Circulating microRNAs Correlate with Multiple Myeloma and Skeletal Osteolytic Lesions

by
Sara Reis Moura
1,2,3,†,
Hugo Abreu
1,2,†,
Carla Cunha
1,2,
Cláudia Ribeiro-Machado
1,2,
Carla Oliveira
1,2,
Mario Adolfo Barbosa
1,2,3,
Herlander Marques
4,5,6 and
Maria Inês Almeida
1,2,3,*
1
Instituto de Investigação e Inovação em Saúde (i3S), University of Porto, 4200-135 Porto, Portugal
2
Instituto de Engenharia Biomédica (INEB), University of Porto, 4200-135 Porto, Portugal
3
Instituto de Ciências Biomédicas Abel Salazar (ICBAS), University of Porto, 4050-313 Porto, Portugal
4
Hospital de Braga, 4710-243 Braga, Portugal
5
Centro Clínico Académico 2CABraga/ICVS, University of Minho, 4710-243 Braga, Portugal
6
Centro de Investigação em Tecnologias e Serviços de Saúde (CINTESIS), Faculty of Medicine, University of Porto, 4169-007 Porto, Portugal
*
Author to whom correspondence should be addressed.
†
These authors equally contributed to this work.
Cancers 2021, 13(21), 5258; https://doi.org/10.3390/cancers13215258
Submission received: 30 July 2021 / Revised: 1 October 2021 / Accepted: 12 October 2021 / Published: 20 October 2021
(This article belongs to the Special Issue Circulating miRNAs as Tumor Biomarkers)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Molecular biomarkers for the diagnosis of multiple myeloma and for the early detection of the associated osteolytic lesions are needed. MicroRNAs are a class of small non-coding RNAs that regulate gene expression post-transcriptionally and have been explored as circulating (extracellular) biomarkers for distinct diseases. Results show that miR-16-5p, miR-20a-5p, and miR-21-5p levels are differently expressed in the plasma of multiple myeloma patients compared with the control group and suggest that their combined expression could be used as a potential circulating biomarker. Furthermore, the expression of plasma microRNAs significantly correlates with myeloma bone disease and with bone lesions in the spine.

Abstract

Multiple myeloma (MM) is the second most frequent hematological disease and can cause skeletal osteolytic lesions. This study aims to evaluate the expression of circulating microRNAs (miRNAs) in MM patients and to correlate those levels with clinicopathological features, including bone lesions. A panel of miRNAs associated with MM onset and progression, or with bone remodeling, was analyzed in the plasma of 82 subjects (47 MM patients; 35 healthy controls). Results show that miR-16-5p, miR-20a-5p, and miR-21-5p are differently expressed between MM patients and healthy controls. Receiver operating characteristic analyses indicate that their combined expression has potential as a molecular marker (Area Under the Curve, AUC of 0.8249). Furthermore, significant correlations were found between the analyzed miRNAs and disease stage, treatment, β2 microglobulin, serum albumin and creatinine levels, but not with calcium levels or genetic alterations. In this cohort, 65.96% of MM patients had bone lesions, the majority of which were in the vertebrae. Additionally, miR-29c-3p was decreased in patients with osteolytic lesions compared with patients without bone disease. Interestingly, circulating levels of miR-29b-3p correlated with cervical and thoracic vertebral lesions, while miR-195-5p correlated with thoracic lesions. Our findings suggest circulating miRNAs can be promising biomarkers for MM diagnosis and that their levels correlate with myeloma bone disease and osteolytic lesions.

Graphical Abstract

1. Introduction

Multiple myeloma (MM) is a malignant hematological disease characterized by excessive proliferation of abnormal plasma cells (PCs) in the bone marrow [1]. It accounts for more than 10% of all hematologic cancers and it remains incurable in the majority of cases, with a 45% 5-year survival rate [2,3]. Up to 80% of the newly diagnosed MM patients have osteolytic lesions that develop as a consequence of the exacerbated proliferation of malignant PCs in the bone marrow, causing an unbalance between the activity of bone-forming cells (osteoblasts) and bone-resorbing cells (osteoclasts) [1,4]. This leads to the disruption of the normal bone remodeling process and to the development of myeloma bone disease (MBD) [5,6]. Consequently, MM patients are at high risk of skeletal-related incidents such as pathological fractures, spinal cord compression, vertebral collapse, severe bone pain and hypercalcemia [7]. Notably, the spine is the most frequent location of the bone lesions caused by MM and vertebral crush fractures are common [8]. Bisphosphonates are pyrophosphate analogues that cause apoptosis of osteoclasts and are currently the gold standard pharmacological approach to prevent and treat MBD [7,9]. More recently, denosumab, a humanized monoclonal antibody against RANKL, has also been used to decrease osteoclastogenesis and resorption by osteoclasts [7]. Imaging tools, such as whole-body low-dose computed tomography (CT), positron emission tomography-CT, and whole-body magnetic resonance imaging, are used to detect osteolytic lesions [10], and there are no currently available molecular biomarkers to detect MBD. Therefore, new tools for the early diagnosis of MM, as well as for the rapid identification of MBD, are needed. This will contribute to initiate the treatment at an early stage, avoiding further complications.
MicroRNAs (miRNAs) are a class of small RNA molecules that is transcribed from the genome but do not encode proteins [11]. In their mature form, miRNAs are about 19–22 nucleotides in length and are able to regulate the expression of messenger RNAs [11] and long non-coding RNAs [12]. Importantly, miRNA expression is disrupted in a vast number of diseases and hematological malignancies, such as MM [13,14]. Specifically, miRNAs are involved in MM onset and progression, which triggered their study as therapeutic targets [14]. Since the last two decades, circulating (extracellular) miRNAs have emerged as potential biomarkers for a wide range of diseases [11,15], including for MM [13]. Considering miRNAs can be detected in human blood serum/plasma, urine, saliva and other body fluids [11], these molecules are promising diagnosis tools through non-invasive or minimally invasive methods [15], leading to the development of several miRNA-based products for disease diagnosis. The stability of endogenous miRNAs in plasma samples is an important prerequisite for biomarkers [16]. It has been shown that miRNAs can be packed into extracellular vesicles [17] or can be associated with protein complexes [18], namely AGO proteins, protecting them from degradation by nucleases. Still, one of the main current challenges is to increase the sensibility and specificity of these circulating biomarkers [15].
The aim of this study is to identify circulating (extracellular) miRNAs to be used as biomarkers for the diagnosis of MM or MBD. For this purpose, a panel of miRNAs associated with MM onset and progression or with MBD pathophysiology, including regulation of osteoblasts and osteoclasts-related mechanisms, was evaluated. Levels of circulating miRNAs were analyzed in plasma samples and correlated with MM clinico-pathological characteristics and with bone lesions, including in the spine. Results show that circulating miRNAs are promising candidates for the diagnosis of MM and for the MM-associated skeletal osteolytic lesions.

2. Materials and Methods

2.1. Patients and Plasma Samples

Human peripheral blood from patients diagnosed with MM (n = 47) was collected at Oncology Department at Hospital de Braga (Braga, Portugal) between July 2018 and January 2021. Patients with other tumor types or inflammatory diseases were excluded from the study. Peripheral blood from healthy donors (n = 35) was collected at Immunohemotherapy Department of Centro Hospitalar de São João (CHSJ, Porto, Portugal). The blood donors used in this study do not have history of cancer or inflammatory diseases, had not been diagnosed with infectious diseases and were not under the effect of anti-inflammatory or anti-depressive drugs. A written informed consent was obtained from all subjects before blood collection. The patients’ clinical data was pseudonymized. The study was approved by Hospital de Braga and CHSJ Ethics Committee for Health, and the protocol conforms to the declaration of Helsinki. Clinical data was collected from medical records. Stage of the disease was determined following the International Staging System (ISS) for multiple myeloma, genetic alterations were assessed by cytogenetics/fluorescence in situ hybridization (FISH), and bone osteolytic lesions were assessed by magnetic resonance imaging, computerized tomography or conventional X-ray.
The blood samples were transferred from the EDTA collection tubes to 15 mL tubes, centrifuged at 1200× g, for 20 min, at room temperature (RT), and the plasma was collected. All samples were processed within three hours following blood collection to prevent hemolysis. No signs of colors from faint pink to bright red were detected in plasma samples following blood centrifugation. Plasma was further centrifuged at 2500× g, for 10 min, at 4 °C, to remove cell debris and stored at −80 °C in microtubes until further use.

2.2. Isolation and Purification of Circulating microRNAs

MiRNAs present in plasma samples were isolated and purified with miRNeasy Serum/Plasma Kit (Qiagen), according to the manufacturer’s instructions. Briefly, 1 mL of QIAzol Lysis Reagent was added to 200 µL of plasma. After vortexing and incubating at RT for 5 min, 2 µL (1.5 nM) of spike-in controls (cel-miR-39-3p and cel-miR-54-3p) were added. Then, 200 µL of chloroform was added, samples vigorously mixed, centrifuged at 12,000× g, for 15 min, at 4 °C, and the upper aqueous phase was collected to a new tube, with addition of 900 µL of ice-cold pure ethanol. Samples were transferred to the RNeasy MinElute spin columns, washed, and RNA was eluted in 15 µL of RNase-free water. Finally, RNA concentration was determined by measuring the absorbance in a NanoDrop Spectrophotometer ND-1000 (ThermoFisher Scientific, CA, USA).

2.3. Selection of the miRNA Panel, Reverse Transcription and Real Time Quantitative Polymerase Chain Reaction (RT-qPCR)

The expression levels of 10 miRNAs associated with MM onset and progression or with bone biological processes (miR-16-5p; miR-20a-5p; miR-21-5p; miR-29a-3p; miR-29b-3p; miR-29c-3p; miR-93-5p; miR-99a-5p; miR-146a-5p and miR-195-5p) were tested. In detail, selection of this panel was based on the following criteria: (1) miRNAs involvement in osteoblasts (bone forming cells) differentiation (miR-20a-5p, miR-29b-3p, miR-29c-3p, miR-146a-5p, miR-195-5p [19], miR-99a-5p [20], miR-93-5p [21], miR-29a-3p [22]), (2) miRNAs in multiple myeloma that mediate communication with cells in the tumor microenvironment (miR-16-5p [23], miR-21-5p [24]), (3) miRNAs dysregulated in MM cells or cell lines [25]. Importantly, several of the above mentioned miRNAs are concomitantly involved in two or more of these processes.
RNA was converted to cDNA through TaqMan® Advanced miRNA cDNA Synthesis Kit (Applied Biosystems), according to the manufacturer´s protocol. Firstly, for the poly(A) tailing reaction, RNA samples were mixed with Poly(A) Buffer, ATP, Poly(A) Enzyme and RNase-free water, briefly vortexed and incubated in a thermocycler (MyCycler Thermo Cycler, Bio-Rad, Hercules, CA, USA) at 37 °C for 45 min followed by 10 min at 65 °C. For the adaptor ligation reaction, DNA Ligase Buffer, 50% PEG 8000, Ligation Adaptor, RNA Ligase and RNase-free water were mixed with the previous reaction and incubated at 16 °C for 1 h. Then, the reverse transcription was performed by adding RT Buffer, dNTP, Universal RT Primers, RT Enzyme Mix and RNase-free water and incubated at 42 °C for 15 min and 5 min at 85 °C. Finally, for cDNA amplification, miR-Amp Master Mix, miR-Amp Primer Mix and RNase-free water was added to 5 µL of the RT reaction product and incubated with the following settings: enzyme activation for 5 min, at 95 °C; denaturation at 95 °C for 3 s and annealing/extension at 60 °C for 30 s, for 14 cycles; stop reaction at 99 °C for 10 min. Then, cDNA was diluted in 50 µL of RNase-free water and stored at −20 °C until further use.
For the RT-qPCR reaction, cDNA was mixed with TaqMan Fast Advanced Master Mix (Applied Biosystems), TaqMan Advanced miRNA Assay for miR-16-5p, miR-20a-5p, miR-21-5p, miR-29a-3p, miR-29b-3p, miR-29c-3p, miR-93-5p, miR-99a-5p, miR-146a-5p or miR-195-5p, and nuclease-free water in a 7500 Fast Real-Time PCR System (Applied Biosystems) and incubated as follows: 95 °C for 20 s; 40 cycles of 95 °C for 3 s and 60 °C for 30 s. RT-qPCR reactions were performed in duplicates. The results were analyzed in the 7500 Software v2.0.6 (Applied Biosystems). Exogenous cel-miR-39-3p and cel-miR-54-3p were used as normalizers. The relative gene expression was calculated using the ΔCq method, in accordance with MIQE guidelines [26].

2.4. Statistical Analysis

Statistical data analysis was performed using the GraphPad Prism v8.1.2 software (GraphPad Software Inc, San Diego, CA, USA) and SPSS Statistics v26 (IBM Corp, Armonk, NY, USA). D’Agostino and Pearson omnibus normality tests were performed to evaluate if data followed a Gaussian distribution. Considering that data did not follow a normal distribution, the non-parametric Mann–Whitney U test was used to evaluate significant differences between groups. Categorical data was tested using the Fisher’s exact test.
Relative expression of miRNAs was compared across groups using a linear regression, adjusted for the baseline characteristic identified as a possible confounder (age). The diagnostic value of the circulating miRNAs in plasma samples was calculated using Receiver Operating Characteristics (ROC) curve. To obtain the value of combined miRNAs expression a binary logistic regression was performed, considering the disease as the dependent variable and the miRNAs levels as covariates. For the ROC curve, the p-value tests the null hypothesis that the area under the curve equals 0.5. The cut-off points with the highest sensitivity and specificity were determined. Considering the data follows a non-parametric distribution, Spearman’s rank correlation was used to calculate correlation coefficient and significance between miRNA expression and clinicopathological parameters. Recursive partitioning (RP) was used to split the dataset iteratively based on certain statistical criteria. The RP classification trees were performed using rpart package from R software and Gini index criteria. Statistical significance was considered for p < 0.05.

3. Results

3.1. Patients Characteristics

Demographic, clinical and biochemical characteristics, treatment regimen, and genetic alterations in the patients’ cohort are shown in Table 1. Patients with MM were significantly older than healthy donors (p < 0.001), but no differences were found between groups regarding gender. The majority of the MM patients (53%) are at stage II, according to the International Staging System (ISS). In addition, IgG was the immunoglobulin most commonly produced by these patients (70%). Genetic alterations are present in MM patients, mainly P53 and RB1 deletion in 23% and 35% of the patients, respectively (Table 1).

3.2. Differences in the Expression of Circulating miRNAs between Multiple Myeloma Patients and Control Group

The expression of the selected 10 miRNAs was evaluated in the plasma of 47 MM patients and 35 healthy donors. Results show that miR-16-5p and miR-20a-5p were significantly downregulated, while levels of miR-21-5p and miR-29c-3p were significantly upregulated, in MM patients compared with the control group (Figure 1). Other members of the miR-29 family, namely miR-29a-3p and miR-29b-3p, were not statistically different between groups (Figure 1). Levels of four additional miRNAs, namely miR-93-3p, miR-99a-5p, miR-146a-5p and miR-195-5p, were also tested but no statistically significant differences were found between groups (Figure 1). Notably, miR-16-5p and miR-93-3p were among the most highly expressed miRNAs in plasma samples with median CT values of 13.96 and 19.84 cycles, respectively. However, miR-93-3p shows high variability, both in the MM and in the control group. On the other hand, miR-146a-5p was the miRNA with the lowest expression in the plasma, with a median Ct value of 28.77 cycles. Importantly, all tested miRNAs show Ct values below the threshold of 35 cycles.
Considering age as a variable that is significantly different between MM samples and the control group (Table 1), miRNAs miR-16-5p, miR-20a-5p, miR-21-5p and miR-29c-3p levels were tested again, adjusted for the age variable. The results from the multivariate linear regression analysis show that miR-16-5p, miR-20a-5p, and miR-21-5p maintain statistical significance, while miR-29c-3p differences could not be proven as age-independent (Table 2). Moreover, results from RP analysis show that the worst prognostic group refers to patients who have low expression of miR-20a-5p (miR-20a-5p < 0.011). The RP was subsequently used to maximize the variance across different nodes and minimize the variance within each node, allowing the identification of potential interactions between miRNAs biomarkers and MM disease. Using combined data, RP analysis identified different risk groups for MM disease (Figure 2). Only age-independent miRNAs differentially expressed between healthy individuals and MM patients were used for this analysis. The worst prognostic group referred to individuals who had low values of miR-20a-5p (94% out of miR-20a-5p < 0.011 individuals have MM disease). On the other hand, individuals with high levels of miR-20a-5p (miR-20a-5p ≥ 0.011), low levels of miR-21-5p (miR-21-5p < 0.0023) and high levels of miR-16-5p (miR-16-5p ≥ 4) or with high levels of miR-20a-5p (miR-20a-5p ≥ 0.011), low levels of miR-21-5p (miR-21-5p < 0.0023), low levels of miR-16-5p (miR-16-5p < 4) and low levels of miR-21-5p (miR-21-5p < 970 × 10−6) had the best prognosis (90% of these individuals are healthy) (Figure 2). Herein, among the tested miRNAs, miR-16-5p, miR-20a-5p, and miR-21-5p were identified as differently expressed in the plasma of MM patients compared with healthy donors.

3.3. Specificity and Sensibility of Circulting microRNAs as Multiple Myeloma Biomarkers

Next, attempting to find a miRNA or a combination of miRNAs that could successfully be used as a biomarker for MM, ROC analysis was performed for the miRNAs that were previously found differentially expressed between the MM group and the control group.
ROC curve analysis of miR-16-5p expression showed 78.72% sensitivity and 51.43% specificity in discriminating between MM and healthy donors, corresponding to an Area Under the Curve (AUC) of 0.665 (p < 0.05, 95% CI: 0.5458–0.7843). The miR-21-5p has an AUC of 0.654 (p < 0.05; sensitivity = 57.45%; specificity = 71.43%, 95% CI: 0.5365–0.7717). Considering the analysis of individual miRNAs, the highest AUC value was found for miR-20a-5p expression (0.707, p < 0.01; sensitivity = 36.17%; specificity = 97.14; 95% CI: 0.5961–0.8179). Nevertheless, combination of these three miRNAs improved the efficacy of the predictions, increasing both AUC and statistical significance. Specifically, combining miR-16-5p, miR-20a-5p and miR-21-5p improves AUC to 0.825 (p < 0.0001, sensitivity = 78.72%; specificity = 80.00%, 95% CI: 0.7333–0.9166) (Figure 3, Supplemental Table S1).

3.4. Correlation between Circulating microRNAs and Clinicopathological Variables in Multiple Myeloma Patients

Correlations of the miRNA levels in MM plasma samples with clinical and biochemical properties, treatment or chromosomal abnormalities was also tested. The miRNAs miR-20a-5p, miR-99a-5p and miR-195-5p negatively correlate with stage of the disease, following ISS classification (Table 3). Moreover, negative correlations were found between circulating miRNA levels and serum proteins. Among those, miR-16-5p, miR-20a-5p, miR-29b-3p, miR-29c-3p, miR-99a-5p and miR-195-5p were negatively correlated with β2M levels, a serum marker of tumor burden in hematologic malignancies [27,28], while miR-29a-3p was the only extracellular miRNA that negatively correlated with serum albumin (Table 3), a protein with prognostic value in MM patients [29,30]. Interestingly, besides miR-20a-5p, the three members of the miR-29 family showed a negative correlation with creatinine serum levels (Table 3), commonly used to measure the degree of insufficiency or renal failure [31]. No associations were found between the expression of circulating miRNAs and calcium levels (Table 3).
Regarding treatment of MM patients, a negative correlation was found between two members of the miR-29 family (namely miR-29b-3p and miR-29c-3p) and the use of chemotherapy or VTD regimen, while miR-93-5p and miR-146a-5p were negatively correlated with bone therapy (Table 4). Regarding genetic alterations, specifically del17p, del13q, translocation t(11;14), and translocation t(4;14), no significant correlations were found for any of the miRNAs (Table 4).

3.5. Correlation between Circulating microRNAs and Skeletal Osteolytic Lesions

Thirty one, out of the 47 MM patients included in the cohort, have bone lesions (65.96%). Of those, 55% have osteolytic lesions in the vertebrae, with or without lesions in other bone regions of the skeleton (Figure 4A). The most affected regions are the thoracic (80.77%; T1-T12) and the lumbar vertebrae (69.23%; L1-L5), while sacrum and cervical lesions account for 30.77% and 19.23%, respectively (Figure 4B). Considering lesions in spinal regions, 61.52% of the patients had lesions in at least two spinal regions, while only 7.69% had injuries in all the regions (Figure 4C).
Our results show that circulating miR-29c-3p was the only miRNA significantly different between MM patients with and without bone lesions. Patients with bone lesions have a downregulation of about 36% in the expression of miR-29c-3p compared with patients without bone lesions, independent of the location of the lesions (Figure 4D, Supplemental Table S2). In agreement, when analyzing all the miRNAs of the panel by RP, only miRNA-29c-3p was identified and the results show that the worst prognostic group refers to patients who have low values of miR-29c-3p (miR-29c-3p < 0.0074) (Figure 4E). Importantly, age and gender are not significantly different between these two groups, miR-29c-3p being an independent marker. Furthermore, ROC curve analysis revealed that miR-29c-3p has discriminatory ability as diagnostic biomarker (AUC = 0.695; p < 0.05) (Figure 4F). In agreement, miR-29c-3p is statistically correlated with bone lesions (r = −0.308; p = 0.042) (Table 5). When considering lesions in specific spinal regions, miR-29b-3p correlates with cervical and thoracic vertebral lesions, while miR-195-5p correlates with thoracic lesions only (Table 5). No correlations were found between miRNAs expression and lesions in the lumbar or sacrum regions.

4. Discussion

Early detection of MM and of MM-related bone disease may lead to faster interventions and more effective treatments. Circulating miRNAs are excellent candidates to be used as biomarkers for the early diagnosis of diseases through non-invasive or minimally invasive methods [11]. In this study, plasma samples were screened for the expression of 10 miRNAs that had been selected due to their role in MM onset/progression or in the bone remodeling process. Herein, miR-16-5p, miR-20a-5p and miR-21-5p were identified as potential biomarkers for MM, particularly when combined together. Interestingly, Rocci et al. identified miR-16 and miR-21 as consistently expressed in at least 20% of MM samples, and further associated low miR-16 expression with poor overall survival [32], while Qu et al. showed miR-20a as downregulated in 33 MM patients versus 20 control samples in a Chinese cohort [33]. An advantage of using a panel of miRNAs compared with individual miRNA levels is that it increases accuracy and diagnosis value, as shown for other tumor types, including breast [34,35], lung [36], and bladder [37], among others.
Importantly, differences in circulating miRNAs between MM patients and healthy controls do not necessarily translate the differential expression found in the malignant PCs versus normal PCs. Interestingly, the secretion of miRNAs to the plasma may be caused by the tumor cells themselves, by the cells in the tumor microenvironment, including immune cells, or by cells that have their function indirectly dysregulated by the tumor. The miR-16-5p is a tumor suppressor downregulated in MM-PCs, while miR-21-5p is an oncogene overexpressed in malignant PCs [25]. The same expression profile was identified in our cohort of plasma/blood samples from MM patients, compared with healthy controls, for these two miRNAs. However, this is not the case for miR-20a-5p. In vitro and in vivo studies showed that miR-20a-5p has an oncogenic role, being upregulated in MM-PCs [38,39]. Nevertheless, levels of this transcript were found to be decreased in the plasma samples from MM patients, suggesting alterations in miR-20a-5p secretion to the plasma may not be directly related with the changes in the malignant cells themselves, but rather related with changes in cells from the tumor microenvironment. This is in line with previous findings that miR-20a-5p levels are decreased in both the bone marrow microenvironment and in the peripheral blood of MM patients [40].
In this study, the association between circulating miRNAs and osteolytic lesions was also investigated, with a particular focus in the spine. The extracellular levels of miR-29c-3p were decreased in MM patients with bone disease compared with those without bone lesions. Recently, Papanota et al. identified let-7b-5p, miR-143-3p, miR-17-5p, miR-214-3p and miR-335-5p as associated with osteolytic bone disease in MM. Both studies are complementary as they use distinct miRNA panels, and contribute to reach the goal of finding novel non-invasive markers for osteolytic lesions [41]. Interestingly, Kapinas et al. showed that miR-29c-3p is highly expressed in fully mature and mineralized mouse-derived osteoblasts [42]. This miRNA is a member of the miR-29 family, known to directly target 16 extracellular matrix genes [43], including genes expressed in the bone. Although miR-29c-3p correlates with bone disease, no significant correlations were found when considering the distinct spinal regions. On the other hand, miR-29b-3p negatively correlated with vertebral lesions in the cervical and thoracic regions. This miRNA was shown to promote osteoblasts differentiation at early stages by directly targeting inhibitors of osteogenesis, while it reduces collagen levels at late stages to facilitate maturation of the collagen fibrillar matrix and, consequently, mineral deposition [44]. Circulating levels of miR-29b-3p were previously found to correlate with histomorphometric parameters in the bone [45]. However, to the best of our knowledge, this is the first study correlating extracellular miRNA levels and osteolytic lesions in the spine caused by MM.

5. Conclusions

Globally, we identified a promising combination of miRNAs in the plasma as biomarkers for MM. In addition, we specifically correlated circulating miRNA levels with the bone disease caused by MM, particularly with spinal osteolytic lesions. Future studies should validate these results in a larger multi-centric cohort of patients.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13215258/s1, Table S1: Receiver Operating Characteristic (ROC) curve analysis for circulating miR-16-5p, miR-20a-5p and miR-21-5p, or their combination (miRNA panel), when analyzing MM patients versus healthy controls. Table S2: Statistical differences in miRNA levels between multiple myeloma patients with and without bone lesions.

Author Contributions

Conceptualization, M.A.B., H.M. and M.I.A.; methodology, S.R.M., H.A. and C.C.; data collection, S.R.M., H.A., H.M. and M.I.A.; statistical analysis, S.R.M., C.R.-M. and C.O.; writing—original draft preparation, S.R.M., H.A. and M.I.A.; writing—review and editing, S.R.M., H.A., C.C., C.R.-M., C.O., M.A.B., H.M. and M.I.A.; funding acquisition, M.I.A. All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by AO Spine-ESA Grant Award 2018 (AO Foundation); Associação Portuguesa Contra a Leucemia, Sociedade Portuguesa de Hematologia, AMGEN; and by Portuguese funds through FCT-Fundação para a Ciência e a Tecnologia (FCT)/Ministério da Ciência, Tecnologia e Ensino Superior in the framework of the project POCI-01-0145-FEDER-031402- R2Bone (FEDER-Fundo Europeu de Desenvolvimento Regional funds through the COMPETE 2020-Operacional Programme for Competitiveness and Internationalisation—POCI, Portugal 2020).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Ethics Committee of Centro Hospitalar São João and Hospital de Braga, Portugal (Ref. 189-18, 89-21).

Informed Consent Statement

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

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical reasons.

Acknowledgments

The authors would like to thank Hospital de Dia de Oncologia, at Hospital de Braga, for kindly donating the samples of multiple myeloma patients used to conduct this study; Serviço de Imunohemoterapia, at Centro Hospitalar São João (CHSJ), for kindly providing the samples of healthy donors; IPO Porto, particularly to Rui Medeiros and Mara Fernandes for assistance with samples collection; the i3s Genomics service; and to Porto.Comprehensive Cancer Centre. SRM, CC and MIA are supported by FCT (SFRH/BD/147229/2019 and BiotechHealth Program; CEECIND/00184/2017/CP1392/CT0001; CEECINST/00091/2018/CP1500/CT0011, respectively).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. García-Ortiz, A.; Rodríguez-García, Y.; Encinas, J.; Maroto-Martín, E.; Castellano, E.; Teixidó, J.; Martínez-López, J. The Role of Tumor Microenvironment in Multiple Myeloma Development and Progression. Cancers 2021, 13, 217. [Google Scholar] [CrossRef]
  2. Oldham, R.A.A.; Faber, M.L.; Keppel, T.R.; Buchberger, A.R.; Waas, M.; Hari, P.; Gundry, R.L.; Medin, J.A. Discovery and validation of surface N-glycoproteins in MM cell lines and patient samples uncovers immunotherapy targets. J. Immunother. Cancer 2020, 8, e000915. [Google Scholar] [CrossRef] [PubMed]
  3. Kazandjian, D.; Branch, L.M. Multiple Myeloma Epidemiology and Survival, a Unique Malignancy. Semin. Oncol. 2017, 43, 676–681. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Terpos, E.; Ntanasis-Stathopoulos, I.; Gavriatopoulou, M.; Dimopoulos, M.A. Pathogenesis of bone disease in multiple myeloma: From bench to bedside. Blood Cancer J. 2018, 8, 1–12. [Google Scholar] [CrossRef] [Green Version]
  5. Terpos, E.; Dimopoulos, M.-A. Myeloma bone disease: Pathophysiology and management. Ann. Oncol. 2005, 16, 1223–1231. [Google Scholar] [CrossRef] [PubMed]
  6. Visconti, R.J.; Kolaja, K.; Cottrell, J.A. A functional three-dimensional microphysiological human model of myeloma bone disease. J. Bone Miner. Res. 2021. [Google Scholar] [CrossRef]
  7. Parrondo, R.D.; Sher, T. Prevention Of Skeletal Related Events In Multiple Myeloma: Focus On The RANK-L Pathway In The Treatment Of Multiple Myeloma. OncoTargets Ther. 2019, 12, 8467–8478. [Google Scholar] [CrossRef] [Green Version]
  8. Lasocki, A.; Gaillard, F.; Harrison, S.J. Multiple myeloma of the spine. Neuroradiol. J. 2017, 30, 259–268. [Google Scholar] [CrossRef]
  9. Terpos, E.; Sezer, O.; Croucher, P.I.; García-Sanz, R.; Boccadoro, M.; Miguel, J.S.; Ashcroft, J.; Bladé, J.; Cavo, M.; Delforge, M.; et al. The use of bisphosphonates in multiple myeloma: Recommendations of an expert panel on behalf of the European Myeloma Network. Ann. Oncol. 2009, 20, 1303–1317. [Google Scholar] [CrossRef]
  10. Zamagni, E.; Tacchetti, P.; Cavo, M. Imaging in multiple myeloma: How? When? Blood 2019, 133, 644–651. [Google Scholar] [CrossRef] [Green Version]
  11. Almeida, M.I.; Reis, R.M.; Calin, G.A. MicroRNA history: Discovery, recent applications, and next frontiers. Mutat. Res. Mol. Mech. Mutagen. 2011, 717, 1–8. [Google Scholar] [CrossRef] [Green Version]
  12. Almeida, M.I.; Reis, R.M.; Calin, G.A. Decoy activity through microRNAs: The therapeutic implications. Expert Opin. Biol. Ther. 2012, 12, 1153–1159. [Google Scholar] [CrossRef] [Green Version]
  13. Soliman, A.M.; Lin, T.S.; Mahakkanukrauh, P.; Das, S. Role of microRNAs in Diagnosis, Prognosis and Management of Multiple Myeloma. Int. J. Mol. Sci. 2020, 21, 7539. [Google Scholar] [CrossRef] [PubMed]
  14. Handa, H.; Murakami, Y.; Ishihara, R.; Kimura-Masuda, K.; Masuda, Y. The Role and Function of microRNA in the Pathogenesis of Multiple Myeloma. Cancers 2019, 11, 1738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Mori, M.A.; Ludwig, R.G.; Garcia-Martin, R.; Brandão, B.B.; Kahn, C.R. Extracellular miRNAs: From Biomarkers to Mediators of Physiology and Disease. Cell Metab. 2019, 30, 656–673. [Google Scholar] [CrossRef] [PubMed]
  16. Mitchell, P.S.; Parkin, R.K.; Kroh, E.M.; Fritz, B.R.; Wyman, S.K.; Pogosova-Agadjanyan, E.L.; Peterson, A.; Noteboom, J.; O’Briant, K.C.; Allen, A.; et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc. Natl. Acad. Sci. USA 2008, 105, 10513–10518. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Silva, A.M.; Teixeira, J.H.; Almeida, M.I.; Gonçalves, R.M.; Barbosa, M.A.; Santos, S.G. Extracellular Vesicles: Immunomodulatory messengers in the context of tissue repair/regeneration. Eur. J. Pharm. Sci. 2017, 98, 86–95. [Google Scholar] [CrossRef] [PubMed]
  18. Turchinovich, A.; Weiz, L.; Langheinz, A.; Burwinkel, B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011, 39, 7223–7233. [Google Scholar] [CrossRef] [PubMed]
  19. Almeida, M.I.; Silva, A.M.; Vasconcelos, D.M.; Almeida, C.R.; Caires, H.R.; Pinto, M.T.; Calin, G.A.; Santos, S.G.; Barbosa, M.A. miR-195 in human primary mesenchymal stromal/stem cells regulates proliferation, osteogenesis and paracrine effect on angiogenesis. Oncotarget 2016, 7, 7–22. [Google Scholar] [CrossRef] [Green Version]
  20. Moura, S.R.; Brás, J.P.; Freitas, J.; Osório, H.; Barbosa, M.A.; Santos, S.G.; Almeida, M.I. miR-99a in bone homeostasis: Regulating osteogenic lineage commitment and osteoclast differentiation. Bone 2020, 134, 115303. [Google Scholar] [CrossRef]
  21. Yang, L.; Cheng, P.; Chen, C.; He, H.-B.; Xie, G.-Q.; Zhou, H.-D.; Xie, H.; Wu, X.-P.; Luo, X.-H. miR-93/Sp7 function loop mediates osteoblast mineralization. J. Bone Miner. Res. 2012, 27, 1598–1606. [Google Scholar] [CrossRef] [PubMed]
  22. Tan, K.; Peng, Y.-T.; Guo, P. MiR-29a promotes osteogenic differentiation of mesenchymal stem cells via targeting HDAC4. Eur. Rev. Med. Pharmacol. Sci. 2018, 22, 3318–3326. [Google Scholar] [PubMed]
  23. Khalife, J.; Ghose, J.; Martella, M.; Viola, D.; Rocci, A.; Troadec, E.; Terrazas, C.; Satoskar, A.R.; Gunes, E.G.; Dona, A.; et al. MiR-16 regulates crosstalk in NF-κB tolerogenic inflammatory signaling between myeloma cells and bone marrow macrophages. JCI Insight 2019, 4. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, J.-H.; Zheng, W.-W.; Cheng, S.-T.; Liu, B.-X.; Liu, F.-R.; Song, J.-Q. Correlation between microRNA-21 and sprouty homolog 2 gene expression in multiple myeloma. Mol. Med. Rep. 2015, 11, 4220–4224. [Google Scholar] [CrossRef] [Green Version]
  25. Federico, C.; Sacco, A.; Belotti, A.; Ribolla, R.; Cancelli, V.; Giacomini, A.; Ronca, R.; Chiarini, M.; Imberti, L.; Marini, M.; et al. Circulating microRNAs and Their Role in Multiple Myeloma. Non-Coding RNA 2019, 5, 37. [Google Scholar] [CrossRef] [Green Version]
  26. Bustin, S.A.; Benes, V.; Garson, J.A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M.W.; Shipley, G.L.; et al. The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments. Clin. Chem. 2009, 55, 611–622. [Google Scholar] [CrossRef] [Green Version]
  27. Tsimberidou, A.M.; Kantarjian, H.M.; Wen, S.; O’Brien, S.; Cortes, J.; Wierda, W.G.; Koller, C.; Pierce, S.; Brandt, M.; Freireich, E.J.; et al. The Prognostic Significance of Serum β2 Microglobulin Levels in Acute Myeloid Leukemia and Prognostic Scores Predicting Survival: Analysis of 1,180 Patients. Clin. Cancer Res. 2008, 14, 721–730. [Google Scholar] [CrossRef] [Green Version]
  28. Litam, P.; Swan, F.; Cabanillas, F.; Tucker, S.L.; McLaughlin, P.; Hagemeister, F.B.; Rodríguez, M.A.; Velasquez, W.S. Prognostic Value of Serum ?-2 Microglobulin in Low-Grade Lymphoma. Ann. Intern. Med. 1991, 114, 855. [Google Scholar] [CrossRef]
  29. Kim, J.E.; Yoo, C.; Lee, D.H.; Kim, S.-W.; Lee, J.-S.; Suh, C. Serum albumin level is a significant prognostic factor reflecting disease severity in symptomatic multiple myeloma. Ann. Hematol. 2010, 89, 391–397. [Google Scholar] [CrossRef]
  30. Chen, Y.H. Hypoalbuminemia in patients with multiple myeloma. Arch. Intern. Med. 1990, 150, 605–610. [Google Scholar] [CrossRef]
  31. Dimopoulos, M.A.; Kastritis, E.; Rosiñol, L.; Blade, J.; Ludwig, H.; Blad, J. Pathogenesis and treatment of renal failure in multiple myeloma. Leukemia 2008, 22, 1485–1493. [Google Scholar] [CrossRef] [Green Version]
  32. Rocci, A.; Hofmeister, C.C.; Geyer, S.; Stiff, A.; Gambella, M.; Cascione, L.; Guan, J.; Benson, D.M.; Efebera, Y.A.; Talabere, T.; et al. Circulating miRNA markers show promise as new prognosticators for multiple myeloma. Leukemia 2014, 28, 1922–1926. [Google Scholar] [CrossRef] [Green Version]
  33. Qu, X.; Zhao, M.; Wu, S.; Yu, W.; Xu, J.; Xu, J.; Li, J.; Chen, L. Circulating microRNA 483-5p as a novel biomarker for diagnosis survival prediction in multiple myeloma. Med. Oncol. 2014, 31, 219. [Google Scholar] [CrossRef]
  34. Li, M.; Zou, X.; Xia, T.; Wang, T.; Liu, P.; Zhou, X.; Wang, S.; Zhu, W. A five-miRNA panel in plasma was identified for breast cancer diagnosis. Cancer Med. 2019, 8, 7006–7017. [Google Scholar] [CrossRef] [Green Version]
  35. Adam-Artigues, A.; Garrido-Cano, I.; Carbonell-Asins, J.A.; Lameirinhas, A.; Simón, S.; Ortega-Morillo, B.; Martínez, M.T.; Hernando, C.; Constâncio, V.; Burgues, O.; et al. Identification of a Two-MicroRNA Signature in Plasma as a Novel Biomarker for Very Early Diagnosis of Breast Cancer. Cancers 2021, 13, 2848. [Google Scholar] [CrossRef]
  36. Zhou, X.; Wen, W.; Shan, X.; Zhu, W.; Xu, J.; Guo, R.; Cheng, W.; Wang, F.; Qi, L.-W.; Chen, Y.; et al. A six-microRNA panel in plasma was identified as a potential biomarker for lung adenocarcinoma diagnosis. Oncotarget 2017, 8, 6513–6525. [Google Scholar] [CrossRef]
  37. Usuba, W.; Urabe, F.; Yamamoto, Y.; Matsuzaki, J.; Sasaki, H.; Ichikawa, M.; Takizawa, S.; Aoki, Y.; Niida, S.; Kato, K.; et al. Circulating miRNA panels for specific and early detection in bladder cancer. Cancer Sci. 2018, 110, 408–419. [Google Scholar] [CrossRef]
  38. Wang, T.; Tao, W.; Zhang, L.; Li, S. Oncogenic role of microRNA-20a in human multiple myeloma. OncoTargets Ther. 2017, 10, 4465–4474. [Google Scholar] [CrossRef] [Green Version]
  39. Jiang, Y.; Chang, H.; Chen, G. Effects of microRNA-20a on the proliferation, migration and apoptosis of multiple myeloma via the PTEN/PI3K/AKT signaling pathway. Oncol. Lett. 2018, 15, 10001–10007. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, W.; Corrigan-Cummins, M.; Barber, E.A.; Saleh, L.M.; Zingone, A.; Ghafoor, A.; Costello, R.; Zhang, Y.; Kurlander, R.J.; Korde, N.; et al. Aberrant Levels of miRNAs in Bone Marrow Microenvironment and Peripheral Blood of Myeloma Patients and Disease Progression. J. Mol. Diagn. 2015, 17, 669–678. [Google Scholar] [CrossRef]
  41. Papanota, A.-M.; Tsiakanikas, P.; Kontos, C.; Malandrakis, P.; Liacos, C.-I.; Ntanasis-Stathopoulos, I.; Kanellias, N.; Gavriatopoulou, M.; Kastritis, E.; Avgeris, M.; et al. A Molecular Signature of Circulating MicroRNA Can Predict Osteolytic Bone Disease in Multiple Myeloma. Cancers 2021, 13, 3877. [Google Scholar] [CrossRef] [PubMed]
  42. Kapinas, K.; Kessler, C.B.; Delany, A.M. miR-29 suppression of osteonectin in osteoblasts: Regulation during differentiation and by canonical Wnt signaling. J. Cell. Biochem. 2009, 108, 216–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kriegel, A.J.; Liu, Y.; Fang, Y.; Ding, X.; Liang, M. The miR-29 family: Genomics, cell biology, and relevance to renal and cardiovascular injury. Physiol. Genom. 2012, 44, 237–244. [Google Scholar] [CrossRef] [PubMed]
  44. Li, Z.; Hassan, M.Q.; Jafferji, M.; Aqeilan, R.I.; Garzon, R.; Croce, C.M.; van Wijnen, A.J.; Stein, J.L.; Stein, G.S.; Lian, J.B. Biological Functions of miR-29b Contribute to Positive Regulation of Osteoblast Differentiation. J. Biol. Chem. 2009, 284, 15676–15684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Feichtinger, X.; Muschitz, C.; Heimel, P.; Baierl, A.; Fahrleitner-Pammer, A.; Redl, H.; Resch, H.; Geiger, E.; Skalicky, S.; Dormann, R.; et al. Bone-related Circulating MicroRNAs miR-29b-3p, miR-550a-3p, and miR-324-3p and their Association to Bone Microstructure and Histomorphometry. Sci. Rep. 2018, 8, 1–11. [Google Scholar] [CrossRef] [PubMed]
Cancers 13 05258 g001 550
Figure 1. Expression levels of circulating microRNAs in plasma samples from multiple myeloma patients (MM, n = 47) and healthy donors (HD, n = 35) (median, interquartile range; Mann–Whitney test, * p < 0.05 and ** p < 0.01).
Figure 1. Expression levels of circulating microRNAs in plasma samples from multiple myeloma patients (MM, n = 47) and healthy donors (HD, n = 35) (median, interquartile range; Mann–Whitney test, * p < 0.05 and ** p < 0.01).
Cancers 13 05258 g001
Cancers 13 05258 g002 550
Figure 2. Decision tree created with 4 splits of the data. The root node (Node 1) contains all the data in the dataset. The algorithm splits the data based on a defined statistic, creating two new nodes (Node 2 and Node 3). Using the same statistic, it splits the data again at Node 2, creating two more leaf nodes (Nodes 4 and 5). For instance, Node 4 splits into Node 8 and 9; and Node 5 splits in two leaf nodes (Node 10 and 11). Lastly, it splits the data at Node 9, creating two more nodes (Nodes 18 and 19). The decision tree makes its prediction for each data row by traversing to the leaf nodes (one of the terminal nodes: Node 8, 18, 19, 10, 11 or 3).
Figure 2. Decision tree created with 4 splits of the data. The root node (Node 1) contains all the data in the dataset. The algorithm splits the data based on a defined statistic, creating two new nodes (Node 2 and Node 3). Using the same statistic, it splits the data again at Node 2, creating two more leaf nodes (Nodes 4 and 5). For instance, Node 4 splits into Node 8 and 9; and Node 5 splits in two leaf nodes (Node 10 and 11). Lastly, it splits the data at Node 9, creating two more nodes (Nodes 18 and 19). The decision tree makes its prediction for each data row by traversing to the leaf nodes (one of the terminal nodes: Node 8, 18, 19, 10, 11 or 3).
Cancers 13 05258 g002
Cancers 13 05258 g003 550
Figure 3. Receiver Operating Characteristic (ROC) curve of miR-16-5p, miR-20a-5p, miR-21-5p or combined panel of these 3 miRNAs between multiple myeloma patients and healthy controls (AUC: area under the ROC curve).
Figure 3. Receiver Operating Characteristic (ROC) curve of miR-16-5p, miR-20a-5p, miR-21-5p or combined panel of these 3 miRNAs between multiple myeloma patients and healthy controls (AUC: area under the ROC curve).
Cancers 13 05258 g003
Cancers 13 05258 g004a 550Cancers 13 05258 g004b 550
Figure 4. Bone lesions in multiple myeloma patients and miR-29c-3p expression. (A) Bone injury sites in multiple myeloma patients; (B) percentage of multiple myeloma patients with spinal lesions per region (cervical, thoracic, lumbar or sacrum); (C) percentage of multiple myeloma patients with spinal lesions in more than one region; (D) expression of circulating miR-29c-3p in multiple myeloma patients with and without bone lesions (median, interquartile range; Mann-Whitney test, * p < 0.05); (E) decision tree created with 1 split of the data. The root node (Node 1) contains all the data in the dataset. The algorithm splits the data based on a defined statistic, creating two new nodes (Node 2 and Node 3). MBD—myeloma bone disease, and MWBD—myeloma without bone disease; (F) Receiver Operating Characteristic curve of miR-29c-3p between patients with and without bone lesions (AUC: area under the ROC curve).
Figure 4. Bone lesions in multiple myeloma patients and miR-29c-3p expression. (A) Bone injury sites in multiple myeloma patients; (B) percentage of multiple myeloma patients with spinal lesions per region (cervical, thoracic, lumbar or sacrum); (C) percentage of multiple myeloma patients with spinal lesions in more than one region; (D) expression of circulating miR-29c-3p in multiple myeloma patients with and without bone lesions (median, interquartile range; Mann-Whitney test, * p < 0.05); (E) decision tree created with 1 split of the data. The root node (Node 1) contains all the data in the dataset. The algorithm splits the data based on a defined statistic, creating two new nodes (Node 2 and Node 3). MBD—myeloma bone disease, and MWBD—myeloma without bone disease; (F) Receiver Operating Characteristic curve of miR-29c-3p between patients with and without bone lesions (AUC: area under the ROC curve).
Cancers 13 05258 g004aCancers 13 05258 g004b
Table
Table 1. Characterization of multiple myeloma patients (MM) and healthy donors.
Table 1. Characterization of multiple myeloma patients (MM) and healthy donors.
Healthy DonorsMMp-Value
Demographic characteristics
Age (years)39.00 (24–60)71.00 (47–89)<0.001
Gender (Male-Female) (%)45.7–54.366.0–34.00.08
Clinical and biochemical characteristics
ISS stage (I-II-III) (%)N/A23.53–52.94–23.53
Ig isotype (IgA-IgD-IgE-IgG-IgM) (%)23.91-0-0-69.57-0
Light chains (%)6.52
β2-microglobulin (ng/mL)4389(1386–29094)
Albumin (g/dL)3.35 (2.3–4.2)
Creatinine (mg/dL)1.1 (0.5–6.3)
Calcium (mg/dL)8.8 (6.4–11.8)
Treatment
Chemotherapy (%)N/A38.10
Radiotherapy (%)13.64
Bone therapeutic (%)28.57
Genetic alterations
del(17)(p13.1) (%)N/A23.26
del(13)(q14.3) (%)34.88
t(11;14) (%)13.64
t(4;14) (%)2.33
For continuous variables: values represent median (min value–max value) and differences were assessed using non-parametric Mann-Whitney test; for categorial variables: values represent percentages and differences were assessed using Fisher's exact test; N/A: non-applicable.
Table
Table 2. Statistical differences in microRNA levels between multiple myeloma patients and healthy donors, when adjusted to age.
Table 2. Statistical differences in microRNA levels between multiple myeloma patients and healthy donors, when adjusted to age.
miRNAp-Value (Adjusted to Age)
miR-16-5p0.027
miR-20a-5p0.024
miR-21-5p0.031
miR-29c-3p0.334
Table
Table 3. Correlation between circulating microRNAs and stage of the disease or serum protein levels in multiple myeloma patients. For each microRNA, the Spearman’s correlation coefficient (r) and significance level (p-value) are shown. Statistically significant correlations are highlighted in blue.
Table 3. Correlation between circulating microRNAs and stage of the disease or serum protein levels in multiple myeloma patients. For each microRNA, the Spearman’s correlation coefficient (r) and significance level (p-value) are shown. Statistically significant correlations are highlighted in blue.
Stage (ISS)Serum Proteins
IIIIIIβ2M
(ng/mL)		Albumin (g/dL)Creatinine (mg/dL)Total Calcium (mg/dL)
miR-16-5pr0.1550.102−0.276−0.414−0.178−0.2870.053
p0.3800.5660.1150.0130.2860.0770.756
miR-20a-5pr0.2830.078−0.375−0.565−0.183−0.4610.043
p0.1050.6610.0290.0000.2710.0030.798
miR-21-5pr−0.1700.0960.0570.146−0.296−0.150−0.175
p0.3380.5890.7510.4020.0710.3620.300
miR-29a-3pr−0.0070.156−0.177−0.315−0.349−0.3280.042
p0.9680.3780.3170.0650.0320.0420.804
miR-29b-3pr0.2330.066−0.311−0.499−0.124−0.4140.054
p0.1840.7100.0730.0020.4640.0100.754
miR-29c-3pr0.198−0.060-0.127−0.419−0.157−0.326−0.060
p0.2620.7360.4730.0120.3470.0430.726
miR-93-5pr0.2830.018−0.304−0.275−0.041−0.039−0.123
p0.1050.9190.0810.1110.8050.8150.469
miR-99a-5pr0.0000.312−0.368−0.357−0.197−0.0070.013
p1.0000.0720.0320.0350.2350.9660.938
miR-146a-5pr−0.2050.240−0.0780.082−0.196−0.015−0.085
p0.2450.1710.6620.6420.2390.9270.618
miR-195-5pr0.2540.078−0.346−0.344−0.109−0.238−0.194
p0.1460.6610.0450.0430.5130.1440.250
Table
Table 4. Correlation between circulating miRNAs and treatment regimens or chromosomal abnormalities in multiple myeloma patients. For each miRNA, the Spearman’s correlation coefficient (r) and significance level (p-value) are shown. Statistically significant correlations are highlighted in blue.
Table 4. Correlation between circulating miRNAs and treatment regimens or chromosomal abnormalities in multiple myeloma patients. For each miRNA, the Spearman’s correlation coefficient (r) and significance level (p-value) are shown. Statistically significant correlations are highlighted in blue.
TreatmentGenetic Alterations
Chemotherapy or
VTd Regimen		RadiotherapyBone
Therapy		del(17)(p13.1)del(13)(q14.3)t(11;14)
Translocation		t(4;14)
Translocation
miR-16-5pr−0.150−0.1620.2160.2530.016−0.114−0.114
p0.3440.2940.1640.1020.9200.4680.461
miR-20a-5pr−0.138−0.0890.1140.271−0.063−0.081−0.078
p0.3850.5670.4660.0790.6890.6050.615
miR-21-5pr−0.121−0.063−0.155−0.200−0.0940.027−0.162
p0.4440.6870.3210.1990.5470.8630.293
miR-29a-3pr−0.129−0.047−0.0080.2220.031−0.027−0.210
p0.4140.7620.9590.1530.8410.8630.171
miR-29b-3pr−0.3850.0380.0790.065−0.163−0.056−0.149
p0.0130.8100.6210.6850.3040.7240.340
miR-29c-3pr−0.360−0.0990.0530.1150.0470.011−0.186
p0.0190.5220.7350.4610.7640.9450.226
miR-93-5pr−0.174−0.188−0.461−0.146−0.1340.0810.174
p0.2710.2220.0020.3490.3930.6050.258
miR-99a-5pr−0.231−0.016−0.261−0.031−0.181−0.0220.066
p0.1420.9200.0910.8430.2460.8900.670
miR-146a-5pr−0.170−0.130−0.510−0.0750.035−0.0430.030
p0.2820.399<0.0010.6310.8220.7830.847
miR-195-5pr−0.162−0.005−0.0730.129−0.295−0.011−0.078
p0.3060.9730.6400.4110.0550.9450.615
Table
Table 5. Correlation between circulating microRNAs and bone disease or spinal osteolytic lesions in multiple myeloma patients. For each miRNA, the Spearman’s correlation coefficient (r) and significance level (p-value) are shown. Statistically significant correlations are highlighted in blue.
Table 5. Correlation between circulating microRNAs and bone disease or spinal osteolytic lesions in multiple myeloma patients. For each miRNA, the Spearman’s correlation coefficient (r) and significance level (p-value) are shown. Statistically significant correlations are highlighted in blue.
Bone
Disease		Regions of the Spine
CervicalThoracicLumbarSacrum
miR-16-5pr−0.1470.1180.127−0.1630.175
p0.3410.4560.4150.3030.267
miR-20a-5pr−0.0730.0760.146−0.1190.155
p0.6400.6330.3490.4530.327
miR-21-5pr0.057−0.270−0.0070.2020.060
p0.7140.0840.9620.1990.706
miR-29a-3pr−0.041−0.112−0.022−0.1030.270
p0.7910.4790.8860.5150.084
miR-29b-3pr−0.135−0.340−0.304−0.2050.153
p0.3890.0300.0500.1980.338
miR-29c-3pr−0.308−0.215−0.191−0.226−0.125
p0.0420.1710.2190.1500.430
miR-93-5pr−0.100−0.185−0.0790.0870.115
p0.5180.2410.6160.5820.468
miR-99a-5pr−0.151−0.215−0.120−0.0080.285
p0.3280.1710.4440.9600.067
miR-146a-5pr−0.226−0.112−0.0150.032−0.025
p0.1410.4790.9240.8420.875
miR-195-5pr−0.080−0.143−0.3070.0000.095
p0.6040.3680.0451.0000.549
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Moura, S.R.; Abreu, H.; Cunha, C.; Ribeiro-Machado, C.; Oliveira, C.; Barbosa, M.A.; Marques, H.; Almeida, M.I. Circulating microRNAs Correlate with Multiple Myeloma and Skeletal Osteolytic Lesions. Cancers 2021, 13, 5258. https://doi.org/10.3390/cancers13215258

AMA Style

Moura SR, Abreu H, Cunha C, Ribeiro-Machado C, Oliveira C, Barbosa MA, Marques H, Almeida MI. Circulating microRNAs Correlate with Multiple Myeloma and Skeletal Osteolytic Lesions. Cancers. 2021; 13(21):5258. https://doi.org/10.3390/cancers13215258

Chicago/Turabian Style

Moura, Sara Reis, Hugo Abreu, Carla Cunha, Cláudia Ribeiro-Machado, Carla Oliveira, Mario Adolfo Barbosa, Herlander Marques, and Maria Inês Almeida. 2021. "Circulating microRNAs Correlate with Multiple Myeloma and Skeletal Osteolytic Lesions" Cancers 13, no. 21: 5258. https://doi.org/10.3390/cancers13215258

APA Style

Moura, S. R., Abreu, H., Cunha, C., Ribeiro-Machado, C., Oliveira, C., Barbosa, M. A., Marques, H., & Almeida, M. I. (2021). Circulating microRNAs Correlate with Multiple Myeloma and Skeletal Osteolytic Lesions. Cancers, 13(21), 5258. https://doi.org/10.3390/cancers13215258

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