Skip to main content

In non-transplant patients with multiple myeloma, the pre-treatment level of clonotypic cells predicts event-free survival



In multiple myeloma (MM), the immunoglobulin heavy chain VDJ gene rearrangement is a unique clonotypic signature that identifies all members of the myeloma clone independent of morphology or phenotype. Each clonotypic MM cell has only one genomic copy of the rearranged IgH VDJ.


Pre-treatment bone marrow aspirates from myeloma patients at diagnosis or in relapse were evaluated for the number of clonotypic cells using real time quantitative PCR (RPCR). RPCR measured the level of clonal cells, termed VDJ%, in 139 diagnosis and relapse BM aspirates from MM patients.


Patients with a VDJ% below the median had a significantly longer event free survival (EFS) then those with a VDJ% higher than the median (p=0.0077, HR=0.57). Further, although the VDJ% from non-transplant patients predicted EFS (p=0.0093), VDJ% failed to predict outcome after autologous stem cell transplant (p=0.53).


Our results suggest that for non-transplant patients, the tumor burden before treatment, perhaps reflecting cancer stem cell progeny/output, is an indirect measure that may indicate the number of MM cancer stem cells and hence event free survival.


Multiple myeloma (MM) is a clonal B-cell malignancy with clinical identification of plasma cells (PCs) that accumulate in the bone marrow (BM). Newer treatments continue to increase the incidence of remission and decrease side effects, but relapse invariably occurs. Although the MM clone is known to be heterogeneous, including earlier stage B cells[16] and putative cancer stem cells[47] as well as PC, current diagnostic and monitoring criteria are based mainly on BM plasmacytosis (BMPC) and serum M-protein levels (Mpr) as measured in the clinical laboratory. Although measures of BMPC morphology and Mpr are clinically informative[812], they are not among the parameters recommended for risk stratification[13, 14]. The IgH VDJ gene rearrangement that characterizes the MM clone in each patient, termed clonotypic, is a molecular signature that identifies all malignant components regardless of differentiation stage, phenotype, morphology or other molecular markers[2, 15]. Even though abnormal circulating or BM PC phenotypes[12, 1619] provide informative prognostic measures, particularly for autologous stem cell transplantation, it would be of value to confirm that the malignant cells themselves express these markers, as well as providing additional predictive value[2025]. Molecular analysis may also be important for evaluating risk stratification in patients receiving novel therapies[2628].

Because each malignant cell harbors only one copy of the clonotypic IgH VDJ, using quantitative PCR the malignant clone can be enumerated in genomic DNA from a heterogeneous population of cells. Here we used SYBR Green real-time quantitative PCR (RPCR) with patient specific CDR2/CDR3 primers to measure the fraction of clonal cells among total BM cells, termed VDJ%, in 139 diagnosis and relapse BM aspirates from MM patients. The results presented below suggest that the VDJ% from the diagnosis and relapse BM prior to initiating a given treatment is an independent measure of the remission length in non-transplant patients, unlike BMPC or Mpr. While standard clinical tests failed to predict remission duration, the percent of clonal cells as measured in BM aspirates with RPCR stratified patients into high risk and low risk categories. The VDJ% enumerates all cells containing the clonal sequence and thus may include relevant cell types not included among BMPC or contributing to Mpr levels. Pre-treatment VDJ% was independent of disease response and drug choice suggesting that absent autologous stem cell transplant, it may be the specific biology of each tumor that determines remission length and not the specific treatment.


Although it has been well documented that clinically-measured levels of residual disease during remissions of hematological cancers can be significantly associated with event free survival, there has been little documentation of the impact on outcome of total tumor burden prior to treatment. This study utilizes patient-specific RPCR to molecularly estimate the total tumor burden, as defined by the IgH VDJ molecular signature, in the BM of MM patients. Patient-specific RPCR quantifies all cell types harboring the clonotypic signature, whether or not they are clinically detectable.

Sensitivity of the RPCR assay

To determine the lowest level of detection, DNA from a MM PC cell line (LP-1) with a known clonotypic IgH VDJ rearrangement was diluted in a constant amount of normal PBMC DNA to simulate patient samples with different clonal percentages. The average molecule count for each dilution, determined using a standard curve for the VDJ and β2m PCRs, was used to calculate the percent of clonal cells (Figure1). The lowest dilution with detectable VDJ signal was 0.001%. Thus, patient samples that fell below this threshold value were arbitrarily assigned a value of 0.001% for inclusion in the analysis.

Figure 1
figure 1

Sensitivity of patient-specific RPCR. The bars represent the mean of 3 replicate PCR reactions to determine the VDJ% for each dilution noted by comparing the VDJ PCR molecule count plotted from the VDJ control curve to the β2m molecule count using half the DNA. The absence of bars at the 0.0001%, 0.00001% and 0.000001% levels indicate the absence of amplification.

A limiting dilution assay was used to confirm the RPCR results from a BM sample. The VDJ% was calculated using the cell equivalents added and the PCR results for each dilution. For the same BM sample, the clonal percent estimate was 7.97% based on limiting dilution and 7.93% based on RPCR.

Pre-treatment VDJ% correlates with EFS

