Extranodal natural killer/T-cell lymphoma (ENKTL) is a highly aggressive non-Hodgkin lymphoma strongly associated with Epstein-Barr virus (EBV) infection and is characterized by a high relapse rate.1 Recently, the 5-year overall survival (OS) of limited-stage ENKTL has increased to 72-74% as the result of the introduction of a novel strategy of concurrent chemoradiotherapy.2 However, the 5-year OS in advanced-stage disease remains around 15-25%.3 Since risk-adapted therapy plays a pivotal role in improving the survival of newly diagnosed patients,4 a better risk classification model could assist in precisely stratifying patients with ENKTL into risk groups and formulating appropriate individualized treatments to improve prognosis.
Several risk scoring systems for newly diagnosed ENKTL patients are currently available in clinical practice, including the International Prognostic Index (IPI), the Korean Prognostic Index (KPI), the prognostic index for natural killer lymphoma with or without EBV-DNA (PINK/PINK-E), and the nomogram-revised risk index (NRI).5-8 Among these systems, the PINK-E is the only scoring system that includes data regarding EBV-DNA and is widely used in clinical practice. Several studies have demonstrated the good prognostic value of EBV-DNA and maximum standardized uptake value (SUVmax) in patients with ENKTL.9-12 However, no risk classification model has used SUVmax in 18F-fluorodeoxyglucose positron emission tomography/ computed tomography (18F-FDG PET/CT) and the quantitative value of EBV-DNA. In previous work, we integrated circulating tumor DNA (ctDNA) into the PINK-E to construct the PINK-EC model, which could overcome the poor discrimination efficiency of PINK-E for patients in low-risk and intermediate-risk groups.13 However, compared with the classic PINK-E system, our new PINK-EC model showed only slight improvement in terms of Harrell’s C-index for OS, which might not satisfy personalized clinical demand. Therefore, a more accurate and precise risk classification model is urgently needed. Here, we developed a nomogram model (we named it SEC, inspired by the initials of “SUVmax”, “EBV-DNA” and “ctDNA”), which included the semiquantitative radiomic parameter SUVmax and quantitative molecular parameters such as EBV-DNA and ctDNA to accurately stratify newly diagnosed ENKTL patients for optimal personalized treatment and management.
In this study, 91 patients newly diagnosed with ENKTL were enrolled at Xinqiao Hospital between February 2017 and November 2023 (ClinicalTrials identifier: ChiCTR1800014813).
The data collected at diagnosis, as previously described,13,14 included age, gender, B-symptoms, Eastern Cooperative Oncology Group performance status (ECOG PS), primary site, regional lymph node involvement, distant lymph node involvement, numbers of extranodal sites, serum lactate dehydrogenase (LDH) level, whole blood EBV-DNA copy number, 18FDG PET/CT SUVmax value and ctDNA concentration. In this study, we extended the follow-up time to the date of this analysis. The methodology of ctDNA measurements was reported previously.13 This study was approved by the China Ethics Committee of Registering Clinical Trials (ChiECRCT-20,180,005).
Based on the conventional cutoff value, continuous variables such as age and LDH concentration were divided into two categories. The ideal cutoff values for SUVmax, ctDNA, and EBV-DNA for survival prediction were determined by the area under the curve (AUC) of the receiver operating characteristic (ROC) curve. Progression-free survival (PFS) was calculated from diagnosis to disease progression, death from any cause, or the date of last follow-up. OS was measured from the date of diagnosis to the date of death due to any cause or the date of the last follow-up. Survival time was estimated using Kaplan-Meier survival curves and compared by log-rank tests. Univariate and multivariate analyses and calculations of hazard ratios (HR) with 95% confidence intervals (95% CI) were performed using Cox regression models. A nomogram was generated based on the independent predictors of survival outcomes determined by univariate and multivariate analyses. A calibration curve (1,000 bootstrap resamples) was constructed to assess the consistency between the predicted and observed survival. The discriminatory ability of the model was evaluated by Harrell’s C-index. A time-dependent ROC curve was used for the comparison of risk stratification for these prognostic models. All the statistical analyses were performed with IBM SPSS statistical software (version 25.0; IBM Inc., NY, USA) and R version 4.3.1 (R Foundation for Statistical Computing, Vienna, Austria). A two-sided P<0.05 was considered statistically significant.
The baseline characteristics of the whole cohort are listed in Online Supplementary Table S1. The optimal cutoff values determined by the ROC analysis for SUV-max, EBV-DNA, and ctDNA were 9.950, 1.4×104 copies/ mL, and 4.026 haploid genome equivalents per milliliter (hGE/mL), respectively (Figure 1).
According to the univariate and multivariate Cox regression analyses, the outcomes showed that only the SUVmax, EBV-DNA, and ctDNA concentration were independent prognostic factors for OS (Figure 2A). We did not include LDH in construction of the nomogram considering the inconsistent identification of LDH >245 U/L as a risk factor in univariate and multivariate Cox regression analyses. We created a prognostic nomogram, SEC, to predict the 3- and 5-year OS of ENKTL patients based on the SUVmax, ctDNA, and EBV-DNA level (Figure 2B). To clarify the predictive ability of the SEC score, calibration plots were generated, and the results demonstrated satisfactory consistency between the nomogram prediction for OS values and the actual observation (Figure 3A, B). Harrell’s C-index of the SEC for PFS and OS prediction were 0.771 (95% CI: 0.710-0.831) and 0.817 (95% CI: 0.768-0.866), respectively, which were better than those of the IPI, KPI, PINK, PINK-E and PINK-EC (Online Supplementary Table S2). These results suggest that the SEC is a more accurate and powerful tool for the prediction of PFS and OS in patients with ENKTL.
To visually display the stratification power of the SEC in ENKTL patients, Kaplan-Meier analysis was performed to evaluate survival outcomes. According to the nomogram, the largest contributing component was less than twice the size of the smallest and we, therefore, assigned equal weights, namely 1 point for each risk factor. Patients were stratified into three risk groups based on the SEC score: low risk (0), intermediate risk (1-2), and high risk (3). The 3-year PFS rates of patients in the low-, intermediate- and high-risk groups were 100% (95% CI: not applicable [NA]), 52.6% (95% CI: 39.4-70.3), and 4.2% (95% CI: 0.6-28.6), respectively, and the 3-year OS rates of patients in the three groups were 100% (95% CI: NA), 62.2% (95% CI: 48.6-79.6) and 4.5% (95% CI: 0.7-30.4), respectively. Furthermore, a pairwise comparison analysis showed that the SEC could discriminate the intermediate-risk group from the low-risk group (PFS: HR=NA; P=0.003; OS: HR=NA, P=0.015) and the high-risk group (PFS: HR=3.424, 95% CI: 1.718-6.824; P<0.001; OS: HR=4.965, 95% CI: 2.287-10.777; P<0.001) (Figure 3C, D).
To further evaluate the risk stratification power of the SEC with current prognostic scoring systems in ENKTL patients, the IPI, KPI, PINK, PINK-E, and PINK-EC also divided patients into three risk groups (Online Supplementary Table S3). Next, time-dependent ROC curves were plotted, and the corresponding AUC were calculated to compare the predictive accuracy of SEC with those of the IPI, KPI, PINK, PINK-E, and PINK-EC. The AUC of PFS and OS for SEC were greater and more stable than those of the other five prognostic models (Figure 3E, F).
The SEC nomogram is the first to use three semiquantitative or quantitative parameters with specific numerical cutoff values and could distinguish low-, intermediate-, and high-risk groups well. The good survival of low-risk patients classified by the SEC model proved that they benefited significantly from the available therapies. The prognosis of high-risk patients distinguished by the SEC model was extremely poor, with 3-year PFS and OS values of only 4.2% and 4.5%, respectively, which indicated that we need to develop more effective treatment options. Compared with IPI, KPI, PINK/PINK-E, or PINK-EC, the SEC nomogram had better prognostic prediction, risk stratification, and clinical use. We are currently conducting further validation of the model in a multicenter trial for the rare morbidity of ENKTL.
Footnotes
- Received February 26, 2024
- Accepted August 9, 2024
Correspondence
Disclosures
No conflicts of interest to disclose.
Funding
Acknowledgments
Our deepest gratitude goes to all the ENKTL patients enrolled in this study.
References
- Fox CP, Ahearne MJ, Pettengell R. Guidelines for the management of mature T- and natural killer-cell lymphomas (excluding cutaneous T-cell lymphoma): a British Society for Haematology guideline. Br J Haematol. 2022; 196(3):507-522. Google Scholar
- Wang X, Zhang L, Liu X. Efficacy and safety of a pegasparaginase-based chemotherapy regimen vs an L-asparaginase-based chemotherapy regimen for newly diagnosed advanced extranodal natural killer/T-cell lymphoma: a randomized clinical trial. JAMA Oncol. 2022; 8(7):1035-1041. Google Scholar
- Sun K-HM, Wong Y-TH, Cheung K-MC. Update on molecular diagnosis in extranodal NK/T-cell lymphoma and its role in the era of personalized medicine. Diagnostics (Basel). 2022; 12(2):409. Google Scholar
- Huang D, Ma N, Li X. Advances in single-cell RNA sequencing and its applications in cancer research. J Hematol Oncol. 2023; 16(1):98. Google Scholar
- Chim C-S, Ma S-Y, Au W-Y. Primary nasal natural killer cell lymphoma: long-term treatment outcome and relationship with the International Prognostic Index. Blood. 2004; 103(1):216-221. Google Scholar
- Lee J, Suh C, Park YH. Extranodal natural killer T-cell lymphoma, nasal-type: a prognostic model from a retrospective multicenter study. J Clin Oncol. 2006; 24(4):612-618. Google Scholar
- Kim SJ, Yoon DH, Jaccard A. A prognostic index for natural killer cell lymphoma after non-anthracycline-based treatment: a multicentre, retrospective analysis. Lancet Oncol. 2016; 17(3):389-400. Google Scholar
- Chen S-Y, Yang Y, Qi S-N. Validation of nomogram-revised risk index and comparison with other models for extranodal nasal-type NK/T-cell lymphoma in the modern chemotherapy era: indication for prognostication and clinical decision-making. Leukemia. 2021; 35(1):130-142. Google Scholar
- Suzuki R, Yamaguchi M, Izutsu K. Prospective measurement of Epstein-Barr virus-DNA in plasma and peripheral blood mononuclear cells of extranodal NK/T-cell lymphoma, nasal type. Blood. 2011; 118(23):6018-6022. Google Scholar
- Kwong Y-L, Pang AWK, Leung AYH, Chim C-S, Tse E. Quantification of circulating Epstein-Barr virus DNA in NK/T-cell lymphoma treated with the SMILE protocol: diagnostic and prognostic significance. Leukemia. 2014; 28(4):865-870. Google Scholar
- Pak K, Kim BS, Kim K. Prognostic significance of standardized uptake value on F18-FDG PET/CT in patients with extranodal nasal type NK/T cell lymphoma: a multicenter, retrospective analysis. Am J Otolaryngol. 2018; 39(1):1-5. Google Scholar
- Qin C, Yang S, Sun X, Xia X, Li C, Lan X. 18F-FDG PET/CT for prognostic stratification of patients with extranodal natural killer/T-cell lymphoma. Clin Nucl Med. 2019; 44(3):201-208. Google Scholar
- Huang D, Li Q, Li X. The novel prognostic index model of combining circulating tumor DNA and PINK-E predicts the clinical outcomes for newly diagnosed extranodal NK/T-cell lymphoma. Hemasphere. 2023; 7(1):e822. Google Scholar
- Li Q, Zhang W, Li J. Plasma circulating tumor DNA assessment reveals KMT2D as a potential poor prognostic factor in extranodal NK/T-cell lymphoma. Biomark Res. 2020; 8:27. Google Scholar
Data Supplements
Figures & Tables
Article Information
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.