Profound impact of sample processing delay on gene expression of multiple myeloma plasma cells

Tobias Meißner, Anja Seckinger, Kari Hemminki, Uta Bertsch, Asta Foersti, Mathias Haenel, Jan Duering, Hans Salwender, Hartmut Goldschmidt, Gareth J. Morgan, Dirk Hose, Niels Weinhold

Gene expression profiling (GEP) has significantly contributed to the elucidation of the molecular heterogeneity of multiple myeloma plasma cells (MMPC) and only recently it has been recommended for risk stratification. Prior to GEP MMPC need to be enriched resulting in an inability to immediately freeze bone marrow aspirates or use RNA stabilization reagents. As a result in multi-center MM trials sample processing delay due to shipping may be an important confounder of molecular analyses and risk stratification based on GEP data. We compared GEP data of 145 in-house and 246 shipped samples and detected 3301 down-regulated and 3501 up-regulated genes in shipped samples. For 3994 genes we confirmed differential expression in an independent set of 85 in-house and 97 shipped samples. Differentially expressed genes were enriched in processes like ribosome biogenesis, cell cycle, and apoptosis. Among GEP based risk predictors the IFM-15 seemed to underestimate high risk in shipped samples, whereas the GEP70 and the EMC-92 gene signatures were more robust. In order to provide a tool to assess the "shipping effect" in public repositories, we generated a 17-gene predictor for shipped samples with a 10-fold cross validation error rate of 0.06 for the training set and an error rate of 0.15 for the validation set. Sample processing delay significantly influences GEP of MMPC, implying it should be avoided if samples were used for risk stratification.