One hundred and thirty nine BM aspirate samples from MM patients were collected just before treatment, 89 of which were from previously untreated MMs and 50 were relapsed/refractory (RR) MMs starting a new line of therapy (Table1). BM aspirates should be understood as sampling a geographically distinct location, a caveat that applies to all aspirate-based studies; the generally held but unproven working assumption is that the aspirate is representative of the BM as a whole. Comparison with standard curves gave a molecule count for VDJ targets and β2m targets. A percent of clonal cells in BMMC, termed VDJ% was then determined. Based on the VDJ%, the patient cohort was dichotomized into two groups: above or below the median of 18.2%. Event free survival (EFS), was significantly different between them, as shown by Kaplan Meier analysis (p=0.0269). Those patient samples with a VDJ% falling below the median had a longer EFS than those above the median (Figure2A). The hazard ratio (HR) for patients below the median VDJ% compared to patients above the median VDJ% is 0.6202. The median EFS time for above and below the median VDJ% was 490 days and 961 days, respectively.

Table 1 Clinical characteristics: the high VDJ% group includes more patients with ASCT than with conventional chemotherapy
Figure 2
figure 2

Kaplan-Meier analysis of pre-treatment BM VDJ% and EFS. The VDJ%, dichotomized by the median, shows a significant correlation with EFS (A) - those who fall below the median (light line) have longer EFSs than those who fall above the median (dark line). Statistical determination of the optimal cut point for dichotomizing into two groups gave a value nearly identical to the median (18.6%). There were 39 relapses out of 69 patient samples below the median and 43 of 70 above the median. When patient samples were stratified to exclude those receiving an ASCT (the group of non-transplant patients), the significance remains (B). When the remaining patients who did receive an ASCT were analyzed, the significant correlation was lost (C). D represents the Kaplan-Meier analyses for three significantly different EFS categories: lightest line=all ASCT treated patients, medium grey line= non-transplant patients below the VDJ% median, and the dark line= non-transplant patients above the VDJ% median, with respectively 20 relapses out of 50 samples, 30 out of 45 and 32 out of 44. E summarizes the logrank p-value, hazard ratio and median survivals for the above plots.

The VDJ% prior to therapy was not associated with overall survival (not shown). A significant association was found between pre-treatment VDJ% and β2m, LDH, hemoglobin or creatinine, but not with any other factors noted in Table2.

Table 2 Clinical characteristics of patients with high or low VDJ%

Sub-dividing the patient cohort to include only non- autologous stem cell transplant (ASCT) patients showed a significant association between VDJ% and EFS in a Kaplan Meier analysis (p=0.0081, HR=0.5252) (Figure2B). In contrast, the VDJ% from the subset of patients who received an ASCT as part of their treatment regimen failed to show a significant association with EFS by Kaplan Meier analysis (p=0.5306, HR=0.7569) (not shown), despite the fact that a large proportion of patients who received ASCT fell into the high VDJ% group (Table1). Significantly different EFS categories were observed between ASCT patients (1012 days), non-ASCT treated patients with a VDJ% below median (740 days) and non-ASCT treated patients with a VDJ% above the median (319 days) (Figure2D, E). The prolonged survival for ASCT patients occurred even though a large proportion had a high VDJ% (Table1).

In a Cox multivariate analysis adjusted for International Staging System (ISS), patients with a VDJ% higher than the median have a shorter EFS than do patients with a low VDJ% (p = .002 for non-ASCT and a trend for ASCT patients p=.11). In a Cox regression univariate model of the combined patient groups adjusted for ISS, patients with high VDJ% also have shorter EFS (p=.005).

VDJ% but not BMPC significantly associates with EFS

Part of the clinical evaluation of MM is to determine the % PC in the BM aspirate, termed BM plasmacytosis (BMPC), a measure not thought to be associated with outcome. RPCR utilized DNA from pre-treatment BM samples for which BMPC values were available from the routine clinical analysis of BM aspirates. For the cohort analyzed here in a Kaplan Meier analysis where the median dichotomized the group, the BMPC failed to show a significant association with EFS (Figure3A). Furthermore, no significant correlations were seen with any clinical outcome when monoclonal protein levels (Mpr) were evaluated (not shown). In contrast, when the VDJ% dichotomized the full cohort for a Kaplan Meier analysis, a significant correlation was observed. If only non-ASCT patients were included in the Kaplan-Meier analysis of BMPC and EFS, there continued to be no significant association (Figure3C). Again in contrast, for VDJ%, a significant correlation with EFS was seen for the non-ASCT cohort (Figures3B and D).

Figure 3
figure 3

Kaplan Meier Curves of BMPC or Matched VDJ% to EFS. The median BMPC percentage dichotomized the group but did not indicate significantly different EFSs (A). The same patient cohort, categorized by VDJ% above and below the median, significantly correlates with EFS (B). If only non-transplant patient samples were included, the EFS compared to the BMPC remains insignificant (C) and the EFS compared to VDJ% remains significant (D).

For a subset of 47 patients, we evaluated the relationship between %PC (as measured by flow cytometry) and VDJ% for the same sample of purified BMMC (data not shown). Given that the majority of clonotypic cells in untreated or relapsed/refractory MM BMMC are likely to be PC, the expected significant correlation was observed between %PC and VDJ% (p<.0001). This however does not allow us to determine the extent to which non-PC contribute to the VDJ%, as implied by the fact that VDJ% but not BMPC predicts EFS (Figure3).

