Mathematical modelling of stem and progenitor cell dynamics during ruxolitinib treatment of patients with myeloproliferative neoplasms.

Introduction: The Philadelphia chromosome-negative myeloproliferative neoplasms are a group of slowly progressing haematological malignancies primarily characterised by an overproduction of myeloid blood cells. Patients are treated with various drugs, including the JAK1/2 inhibitor ruxolitinib. Mathematical modelling can help propose and test hypotheses of how the treatment works. Materials and methods: We present an extension of the Cancitis model, which describes the development of myeloproliferative neoplasms and their interactions with inflammation, that explicitly models progenitor cells and can account for treatment with ruxolitinib through effects on the malignant stem cell response to cytokine signalling and the death rate of malignant progenitor cells. The model has been fitted to individual patients' data for the JAK2 V617F variant allele frequency from the COMFORT-II and RESPONSE studies for patients who had substantial reductions (20 percentage points or 90% of the baseline value) in their JAK2 V617F variant allele frequency (n = 24 in total). Results: The model fits very well to the patient data with an average root mean square error of 0.0249 (2.49%) when allowing ruxolitinib treatment to affect both malignant stem and progenitor cells. This average root mean square error is much lower than if allowing ruxolitinib treatment to affect only malignant stem or only malignant progenitor cells (average root mean square errors of 0.138 (13.8%) and 0.0874 (8.74%), respectively). Discussion: Systematic simulation studies and fitting of the model to the patient data suggest that an initial reduction of the malignant cell burden followed by a monotonic increase can be recapitulated by the model assuming that ruxolitinib affects only the death rate of malignant progenitor cells. For patients exhibiting a long-term reduction of the malignant cells, the model predicts that ruxolitinib also affects stem cell parameters, such as the malignant stem cells' response to cytokine signalling.

Case Report: First longitudinal study of a patient with CALR positive clonal hematopoiesis of indeterminate potential developing into pre-fibrotic myelofibrosis.

Initial diagnosis of overt myeloproliferative neoplasms (MPNs) represents the juncture during clonal evolution when symptoms or complications prompt an afflicted individual to seek medical attention. In 30-40% of the MPN subgroups essential thrombocythemia (ET) and myelofibrosis (MF), somatic mutations in the calreticulin gene (CALR) are drivers of the disease resulting in constitutive activation of the thrombopoietin receptor (MPL). In the current study, we describe a healthy CALR mutated individual during a 12 year follow-up from initial identification of CALR clonal hematopoiesis of indeterminate potential (CHIP) to the diagnosis of pre-MF. The pre-diagnostic exponential development dynamics of the malignant clone demonstrated close correlation with the platelet counts, neutrophil-to-lymphocyte (NLR) ratio, and inversely correlated to hemoglobin and erythrocyte counts. Backward extrapolation of the growth rate indicated the potential for discovery of the malignant clone many years prior to presentation of overt disease, opening a window of opportunity for early treatment intervention. We did not find any additional mutations associated with MPNs and the current case report provides novel information regarding the development of a driver mutation and the association with blood cell counts prior to clinical manifestation of symptoms suggesting that pre-diagnostic dynamics may supplement future diagnostic criteria for early diagnosis and intervention in MPN patients.