Proteomic analysis of the bone marrow microenvironment in extramedullary multiple myeloma patients.

Multiple myeloma (MM) is a heterogeneous hematological malignancy characterized by the uncontrolled clonal proliferation of bone marrow (BM) plasma cells. The poor prognosis of patients is associated with the presence of extramedullary disease (EMD). Previously, different mechanisms involved in the colonization of BM niches by MM cells and their escape during EMD have been described. Thus, we aimed to investigate the expression of selected cytokines in the BM plasma of MM patients as well as EMD patients to reveal novel molecules involved in EMD pathogenesis. Expression of 120 different cytokines was measured in BM plasma of 13 MM and 11 EMD patients using Proteome Profiler Antibody Arrays. The correlation between statistically significant cytokines and clinicopathological parameters of patients was determined using the Spearman correlation analysis. Finally, protein-protein interactions were analyzed, and GO and KEGG pathways enrichment analysis was performed. In total, 27 cytokines were found to be differently expressed between MM and EMD patients. After the Benjamini-Hochberg correction for multiple testing, the statistical significance of two cytokines downregulated in EMD (EGF, BDNF) and six cytokines upregulated in EMD (NAP-2, ADIPOQ, CRP, MIG, BAFF, and THBS1) was maintained. Correlation analysis proved a significant association between the expression of these molecules and selected clinical-pathological features of MM/EMD patients. Protein association network analysis revealed important protein-protein interactions between THBS1/EGF, MIG/NAP-2, THBS1/NAP-2, EGF/NAP-2, and ADIPOQ/CRP. Finally, identified cytokines were proved to be significantly involved in focal adhesion, PI3K/AKT, and MAPK signaling pathways, and regulation of cell development, localization, proliferation, migration, differentiation, immune system processes, and stress response. Obtained results confirm the key function of the BM microenvironment in the pathogenesis of MM and indicate the essential role of numerous cytokines in disease progression and EMD development. However, the exact mechanisms need to be further clarified.

Involvement of Small Non-Coding RNA and Cell Antigens in Pathogenesis of Extramedullary Multiple Myeloma.

Extramedullary multiple myeloma (EMD) is an aggressive disease; malignant plasma cells lose their dependence in the bone marrow microenvironment and migrate into tissues. EMD is a negative prognostic factor of survival. Using flow cytometry and next-generation sequencing, we aimed to identify antigens and microRNAs (miRNAs) involved in EMD pathogenesis. Flow cytometry analysis revealed significant differences in the level of clonal plasma cells between MM and EMD patients, while the expression of CD markers was comparable between these two groups. Further, miR-26a-5p and miR-30e-5p were found to be significantly down-regulated in EMD compared to MM. Based on the expression of miR-26a-5p, we were able to distinguish these two groups of patients with high sensitivity and specificity. In addition, the involvement of deregulated miRNAs in cell cycle regulation, ubiquitin-mediated proteolysis and signaling pathways associated with infections or neurological disorders was observed using GO and KEGG pathways enrichment analysis. Subsequently, a correlation between the expression of analyzed miRNAs and the levels of CD molecules was observed. Finally, clinicopathological characteristics as well as CD antigens associated with the prognosis of MM and EMD patients were identified. Altogether, we identified several molecules possibly involved in the transformation of MM into EMD.

Necessity of flow cytometry assessment of circulating plasma cells and its connection with clinical characteristics of primary and secondary plasma cell leukaemia.

Plasma cell leukaemia (PCL) is a rare and very aggressive plasma cell disorder. Preventing a dismal outcome of PCL requires early diagnosis with appropriate analytical tools. Therefore, the investigation of 33 patients with primary and secondary PCL was done when the quantity of circulating plasma cells (PCs) using flow cytometry (FC) and morphology assessment was evaluated. The phenotypic profile of the PCs was also analysed to determine if there is an association with clinical outcomes and to evaluate the prognostic value of analysed markers. Our results revealed that FC is an excellent method for identifying circulating PCs as a significantly higher number was identified by FC than by morphology (26·7% vs. 13·5%, P = 0·02). None of secondary PCL cases expressed CD19 or CD20. A low level of expression with similar positivity of CD27, CD28, CD81 and CD117 was found in both PCL groups. A decrease of CD44 expression was detected only in secondary PCL. Expression of CD56 was present in more than half of PCL cases as well as cytoplasmic nestin. A decreased level of platelets, Eastern Cooperative Oncology Group score of 2-3 and lack of CD20+ PC were associated with a higher risk of death. FC could be incorporated in PCL diagnostics not only to determine the number of circulating PCs, but also to assess their phenotype profile and this information should be useful in patients' diagnosis and possible prognosis.