The VDJ% from remission BMs may not associate with outcome

Little research has addressed the VDJ% for patients in remission after treatment for relapsed MM. In a preliminary analysis, for a small set of 32 relapsed patients who were subsequently treated and later achieved remission, the tumor burden in remission BM samples (post-treatment) was quantified using patient-specific RPCR. The samples were from different patients with different treatment courses, including frontline VAD/ASCT and bortezomib/ASCT (20/32 patients), or relapse samples treated with dexamethasone, bortezomib or lenalidomide (12/32 patients). Some remission samples did not have detectable VDJ%: these were arbitrarily assigned a value of 0.001%.

In a Kaplan Meier analysis, the remission VDJ% did not show a significant correlation with EFS (p=0.6934, HR=0.8026). Of the 25 patients with matched pre-treatment and remission samples, the relative reduction of VDJ% was calculated. This ‘relative reduction’ was not associated with EFS in a Kaplan-Meier analysis (p=0.4620, HR=0.6041).

As the remission BM samples were taken at varying times in the treatment course, time to progression (TTP) provides an alternate measure of outcome. TTP is calculated as the number of days between aspirating the remission BM and the time of clinical relapse. Neither the remission VDJ% or the relative reduction were significantly associated with TTP (p=0.4407 and p=0.5430, respectively) (data not shown).

The remission VDJ% values and relative reduction were compared for VAD and bortezomib as frontline therapies (Figure4A and B), and for dexamethasone, bortezomib or lenalidomide (Figure4C). The number of patients is low, but no significant differences were found among these groups.

Figure 4
figure 4

Remission VDJ% and relative reduction of remission BM samples. The post treatment VDJ% from a front-line trial of relapsed MM patients, comparing standard VAD to bortezomib just prior to ASCT (A). The times reduction of the VDJ% in the same samples are plotted in B. In relapsed patients, the VDJ% is noted in C after having been treated with dexamethasone, lenalidomide or bortezomib. The line through each column represents the mean. Dex=dexamethasone.


This study demonstrates that pre-treatment BM VDJ% has an inverse correlation with event free survival. The percent of clonal cells measured by RPCR in aspirated pre-treatment BM cells significantly correlated with EFS in a Kaplan-Meier analysis. Patients were treated with a variety of therapies for front-line and relapsed/refractory MM. VDJ% in non-transplant patients receiving chemotherapy or novel therapies (e.g. lenalidomide or bortezomib) correlated inversely with EFS, but VDJ% was not predictive for patients receiving an ASCT. Subsequently, three significantly different remission durations were identified. Patients receiving an ASCT had the longest remission duration. Patients not eligible for an ACST and having a VDJ% below the median (17% in this data set) had a longer remission duration than did those having a VDJ% above the median.

Others have shown that patients receiving high dose cytoreduction with auto- or allo-graft rescue are more likely to achieve molecular complete remissions and a subset of these patients have improved progression-free survival[29, 30]. It appears that the benefits of complete remissions may occur mainly in the context of intensive therapy and stem cell rescue[11, 21, 31, 32]. Consistent with these observations, the loss of correlation between pre-treatment VDJ% and outcome, compared to the non-transplant group, may further highlight the efficacy of high dose therapy followed by auto- or allo- stem cell transplant. It appears that an ASCT overcomes the contribution to outcome of pre-treatment tumor burden, even when that burden is high. Others have reported the pre-treatment levels of clonal cells[20, 25, 33] but did not evaluate associations between tumor burden and outcome. Most studies evaluated cohorts of patients who had received ASCT. Consistent with this, our work shows an absence of prognostic significance for the pre-treatment VDJ% in MM patients who ultimately received an ASCT.

It is important to note that most studies identifying EFS prognostic indicators, including serum β2m and albumin levels[34, 35], chromosome 13 monosomy[36, 37], and t(4;14) translocations[38], have evaluated samples from previously untreated patients. However, patients can only be ‘previously untreated’ once, and for the purposes of defining relapse risks for each cycle of therapy, the influence of these factors is less defined. Our study considers all the treatment cycles. VDJ% is predictive of EFS in untreated and relapsed MM, with the caveat that this correlation does not hold for patients who were ultimately treated with ASCT. As part of clinical testing, VDJ% may indicate to a treating clinician the current likelihood of early relapse.

In this study, the VDJ% from remission/post treatment BM samples could not predict outcome, albeit with only small numbers of patients. This is in contrast to previously published work, also on relatively small sample sets, showing that the frequency of clonal cells does correlate with EFS[29, 32, 39]. However, the small cohorts in previous studies were from uniformly treated patients, often collected at similar time points. In contrast, the remission BM samples used in this study were from patients treated with a variety of therapeutic regimens, although the majority (20/32) were in remission after relapsing from intensive therapy and ASCT. Samples were also collected from patients receiving lenalidomide, bortezomib, or dexamethasone for relapsed disease (12/32), taken at varying time points during remission. As no significant correlation was observed between the remission VDJ% or relative reduction in VDJ% after treatment and EFS from these unequal time points, the TTP was used to standardize the time points of BM sampling. Since the remission samples studied here were obtained from heterogeneous patient populations, the results may be broadly applicable. Analysis of a much larger cohort will be needed to determine whether or not these are equally effective treatments for reducing the disease burden as measured molecularly.

