MicroRNA expression in multiple myeloma is associated with genetic subtype, isotype and survival
© Chi et al; licensee BioMed Central Ltd. 2011
Received: 20 January 2011
Accepted: 18 May 2011
Published: 18 May 2011
MicroRNAs are small RNA species that regulate gene expression post-transcriptionally and are aberrantly expressed in many cancers including hematological malignancies. However, the role of microRNAs in the pathogenesis of multiple myeloma (MM) is only poorly understood. We therefore used microarray analysis to elucidate the complete miRNome (miRBase version 13.0) of purified tumor (CD138+) cells from 33 patients with MM, 5 patients with monoclonal gammopathy of undetermined significance (MGUS) and 9 controls.
Unsupervised cluster analysis revealed that MM and MGUS samples have a distinct microRNA expression profile from control CD138+ cells. The majority of microRNAs aberrantly expressed in MM (109/129) were up-regulated. A comparison of these microRNAs with those aberrantly expressed in other B-cell and T-cell malignancies revealed a surprising degree of similarity (~40%) suggesting the existence of a common lymphoma microRNA signature. We identified 39 microRNAs associated with the pre-malignant condition MGUS. Twenty-three (59%) of these were also aberrantly expressed in MM suggesting common microRNA expression events in MM progression. MM is characterized by multiple chromosomal abnormalities of varying prognostic significance. We identified specific microRNA signatures associated with the most common IgH translocations (t(4;14) and t(11;14)) and del(13q). Expression levels of these microRNAs were distinct between the genetic subtypes (by cluster analysis) and correctly predicted these abnormalities in > 85% of cases using the support vector machine algorithm. Additionally, we identified microRNAs associated with light chain only myeloma, as well as IgG and IgA-type MM. Finally, we identified 32 microRNAs associated with event-free survival (EFS) in MM, ten of which were significant by univariate (logrank) survival analysis.
In summary, this work has identified aberrantly expressed microRNAs associated with the diagnosis, pathogenesis and prognosis of MM, data which will prove an invaluable resource for understanding the role of microRNAs in this devastating disease.
This article was reviewed by Prof. Neil Smalheiser, Prof. Yuriy Gusev, and an unknown reviewer.
Multiple myeloma (MM) is a plasma cell (PC) malignancy with an annual incidence of over 14,000 cases in the US alone. MM is essentially an incurable disease with a median survival of ~3 years that accounts for nearly 2% of deaths from cancer and about 20% of deaths from hematological cancers . Newer therapies, however, are resulting in improvements in the median survival . Recent advances in molecular and genetic research have lead to the realization that MM, although defined histologically as a single entity, encompasses a wide range and frequently complex mixture of genomic abnormalities which differ in both their molecular pathogenesis and prognostic significance. The recent discovery of short non-coding RNA molecules that regulate gene expression post-transcriptionally, known as microRNAs, represent yet another level of complexity in our understanding of gene regulation and as such could further our understanding of the pathogenesis of MM.
MicroRNAs have been demonstrated to have diagnostic and prognostic potential in cancer [3–7] and it has been suggested that microRNA expression profiling can distinguish cancers according to both the cellular nature and the developmental stage of the tumor with a greater degree of accuracy than traditional gene expression analysis . There is increasingly strong evidence that microRNAs are involved in the pathogenesis of many cancers including B and T cell lymphomas [9, 10]. There is however little known about the role that microRNAs play in MM.
Therefore we undertook a comprehensive study using microarray technology to elucidate the complete miRNome (miRBase version 13.0) of purified tumor (CD138+) cells from the bone marrow of 33 MM and 5 MGUS patients (and 9 controls). In order to investigate microRNA expression in different genetic subtypes MM cases were classified cytogenetically by FISH. These data were then correlated with genetic subtype and clinical parameters.
Most aberrantly expressed microRNAs associated with MM are up-regulated
MicroRNAs differentially expressed (P < 0.05) between MM (n = 33) and controls (n = 9) depicting chromosomal location and where relevant their expression within a microRNA cluster.
adj. P value
Expression levels of the 129 MM-associated microRNAs were distinct between MM and control samples (Additional file 1 Figure S1). Indeed using only expression levels of the ten most discriminatory (by P-value) up- and down-regulated microRNAs distinguished samples by cluster analysis (Figure 1B) and correctly predicted samples as either MM or control for 36/41 (88%) samples using the leave-one-out cross-validation support vector machine (SVM) algorithm.
The microRNA profile of MGUS patient samples is distinct from MM and controls
MGUS is generally believed to represent a pre-malignant form of MM . To examine microRNA expression in this condition we initially compared the profile of PC cells from MGUS patients with that of controls. Thirty-nine microRNAs were identified as being aberrantly expressed in MGUS samples (P < 0.05; Additional file 1 Table S3); 28 microRNAs were up-regulated and 11 down-regulated. Nine of these were consistent with a previous study , including up-regulation of miR-21, the most up-regulated microRNA identified (10.3 fold change; Figure 2A), miR-222 and miR-342 (Figure 2E, H respectively), and down-regulation of miR-200b (Figure 2F). Expression levels of MGUS-associated microRNAs distinguished between MGUS and controls by cluster analysis (Additional file 1 Figure S2), and correctly predicted 12/13 (92%) the presence of this disorder by SVM analysis.
Different genetic subtypes of MM have distinct microRNA expression profiles
Recently it has become increasingly clear that MM displays an enormous genomic complexity that underlies the clinical heterogeneity observed with this disease. The presence of recurrent chromosomal abnormalities such as IgH translocations and del(13q) have been shown to have prognostic significance in multiple studies and as a consequence cytogenetic classification is increasingly being used to stratify MM patients in clinical practice although not yet routinely [15, 16]. In order to investigate microRNA expression in the different genetic subtypes of MM, twenty-six cases (for which there was sufficient material) were cytogenetically classified using the micro-FISH technique . Individual patient karyotype data are summarized in Additional file 1 Table S1. Fourteen of 26 (54%) of MM cases harbored an IgH translocation and 7/26 (27%) exhibited the 13q (RB) deletion. Of the translocation group eight (57%) displayed the t(11;14) translocation, three (21%) cases the t(4;14) translocation, and three cases had an IgH translocation with an un-known (tested) fusion partner. Ten (38%) MM patients had no detectable chromosomal abnormality using our FISH probe sets, although it is probable that these cases had undetected abnormalities.
MicroRNAs associated (P < 0.05) with different genetic subtypes depicting chromosomal location.
IgH trans (n = 14)
t(11;14) (n = 8)
t(4;14) (n = 3)
13q(del) (n = 6)
A comparison between microRNA expression in del(13q) MM patients with patients not harboring this deletion identified 28 microRNAs by ANOVA (P < 0.05; Table 2). Expression levels of these microRNAs distinguished cases that were del(13q) by cluster analysis (Figure 3D), and correctly predicted its presence in 22/26 (85%) of cases by SVM analysis. Eight of these microRNAs were consistent with research by Gutiérrez and colleagues , including down-regulation of miR-15a and miR-16, encoded within the 13q14 locus (Figure 2K and 2L respectively).
MicroRNA expression is associated with clinical parameters and survival
The most common paraprotein isotypes associated with MM are IgG (60%) and IgA (24%) . It has been reported that IgA myeloma has a poorer prognostic outcome than patients with IgG myeloma , although it should be noted that we (data not shown) and others  did not find any survival differences between isotypes. Consequently we compared expression levels between thirteen cases of IgG-type MM with eight cases of IgA-type. This resulted in the identification of twenty-one differentially expressed microRNAs (P < 0.05; Additional file 1 Table S5). Expression levels of these microRNAs were distinct between IgA and IgG MM cases (Additional file 1 Figure S4), and correctly predicted isotype for 20/21 (95%) of cases by SVM analysis.
Approximately 10% of MM patients secrete only light chains instead of the complete immunoglobulin molecule . So called "light chain only myeloma" or "Bence Jones myeloma" had been reported to have a poorer prognostic outcome [22, 24] although a more recent (and comprehensive) study found no such difference . In our cohort eight MM cases had light chain only myeloma. Twenty-seven microRNAs were identified as being differentially expressed compared to non-light chain only MM cases (Additional file 1 Table S6). Expression levels of these microRNAs correctly predicted light-chain only myeloma in all cases.
In this study we used microarray technology to elucidate the complete miRNome (miRBase version 13.0) of purified PCs from 33 MM and 5 MGUS patients and compared their expression with counterpart normal PCs. One hundred and twenty-nine microRNAs were identified as being aberrantly expressed in MM. Thirty-nine (30%) of these were previously reported to be aberrantly expressed in MM by Pichiorri et al  and/or Zhou et al . The reason for these discrepancies probably reflects the use of differing statistical, microarray and sampling methods. Of note the former study only examined 10 MM and 4 controls, whilst the latter study included only 2 controls for analysis. Consistent with the findings of Pichiorri and colleagues we found that miR-21 was up-regulated in both MM and MGUS as were all seven microRNAs encoded by miR-17-92 cluster (average fold-change = 2.1; range 1.75-2.78). However unlike that study, although miR-32 and miR-181a were identified as being up-regulated by microarray analysis, we failed to find any significant difference in expression levels by qRT-PCR (data not shown). Surprisingly 10 of the 20 down-regulated microRNAs identified in the current study are encoded by the miR-371-1323 cluster located at 19q13.42. This cluster consists of 49 microRNAs, 39 of which were present on our microarray. All members of this cluster were down-regulated (average 1.6-fold; range 1.1-2.5). Interestingly, members of this cluster (miR-373 and miR-520c) have been demonstrated to be associated with tumor invasion and metastasis .
The vast majority (109/129 (85%)) of microRNA aberrantly expressed in MM were up-regulated, consistent with the findings of Zhou et al . We [7, 26] (and others) found the same pattern in other B-cell malignancies. In contrast, we  (and others ) recently demonstrated that T-cell malignancies are associated with a global decrease in microRNA expression. This suggests a fundamental difference in microRNA biogenesis between B and T-cell malignancies. A possible explanation for this distinction comes from our previous observation that components of the microRNA biosynthetic pathway (i.e. the microprocessor complex (DGCR8 and Drosha)) are up-regulated in B-cell malignancies (including MM) but down-regulated in T-cell malignancies .
MGUS has been shown to consistently proceed the development of MM  and as such is considered to represent a pre-malignant model for MM pathogenesis. Twenty-three of the 39 (59%) microRNAs aberrantly expressed in MGUS were also aberrantly expressed in MM (Additional file 1 Table S3) suggesting common microRNA expression events in MM progression. Thirty-seven microRNAs were identified as being differentially expressed between MGUS and MM (Additional file 1 Table S4), only 9 of these were also differentially expressed in MGUS compared to controls (i.e. specific for MGUS) including miR-21 as previously identified . The relationship between microRNAs dysfunctionally expressed in MGUS and MM are shown in Additional file 1 Figure S3. miR-21 over-expression has been linked to an anti-apoptotic phenotype and enhanced tumor growth [38, 39] and as such may represent an important aberration in the early pathogenesis of MM.
Next we investigated the identity of microRNAs associated with the major genetic subtypes of MM. Two previous studies by Lionetti et al  and Gutiérrez et al  have also examined microRNA expression associated with recurrent chromosomal abnormalities in MM. Generally our findings were consistent with these studies although microRNAs associated with specific abnormalities did not always completely overlap. The cause for discrepancies between studies, including lack of concordance between the two previously published studies, again probably represents differences between microarray platforms, statistical techniques and samples used for analysis. For example the study of Gutiérrez et al considered cases in which del(13q) was the sole chromosomal abnormality and deregulated microRNAs were identified on the basis of their expression relative to control PCs, whereas our study and that of Lionetti et al considered the difference between MM cases with and without specific chromosomal abnormalities [19, 20]. The latter study considered genetic abnormalities in the context of a gene-expression defined molecular translocation/cyclin (TC) classification system . Such differences between studies are a reflection of the relative immaturity of the microRNA field that has only formally existed for the last 9 years, and is something that will surely be resolved in future years.
Interestingly, from the 28 microRNAs associated with the presence of del(13q), all seven microRNAs encoded by the miR-17-92 cluster (encoded at 13q31) were down-regulated (average 3.1 fold; range 1.2-4.1). Recent data suggests that unlike the characteristic del(13q) karyotype observed in chronic lymphocytic leukemia (CLL) defined by loss of the 13q14 locus, the minimal deleted region of 13q loss in MM extends from 13q14 to the 13q34 locus , consistent with down-regulation of the miR-17-92 cluster observed in these cases. Although down-regulation of the miR-17-92 cluster was also reported by Gutiérrez et al , this finding is counter-intuitive as members of this cluster are widely believed to function as oncogenes  whilst the presence of del(13q) is associated with an adverse prognostic outcome in MM . A possible explanation for this apparent paradox is that although widely up-regulated in a wide range of cancers [9, 26, 43], the miR-17-92 cluster has been proposed to act as both tumor-suppressor or oncogene depending upon the cellular context [44–46]. Recently we reported that this cluster is down-regulated in T-cell lymphoma where it was shown to have tumor-suppressor properties . Consequently, it is possible that the cluster is similarly acting as a tumor-suppressor in del(13q) MM cases although it should be noted that up-regulation of this cluster is associated with MM diagnosis (Table 1) and at least one member of this cluster (miR-19) was found to antagonize MM tumor development in a murine model .
We also report for the first time a correlation between microRNA expression and isotype in MM including the presence of light-chain only myeloma. Such information might be useful in providing an alternative clinical classification of new MM patients. Finally, we examined how the expression levels of individual microRNAs might influence clinical outcome of MM patients by correlating expression values with adverse effects using Kaplan-Meier survival analysis. We found that high expression of miR-153, miR-490, miR-455, miR-642, miR-500, and miR-296 were associated with better EFS in our cohort of patients whilst the opposite was true for miR-548d, miR-373, miR-554 and miR-888. Interestingly, the most discriminatory of these microRNAs (i.e. miR-153) has been shown to target BCL2 and MCL1  and that over-expression of both of these are associated with poor prognostic outcome in MM . Low expression of miR-455 has also been associated poor overall survival in endometrial serous adenocarcinoma , and low levels of miR-296 shown to be associated with increased prostate tumor growth and invasion properties via targeting of HMGA . It should be noted however that these data need to be tested in an independent test cohort before any firm conclusions can be drawn about the prognostic ability of these microRNAs in MM patient survival, whilst the current results provide the experimental basis for future validation experiments.
In summary, these results indicate that aberrant expression of microRNAs is a common feature of MM as well as pre-malignant MGUS, and that this phenomenon is associated with genetic and clinical subtype as well as clinical outcome. This data reinforces the complex nature of MM and hopefully will provide useful information in clarifying the molecular basis of this disease in future studies.
Patient samples and cell lines
Bone marrow samples were obtained from patients attending the Department of Haematology, John Radcliffe Hospital, Oxford, UK or Dubrava University hospital Zagreb, Croatia. Patients were diagnosed as MM (n = 33) or MGUS (n = 5) according to criteria described in the myeloma management guidelines . The median age of patients was 67 yrs (range 43-89 yrs). Individual patient details are shown in Additional file 1 Table S1. Plasma cells (PC) were purified from bone marrow aspirates by positive immunomagnetic selection (CD138+) according to manufacturer's instructions (Miltenyi Biotec, Bisley, UK). PCs obtained from the bone marrow of nine healthy individuals were used as controls. The purity of PC samples was > 90% CD138+ as measured by immunohistochemical staining.
Patient samples were collected in accordance with the Declaration of Helsinki and approved by the local ethics committee (Oxforshire Regional Ethics Committee; Ref no. 07/Q1606/25). Patients gave written informed consent for the sample collection.
MM cell lines, NCI-H929, JJN3, RPMI8226 and Thiel were routinely cultured in RPMI 1640 containing 10% fetal calf serum (Invitrogen, Paisley, UK). All cell lines were obtained from the DSMZ cell collection (Braunschweig, Germany) except Thiel which was a kind gift from Prof. Diehl (University of Cologne, Germany).
Cytogenetic classification of myeloma samples
Purified MM patient PCs (~105 cells) were classified according to cytogenetic criteria using the micro-FISH technique as previously described . The following probe sets obtained from Vysis (Abbott Diagnostics, Maidenhead, UK) were used for classification; D13S319 Spectrum-Orange/LSI-13q34-Spectrum-Green (del(13q) probe), IgH Spectrum-Green/FGFR3 Spectrum-Orange (t(4;14) fusion probe), IgH Spectrum-Green/CCND1 Spectrum-Orange (t(11;14) fusion probe), and IgH dual-color (Spectrum-Green/Spectrum-Orange) break-apart probe . At least 100 nuclei were scored for each probe. A 20% cut-off was used for numerical abnormalities and 10% for fusions and break-apart probes as recommended by the European Myeloma Network FISH workshop guidelines http://www.myeloma-europe.org.
RNA purification and microarray analysis
MicroRNA was isolated from samples (~106 cells) using Trizol (Invitrogen, Paisley, UK) and RNeasy columns as described by the manufacturer (Qiagen, Crawley, UK). MicroRNA (~500 ng) were labeled and hybridized to μRNA microarrays as previously described  using tonsillar material (pooled from twelve healthy individuals) as a common reference in a dye-balanced design. The arrays contained 655 human probes (miRBase v.10.1). Probe details can be found at http://www.microRNAworld.com.
Image analysis was carried out with BlueFuse software (BlueGnome, Cambridge, UK). Raw image data were global loess-normalized within arrays and normalized between arrays using the LIMMA package . The normalized log ratios (average of four replicates per probe) were used for subsequent analysis in Genespring 7.2 (Agilent Technologies, CA, US). ANOVA analysis was used to identify microRNAs differentially expressed between sample types and P values were adjusted using the Benjamini-Hogberg correction method. Differentially expressed genes were tested for their ability to predict sample class using the leave-one-out cross-validation support vector machine (SVM) function in Genespring (Polynomial Dot Product (Order 1) Kernel Function. Diagonal Scaling Factor: 0). All microarray data was MIAME compliant and raw data has been deposited in the GEO database (Accession series # GSE243371).
Quantitative RT-PCR (qRT-PCR)
MicroRNA qRT-PCR was carried out using Taqman probes as described by the manufacturer (Applied Biosystems, Warrington, UK) using 20 ng of microRNA per reaction in a Roche LightCycler 480 machine. Triplicate samples were used throughout. The mean Ct value of each triplicate was used for analysis, by the ΔCt method (ΔCt = mean Ct of microRNA of interest-mean Ct of U6). Expression levels were compared using Mann-Whitney independent t-test (Graphpad Prism v.4.0, La Jolla, CA).
Kaplan-Meier survival analysis was carried out on event free survival (EFS) times of MM patients as a function of microRNA expression, using the median value as cutoff. EFS was calculated as the time of diagnosis to the date of clinical relapse, death or last contact. Patients who were relapse-free at time of last contact were censored for analysis. Mean follow-up time was 20 months (range 1-60 months). Curves were compared by univariate (logrank) analysis using GraphPad Prism version 4.00 (La Jolla, CA).
Reviewer #1: Prof. Neil Smalheiser, Department of Psychiatry, University of Illinois, Chicago, IL, US.
This is a nice paper and I only have relatively minor comments, which can be addressed by page number:
p. 2. It would not hurt to have some discussion of the miR expression in the control group, in terms of the cell type profile, heterogeneity across individuals, etc. If you included that, I missed it.
Whilst the reviewer has an interesting suggestion i.e. to look at variation in miRNA expression across the control group, the main theme of this research was to look at expression in multiple myeloma, therefore the number of controls samples used (n = 9) does not have sufficient power in order to address that specific question. However, our data does suggest that the control samples varied little in their miRNA expression profile as shown by unsupervised cluster analysis (Figure 1).
p. 2. You speak of altered miR expression in tumor cells as "aberrant" throughout the paper. Yet not every difference is necessarily aberrant. Some changes may be compensatory; some changes may simply reflect the physiological state of the tumor cells without being linked to tumorigenicity per se (e.g., they may differ in O2 tension or pH).
We whole-heartedly agree that not all differentially expressed miRNAs identified in MM tumor cells will be directly linked with their neoplastic nature but instead may reflect indirect effects such as their abnormal microenvironment and/or interaction with non-tumor cells for example. We believe however that as these effects, direct and indirect, cumulatively represent tumor-associated changes in the miRNA expression profile of cells, that the use of the adjective 'aberrant' in this context is an accurate description of the expression behavior of miRNAs that are 'diverging from the normal' (Oxford Concise Dictionary).
p. 10 and Methods. It is worthwhile to describe the method of correcting p-values for multiple-testing in more detail.
As stated in the methods section p-values were corrected using the Benjamini-Hogberg multiple testing correction algorithm as implemented in the ANOVA function of the Genespring software package. This algorithm is based on false discovery rate error type and further details can be found in the original publication of this method (Benjamini and Hochberg (1995). Journal of the Royal Statistical Society B, 57, 289-300)). As this method is widely used in microarray analyses we believe that a more in-depth and mathematical description of this method is not warranted for this manuscript.
p. 10 and throughout. The paper has a nice balance between considering miRs as a clinical tool and as an indication of tumor biology. However, you can go still further to analyze in biological terms WHY the pattern of altered miR expression is the way it is, or (stated a different way) to analyze what dimensions are most important for SVM prediction. For example, do the altered miRs all share the same TFs driving them?
The reviewer has raised a very important point about addressing the reason for aberrant miRNA expression in MM. However, the causal mechanisms behind miRNA dysregulation in any cancer are complex, although we did try to proffer some insight in the context of MM in both results and discussion sections of this manuscript. Possible mechanisms include chromosomal lesions at miRNA-encoding regions and defects in miRNA biosynthetic machinery, both of which were addressed in this manuscript. Further mechanisms include epigenetic regulation, changes in expression of genes such as E2F that can bind to promoter regions of pri-miRNA sequences, as well as any factor that can regulate levels of the protein-encoding genes that encode for >40% of human miRNAs. The problem with identifying genes such as TFs that are potentially targeted by commonly differentially expressed miRNAs, is that currently the predictive algorithms to carry out this type of analysis perform very poorly (Sethupathy, et al (2006) Nat Methods, 3 , 881-886). Therefore we believe that the robustness and hence usefulness of such an in silico approach is limited. Consequently we decided to restrict our discussion to those target genes that have been independently experimentally validated.
p. 13. You mention that dicer and drosha tend to be up-regulated in B-cell malignancies. It would be desirable to have measured these in the SAME samples reported here and see if dicer and/or drosha levels are tightly correlated with miR up-regulation.
This would be the ideal situation, and in the case of the DLBCL case this comparison was carried out (Lawrie et al (2009) Br. J. Haematol., 145 , 545-548). Unfortunately there was insufficient material left from the current study to carry out a similar analysis on the MM samples.
p. 19, Methods. Were control and tumor samples collected during the same period? stored for similar amounts of time? Assayed in parallel? If not, this could create confounds.
Due to the logistics of obtaining bone marrow samples from ~50 patients/controls it would have been extremely difficult to collect samples over the same period- samples were collected over a ~2 year period from two centers (Oxford and Zagreb). Patient/control samples were assayed in parallel and randomized/anonymized to the researchers processing them in order to minimize confounding variables.
p. 20, Methods. Was U6 expression the same in control and tumor samples? This should be stated explicitly.
U6 levels were measured in controls and tumor in parallel using the same assay system (i.e. Taqman probe qRT-PCR). There was no significant difference in absolute U6 expression levels (i.e. Ct values) between controls and tumor samples.
p. 29, Table 1. You have mir-200a listed twice with two different values?! Please check this (very long) table.
This was a mistake. Position 25 in Table 1should have read miR-200a* and not miR-200a, this has been corrected.
p. 35, Figure 3 Legend. Please spell out DL, MM and SZ
This issue has been addressed.
Reviewer #2: Prof. Yuriy Gusev, Department of Surgery, University of Oklahoma, OK, USA.
This reviewer provided no comments for publication.
Reviewer #3: Unknown reviewer
This manuscript includes a description of microarray analysis for the expression of microRNA in purified CD138+ cells from 33 patients with multiple myeloma (MM), 5 patients with monoclonal gammopathy of undetermined significance (MGUS), and 9 healthy volunteers. The authors found distinct microRNA expression profiles from each group. Similar studies from 2 different groups have been reported. However, in one study only 16 samples from patients with MM, 6 from patients with MGUS, and 6 controls were included and in the another study samples from 52 patients but only 2 controls were studied. The manuscript is clear and concise. The data are convincing and the area of investigation is important. This manuscript is a valuable contribution. Specific comments follow:
1. The authors state several times that they have performed a comprehensive study to elucidate the complete miRNome of the CD138+ cells; however, I could not find how "complete" was defined. The manuscript would be improved by including a consideration about the completeness of the analysis performed. How the microarray analysis was defined as a complete analysis of miRNA and citations would enhance the manuscript.
The 'complete miRNome' refers to the microarrays (μRNA microarrays) used for these experiments which contain probes against every known miRNA (the miRNome) as defined by the reference depository for miRNAs (miRBase- http://www.mirbase.org) (Griffiths-Jones et al (2006), Nucleic Acids Research, 34 , D140-144). At the time when these experiments were started (June 2009), the miRBase version was 13.0 which included 655 human sequences. The authors readily accept that the miRNA database has grown since this time and that the current version of the miRBase miRNome is version 16 and contains 1048 human sequences. This assertion has been modified in the text to read 'complete miRNome (miRBase version 13.0)' in the text in order to avoid confusion.
2. The "leave-one-out cross-validation support vector machine algorithm" (mentioned on page 7) is not defined or described in the methods section. The support vector machine algorithm is also mentioned in the abstract but the manuscript would be improved by describing and referencing this technique.
The support vector machine (SVM) algorithm is a supervised learning method commonly used for classification in microarray experiments (along with k-nearest neighbor algorithm) and is based on non-probabilistic binary linear classification (Cortes & Vapnik (1995) Machine learning, 20 ). We believe that as the parameters used for its implementation within the Genespring package as described in the materials and methods section are already given, that this would be sufficient information for others to carry out the same analysis.
3. In the results the authors state that 8 miRNA were chosen to be measured by qRT-PCR to validate the microarray data. They state, "These data were consistent with the microarray results". In the discussion the authors indicate that "miR-32 and miR-181a were identified as being up-regulated by microarray analysis" by qRT-PCR failed to show any significant differences in expression levels. Were there other miRNA studied by qRT-PCR? Why weren't miR-32 and miR181a included in Figure 2?
The eight miRNAs that were tested by qRT-PCR shown in Figure 2were chosen on the basis that they were differentially expressed in our microarray analysis and had previously been shown to be dysregulated in MM or other hematological malignancies. As all of these miRNAs were also significantly differentially expressed when measured by qRT-PCR they were described as being 'consistent with the microarray data'. The only other miRNAs that we measured (miR-32 and miR-181a) by qRT-PCR which were also chosen on the basis shown to be dysregulated in previous studies but we did not find them to be significantly different by qRT-PCR which is why the data was not included in Figure 2.
4. In the discussion the authors state that the investigators of reference 17 "only examined 10 MM and 4 controls". In the abstract of this manuscript, Pichiorri et al. indicate that 16 samples from patients with MM were studied along with 6 with MGUS and 6 normal donors. Also, in this same spot of the discussion, the authors indicate that the investigators of reference 16 only studied 2 controls but the 52 patient samples analyzed was not mentioned. Moreover, the purpose of that study was to risk stratification.
These comments in the discussion were made in the context of the identification of aberrantly expressed miRNA associated with MM diagnosis. The Pichiorri study states in the results section (pg 12886) that their MM-associated miRNA signature was derived from the array data of 41 MM cell lines, 10 MM patients and 4 controls. Whereas for the Zhou et al study (pg 7904) the MM miRNA signature was derived from a comparison of array data from 52 MM cases with samples from 2 healthy donors. Whilst it is true that the Zhou et al study mainly focused on the identification of risk groups within the 52 MM cases, a significant part of the results section of this paper refers to the identity of miRNAs that are differentially expressed between controls and MM cases. The point that we were making in the discussion is that insufficient numbers of controls were used in both of these previous studies, and that our study was therefore an improvement on these studies, in the context of the identification of aberrantly expressed miRNAs in MM tumor cells.
Acknowledgements and Funding
This work was funded by grants from Leukaemia and Lymphoma Research (JC, EB, X-HC, DT, JB and JSW) and the Julian Starmer-Smith Memorial Fund (CHL). The authors acknowledge financial support from the Department of Health via the National Institute for Health Research (NIHR) comprehensive Biomedical Research Centre award to Oxford Radcliffe NHS Trust.
The authors wish to thank Maite Cabes, Graham Collins, Narendra Kaushik, Andy Campbell and Pamela Roberts (John Radcliffe Hospital) for their help in collecting patient material and collating clinical data. We are indebted to Fiona Ross and her colleagues at the LLR Myeloma Cytogenetic Database, Wessex Regional Genetic Laboratory, Salisbury, UK for their help in implementing and interpreting the microFISH technique as well as their input in the writing this manuscript.
List of abbreviations
analysis of variance
- BCL2 :
- B-cell CLL/lymphoma 2
- CCND1 :
- cyclin D1
- DGCR8 :
- DiGeorge syndrome critical region gene 8
- FGFR3 :
- fibroblast growth factor receptor 3
- HGMA1 :
- Homogentisate 1,2-dioxygenase
immunoglobulin (heavy chain)
- LMO2 :
- LIM domain only 2 (rhombotin-like 1)
- MCL1 :
- myeloid cell leukemia sequence 1 (BCL2-related)
monoclonal gammopathy of undetermined significance
Minimum Information About a Microarray Experiment
quantitative reverse transcriptase polymerase chain reaction
- RB1 :
- retinoblastoma 1
support vector machine
- Ries L, Eisner M, Kosary C: SEER Cancer Statistics Review, 1975-2001. 2004, Bethesda: National Cancer InstituteGoogle Scholar
- Mateos MV, Richardson PG, Schlag R, Khuageva NK, Dimopoulos MA, Shpilberg O, Kropff M, Spicka I, Petrucci MT, Palumbo A, et al: Bortezomib plus melphalan and prednisone compared with melphalan and prednisone in previously untreated multiple myeloma: updated follow-up and impact of subsequent therapy in the phase III VISTA trial. J Clin Oncol. 2010, 28: 2259-2266. 10.1200/JCO.2009.26.0638.PubMedView ArticleGoogle Scholar
- Calin GA, Ferracin M, Cimmino A, Di Leva G, Shimizu M, Wojcik SE, Iorio MV, Visone R, Sever NI, Fabbri M, et al: A MicroRNA signature associated with prognosis and progression in chronic lymphocytic leukemia. N Engl J Med. 2005, 353: 1793-1801. 10.1056/NEJMoa050995.PubMedView ArticleGoogle Scholar
- Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M, et al: MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005, 65: 7065-7070. 10.1158/0008-5472.CAN-05-1783.PubMedView ArticleGoogle Scholar
- Lawrie CH, Soneji S, Marafioti T, Cooper CD, Palazzo S, Paterson JC, Cattan H, Enver T, Mager R, Boultwood J, et al: MicroRNA expression distinguishes between germinal center B cell-like and activated B cell-like subtypes of diffuse large B cell lymphoma. Int J Cancer. 2007, 121: 1156-1161. 10.1002/ijc.22800.PubMedView ArticleGoogle Scholar
- Lawrie CH, Gal S, Dunlop HM, Pushkaran B, Liggins AP, Pulford K, Banham AH, Pezzella F, Boultwood J, Wainscoat JS, et al: Detection of elevated levels of tumour-associated microRNAs in serum of patients with diffuse large B-cell lymphoma. Br J Haematol. 2008, 141: 672-675. 10.1111/j.1365-2141.2008.07077.x.PubMedView ArticleGoogle Scholar
- Lawrie CH, Chi J, Taylor S, Tramonti D, Ballabio E, Palazzo S, Saunders NJ, Pezzella F, Boultwood J, Wainscoat JS, Hatton CS: Expression of microRNAs in diffuse large B cell lymphoma is associated with immunophenotype, survival and transformation from follicular lymphoma. J Cell Mol Med. 2009, 13: 1248-1260. 10.1111/j.1582-4934.2008.00628.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838. 10.1038/nature03702.PubMedView ArticleGoogle Scholar
- Lawrie CH: microRNA expression in lymphoid malignancies: new hope for diagnosis and therapy?. J Cell Mol Med. 2008, 12: 1432-1444. 10.1111/j.1582-4934.2008.00399.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Ballabio E, Mitchell T, van Kester MS, Taylor S, Dunlop HM, Chi J, Tosi I, Vermeer MH, Tramonti D, Saunders NJ, et al: MicroRNA expression in Sezary syndrome: identification, function, and diagnostic potential. Blood. 2010, 116: 1105-1113. 10.1182/blood-2009-12-256719.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhou Y, Chen L, Barlogie B, Stephens O, Wu X, Williams DR, Cartron MA, van Rhee F, Nair B, Waheed S, et al: High-risk myeloma is associated with global elevation of miRNAs and overexpression of EIF2C2/AGO2. Proc Natl Acad Sci USA. 107: 7904-7909.
- Pichiorri F, Suh SS, Ladetto M, Kuehl M, Palumbo T, Drandi D, Taccioli C, Zanesi N, Alder H, Hagan JP, et al: MicroRNAs regulate critical genes associated with multiple myeloma pathogenesis. Proc Natl Acad Sci USA. 2008, 105: 12885-12890. 10.1073/pnas.0806202105.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee Y, Jeon K, Lee JT, Kim S, Kim VN: MicroRNA maturation: stepwise processing and subcellular localization. EMBO J. 2002, 21: 4663-4670. 10.1093/emboj/cdf476.PubMedPubMed CentralView ArticleGoogle Scholar
- Kyle RA, Therneau TM, Rajkumar SV, Larson DR, Plevak MF, Offord JR, Dispenzieri A, Katzmann JA, Melton LJ: Prevalence of monoclonal gammopathy of undetermined significance. N Engl J Med. 2006, 354: 1362-1369. 10.1056/NEJMoa054494. 3PubMedView ArticleGoogle Scholar
- Zhan F, Huang Y, Colla S, Stewart JP, Hanamura I, Gupta S, Epstein J, Yaccoby S, Sawyer J, Burington B, et al: The molecular classification of multiple myeloma. Blood. 2006, 108: 2020-2028. 10.1182/blood-2005-11-013458.PubMedPubMed CentralView ArticleGoogle Scholar
- Fonseca R, Bergsagel PL, Drach J, Shaughnessy J, Gutierrez N, Stewart AK, Morgan G, Van Ness B, Chesi M, Minvielle S, et al: International Myeloma Working Group molecular classification of multiple myeloma: spotlight review. Leukemia. 2009, 23: 2210-2221. 10.1038/leu.2009.174.PubMedPubMed CentralView ArticleGoogle Scholar
- Ross FM, Ibrahim AH, Vilain-Holmes A, Winfield MO, Chiecchio L, Protheroe RK, Strike P, Gunasekera JL, Jones A, Harrison CJ, et al: Age has a profound effect on the incidence and significance of chromosome abnormalities in myeloma. Leukemia. 2005, 19: 1634-1642. 10.1038/sj.leu.2403857.PubMedView ArticleGoogle Scholar
- Kuehl WM, Bergsagel PL: Early genetic events provide the basis for a clinical classification of multiple myeloma. Hematology (Am Soc Hematol Educ Program). 2005, 346-352.Google Scholar
- Gutierrez NC, Sarasquete ME, Misiewicz-Krzeminska I, Delgado M, De Las Rivas J, Ticona FV, Ferminan E, Martin-Jimenez P, Chillon C, Risueno A, et al: Deregulation of microRNA expression in the different genetic subtypes of multiple myeloma and correlation with gene expression profiling. Leukemia. 2010, 24: 629-637. 10.1038/leu.2009.274.PubMedView ArticleGoogle Scholar
- Lionetti M, Biasiolo M, Agnelli L, Todoerti K, Mosca L, Fabris S, Sales G, Deliliers GL, Bicciato S, Lombardi L, et al: Identification of microRNA expression patterns and definition of a microRNA/mRNA regulatory network in distinct molecular groups of multiple myeloma. Blood. 2009, 114: e20-26. 10.1182/blood-2009-08-237495.PubMedView ArticleGoogle Scholar
- Greipp PR, San Miguel J, Durie BG, Crowley JJ, Barlogie B, Blade J, Boccadoro M, Child JA, Avet-Loiseau H, Kyle RA, et al: International staging system for multiple myeloma. J Clin Oncol. 2005, 23: 3412-3420. 10.1200/JCO.2005.04.242.PubMedView ArticleGoogle Scholar
- Sirohi B, Powles R, Kulkarni S, Rudin C, Saso R, Lal R, Singhal S, Mehta J, Horton C, Treleaven J: Comparison of new patients with Bence-Jones, IgG and IgA myeloma receiving sequential therapy: the need to regard these immunologic subtypes as separate disease entities with specific prognostic criteria. Bone Marrow Transplant. 2001, 28: 29-37. 10.1038/sj.bmt.1703093.PubMedView ArticleGoogle Scholar
- Kyle RA, Gertz MA, Witzig TE, Lust JA, Lacy MQ, Dispenzieri A, Fonseca R, Rajkumar SV, Offord JR, Larson DR, et al: Review of 1027 patients with newly diagnosed multiple myeloma. Mayo Clin Proc. 2003, 78: 21-33. 10.4065/78.1.21.PubMedView ArticleGoogle Scholar
- Shustik C, Bergsagel DE, Pruzanski W: Kappa and lambda light chain disease: survival rates and clinical manifestations. Blood. 1976, 48: 41-51.PubMedGoogle Scholar
- Huang Q, Gumireddy K, Schrier M, le Sage C, Nagel R, Nair S, Egan DA, Li A, Huang G, Klein-Szanto AJ, et al: The microRNAs miR-373 and miR-520c promote tumour invasion and metastasis. Nat Cell Biol. 2008, 10: 202-210. 10.1038/ncb1681.PubMedView ArticleGoogle Scholar
- Lawrie CH, Saunders NJ, Soneji S, Palazzo S, Dunlop HM, Cooper CD, Brown PJ, Troussard X, Mossafa H, Enver T, et al: MicroRNA expression in lymphocyte development and malignancy. Leukemia. 2008, 22: 1440-1446. 10.1038/sj.leu.2405083.PubMedView ArticleGoogle Scholar
- Zhang J, Jima DD, Jacobs C, Fischer R, Gottwein E, Huang G, Lugar PL, Lagoo AS, Rizzieri DA, Friedman DR, et al: Patterns of microRNA expression characterize stages of human B cell differentiation. Blood. 2009, 113: 4586-4594. 10.1182/blood-2008-09-178186.PubMedPubMed CentralView ArticleGoogle Scholar
- Merkel O, Hamacher F, Laimer D, Sifft E, Trajanoski Z, Scheideler M, Egger G, Hassler MR, Thallinger C, Schmatz A, et al: Identification of differential and functionally active miRNAs in both anaplastic lymphoma kinase (ALK)+ and ALK- anaplastic large-cell lymphoma. Proc Natl Acad Sci USA. 2010, 107: 16228-16233. 10.1073/pnas.1009719107.PubMedPubMed CentralView ArticleGoogle Scholar
- Lawrie CH, Cooper CD, Ballabio E, Chi J, Tramonti D, Hatton CS: Aberrant expression of microRNA biosynthetic pathway components is a common feature of haematological malignancy. Br J Haematol. 2009, 145: 545-548. 10.1111/j.1365-2141.2009.07642.x.PubMedView ArticleGoogle Scholar
- Wong QW, Lung RW, Law PT, Lai PB, Chan KY, To KF, Wong N: MicroRNA-223 is commonly repressed in hepatocellular carcinoma and potentiates expression of Stathmin1. Gastroenterology. 2008, 135: 257-269. 10.1053/j.gastro.2008.04.003.PubMedView ArticleGoogle Scholar
- Fazi F, Racanicchi S, Zardo G, Starnes LM, Mancini M, Travaglini L, Diverio D, Ammatuna E, Cimino G, Lo-Coco F, et al: Epigenetic silencing of the myelopoiesis regulator microRNA-223 by the AML1/ETO oncoprotein. Cancer Cell. 2007, 12: 457-466. 10.1016/j.ccr.2007.09.020.PubMedView ArticleGoogle Scholar
- Stamatopoulos B, Meuleman N, Haibe-Kains B, Saussoy P, Van Den Neste E, Michaux L, Heimann P, Martiat P, Bron D, Lagneaux L: microRNA-29c and microRNA-223 down-regulation has in vivo significance in chronic lymphocytic leukemia and improves disease risk stratification. Blood. 2009, 113: 5237-5245. 10.1182/blood-2008-11-189407.PubMedView ArticleGoogle Scholar
- Lawrie CH, Ballabio E, Dyar OJ, Jones M, Ventura R, Chi J, Tramonti D, Gooding S, Boultwood J, Wainscoat JS, et al: MicroRNA expression in chronic lymphocytic leukaemia. Br J Haematol. 2009, 147: 398-402. 10.1111/j.1365-2141.2009.07857.x.PubMedView ArticleGoogle Scholar
- Felli N, Pedini F, Romania P, Biffoni M, Morsilli O, Castelli G, Santoro S, Chicarella S, Sorrentino A, Peschle C, Marziali G: MicroRNA 223-dependent expression of LMO2 regulates normal erythropoiesis. Haematologica. 2009, 94: 479-486. 10.3324/haematol.2008.002345.PubMedPubMed CentralView ArticleGoogle Scholar
- Natkunam Y, Zhao S, Mason DY, Chen J, Taidi B, Jones M, Hammer AS, Hamilton Dutoit S, Lossos IS, Levy R: The oncoprotein LMO2 is expressed in normal germinal-center B cells and in human B-cell lymphomas. Blood. 2007, 109: 1636-1642. 10.1182/blood-2006-08-039024.PubMedPubMed CentralView ArticleGoogle Scholar
- Akao Y, Nakagawa Y, Kitade Y, Kinoshita T, Naoe T: Downregulation of microRNAs-143 and -145 in B-cell malignancies. Cancer Sci. 2007, 98: 1914-1920. 10.1111/j.1349-7006.2007.00618.x.PubMedView ArticleGoogle Scholar
- Landgren O, Kyle RA, Pfeiffer RM, Katzmann JA, Caporaso NE, Hayes RB, Dispenzieri A, Kumar S, Clark RJ, Baris D, et al: Monoclonal gammopathy of undetermined significance (MGUS) consistently precedes multiple myeloma: a prospective study. Blood. 2009, 113: 5412-5417. 10.1182/blood-2008-12-194241.PubMedPubMed CentralView ArticleGoogle Scholar
- Chan JA, Krichevsky AM, Kosik KS: MicroRNA-21 is an antiapoptotic factor in human glioblastoma cells. Cancer Res. 2005, 65: 6029-6033. 10.1158/0008-5472.CAN-05-0137.PubMedView ArticleGoogle Scholar
- Si ML, Zhu S, Wu H, Lu Z, Wu F, Mo YY: miR-21-mediated tumor growth. Oncogene. 2006, 26: 2799-2803.PubMedView ArticleGoogle Scholar
- Fonseca R, Oken MM, Harrington D, Bailey RJ, Van Wier SA, Henderson KJ, Kay NE, Van Ness B, Greipp PR, Dewald GW: Deletions of chromosome 13 in multiple myeloma identified by interphase FISH usually denote large deletions of the q arm or monosomy. Leukemia. 2001, 15: 981-986. 10.1038/sj.leu.2402125.PubMedView ArticleGoogle Scholar
- He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ, Hammond SM: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435: 828-833. 10.1038/nature03552.PubMedPubMed CentralView ArticleGoogle Scholar
- Zojer N, Konigsberg R, Ackermann J, Fritz E, Dallinger S, Kromer E, Kaufmann H, Riedl L, Gisslinger H, Schreiber S, et al: Deletion of 13q14 remains an independent adverse prognostic variable in multiple myeloma despite its frequent detection by interphase fluorescence in situ hybridization. Blood. 2000, 95: 1925-1930.PubMedGoogle Scholar
- Lawrie CH: MicroRNAs and haematology: small molecules, big function. Br J Haematol. 2007, 137: 503-512. 10.1111/j.1365-2141.2007.06611.x.PubMedView ArticleGoogle Scholar
- Cloonan N, Brown MK, Steptoe AL, Wani S, Chan WL, Forrest AR, Kolle G, Gabrielli B, Grimmond SM: The miR-17-5p microRNA is a key regulator of the G1/S phase cell cycle transition. Genome Biol. 2008, 9: R127-10.1186/gb-2008-9-8-r127.PubMedPubMed CentralView ArticleGoogle Scholar
- Aguda BD, Kim Y, Piper-Hunter MG, Friedman A, Marsh CB: MicroRNA regulation of a cancer network: consequences of the feedback loops involving miR-17-92, E2F, and Myc. Proc Natl Acad Sci USA. 2008, 105: 19678-19683. 10.1073/pnas.0811166106.PubMedPubMed CentralView ArticleGoogle Scholar
- Hossain A, Kuo MT, Saunders GF: Mir-17-5p regulates breast cancer cell proliferation by inhibiting translation of AIB1 mRNA. Mol Cell Biol. 2006, 26: 8191-8201. 10.1128/MCB.00242-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu J, Liao X, Wong C: Downregulations of B-cell lymphoma 2 and myeloid cell leukemia sequence 1 by microRNA 153 induce apoptosis in a glioblastoma cell line DBTRG-05MG. Int J Cancer. 2010, 126: 1029-1035.PubMedGoogle Scholar
- Wuilleme-Toumi S, Robillard N, Gomez P, Moreau P, Le Gouill S, Avet-Loiseau H, Harousseau JL, Amiot M, Bataille R: Mcl-1 is overexpressed in multiple myeloma and associated with relapse and shorter survival. Leukemia. 2005, 19: 1248-1252. 10.1038/sj.leu.2403784.PubMedView ArticleGoogle Scholar
- Hiroki E, Akahira J, Suzuki F, Nagase S, Ito K, Suzuki T, Sasano H, Yaegashi N: Changes in microRNA expression levels correlate with clinicopathological features and prognoses in endometrial serous adenocarcinomas. Cancer Sci. 2010, 101: 241-249. 10.1111/j.1349-7006.2009.01385.x.PubMedView ArticleGoogle Scholar
- Wei JJ, Wu X, Peng Y, Shi G, Olca B, Yang X, Daniels G, Osman I, Ouyang J, Hernando E, et al: Regulation of HMGA1 expression by microRNA296 affects prostate cancer growth and invasion. Clin Cancer Res. 2010Google Scholar
- Durie BG, Kyle RA, Belch A, Bensinger W, Blade J, Boccadoro M, Child JA, Comenzo R, Djulbegovic B, Fantl D, et al: Myeloma management guidelines: a consensus report from the Scientific Advisors of the International Myeloma Foundation. Hematol J. 2003, 4: 379-398. 10.1038/sj.thj.6200312.PubMedView ArticleGoogle Scholar
- Wuilleme S, Robillard N, Lode L, Magrangeas F, Beris H, Harousseau JL, Proffitt J, Minvielle S, Avet-Loiseau H: Ploidy, as detected by fluorescence in situ hybridization, defines different subgroups in multiple myeloma. Leukemia. 2005, 19: 275-278. 10.1038/sj.leu.2403586.PubMedView ArticleGoogle Scholar
- Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.