Dynamics of tumor-specific cfDNA in response to therapy in multiple myeloma patients.

OBJECTIVES: Progress in multiple myeloma treatment allows patients to achieve deeper responses, for which the assessment of minimal residual disease (MRD) is critical. Typically, bone marrow samples are used for this purpose; however, this approach is site-limited. Liquid biopsy represents a minimally invasive and more comprehensive technique that is not site-limited, but equally challenging. METHODS: While majority of current data comes from short-term studies, we present a long-term study on blood-based MRD monitoring using tumor-specific cell-free DNA detection by ASO-qPCR. One hundred and twelve patients were enrolled into the study, but long-term sampling and analysis were feasible only in 45 patients. RESULTS: We found a significant correlation of quantity of tumor-specific cell-free DNA levels with clinically meaningful events [induction therapy (P = .004); ASCT (P = .012)]. Moreover, length of cfDNA fragments is associated with better treatment response of patients. CONCLUSIONS: These results support the concept of tumor-specific cell-free DNA as a prognostic marker.

Multiplex Assay for Quantification of Acute Phase Proteins and Immunoglobulin A in Dried Blood Spots.

Inflammation is the first line defense mechanism against infection, tissue damage, or cancer development. However, inappropriate inflammatory response may also trigger diseases. The quantification of inflammatory proteins is essential to distinguish between harmful and beneficial immune response. Currently used immunoanalytical assays may suffer specificity issues due to antigen-antibody interaction and possible cross-reactivity of antibody with other protein species. In addition, immunoanalytical assays typically require invasive blood sampling and additional logistics; they are relatively costly and highly challenging to multiplex. We present a multiplex assay based on selected reaction monitoring (SRM) for quantification of seven acute-phase proteins (i.e., SAA1, SAA2-isoform1, SAA4, CRP, A1AT-isoform1, A1AG1, A1AG2) and the adaptive immunity effector IGHA1 in dried blood spots. This type of sample is readily available from all human subjects including newborns. The study utilizes proteotypic isotopically labeled peptides with trypsin-cleavable tag and presents optimized and reproducible workflow and several important practical remarks regarding quantitative SRM assays development. The panel of inflammatory proteins was quantified with sequence specificity capable to differentiate protein isoforms with intra- and interday precision (<16.4% coefficient of variation (CV) and <14.3% CV, respectively). Quantitative results were correlated with immuno-nephelometric assay (typically greater than 0.9 Pearson's R).

Circulating exosomal long noncoding RNA PRINS-First findings in monoclonal gammopathies.

Multiple myeloma is the second most common hematological malignancy characterized by focal lesions of malignant plasma cells in the bone marrow. These lesions contain subclones that directly influence survival of patients. Bone marrow biopsies are single-site biopsies and thus cannot contain all information about the tumor. In contrast, liquid biopsies analyze circulating cells and molecules that are secreted from all sites of the tumor. Long noncoding RNA molecules are one class of these molecules. We performed a two-phase biomarker study investigating lncRNA expression profiles in exosomes of peripheral blood serum of newly diagnosed multiple myeloma (MM) patients, monoclonal gammopathy of undetermined significance (MGUS) patients in comparison with healthy donors (HD). Surprisingly, this analysis revealed dysregulation of only one exosomal lncRNA PRINS in MM vs HD. Overall, MM and MGUS patients were distinguished from HD with sensitivity of 84.9% and specificity of 83.3%. Our study suggests a possible diagnostic role for exosomal lncRNA PRINS in monoclonal gammopathies patients.

Selected Genetic Polymorphisms Associated with Hypoxia and Multidrug Resistance in Monoclonal Gammopathies Patients.

BACKGROUND: Adaptive response to hypoxia is regulated by several mechanisms and transcription factors, including hypoxia-inducible factors (HIFs). Activation of HIF-1α is associated with increased expression of P-glycoprotein and multidrug resistance in cancer cells. In this retrospective study, we analyzed candidate single-nucleotide polymorphisms (SNPs) in HIF-1α and HIF-1ß associated with risk of monoclonal gammopathy of undetermined significance (MGUS) or multiple myeloma (MM). PATIENTS AND METHODS: Genotypes of SNPs associated with hypoxia were determined in an independent cohort of monoclonal gammopathies (MG) (275 MM and 228 MGUS patients) and in 219 cancer-free controls by real time polymerase chain reaction allelic discrimination. RESULTS: When MM patients were compared to controls, protective role of CG genotype compared to CC in HIF-1ß (rs2228099) for MM development was observed (OR = 0.65; CI 0.45-0.95; p = 0.026). Even after adjustment for patients' age and body mass index (BMI), there were significantly lower odds (OR = 0.55; p = 0.045) of developing MM patients of CG genotype in comparison to CC genotype. Log-rank test confirmed association of GT haplotype (rs11549467, rs2057482) in HIF-1α with better overall survival (median 41.8 months; (CI 35.1-48.5)) for "none GT" and median 93.8 months (CI 31.3-156.4) for "at least one GT" haplotype (p = 0.0500). Further, significant associations between SNPs in MDR1 and outcome of MM were found in 110 MM patients that underwent bortezomib-based treatment. CONCLUSION: Our study showed a genetic predisposition for risk of MG development and/or outcome of MM patients; nevertheless, further studies are needed to confirm our initial analysis.

The spectrum of somatic mutations in monoclonal gammopathy of undetermined significance indicates a less complex genomic landscape than that in multiple myeloma.

Monoclonal gammopathy of undetermined significance is a pre-malignant precursor of multiple myeloma with a 1% risk of progression per year. Although targeted analyses have shown the presence of specific genetic abnormalities such as IGH translocations, RB1 deletion, 1q gain, hyperdiploidy or RAS gene mutations, little is known about the molecular mechanism of malignant transformation. We performed whole exome sequencing together with comparative genomic hybridization plus single nucleotide polymorphism array analysis in 33 flow-cytometry-separated abnormal plasma cell samples from patients with monoclonal gammopathy of undetermined significance to describe somatic gene mutations and chromosome changes at the genome-wide level. Non-synonymous mutations and copy-number alterations were present in 97.0% and in 60.6% of cases, respectively. Importantly, the number of somatic mutations was significantly lower in monoclonal gammopathy of undetermined significance than in myeloma (P<10-4) and we identified six genes that were significantly mutated in myeloma (KRAS, NRAS, DIS3, HIST1H1E, EGR1 and LTB) within the monoclonal gammopathy of undetermined significance dataset. We also found a positive correlation with increasing chromosome changes and somatic gene mutations. IGH translocations, comprising t(4;14), t(11;14), t(14;16) and t(14;20), were present in 27.3% of cases and in a similar frequency to myeloma, consistent with the primary lesion hypothesis. MYC translocations and TP53 deletions or mutations were not detected in samples from patients with monoclonal gammopathy of undetermined significance, indicating that they may be drivers of progression to myeloma. Data from this study show that monoclonal gammopathy of undetermined significance is genetically similar to myeloma, however overall genetic abnormalities are present at significantly lower levels in monoclonal gammopathy of undetermined significant than in myeloma.

Activity of aldehyde dehydrogenase in B-cell and plasma cell subsets of monoclonal gammopathy patients and healthy donors.

BACKGROUND: Aldehyde dehydrogenase (ALDH) is highly active in physiological stem cells as well as in tumor-initiating cells of some malignancies including multiple myeloma (MM). Finding higher activity of ALDH in some cell subsets in monoclonal gammopathies (MG) could identify potential source of myeloma-initiating cells (MICs). METHODS: Bone marrow of 12 MM, 9 monoclonal gammopathy of undetermined significance (MGUS), and 10 healthy donors (HD) were analyzed by flow cytometry. ALDH activity of B-cells and plasma cells (PC) was analyzed using Aldefluor. RESULTS: Similar changes of ALDH activity were found during B-cell development in HD and MG. Decreasing of ALDH activity from immature to naïve B-cells was found. In postgerminal stages, the activity started to increase, and in PCs, the ALDH activity was the same as in immature B-cells. Increased ALDH activity of all PC subsets compared to naïve B-cells was found in MM as well as in HD, while in MGUS, only CD19- PCs have higher ALDH activity. In HD, ALDH activity was higher in CD19+ PCs compared with MG. CONCLUSIONS: Our results indicate that changes of ALDH activity are the natural phenomenon in B-cell development; thus, high ALDH activity as a single marker is not appropriate for MICs identification.

Deregulated expression of long non-coding RNA UCA1 in multiple myeloma.

OBJECTIVES: Long non-coding RNAs (lncRNAs) are RNA transcripts longer than 200 nucleotides that are not translated into proteins. They are involved in pathogenesis of many diseases including cancer and have a potential to serve as diagnostic and prognostic markers. We aimed to investigate lncRNA expression profiles in bone marrow plasma cells (BMPCs) of newly diagnosed multiple myeloma (MM) patients in comparison to normal BMPCs of healthy donors (HD) in a three-phase biomarker study. METHODS: Expression profile of 83 lncRNA was performed by RT2 lncRNA PCR Array (Qiagen), followed by quantitative real-time PCR using specific TaqMan non-coding RNA assays analyzing 84 newly diagnosed MM patients and 25 HD. RESULTS: Our analysis revealed dysregulation of two lncRNAs; NEAT1 (sensitivity of 55.0% and specificity of 79.0%) and UCA1 (sensitivity of 85.0% and specificity of 94.7%). UCA1 levels correlated with albumin and monoclonal immunoglobulin serum levels, cytogenetic aberrations, and survival of MM patients. CONCLUSION: Our study suggests a possible prognostic impact of UCA1 expression levels on MM patients.

Cytogenetics in multiple myeloma patients progressing into extramedullary disease.

BACKGROUND: Extramedullary disease in multiple myeloma patients is an uncommon event occurring either at the time of diagnosis, or during disease progression/relapse. This manifestation is frequently associated with poor outcome and resistance to treatment. We evaluated chromosomal alterations of plasma cells of multiple myeloma patients with extramedullary relapse, either in the bone marrow (BM) or at extramedullary sites, and in previous BM collection by interphase fluorescence in situ hybridization. MATERIAL AND METHODS: Thirty-one patients [25 BM plasma cells (BMPCs), 18 extramedullary tumor plasma cells], of which 12 had paired samples of BM and extramedullary plasma cells and 14 had previous collection of BM, were investigated for the presence of chromosomal aberrations (CHAs): del(17)(p13), del(13)(q14), 14q32 disruption, t(4;14)(p16;q32), t(14;16)(q32;q23), gain(1)(q21), and hyperdiploidy status. RESULTS: Overall, in unrelated samples, t(4;14) was more prevalent in extramedullary plasma cells, and hyperdiploidy was more frequent in BMPCs. In paired samples, there was a higher frequency of del(13)(q14) and 14q32 disruption in BMPCs. Frequency of all studied CHAs was higher in BMPCs of extramedullary patients than in their previous sample collection. CONCLUSION: These data show that plasma cells harbor more aberrations during their transformation into extramedullary form.

Genomewide profiling of copy-number alteration in monoclonal gammopathy of undetermined significance.

Monoclonal gammopathy of undetermined significance (MGUS) is a benign condition with an approximate 1% annual risk of symptomatic plasma cell disorder development, mostly to multiple myeloma (MM). We performed genomewide screening of copy-number alterations (CNAs) in 90 MGUS and 33 MM patients using high-density DNA microarrays. We identified CNAs in a smaller proportion of MGUS (65.6%) than in MM (100.0%, P = 1.31 × 10-5 ) and showed median number of CNAs is lower in MGUS (3, range 0-22) than in MM (13, range 4-38, P = 1.82 × 10-10 ). In the MGUS cohort, the most frequent losses were located at 1p (5.6%), 6q (6.7%), 13q (30.0%), 14q (14.4%), 16q (8.9%), 21q (5.6%), and gains at 1q (23.3%), 2p (6.7%), 6p (13.3%), and Xq (7.8%). Hyperdiploidy was detected in 38.9% of MGUS cases, and the most frequent whole chromosome gains were 3 (25.6%), 5 (23.3%), 9 (37.8%), 15 (23.3%), and 19 (32.2%). We also identified CNAs such as 1p, 6q, 8p, 12p, 13q, 16q losses, 1q gain and hypodiploidy, which are potentially associated with an adverse prognosis in MGUS. In summary, we showed that MGUS is similar to MM in that it is a genetically heterogeneous disorder, but overall cytogenetic instability is lower than in MM, which confirms that genetic abnormalities play important role in monoclonal gammopathies.

Improved Screening of Monoclonal Gammopathy Patients by MALDI-TOF Mass Spectrometry.

Monoclonal gammopathies are a group of blood diseases characterized by presence of abnormal immunoglobulins in peripheral blood and/or urine of patients. Multiple myeloma and plasma cell leukemia are monoclonal gammopathies with unclear etiology, caused by malignant transformation of bone marrow plasma cells. Mass spectrometry with matrix-assisted laser desorption/ionization and time-of-flight detection is commonly used for investigation of the peptidome and small proteome of blood plasma with high accuracy, robustness, and cost-effectivity. In addition, mass spectrometry coupled with advanced statistics can be used for molecular profiling, classification, and diagnosis of liquid biopsies and tissue specimens in various malignancies. Despite the fact there have been fully optimized protocols for mass spectrometry of normal blood plasma available for decades, in monoclonal gammopathy patients, the massive alterations of biophysical and biochemical parameters of peripheral blood plasma often limit the mass spectrometry measurements. In this paper, we present a new two-step extraction protocol and demonstrated the enhanced resolution and intensity (>50×) of mass spectra obtained from extracts of peripheral blood plasma from monoclonal gammopathy patients. When coupled with advanced statistics and machine learning, the mass spectra profiles enabled the direct identification, classification, and discrimination of multiple myeloma and plasma cell leukemia patients with high accuracy and precision. A model based on PLS-DA achieved the best performance with 71.5% accuracy (95% confidence interval, CI = 57.1-83.3%) when the 10× repeated 5-fold CV was performed. In summary, the two-step extraction protocol improved the analysis of monoclonal gammopathy peripheral blood plasma samples by mass spectrometry and provided a tool for addressing the complex molecular etiology of monoclonal gammopathies.

LncRNAs LY86-AS1 and VIM-AS1 Distinguish Plasma Cell Leukemia Patients from Multiple Myeloma Patients.

Long non-coding RNAs (lncRNAs) are functional RNAs longer than 200 nucleotides. Due to modern genomic techniques, the involvement of lncRNAs in tumorigenesis has been revealed; however, information concerning lncRNA interplay in multiple myeloma (MM) and plasma cell leukemia (PCL) is virtually absent. Herein, we aimed to identify the lncRNAs involved in MM to PCL progression. We investigated representative datasets of MM and PCL patients using next-generation sequencing. In total, 13 deregulated lncRNAs (p < 0.00025) were identified; four of them were chosen for further validation in an independent set of MM and PCL patients by RT-qPCR. The obtained results proved the significant downregulation of lymphocyte antigen antisense RNA 1 (LY86-AS1) and VIM antisense RNA 1 (VIM-AS1) in PCL compared to MM. Importantly, these two lncRNAs could be involved in the progression of MM into PCL; thus, they could serve as promising novel biomarkers of MM progression.

Levels of angiogenic factors in patients with multiple myeloma correlate with treatment response.

Angiogenesis plays a significant role in the pathogenesis of multiple myeloma (MM). We have measured concentrations of angiogenesis activators, including vascular endothelial growth factor (VEGF), basic fibroblast growth factor, and hepatocyte growth factor (HGF), and inhibitors, including endostatin, thrombospondin-1 (TSP-1), and angiostatin in the peripheral and bone marrow blood of MM patients at diagnosis and after high-dose chemotherapy. We have analyzed 96 patients with secretory MM. Serial measurements of angiogenesis factors/inhibitors were analyzed in the plasma by subgroups based on the best treatment response. Concentrations of angiogenic factors were determined in the peripheral blood and bone marrow plasma. There were significant decreases of VEGF and HGF levels and a significant increase in TSP-1 concentrations in the bone marrow plasma of patients who achieved complete or very good partial response in contrast to those who had partial or no response. VEGF and HGF levels decrease but those of TSP-1 increase after successful treatment for MM, indicating a reduction in the rate of angiogenesis.

Rapid discrimination of multiple myeloma patients by artificial neural networks coupled with mass spectrometry of peripheral blood plasma.

Multiple myeloma (MM) is a highly heterogeneous disease of malignant plasma cells. Diagnosis and monitoring of MM patients is based on bone marrow biopsies and detection of abnormal immunoglobulin in serum and/or urine. However, biopsies have a single-site bias; thus, new diagnostic tests and early detection strategies are needed. Matrix-Assisted Laser Desorption/Ionization Time-of Flight Mass Spectrometry (MALDI-TOF MS) is a powerful method that found its applications in clinical diagnostics. Artificial intelligence approaches, such as Artificial Neural Networks (ANNs), can handle non-linear data and provide prediction and classification of variables in multidimensional datasets. In this study, we used MALDI-TOF MS to acquire low mass profiles of peripheral blood plasma obtained from MM patients and healthy donors. Informative patterns in mass spectra served as inputs for ANN that specifically predicted MM samples with high sensitivity (100%), specificity (95%) and accuracy (98%). Thus, mass spectrometry coupled with ANN can provide a minimally invasive approach for MM diagnostics.

Urinary intermediates of tryptophan as indicators of the gut microbial metabolism.

While over 10% of the human metabolome is directly associated with the gut microbial metabolism, specific metabolites are largely uncharacterized. Therefore, methods for the identification and quantification of microbiota-associated metabolites in biological fluids such as urine or plasma are necessary in order to elucidate the molecular basis of host-microbiota interaction. In this study, we focused on the tryptophan metabolism, employing quantitative assays by ultra-high performance liquid chromatography (UHPLC) and tandem mass spectrometry, specifically selected reaction monitoring (SRM). Metabolite standards were utilized to generate SRM library for 16 intermediates of the tryptophan metabolism which were human endogenous as well as microbiota-associated based on the HMDB classification. Next, the SRM assays were utilized for screening in maternal urine samples and in dried urine specimens from neonates. The approach resulted in the discovery of microbiota-associated metabolites (methyl indole-3-acetate and methyl indol-3-propionate) previously unreported in urine samples and additionally in quantification of 8 intermediates of the tryptophan metabolism. To the best of our knowledge, this study represents the first attempt to explore previously unreported microbial metabolites in urine by UHPLC-SRM and novel methodology for simultaneous determination of microbiota-modulated component of Trp metabolism.

Extramedullary relapse of multiple myeloma defined as the highest risk group based on deregulated gene expression data.

BACKGROUND: Multiple myeloma (MM) is characterized by malignant proliferation of plasma cells (PC) which accumulate in the bone marrow (BM). The advent of new drugs has changed the course of the disease from incurable to treatable, but most patients eventually relapse. One group of MM patients (10-15%) is considered high-risk because they relapse within 24 months. Recently, extramedullary relapse of MM (EM) has been observed more frequently. Due to its aggressivity and shorter survival, EM is also considered high-risk. AIMS: The goal of this study was to determine if the so-called high-risk genes published by the University of Arkansas group (UAMS) are even more deregulated in EM patients than in high-risk MM patients and if these patients may be considered high-risk. METHODS: Nine samples of bone marrow plasma cells from MM patients as well as 9 tumors and 9 bone marrow plasma cells from EM patients were used. Quantitative real-time PCR was used for evaluation of expression of 15 genes connected to the high-risk signature of MM patients. RESULTS: Comparison of high-risk plasma cells vs extramedullary plasma cells revealed 4 significantly deregulated genes (CKS1B, CTBS, NADK, YWHAZ); moreover, comparison of extramedullary plasma cells vs extramedullary tumors revealed significant differences in 9 out of 15 genes. Of these, 6 showed significant changes as described by the UAMS group (ASPM, SLC19A1, NADK, TBRG4, TMPO and LARS2). CONCLUSIONS: Our data suggest that increasing genetic abnormalities as described by the gene expression data are associated with increased risk for EM relapse.

Circulating Serum MicroRNA-130a as a Novel Putative Marker of Extramedullary Myeloma.

Poor outcome of extramedullary disease in multiple myeloma patients and lack of outcome predictors prompt continued search for new markers of the disease. In this report, we show circulating microRNA distinguishing multiple myeloma patients with extramedullary disease from myeloma patients without such manifestation and from healthy donors. MicroRNA-130a was identified by TaqMan Low Density Arrays and verified by quantitative PCR on 144 serum samples (59 multiple myeloma, 55 myeloma with extramedullary disease, 30 healthy donors) in test and validation cohorts as being down-regulated in myeloma patients with extramedullary disease. Circulating microRNA-130a distinguished myeloma patients with extramedullary disease from healthy donors with specificity of 90.0% and sensitivity of 77.1%, patients with extramedullary disease from newly diagnosed multiple myeloma patients with specificity of 77.1% and sensitivity of 34.3% in the test cohort and with specificity of 91.7% and sensitivity of 30.0% in the validation cohort of patients. Circulating microRNA-130a in patients with extramedullary myeloma was associated with bone marrow plasma cells infiltration. Further, microRNA-130a was decreased in bone marrow plasma cells obtained from patients with extramedullary myeloma in comparison to bone marrow plasma cells of myeloma patients without such manifestation, but it was increased in tumor site plasma cells of patients with extramedullary disease compared to bone marrow plasma cells of such patients (p<0.0001). Together, our data suggest connection between lower level of microRNA-130a and extramedullary disease and prompt further work to evaluate this miRNA as a marker of extramedullary disease in multiple myeloma.