Although still very speculative, if the conclusion is correct that the remission BM VDJ% lacks a correlation with outcome, this may provide a new perspective on myeloma growth. The VDJ% enumerates all cells containing the IgH VDJ clonal sequence, regardless of their morphology or differentiation stage. In pre-treatment BM, the number of clonal cells might be expected to have a proportional relationship to the number of progenitor cells/cancer stem cells. A BM aspirate would contain a number of clonal cells that may numerically reflect this generative population. When the tumor is treated, this representative population is disproportionately lost but drug-resistant cancer stem cells persist, such that the aggregate VDJ% no longer reflects the size of the progenitor compartment. The generative compartment escapes therapy and after therapy can replenish the tumor mass. Morphologically defined BMPC on the other hand, explicitly quantify only PC. The difference between VDJ% and BMPC as measures of clonality may reflect a population(s) of cells that are enumerated by RPCR but remain undetected by morphologic inspection. These cells may be responsible for the correlation between EFS and VDJ% as well as the lack of correlation between EFS and BMPC. Furthermore, these non-PC clonal cells may contribute to malignant clonal expansion.

While the VDJ% in pre-treated BMs shows a significant correlation to EFS, BMPC or Mpr levels (as surrogate indicators of PC but not of earlier stage clonotypic cells) do not. We speculate that this may indicate a significant population(s) of cells not detected by routine clinical tests that has a strong influence on the course of disease. Our current methods do not allow us to directly evaluate the nature or number of such non-PC clonotypic cells, but our previous work suggests that they are frequent[2, 6, 15]. Furthermore, for the patients studied, this cell population appears to be depleted by subsequent ASCT but not by other forms of treatment. One would expect that if there were an even distribution of cell types killed by the therapy, then the tumor burden (VDJ%) in a remission BM would be predictive. Given that this may not always apply, it seems likely that therapy might sometimes have a non-random impact by killing some members of the MM clone but not others, because only some compartments respond. Speculatively, bone marrow prior to treatment (Figure5, left) contains a population of cells that proportionately represent the number of putative cancer stem cells[4, 6]. On the right, a remission bone, plasma cells have been depleted by non-transplant therapy, but progenitors remain, so the sampling no longer predicts the number of cancer stem cells. In this situation, unlike other treatments, ASCT may deplete the progenitor population such that the pre-treatment numbers no longer have predictive value.

Figure 5
figure 5

Hypothesis: a possible perspective on myeloma growth. Green cells represent plasma cells; yellow cells represent cancer stem cells/myeloma progenitors. This model proposes that the number of plasma cells accurately reflects the number of generative cancer stem cells prior to therapy, after diagnosis or in relapse (left panel) but does not bear a numerical relationship to cancer stem cells in remission (right panel). We speculate that this may be due to selective killing of plasma cells by therapy that spares cancer stem cells. In this context, the pre-treatment VDJ% may correlate with EFS because it is an indirect measure of cancer stem cell activity.


Our work can be interpreted in terms of a clonotypic cancer stem cell population in MM BM[4]. Taken together, our results suggest that except for patients treated with ASCT, tumor burden/cancer stem cell output before treatment is an indirect measure that may reflect the number of MM cancer stem cells and hence event free survival.

Material and methods

Patients, clinical data and outcome

BM aspirates were collected from 171 MM patients and the clonotypic IgH VDJ rearrangement was identified in all of them. Written informed consent was obtained, as approved by the University of Alberta and Alberta Health Services Human Ethics Review Committees. One hundred and thirty nine BM aspirate samples (Table1) were utilized in the pre-treatment study of clonal (VDJ) percent. These included 89 previously untreated MM patients and 50 relapsed/refractory MM patients. BM plasmacytosis (BMPC) and monoclonal protein levels were determined during routine clinical testing and the values used for analysis were obtained from patient records. The % BMPC was determined by the clinical laboratory from Wright-Giemsa stained BM particle crush slides. Thirty-two remission BM aspirates were also analyzed. These samples were drawn from patients at varying time points during remission (e.g. after ASCT or at CR). Patients received initial or new rounds of therapy after the BM aspirate was taken. Patient treatments include single agent or multi agent therapy using combinations of vincristine, doxorubicin, melphalan, prednisone, thalidomide, dexamethasone, cyclophosphamide, lenalidomide or bortezomib and some received an ASCT. Patient samples were unfractionated dry cell pellets stored at −80°C, from which DNA was purified. Remission duration, as defined by event-free survival (EFS), was defined as the number of days between the BM sampling (at diagnosis or relapse) and clinical relapse (defined by a rise in peripheral or urine M-protein level 40) or death. Time to progression (TTP) was measured as the number of days between the remission BM sampling and relapse/death. Overall survival (OS) is defined as the number of days from diagnosis until last follow-up or death.

Cell purification and DNA isolation

BM aspirates were processed to obtain mononuclear cells (BMMC) using Ficoll-Paque PLUS (GE Healthcare, Piscataway, NJ, USA) as described previously[15]. BM mononuclear cells (BMMC) were archived at −80°C as pellets for later DNA isolation; because the entire pellet is used for the DNA extraction, there is no loss of cellular material from the thawed pellet. DNAzol Reagent (Invitrogen, Carlsbad, CA, USA) was used to extract the DNA according to the manufacturer’s instructions. The DNA quantification was performed using the NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) to ensure equal loading for the PCR applications.

