Clinical features and survival outcomes in IgD myeloma: a study by Asia Myeloma Network (AMN)

Immunoglobulin D (IgD) myeloma is a rare isotype that comprises 1–2% of multiple myeloma (MM) patients [1–3], which has significantly inferior survival for a median overall survival (OS) between 13 and 21 months [4–6]. Given the lack of large cohort with comprehensive clinical and cytogenetic assessment, knowledge about IgD myeloma is obtained mostly from a limited sample size [7]. Therefore, we carried out a multicenter retrospective study to evaluate the prevalence, clinical features, prognosis, and to develop and validate a prognostic model, including 356 patients with IgD myeloma from 14 centers of Asian Myeloma Network (AMN). Data were collected from China, Korea, and Singapore diagnosed from 2002 to 2019 (Supplementary Table 1). Ethical committee approvals were obtained and study protocol was approved by the Institutional Review Board of each institution. To avoid clinical information leak, and get These authors contributed equally: Jin Liu, Xiaoxia Hu, Yanchun Jia

a real sense of the accurate model's outcomes, we split existing 356 IgD MM to three parts, namely training cohort (one center from Shanghai, n = 212), validation cohort 1 (two centers from Beijing, n = 81), and validation cohort 2 (centers from Korea and Singapore, n = 63). The Least Absolute Shrinkage and Selector Operation (LASSO) Cox regression model to determine prognostic factors from the variables with P < 0.05 in the log-rank tests was performed as described [8,9]. The quality of the prediction model was measured using the concordance index (C-index) and areas under the time-dependent receiver-operating characteristics (ROC) curves (AUCs). A bootstrap with 1000 re-samples was used for internal validation. SAS 9.4 and R 3.5.1 were used for the statistical analysis. A total of 356 patients with IgD myeloma represented 2-8.8% of all myeloma patients, especially over 5% IgD myeloma prevalence in Chinese centers. We compared the clinical characteristics of IgD myeloma with 712 (1:2) non-IgD myeloma patients random selected as control matched for year of diagnosis and systemic therapy from Shanghai Changzheng Hospital. Baseline characteristics of total cohort are listed on Table 1 and different centers are shown in Supplementary Table 2. IgD myeloma patients had a higher frequency in male, younger than 65 years, advanced R-ISS stage III, hypercalcemia, elevated creatinine levels, and elevated LDH. Cytogenetic information was available for 301 patients (84.6%), while the 1q21 probe was only performed in 75.8% patients. Notably, 29.2% frequency of t(11;14) was predominantly higher compared to those in non-IgD subtypes (P < 0.001). Among the 88 IgD patients harbored t (11;14), the most frequent chromosome abnormalities (CA) coupled with t(11;14) were 13q-(31.8%), 1q21 + (30.7%), and followed 17p-(11.4%). 'Double-hit' or 'triple-hit' [10,11] only occurred 5.6% and 0.3% patients, respectively.
And then, we compared IgD myeloma patients with IgG, IgA, and light chain patients random selected as matched control (Supplementary Table 3). The median age was younger in IgD compared with others myeloma subtypes. Notably, the frequency of t(11;14) was significantly higher than non-IgD subtypes (IgD 29.2% vs IgG 10.6% vs IgA 8.4%, P < 0.001), but was a slight higher than light chain subtype (29.2 vs 24.9%). 'Double-hit' phenotype was significant lower in IgD myelomas than others subtypes.
Frontline treatment modalities used are shown in Supplementary Table 4. The overall response rate (ORR) was 88.8%, and very good partial response or better was 58.6% (Supplementary Table 5). After a median follow-up of 8.2, 7.3, and 4.9 years for the three cohorts, the median OS were 36.5 months for the total cohort and 31.2 months in training cohort, 52.2 months in validation cohort 1, and 45.7 months in validation cohort 2 (Supplementary Fig. 1 and Supplementary Table 6). Patients received IMiDs showed a relatively longer median OS than others regimens, however, which untranslated into a significant survival benefit (P = 0.17, Supplementary Fig. 2a), and might be the subgroups limitation. Patients received ASCT had a median OS of 45.7 months, which was a slightly longer than 35 months for non-ASCT patients (P = 0.4, Supplementary Fig. 2b). We subsequently investigated whether cytogenetic aberration was a prognostic factor [12], which showed that CA did not have an impact on OS, suggesting other molecular events overcome initial CA risk features and impact prognosis.
Subsequently, the LASSO Cox regression model to determine prognostic factors from the univariate analysis was performed (Supplementary Table 7). Five clinical parameters with statistically relevant, including lambda light chain, plasma cells in BM ≥ 50%, hemoglobin < 100 g/L, LDH ≥ 245 U/L, and extramedullary plasmacytoma, were integrated into multivariate LASSO regression model (Supplementary Table 8). A nomogram was developed and the risk score was computed as follows: 0.9215 × lambda light chain + 0.6376 × plasma cells in BM (≥50%) + 0.5203 × anemia (<100 g/L) + 0.6864 × LDH (≥245 U/L) + 0.4484 × extramedullary plasmacytoma (variable present = 1, absent = 0, Fig. 1a, b). The predictive accuracy for OS calculated using the C-index was 0.705 (95% CI, 0.663-0.747). In the internal  Fig. 1c-e suggested a good fit for the observed nomogram, when compared with the ideal nomogram. The panel displayed an AUC value at 1-year, 3-year, and 5-year OS, and the validation sets had a similar high AUC values at these timepoints ( Fig. 1f-h).
On the basis of the distribution of the risk scores and the 3year survival probability, two categories of risk were created with the cut-off point at 1.56: standard risk (risk score ≤ 1.56, n = 156) and high-risk subgroup (risk score > 1.56, n = 200). The clinical characteristics between derivations were presented in Supplementary Table 6. Patients with IgD myeloma at standard risk were significantly better than high-risk subgroup (Fig. 1i). Similar results were obtained in training and validation cohorts respectively (Fig. 1j-l). Notably we identified that ASCT showed a survival advantage in highrisk group, while it did not improve the survival in standard risk group (Supplementary Fig. 3). Moreover the prediction value of the model was independent of induction modalities ( Supplementary Fig. 4).
The clinical features of IgD myeloma in this study indicated a potentially high-risk population. We identified a higher prevalence difference may be related to the difference in ethnicity, an efficient algorithm diagnosis, and assembled patients. We found 40.6% of patients had abnormal karyotypes. Translocation (11;14) was predominantly high frequency, which showed the similar to others studies [13,14]. By far, the most clinically relevant feature of myeloma with t (11;14) is increased expression of the anti-apoptotic protein BCL-2 [15], and venetoclax of Bcl-2 inhibitor provides a new option for t (11;14), which may have a special significances for IgD patients with t (11;14).
The median OS with 36.5 months in IgD myeloma showed an inferior survival than patients with more common myeloma subtypes [16,17]. Although the ORR (88.8%) has not improved rapidly in this AMN study, patients received IMiD with lenalidomide showed a better trend median OS than patients received PIs, conventional chemotherapy regimens, even PI + thalidomide as PI + IMiD subgroup, while PI + IMiD overcame improvement due to high peripheral neuropathy to limit the application. However, the patients received IMiD is too small to draw the definite conclusion, which might provide a clue to optimize treatment. Meanwhile, we demonstrated ASCT would benefit for IgD myeloma patients with high-risk score.
Previous studies concentrated on the prognostic value of single variable [1,3,5,18,19], as it is a challenge to construct a comprehensive prognostic model due to the rarity of IgD myeloma. We used 356 IgD myeloma patients to develop a prognostic model based on the 17 variables of the multivariate LASSO model to estimate the 3-year OS. Five clinical parameters were identified as clinically relevant to compute the risk score, which a high-risk score was consistently associated with a poor survival outcome. Moreover, the prediction value of the model was independent of induction modalities. This risk model improves the classification of IgD myeloma and may facilitate the development of risk-adapted treatment strategies. Although this retrospective design, presence of missing data, and heterogeneous treatment could have resulted in the introduction of selection bias and underestimate or overestimate of some of results, we believe that the large sample size and the multicenter study performed in such a rare disease could balance out some of these weaknesses.
In conclusion, we described the clinical features of 356 IgD myeloma patients, as well as developed and validated a predictive model containing five baseline clinical variables that could group the IgD patients into standard risk and high risk. Meanwhile, we demonstrated that IMiDs therapy might be a trend to benefit for the patient's outcome, and ASCT could benefit patient with high risk within the predictive model. These findings may provide guidance for management of IgD myeloma and better prognostic stratification for development of risk-adapted treatment strategies.
Author contributions JL, XH and YJ collected and analyzed the data, and wrote the first draft, and approved the final version of the paper; JL, CS, MKK, FZ, and WC performed patient management and approved the final version of the paper; YJ, AL, LY and ZY, performed patients' follow-up, participated in final data analysis and approval of the final version of the paper; QC analyzed the data and approval of the final version of the paper; JHL, KK, and WJC, performed patient management revised the paper critically, and approved the final version of the paper; SKK and BD revised the paper critically, and approved the final version of the paper; WF, and JH designed the study, performed patient management, and approved the final version of the paper, and; JD designed the study, performed patient management, analyzed the data, wrote the first draft, approved the final version of the paper. Fig. 1 Development of a predictive model and validation. LASSO Cox regression model to determine prognostic factors from the variables with P < 0.05 in the log-rank tests was performed. The prediction model was established on the basis of variables selected from the LASSO Cox regression model and weighted using the Cox regression coefficient. Cross-validation for tuning parameter selection in the LASSO model (a). The nomogram based on data from training patients to predict individual prognosis (b). The calibration curves of an alternative nomogram to predict 3-year OS of IgD myeloma in training patients (c), validation cohort 1 (d) and validation cohort 2 (e). The X-axis represents the predicted survival probability calculated using the nomogram, while the Y-axis represents the actual survival probability for patients. The gold 45-degree line represents the ideal nomogram, while the black line represents the observed nomogram.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/.