Abstract
Multiple myeloma (MM) is an incurable cancer of terminally differentiated plasma cells that proliferate in the bone marrow. miRNAs are promising biomarkers for risk stratification in MM and several miRNAs are shown to have a function in disease pathogenesis. However, to date, surprisingly few miRNA-mRNA interactions have been described for and functionally validated in MM. In this study, we performed miRNA-seq and mRNA-seq on CD138 + cells isolated from bone marrow aspirates of 86 MM patients to identify novel interactions between sRNAs and mRNAs. We detected 9.8% significantly correlated miRNA-mRNA pairs of which 5.17% were positively correlated and 4.65% were negatively correlated. We found that miRNA-mRNA pairs that were predicted by in silico target-prediction algorithms were more negatively correlated than non-target pairs, indicating functional miRNA targeting and that correlation between miRNAs and mRNAs from patients can be used to identify miRNA-targets. mRNAs for negatively correlated miRNA-mRNA target pairs were associated with gene ontology terms such as autophagy, protein degradation and endoplasmic stress response, reflecting important processes in MM. Targets for two specific miRNAs, miR-125b-5p and miR-365b-3p, were functionally validated in MM cell line transfection experiments followed by RNA-sequencing and qPCR. In summary, we identified functional miRNA-mRNA target pairs by correlating miRNA and mRNA data from primary MM cells. We identified several target pairs that are of potential interest for further studies. The data presented here may serve as a hypothesis-generating knowledge base for other researchers in the miRNA/MM field. We also provide an interactive web application that can be used to exploit the miRNA-target interactions as well as clinical parameters associated to these target-pairs.