Patient specific PCR primers

Primers were specific to each patient’s VDJ rearrangement as previously described[40]. Briefly, a panel of 12 V region 5’ primers are used with either an IgA or IgG constant region 3’ primer to PCR amplify any rearranged IgH transcripts, followed by direct DNA sequencing. Primers were designed over the CDR1 or CDR2 region and the CDR3 region. The primers were tested on the unfractionated BMMC sample from the same patient. In parallel with a standard PCR control (water), samples from other patients with similar VH regions controlled for primer specificity. When possible, primers were further validated using single sorted CD38+CD138+ PCs from the BMMC of the same patient. Twenty-four single PCs were tested by PCR using patient specific primers; a success rate of 75% or higher was required to validate identification of the clonotypic IgH VDJ.

Clonal cell enumeration by RPCR

RPCR was performed on genomic DNA from BMMCs using the DNA Engine Opticon 2 (Bio-Rad Laboratories, Hercules, CA, USA). The DyNAmo HS SYBR Green qPCR Kit (Finnzymes, Espoo, Finland) was used for the RPCR reaction as follows: 10 μl of DynAmo Master Mix, 0.25μM of each patient specific primer, genomic DNA and water to 20μl. For the patient specific reactions using BMMC samples, 150ng of DNA was added. The control PCR amplified an intron/exon boundary of the β2m gene (to eliminate contamination by cDNA) and 75ng each BM reaction of genomic DNA was added to the PCR to compensate for the two genomic copies of β2m versus the single rearranged clonotypic VDJ copy. The cycling conditions for all the RPCR reactions were 15min at 95°C, 45 cycles of 15sec at 94°C, 30sec at 60°C and 30sec at 72°C, and a final extension phase of 2min at 72°C. Opticon Monitor 3 software (Bio-Rad Laboratories, Hercules, CA, USA) was used to analyze the RPCR data. The threshold was placed at the midway point of the log amplification curve to generate the C(T) value for each RPCR reaction. The amplicon was verified through melting curve analysis.

A cloned amplicon was used to generate a standard curve to determine the exact molecule count for each reaction. The β2m and patient specific IgH VDJ amplicons were cloned using pGEM-T Easy Vector System I (Promega, Madison, WI, USA) and Subcloning Efficiency DH5α cells (Invitrogen, Carlsbad, CA, USA) per manufacturer’s instructions. Plasmids were isolated using QIAprep Spin Miniprep Kit (Qiagen Inc., Mississauga, Ontario, Canada) and the identity of the insert was confirmed by PCR with the primers used to generate the original amplicon. A 10-fold dilution series was constructed and at least three replicates for each dilution were subjected to RPCR analysis concurrent with the BMMC sample. The software generated a standard curve where the slope defines the relationship between the C(T) value to a molecule count. The percentage of MM clonal cells, termed VDJ%, was calculated based on the molecule count of patient specific VDJ targets as compared to β2m the β2M molecule count. Β2M is encoded on chromosome ch15. A subset of MM patients have PC with trisomy 15 and thus likely 3 copies of β2M. This may result in an undefined extent of underestimate for the VDJ% in such patients[41]. However, each MM BM sample analyzed includes a large proportion of normal cells, all of which are diploid. We find that the median % PC in ficoll-purified BMC from MM patients is about 20%, confirming that the majority of BMMC in any given sample are likely to be normal cells with only 2 copies of the B2M gene. Monosomy or trisomy of chromosome 14 would also influence the calculation of VDJ%, depending on which ch14 has been affected; it is likely these would be equivalently distributed in both groups of patients. Given than only one ch14 has the clonotypic IgH VDJ, there is no way to know if the deleted ch14 in monosomies (10% of patients) or the extra ch14 in trisomies (14% of patients) harbors the clonotypic VDJ or an unrearranged/nonproductive VDJ. On a random basis, the affected chromosome will be the productively rearranged VDJ in about half of structural abnormalities, but the prevalence of M protein suggests that for monosomies, the ch14 with a productive (clonotypic) VDJ is not affected. For ch14 translocations, the fact that nearly all such MM make M protein indicates that the productively rearranged VDJ (clonotypic) is not affected by the translocation. Consideration of chromosomal abnormalities as related to RPCR is relevant only if in any given patient all cells have the same chromosomal count, an assumption that is challenged by the presence of normal cells in BMMC and the extensive intra-clonal heterogeneity in individual MM patients[4245]. These caveats aside, the resulting percentage represents the number of clonal cells present in total BMMC. Where both a pre-treatment and remission VDJ% were available, the relative reduction was calculated by dividing the pre-treatment VDJ% by the remission VDJ%.

Cloned amplicons were generated for only the first 53 patients from the cohort. The C(T) values generated from the individual standard curve reactions were plotted together to form an aggregate regression curve. The resolved slope formula was used to generate new aggregate VDJ% values using the original samples C(T) values where x was the C(T) value and y was the molecule count. A significant correlation between the VDJ% values calculated with the individual standard curves versus those using the aggregate standard curve was found using a Spearman’s correlation (p<0.0001, r=0.5530) and a Fishers Exact test based on the medians (p<0.0001). The VDJ% values calculated using the two methods were plotted in a Kaplan-Meier curve and a logrank test showed that both generate significant associations with EFS. Based on this analysis, for the remaining 118 patients, an aggregate VDJ curve was verified for use in place of a patient specific curve.

