Alapi, Hilel

Abstract Background Red blood cell distribution width (RDW) is an easily accessible diagnostic test that is underutilized as a prognostic marker, despite substantial supporting evidence in cancer, cardiovascular diseases, and helathy individuals (Pilling et al. PLoS One 2028;13:e0203504). The association between increased RDW and incident acute myeloid leukemia (AML) has long been suggested in healthy volunteers (Abelson et al. Nature 2019;559:400). We looked into a similar phenomenon in patients with polycythemia vera (PV) or essential thrombocythemia (ET) in large discovery and validation cohorts. Methods A machine learning (ML) approach was used to identify baseline clinical and laboratory parameters that were predictive of leukemic transformation (LT) and develop a predictive model in a large discovery cohort of PV and ET, obtained from electronic medical records of the Maccabi Healthcare Services (MHS; Tel-Aviv, Israel). Receiver operating characteristic (ROC) curve analysis was utilized to determine the optimal cutoff points for RDW and other continuous variables. Mayo Clinic databases (Rochester, MN, USA) for patients with PV and ET were utilized as validation cohorts. Identification of study patients in the latter was according to the International Consensus Classification (Arber et al. Blood 2022; 140:1200) while in the discovery cohort was based on presence of ICD-9 codes associated with documentation of disease-specific driver mutations or treatment records that are consistent. Statistical analyses were performed using the R statistical software 4.4.2 (November, 2024) (Foundation for Statistical Computing, Vienna, Austria) or JMP Pro 18.0.0 software (SAS Institute, Cary, NC, USA). Results The discovery cohorts from Israel included 4,592 patients with PV (median age 63 yesrs; 53% males) and 5,968 with ET (median age 61 yesrs; males 40%). The validation cohorts from the Mayo Clinic included 633 PV (median 65 years; 51% males) and 634 ET (median 60 years; females 64%) patients that were informative for RDW. In the discovery cohorts, median follow-up (incident AML) was 6.1 years (N=186; 4.1%) in PV and 6.1 years (N=85; 1.4%) in ET. The corresponding values for the validation cohorts were 6.6 years (N=24; 3.8%) in PV and 7.0 years (N=22; 3.5%) in ET. Discovery cohort analyses Covariate-adjusted Cox regression analysis was utilized to estimate associations between LT and several clinical variables followed by machine learning pipeline to predict LT in a given timeframe. Conversion from continuous to binary variables was accompliched by ROC analysis. Shapley Additive exPlanations (SHAP) application enabled selection of top-ranked features that included RDW >15.5 (HR 3.2, 1.1-9.5; p=0.03), age >60 years (HR 2.9, 1.4-5.9; p<0.01), and LDH >1.5 x UNL (HR 4.0 (1.6-10.0; p<0.01). A similar analysis in ET also identified RDW >15.5 (p<0.01), age >60 years (p=0.01) as well as male sex (p=0.03) as independent risk factors for LT. Subseqent risk models based on these risk factors were externally validated and enabled robust risk stratification with 10-20 year AUCs rangening from 0.75 to 0.81. Validation cohorts analyses The indpendent predictive value of increased RDW observed in the discovery cohort was confirmed in the validation cohorts for both PV and ET and further gender-specified: female/male RDW ≥16.4/15.8 in PV (HR 7.9, 1.9-33.8) and ≥16.1/14.4 in ET (HR 5.1, 2.2-11.9). Also confirmed were the additional prognostic contributions of age ≥60 years (HR 3.8, 1.3-10.9) and LDH >1.5 x UNL (HR 3.9, 1.6-9.4), in PV, while no other clinical risk factor showed RDW-independent prognostic relevance in ET. The prognostic impact of RDW in ET was most apparent in females (p<0.01), as opposed to males (p=0.32) while the effect was equally significant in PV. In multivariable analysis that included previously recognized high risk mutations in PV (SRSF2) and ET (SF3B1), both RDW and mutations sustained significance (p<0.05). Conclusion The current study identifies higher RDW as a key risk factor for LT in both PV and ET, independent of other clinical or genetic risk factors. Possible underlying mechanisms include increased RDW as a marker of long-term pre-leukemic clones with sub-clinical ineffective erythropoiesis or an associated pro-leukemic inflammatory state.

ABSTRACT Survival prediction models in essential thrombocythemia (ET) include the International Prognostic Scoring System (IPSET) and the more recently introduced triple‐A (AAA) prognostic score. The latter enlists age and absolute neutrophil (ANC) and lymphocyte (ALC) counts as risk variables. In the current study, a Mayo Clinic discovery cohort of 658 patients with ET was used to identify AAA‐independent risk variables. Accordingly, multivariable analysis‐derived HRs (95% CI) were 15.7 (8.4–29.5) for age > 70 years (8 points); 4.2 (2.3–7.5) for age 50 to 70 years (2 points); 1.8 (1.2–2.5) for ANC ≥ 8 × 10 9 /L (1 point); 1.4 (1.03–1.9) for ALC < 1.7 × 10 9 /L (1 point); 1.8 (1.2–2.6) for absolute monocyte count (AMC) ≥ 0.5 × 10 9 /L (1 point); 1.8 (1.2–2.3) for male sex (1 point); 1.8 (1.3–2.4) for arterial hypertension (1 point); and 1.6 (1.2–2.3) for arterial thrombosis (1 point). HR‐weighted scoring enabled a 4‐tiered risk classification: ultra‐low (0–1 points; N = 94; median survival 42.7 years), low (2–4 points; N = 297; 23 years), intermediate (5 points; N = 66; 17.3 years), and high (6–14 points; N = 201; 10.8 years). Time‐dependent predictive performance at 20/25 years favored AAA + (AUC 0.92/0.91) vs. AAA (0.86/0.86) vs. IPSET (0.81/0.84). The AAA + risk model was subsequently validated by two external cohorts from Israel ( N = 5968) and Italy ( N = 682). In the cohort from Israel, disease‐specific mortality was assessed by comparing observed survival to an age‐ and sex‐matched reference population, which suggested near‐normal life expectancy in ultra‐low risk patients. The current study highlights host‐related factors as the primary determinants of longevity in ET and provides a composite risk score (AAA + ) that is based on complete blood count‐derived parameters and host‐related factors. Predictive performance of the new model was shown to be superior to that of IPSET and AAA.