To determine sensitivity of the RPCR, DNA from the LP-1 cell line was used. DNA was purified from normal peripheral blood mononuclear cells and spiked with LP-1 DNA. To each VDJ reaction, 150ng of DNA from each dilution was added, and to each β2M reaction, 75ng of DNA was added to compensate for two β2M gene loci compared to the single rearranged VDJ locus.

Statistical tests

To test for correlations with outcome, including event free survival and overall survival, Kaplan Meier curves were generated using GraphPad Prism 5 and a logrank statistics generated a p-value. Correlation tests were performed using Pearson’s Correlation, Spearman’s Correlation or Fisher Exact tests wherever applicable. Kaplan-Meier curves were generated for data using timed endpoints and longrank tests determined the statistical significance. A Kruskal-Wallis test or an ANOVA was used to compare the medians of more than two groups of data. These tests were done using Sigma Plot Software.



autologous stem cell transplant


Bone marrow


Bone marrow mononuclear cells


Bone marrow plasmacytosis




Complete remission


Event free survival


Hazard ratio


Immunoglobulin heavy chain


International Staging System


Multiple myeloma


Monoclonal protein


Overall survival


Plasma cells


Polymerase chain reaction




Real-time quantitative PCR


Time to progression


Variable, diversity and joining segments of Ig.


  1. Pilarski LM, Hipperson G, Seeberger K: Myeloma progenitors in the blood of patients with aggressive or minimal disease: engraftment and self-renewal of primary human myeloma in the bone marrow of NOD SCID mice. Blood. 2000, 95: 1056-1065.

    CAS  PubMed  Google Scholar 

  2. Pilarski LM, Baigorri E, Mant MJ: Multiple myeloma includes CD20+ B and plasma cells that persist in patients treated with rituximab. Clinical Medicine: Oncology. 2008, 2: 275-285.

    PubMed Central  CAS  Google Scholar 

  3. Reiman T, Seeberger K, Taylor BJ: Persistent pre-switch myeloma cells correlate with decreased survival: evidence for isotype switching within the myeloma clone. Blood. 2001, 98: 2791-2799. 10.1182/blood.V98.9.2791

    Article  CAS  PubMed  Google Scholar 

  4. Kirshner J, Thulien K, Martin LD: A unique 3-D model for evaluating the impact of therapy on multiple myeloma. Blood. 2008, 112: 2935-2945. 10.1182/blood-2008-02-142430

    Article  CAS  PubMed  Google Scholar 

  5. Matsui W, Huff CA, Wang Q: Characterization of clonogenic multiple myeloma cells. Blood. 2004, 103: 2332-2336. 10.1182/blood-2003-09-3064

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Debes Marun CS, Belch AR, Pilarski LM: In multiple myeloma, bone-marrow lymphocytes harboring the same chromosomal abnormalities as autologous plasma cells predict poor survival. Am J Hematol. 2012, 87: 579-587. 10.1002/ajh.23194

    Article  CAS  PubMed  Google Scholar 

  7. Jakubikova J, Adamia S, Kost-Alimova M: Lenalidomide targets clonogenic side population in multiple myeloma: pathophysiologic and clinical implications. Blood. 2011, 117: 4409-4419. 10.1182/blood-2010-02-267344

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Dimopoulos M, Kyle R, Fermand JP: Consensus recommendations for standard investigative workup: report of the International Myeloma Workshop Consensus Panel 3. Blood. 2011, 117: 4701-4705. 10.1182/blood-2010-10-299529

    Article  CAS  PubMed  Google Scholar 

  9. Paul E, Sutlu T, Deneberg S: Impact of chromosome 13 deletion and plasma cell load on long-term survival of patients with multiple myeloma undergoing autologous transplantation. Oncol Rep. 2009, 22: 137-142.

    CAS  PubMed  Google Scholar 

  10. Hwang IH, Chung JS, Shin HJ: Predictive value of post-transplant bone marrow plasma cell percent in multiple myeloma patients undergone autologous transplantation. Korean J Intern Med. 2011, 26: 76-81. 10.3904/kjim.2011.26.1.76

    Article  PubMed Central  PubMed  Google Scholar 

  11. Harousseau JL, Attal M, Vet-Loiseau H: The role of complete response in multiple myeloma. Blood. 2009, 114: 3139-3146. 10.1182/blood-2009-03-201053

    Article  CAS  PubMed  Google Scholar 

  12. Paiva B, Vidriales MB, Perez JJ: Multiparameter flow cytometry quantification of bone marrow plasma cells at diagnosis provides more prognostic information than morphological assessment in myeloma patients. Haematologica. 2009, 94: 1599-1602. 10.3324/haematol.2009.009100

    Article  PubMed Central  PubMed  Google Scholar 

  13. Munshi NC, Anderson KC, Bergsagel PL: Consensus recommendations for risk stratification in multiple myeloma: report of the International Myeloma Workshop Consensus Panel 2. Blood. 2011, 117: 4696-4700. 10.1182/blood-2010-10-300970

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Kapoor P, Rajkumar SV: Update on risk stratification and treatment of newly diagnosed multiple myeloma. Int J Hematol. 2011, 94: 310-320. 10.1007/s12185-011-0947-z

    Article  CAS  PubMed  Google Scholar 

  15. Szczepek AJ, Seeberger K, Wizniak J: A high frequency of circulating B cells share clonotypic IgH VDJ rearrangements with autologous bone marrow plasma cells in multiple myeloma, as measured by single cell and in situ RT-PCR. Blood. 1998, 92: 2844-2855.

    CAS  PubMed  Google Scholar 

  16. Paiva B, Gutierrez NC, Rosinol L: High-risk cytogenetics and persistent minimal residual disease by multiparameter flow cytometry predict unsustained complete response after autologous stem cell transplantation in multiple myeloma. Blood. 2012, 119: 687-691. 10.1182/blood-2011-07-370460

    Article  CAS  PubMed  Google Scholar 

  17. Peceliunas V, Janiulioniene A, Matuzeviciene R, Zvirblis T, Griskevicius L: Circulating plasma cells predict the outcome of relapsed or refractory multiple myeloma. Leuk Lymphoma. 2012, 53: 641-643. 10.3109/10428194.2011.627481

    Article  CAS  PubMed  Google Scholar 

  18. Paiva B, Vidriales MB, Cervero J: Multiparameter flow cytometric remission is the most relevant prognostic factor for multiple myeloma patients who undergo autologous stem cell transplantation. Blood. 2008, 112: 4017-4023. 10.1182/blood-2008-05-159624

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Mateo G, Montalban MA, Vidriales MB: Prognostic value of immunophenotyping in multiple myeloma: a study by the PETHEMA/GEM cooperative study groups on patients uniformly treated with high-dose therapy. J Clin Oncol. 2008, 26: 2737-2744. 10.1200/JCO.2007.15.4120

    Article  PubMed  Google Scholar 

  20. Cremer FW, Ehrbrecht E, Kiel K: Evaluation of the kinetics of the bone marrow tumor load in the course of sequential high-dose therapy assessed by quantitative PCR as a predictive parameter in patients with multiple myeloma. Bone Marrow Transplant. 2000, 26: 851-858. 10.1038/sj.bmt.1702628

    Article  CAS  PubMed  Google Scholar 

  21. Corradini P, Cavo M, Lokhorst H: Molecular remission after myeloablative allogeneic stem cell transplantation predicts a better relapse-free survival in patients with multiple myeloma. Blood. 2003, 102: 1927-1929. 10.1182/blood-2003-01-0189

    Article  CAS  PubMed  Google Scholar 

  22. Lipinski E, Cremer FW, Ho AD, Goldschmidt H, Moos M: Molecular monitoring of the tumor load predicts progressive disease in patients with multiple myeloma after high-dose therapy with autologous peripheral blood stem cell transplantation. Bone Marrow Transplant. 2001, 28: 957-962. 10.1038/sj.bmt.1703276

    Article  CAS  PubMed  Google Scholar 

  23. Bakkus MH, Bouko Y, Samson D: Post-transplantation tumour load in bone marrow, as assessed by quantitative ASO-PCR, is a prognostic parameter in multiple myeloma. Br J Haematol. 2004, 126: 665-674. 10.1111/j.1365-2141.2004.05120.x

    Article  PubMed  Google Scholar 

  24. Raab MS, Cremer FW, Breitkreutz IN: Molecular monitoring of tumour load kinetics predicts disease progression after non-myeloablative allogeneic stem cell transplantation in multiple myeloma. Ann Oncol. 2005, 16: 611-617. 10.1093/annonc/mdi123

    Article  CAS  PubMed  Google Scholar 

  25. Fenk R, Ak M, Kobbe G: Levels of minimal residual disease detected by quantitative molecular monitoring herald relapse in patients with multiple myeloma. Haematologica. 2004, 89: 557-566.

    CAS  PubMed  Google Scholar 

  26. Laubach JP, Mitsiades CS, Mahindra A: Novel therapies in the treatment of multiple myeloma. J Natl Compr Canc Netw. 2009, 7: 947-960.

    CAS  PubMed  Google Scholar 

  27. Kapoor P, Kumar S, Fonseca R: Impact of risk stratification on outcome among patients with multiple myeloma receiving initial therapy with lenalidomide and dexamethasone. Blood. 2009, 114: 518-521. 10.1182/blood-2009-01-202010

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Kumar SK, Uno H, Jacobus SJ: Impact of gene expression profiling-based risk stratification in patients with myeloma receiving initial therapy with lenalidomide and dexamethasone. Blood. 2011, 118: 4359-4362. 10.1182/blood-2011-03-342089

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Martinelli G, Terragna C, Zamagni E: Molecular remission after allogeneic or autologous transplantation of hematopoietic stem cells for multiple myeloma. J Clin Oncol. 2000, 18: 2273-2281.

    CAS  PubMed  Google Scholar 

  30. Martinez-Sanchez P, Montejano L, Sarasquete ME: Evaluation of minimal residual disease in multiple myeloma patients by fluorescent-polymerase chain reaction: the prognostic impact of achieving molecular response. Br J Haematol. 2008, 142: 766-774. 10.1111/j.1365-2141.2008.07263.x

    Article  CAS  PubMed  Google Scholar 

  31. Kroger N, Badbaran A, Lioznov M: Post-transplant immunotherapy with donor-lymphocyte infusion and novel agents to upgrade partial into complete and molecular remission in allografted patients with multiple myeloma. Exp Hematol. 2009, 37: 791-798. 10.1016/j.exphem.2009.03.008

    Article  PubMed  Google Scholar 

  32. Galimberti S, Benedetti E, Morabito F: Prognostic role of minimal residual disease in multiple myeloma patients after non-myeloablative allogeneic transplantation. Leuk Res. 2005, 29: 961-966. 10.1016/j.leukres.2005.01.017

    Article  CAS  PubMed  Google Scholar 

  33. Rimelen V, Ferrand C, Debecker A: Relative quantification of residual tumor cells by lightcycler real-time IgH PCR in autografted multiple myeloma patients. Leukemia. 2005, 19: 492-495. 10.1038/sj.leu.2403620

    Article  CAS  PubMed  Google Scholar 

  34. Greipp PR, San MJ, Durie BG: International staging system for multiple myeloma. J Clin Oncol. 2005, 23: 3412-3420. 10.1200/JCO.2005.04.242

    Article  PubMed  Google Scholar 

  35. Jacobson JL, Hussein MA, Barlogie B, Durie BG, Crowley JJ: A new staging system for multiple myeloma patients based on the Southwest Oncology Group (SWOG) experience. Br J Haematol. 2003, 122: 441-450. 10.1046/j.1365-2141.2003.04456.x

    Article  PubMed  Google Scholar 

  36. Fassas AB, Tricot G: Chromosome 13 deletion/hypodiploidy and prognosis in multiple myeloma patients. Leuk Lymphoma. 2004, 45: 1083-91. 10.1080/10428190310001642710

    Article  CAS  PubMed  Google Scholar 

  37. Zojer N, Konigsberg R, Ackermann J: 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.

    CAS  PubMed  Google Scholar 

  38. Keats JJ, Reiman T, Maxwell CA: In multiple myeloma, t(4;14)(p16;q32) is an adverse prognostic factor irrespective of FGFR3 expression. Blood. 2003, 101: 1520-1529. 10.1182/blood-2002-06-1675

    Article  CAS  PubMed  Google Scholar 

  39. Sarasquete ME, Garcia-Sanz R, Gonzalez D: Minimal residual disease monitoring in multiple myeloma: a comparison between allelic-specific oligonucleotide real-time quantitative polymerase chain reaction and flow cytometry. Haematologica. 2005, 90: 1365-1372.

    CAS  PubMed  Google Scholar 

  40. Taylor BJ, Kriangkum J, Strachan ER, Wizniak J, Pilarski LM: Identification of clonotypic IgH VDJ sequences in multiple myeloma. Methods Mol Med. 2005, 113: 121-144.

    CAS  PubMed  Google Scholar 

  41. Debes-Marun CS, Dewald GW, Bryant S: Chromosome abnormalities clustering and its implications for pathogenesis and prognosis in myeloma. Leukemia. 2003, 17: 427-36. 10.1038/sj.leu.2402797

    Article  CAS  PubMed  Google Scholar 

  42. Cremer FW, Kartal M, Hose D: High incidence and intraclonal heterogeneity of chromosome 11 aberrations in patients with newly diagnosed multiple myeloma detected by multiprobe interphase FISH. Cancer Genet Cytogenet. 2005, 161: 116-124. 10.1016/j.cancergencyto.2005.02.015

    Article  CAS  PubMed  Google Scholar 

  43. Lopez-Corral L, Gutierrez NC, Vidriales MB: The progression from MGUS to smoldering myeloma and eventually to multiple myeloma involves a clonal expansion of genetically abnormal plasma cells. Clin Cancer Res. 2011, 17: 1692-1700. 10.1158/1078-0432.CCR-10-1066

    Article  CAS  PubMed  Google Scholar 

  44. Drach J, Schuster J, Nowotny H: Multiple myeloma: high incidence of chromosomal aneuploidy as detected by interphase fluorescence in situ hybridization. Cancer Res. 1995, 55: 3854-3859.

    CAS  PubMed  Google Scholar 

  45. Alfaro R, Rosell J, Duran MA: Structural rearrangements of trisomies are a risk marker of clinical progression in hyperdiploid multiple myeloma. Anticancer Res. 2011, 31: 1599-1602.

    CAS  PubMed  Google Scholar 

Download references


This work was supported by a grant from the Canadian Institutes of Health and The Alberta Cancer Foundation. We are grateful for funds from the Myeloma Alberta Support Society. KJT received an AHFMR graduate studentship.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Linda M Pilarski.

Additional information

Competing interests

The authors have no conflicts of interest and no competing interests.

Authors' contributions

KJT designed the experiments, did the work, analyzed the data and wrote the paper. TR collected patient samples, guided the analysis and edited the paper. ARB and LMP directed the work and wrote the paper. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Thulien, K.J., Belch, A.R., Reiman, T. et al. In non-transplant patients with multiple myeloma, the pre-treatment level of clonotypic cells predicts event-free survival. Mol Cancer 11, 78 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: