Quiescence Modulates Stem Cell Maintenance and Regenerative Capacity in the Aging Brain.

The function of somatic stem cells declines with age. Understanding the molecular underpinnings of this decline is key to counteract age-related disease. Here, we report a dramatic drop in the neural stem cells (NSCs) number in the aging murine brain. We find that this smaller stem cell reservoir is protected from full depletion by an increase in quiescence that makes old NSCs more resistant to regenerate the injured brain. Once activated, however, young and old NSCs show similar proliferation and differentiation capacity. Single-cell transcriptomics of NSCs indicate that aging changes NSCs minimally. In the aging brain, niche-derived inflammatory signals and the Wnt antagonist sFRP5 induce quiescence. Indeed, intervention to neutralize them increases activation of old NSCs during homeostasis and following injury. Our study identifies quiescence as a key feature of old NSCs imposed by the niche and uncovers ways to activate NSCs to repair the aging brain.

HSC Niche Dynamics in Regeneration, Pre-malignancy, and Cancer: Insights From Mathematical Modeling.

The hematopoietic stem cell (HSC) niche is a crucial driver of regeneration and malignancy. Its interaction with hematopoietic and malignant stem cells is highly complex and direct experimental observations are challenging. We here develop a mathematical model which helps relate processes in the niche to measurable changes of stem and non-stem cell counts. HSC attached to the niche are assumed to be quiescent. After detachment HSC become activated and divide or differentiate. To maintain their stemness, the progeny originating from division must reattach to the niche. We use mouse data from literature to parametrize the model. By combining mathematical analysis and computer simulations, we systematically investigate the impact of stem cell proliferation, differentiation, niche attachment, and detachment on clinically relevant scenarios. These include bone marrow transplantation, clonal competition, and eradication of malignant cells. According to our model, sampling of blood or bulk marrow provides only limited information about cellular interactions in the niche and the clonal composition of the stem cell population. Furthermore, we investigate how interference with processes in the stem cell niche could help to increase the effect of low-dose chemotherapy or to improve the homing of genetically engineered cells.

Clonal haematopoiesis of indeterminate potential and impaired kidney function-A Danish general population study with 11 years follow-up.

The myeloproliferative neoplasms are associated with chronic kidney disease but whether clonal haematopoiesis of indeterminate potential (CHIP) is associated with impaired kidney function is unknown. In the Danish General Suburban Population Study (N = 19 958) from 2010 to 2013, 645 individuals were positive for JAK2V617F (N = 613) or CALR (N = 32) mutations. Mutation-positive individuals without haematological malignancy were defined as having CHIP (N = 629). We used multiple and inverse probability weighted (IPW)-adjusted linear regression analysis to estimate adjusted mean (95% confidence interval) differences in estimated glomerular filtration rate (eGFR; ml/min/1.73 m2 ) by mutation status, variant allele frequency (VAF%), blood cell counts, and neutrophil-to-lymphocyte ratio (NLR). We performed 11-year longitudinal follow-up of eGFR in all individuals. Compared to CHIP-negative individuals, the mean differences in eGFR were -5.6 (-10.3, -0.8, p = .02) for CALR, -11.9 (-21.4, -2.4, p = 0.01) for CALR type 2, and -10.1 (-18.1, -2.2, p = .01) for CALR with VAF ≥ 1%. The IPW-adjusted linear regression analyses showed similar results. NLR was negatively associated with eGFR. Individuals with CALR type 2 had a worse 11-year longitudinal follow-up on eGFR compared to CHIP-negative individuals (p = .004). In conclusion, individuals with CALR mutations, especially CALR type 2, had impaired kidney function compared to CHIP-negative individuals as measured by a lower eGFR at baseline and during 11-year follow-up.

Mathematical modelling of the hematopoietic stem cell-niche system: Clonal dominance based on stem cell fitness.

Human blood cell production is maintained by hematopoietic stem cells (HSC) which give rise to all types of mature blood cells. Experimental observation of HSC in their physiologic bone-marrow microenvironment, the so-called stem cell niche, is challenging. Therefore, the details of HSC dynamics and the cellular interactions in the stem cell niche remain elusive. Mutations that lead to a competitive advantage are the cause of clinical challenges when treating HSC-derived malignancies such as acute myeloid leukemia or the myeloproliferative neoplasms (MPNs). To investigate the significance of the interaction between the HSC and the stem cell niche in these malignancies, we propose and analyse a mechanism-based mathematical model of HSC dynamics within the bone-marrow microenvironment. The model is based on the central hypothesis that HSC self-renewal depends on the niche. In the model, the interaction of HSC with specific niches located in the bone marrow are key to the indefinite HSC renewal necessary for long-term maintenance of blood cell production. We formulate a general model of n distinct clones that differ with respect to cell properties. We identify an attractive trapping region and compute and classify all steady states. A concept of HSC fitness naturally arises from the model analysis. HSC fitness is found to determine the asymptotic behaviour of the model, as the HSC clone with the highest fitness is related to the unique locally stable steady state. Based on biological assumptions about HSC, we propose two reduced models of different complexity. A thorough mathematical analysis reveals that both reduced models have the same asymptotic behaviour as the full model. We compare the simpler of the two models, a logistic equation of the disease burden, to clinical data of MPN-patients. The reduced model is found to agree well with data and suggests a simple interpretation and possible prediction of patient prognosis. The proposed mathematical model and the reduced forms have the potential to provide insights into the regulation of HSC dynamics and blood cell formation, and ultimately for future advances in treatment of hematologic malignancies.

Dynamical properties of feedback signalling in B lymphopoiesis: A mathematical modelling approach.

Haematopoiesis is the process of generation of blood cells. Lymphopoiesis generates lymphocytes, the cells in charge of the adaptive immune response. Disruptions of this process are associated with diseases like leukaemia, which is especially incident in children. The characteristics of self-regulation of this process make them suitable for a mathematical study. In this paper we develop mathematical models of lymphopoiesis using currently available data. We do this by drawing inspiration from existing structured models of cell lineage development and integrating them with paediatric bone marrow data, with special focus on regulatory mechanisms. A formal analysis of the models is carried out, giving steady states and their stability conditions. We use this analysis to obtain biologically relevant regions of the parameter space and to understand the dynamical behaviour of B-cell renovation. Finally, we use numerical simulations to obtain further insight into the influence of proliferation and maturation rates on the reconstitution of the cells in the B line. We conclude that a model including feedback regulation of cell proliferation represents a biologically plausible depiction for B-cell reconstitution in bone marrow. Research into haematological disorders could benefit from a precise dynamical description of B lymphopoiesis.

Mathematical modeling of plant cell fate transitions controlled by hormonal signals.

Coordination of fate transition and cell division is crucial to maintain the plant architecture and to achieve efficient production of plant organs. In this paper, we analysed the stem cell dynamics at the shoot apical meristem (SAM) that is one of the plant stem cells locations. We designed a mathematical model to elucidate the impact of hormonal signaling on the fate transition rates between different zones corresponding to slowly dividing stem cells and fast dividing transit amplifying cells. The model is based on a simplified two-dimensional disc geometry of the SAM and accounts for a continuous displacement towards the periphery of cells produced in the central zone. Coupling growth and hormonal signaling results in a nonlinear system of reaction-diffusion equations on a growing domain with the growth rate depending on the model components. The model is tested by simulating perturbations in the level of key transcription factors that maintain SAM homeostasis. The model provides new insights on how the transcription factor HECATE is integrated in the regulatory network that governs stem cell differentiation.

New targeted approaches for epigenetic age predictions.

BACKGROUND: Age-associated DNA methylation changes provide a promising biomarker for the aging process. While genome-wide DNA methylation profiles enable robust age-predictors by integration of many age-associated CG dinucleotides (CpGs), there are various alternative approaches for targeted measurements at specific CpGs that better support standardized and cost-effective high-throughput analysis. RESULTS: In this study, we utilized 4647 Illumina BeadChip profiles of blood to select CpG sites that facilitate reliable age-predictions based on pyrosequencing. We demonstrate that the precision of DNA methylation measurements can be further increased with droplet digital PCR (ddPCR). In comparison, bisulfite barcoded amplicon sequencing (BBA-seq) gave slightly lower correlation between chronological age and DNA methylation at individual CpGs, while the age-predictions were overall relatively accurate. Furthermore, BBA-seq data revealed that the correlation of methylation levels with age at neighboring CpG sites follows a bell-shaped curve, often associated with a CTCF binding site. We demonstrate that within individual BBA-seq reads the DNA methylation at neighboring CpGs is not coherently modified, but reveals a stochastic pattern. Based on this, we have developed a new approach for epigenetic age predictions based on the binary sequel of methylated and non-methylated sites in individual reads, which reflects heterogeneity in epigenetic aging within a sample. CONCLUSION: Targeted DNA methylation analysis at few age-associated CpGs by pyrosequencing, BBA-seq, and particularly ddPCR enables high precision of epigenetic age-predictions. Furthermore, we demonstrate that the stochastic evolution of age-associated DNA methylation patterns in BBA-seq data enables epigenetic clocks for individual DNA strands.

Oscillations in a white blood cell production model with multiple differentiation stages.

In this work we prove occurrence of a super-critical Hopf bifurcation in a model of white blood cell formation structured by three maturation stages. We provide an explicit analytical expression for the bifurcation point depending on model parameters. The Hopf bifurcation is a unique feature of the multi-compartment structure as it does not exist in the corresponding two-compartment model. It appears for a parameter set different from the parameters identified for healthy hematopoiesis and requires changes in at least two cell properties. Model analysis allows identifying a range of biologically plausible parameter sets that can explain persistent oscillations of white blood cell counts observed in some hematopoietic diseases. Relating the identified parameter sets to recent experimental and clinical findings provides insights into the pathological mechanisms leading to oscillating blood cell counts.

Using mathematical models to improve risk-scoring in acute myeloid leukemia.

Acute myeloid leukemia (AML) is an aggressive cancer of the blood forming (hematopoietic) system. Due to the high patient variability of disease dynamics, risk-scoring is an important part of its clinical management. AML is characterized by impaired blood cell formation and the accumulation of so-called leukemic blasts in the bone marrow of patients. Recently, it has been proposed to use counts of blood-producing (hematopoietic) stem cells (HSCs) as a biomarker for patient prognosis. In this work, we use a non-linear mathematical model to provide mechanistic evidence for the suitability of HSC counts as a prognostic marker. Using model analysis and computer simulations, we compare different risk-scores involving HSC quantification. We propose and validate a simple approach to improve risk prediction based on HSC and blast counts measured at the time of diagnosis.

A structured population model of clonal selection in acute leukemias with multiple maturation stages.

Recent progress in genetic techniques has shed light on the complex co-evolution of malignant cell clones in leukemias. However, several aspects of clonal selection still remain unclear. In this paper, we present a multi-compartmental continuously structured population model of selection dynamics in acute leukemias, which consists of a system of coupled integro-differential equations. Our model can be analysed in a more efficient way than classical models formulated in terms of ordinary differential equations. Exploiting the analytical tractability of this model, we investigate how clonal selection is shaped by the self-renewal fraction and the proliferation rate of leukemic cells at different maturation stages. We integrate analytical results with numerical solutions of a calibrated version of the model based on real patient data. In summary, our mathematical results formalise the biological notion that clonal selection is driven by the self-renewal fraction of leukemic stem cells and the clones that possess the highest value of this parameter are ultimately selected. Moreover, we demonstrate that the self-renewal fraction and the proliferation rate of non-stem cells do not have a substantial impact on clonal selection. Taken together, our results indicate that interclonal variability in the self-renewal fraction of leukemic stem cells provides the necessary substrate for clonal selection to act upon.

Reduced hematopoietic stem cell frequency predicts outcome in acute myeloid leukemia.

In patients with acute myeloid leukemia and low percentages of aldehyde-dehydrogenase-positive cells, non-leukemic hematopoietic stem cells can be separated from leukemic cells. By relating hematopoietic stem cell frequencies to outcome we detected poor overall- and disease-free survival of patients with low hematopoietic stem cell frequencies. Serial analysis of matched diagnostic and follow-up samples further demonstrated that hematopoietic stem cells increased after chemotherapy in patients who achieved durable remissions. However, in patients who eventually relapsed, hematopoietic stem cell numbers decreased dramatically at the time of molecular relapse demonstrating that hematopoietic stem cell levels represent an indirect marker of minimal residual disease, which heralds leukemic relapse. Upon transplantation in immune-deficient mice cases with low percentages of hematopoietic stem cells of our cohort gave rise to leukemic or no engraftment, whereas cases with normal hematopoietic stem cell levels mostly resulted in multi-lineage engraftment. Based on our experimental data, we propose that leukemic stem cells have increased niche affinity in cases with low percentages of hematopoietic stem cells. To validate this hypothesis, we developed new mathematical models describing the dynamics of healthy and leukemic cells under different regulatory scenarios. These models suggest that the mechanism leading to decreases in hematopoietic stem cell frequencies before leukemic relapse must be based on expansion of leukemic stem cells with high niche affinity and the ability to dislodge hematopoietic stem cells. Thus, our data suggest that decreasing numbers of hematopoietic stem cells indicate leukemic stem cell persistence and the emergence of leukemic relapse.

Feedback signals in myelodysplastic syndromes: increased self-renewal of the malignant clone suppresses normal hematopoiesis.

Myelodysplastic syndromes (MDS) are triggered by an aberrant hematopoietic stem cell (HSC). It is, however, unclear how this clone interferes with physiologic blood formation. In this study, we followed the hypothesis that the MDS clone impinges on feedback signals for self-renewal and differentiation and thereby suppresses normal hematopoiesis. Based on the theory that the MDS clone affects feedback signals for self-renewal and differentiation and hence suppresses normal hematopoiesis, we have developed a mathematical model to simulate different modifications in MDS-initiating cells and systemic feedback signals during disease development. These simulations revealed that the disease initiating cells must have higher self-renewal rates than normal HSCs to outcompete normal hematopoiesis. We assumed that self-renewal is the default pathway of stem and progenitor cells which is down-regulated by an increasing number of primitive cells in the bone marrow niche--including the premature MDS cells. Furthermore, the proliferative signal is up-regulated by cytopenia. Overall, our model is compatible with clinically observed MDS development, even though a single mutation scenario is unlikely for real disease progression which is usually associated with complex clonal hierarchy. For experimental validation of systemic feedback signals, we analyzed the impact of MDS patient derived serum on hematopoietic progenitor cells in vitro: in fact, MDS serum slightly increased proliferation, whereas maintenance of primitive phenotype was reduced. However, MDS serum did not significantly affect colony forming unit (CFU) frequencies indicating that regulation of self-renewal may involve local signals from the niche. Taken together, we suggest that initial mutations in MDS particularly favor aberrant high self-renewal rates. Accumulation of primitive MDS cells in the bone marrow then interferes with feedback signals for normal hematopoiesis--which then results in cytopenia.

Lung-targeted RNA interference against angiopoietin-2 ameliorates multiple organ dysfunction and death in sepsis.

OBJECTIVE: Angiopoietin-2, a protein secreted by stimulated endothelium and an antagonist of the endothelium-stabilizing receptor Tie2, contributes to the pathophysiology of septic multiple organ dysfunction. We tested the therapeutic potential of a pulmonary-endothelium-specific RNA interference-based angiopoietin-2 targeting strategy in sepsis. DESIGN: Laboratory and animal research. SETTINGS: Research laboratories of the Medical School Hannover, Department of Nephrology and Hypertension, Hannover and Silence Therapeutics GmbH, Berlin. SUBJECTS: C57Bl/6 mice. INTERVENTIONS: Lung-endothelium-specific angiopoietin-2 small interfering RNA was administered both before and after sepsis induction (cecal ligation and puncture or lipopolysaccharides) intravenously. MEASUREMENTS AND MAIN RESULTS: Angiopoietin-2 small interfering RNA was highly specific and reduced angiopoietin-2 expression in the septic murine lungs up to 73.8% (p = 0.01) and enhanced the phosphorylation of Tie2 both in control and septic animals. Angiopoietin-2 small interfering RNA reduced pulmonary interleukin-6 transcription, intercellular adhesion molecule expression, neutrophil infiltration, and vascular leakage. Manifestations of sepsis were also attenuated in distant organs, including the kidney, where renal function was improved without affecting local angiopoietin-2 production. Finally, angiopoietin-2 small interfering RNA ameliorated the severity of illness and improved survival in cecal ligation and puncture, both as a pretreatment and as a rescue intervention. CONCLUSION: The Tie2 antagonist angiopoietin-2 represents a promising target against sepsis-associated multiple organ dysfunction. A novel RNA interference therapeutic approach targeting gene expression in the pulmonary endothelium could be a clinically relevant pharmacological strategy to reduce injurious angiopoietin-2 synthesis.

Assessing hematopoietic (stem-) cell behavior during regenerative pressure.

Hematopoiesis is a complex and strongly regulated process. In case of regenerative pressure, efficient recovery of blood cell counts is crucial for survival of an individual. We propose a quantitative mathematical model of white blood cell formation based on the following cell parameters: (1) proliferation rate, (2) self-renewal, and (3) cell death. Simulating this model we assess the change of these parameters under regenerative pressure. The proposed model allows to quantitatively describe the impact of these cell parameters on engraftment time after stem cell transplantation. Results indicate that enhanced self-renewal during the posttransplant period is crucial for efficient regeneration of blood cell counts while constant or reduced self-renewal leads to delayed recovery or graft failure. Increased cell death in the posttransplant period has a similar impact. In contrast, reduced proliferation or pre-homing cell death causes only mild delays in blood cell recovery which can be compensated sufficiently by increasing the dose of transplanted cells.

Modelling post-chemotherapy stem cell dynamics in the bone marrow niche of AML patients.

Acute myeloid leukemia (AML) is a stem cell-driven malignancy of the blood forming (hematopoietic) system. Despite of high dose chemotherapy with toxic side effects, many patients eventually relapse. The "7+3 regimen", which consists of 7 days of cytarabine in combination with daunorubicin during the first 3 days, is a widely used therapy protocol. Since peripheral blood cells are easily accessible to longitudinal sampling, significant research efforts have been undertaken to characterize and reduce adverse effects on circulating blood cells. However, much less is known about the impact of the 7+3 regimen on human hematopoietic stem cells and their physiological micro-environments, the so-called stem cell niches. One reason for this is the technical inability to observe human stem cells in vivo and the discomfort related to bone marrow biopsies. To better understand the treatment effects on human stem cells, we consider a mechanistic mathematical model of the stem cell niche before, during and after chemotherapy. The model accounts for different maturation stages of leukemic and hematopoietic cells and considers key processes such as cell proliferation, self-renewal, differentiation and therapy-induced cell death. In the model, hematopoietic (HSCs) and leukemic stem cells (LSCs) compete for a joint niche and respond to both systemic and niche-derived signals. We relate the model to clinical trial data from literature which longitudinally quantifies the counts of hematopoietic stem like (CD34+CD38-ALDH+) cells at diagnosis and after therapy. The proposed model can capture the clinically observed interindividual heterogeneity and reproduce the non-monotonous dynamics of the hematopoietic stem like cells observed in relapsing patients. Our model allows to simulate different scenarios proposed in literature such as therapy-related impairment of the stem cell niche or niche-mediated resistance. Model simulations suggest that during the post-therapy phase a more than 10-fold increase of hematopoietic stem-like cell proliferation rates is required to recapitulate the measured cell dynamics in patients achieving complete remission. We fit the model to data of 7 individual patients and simulate variations of the treatment protocol. These simulations are in line with the clinical finding that G-CSF priming can improve the treatment outcome. Furthermore, our model suggests that a decline of HSC counts during remission might serve as an indication for salvage therapy in patients lacking MRD (minimal residual disease) markers.

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.

Neutrophil-to-lymphocyte ratio and all-cause mortality with and without myeloproliferative neoplasms-a Danish longitudinal study.

The neutrophil-to-lymphocyte ratio(NLR) is increased in chronic inflammation and myeloproliferative neoplasms (MPN). We hypothesize that NLR is associated with all-cause mortality and mortality by comorbidity burden in the general population and individuals with MPN. We included 835,430 individuals from The Danish General Suburban Population Study, general practitioners, and outpatient clinics. We investigated NLR on mortality stratified by prevalent and incident MPN, essential thrombocythemia (ET), polycythemia vera (PV), myelofibrosis (MF), comorbidity burden (CCI-score), and the Triple-A risk score using hazard ratio (HR) and 95% confidence interval (95%CI). NLR 1-1.9 was the reference level. During a median follow-up of 11.2 years, 197,802 deaths were recorded. All-cause mortality increased for a stepwise increasing NLR with a HR (95%CI) for NLR ≥ 6 of 2.06(2.03-2.09) for the whole population and 2.93(2.44-3.50) in prevalent MPN. ET, PV, and MF had a HR (95%CI) for NLR ≥ 2 of 2.14(1.71-2.69), 2.19(1.89-2.54), and 2.31(1.91-2.80). Results were similar for incident MPN. Mortality was higher for stepwise increasing NLR and CCI-score(pinteraction < 2×10-16), with a HR for NLR ≥ 6 of 2.23(2.17-2.29), 4.10(4.01-4.20), and 7.69(7.50-7.89), for CCI-score 0, 1-2, or ≥3. The Triple-A risk score demonstrated alignment with NLR. Increasing NLR and comorbidity burden were associated with lower survival in individuals without MPN but were even worse in prevalent and incident MPN, ET, PV, and MF.

Transformative Materials for Interfacial Drug Delivery.

Drug delivery systems (DDS) are designed to temporally and spatially control drug availability and activity. They assist in improving the balance between on-target therapeutic efficacy and off-target toxic side effects. DDS aid in overcoming biological barriers encountered by drug molecules upon applying them via various routes of administration. They are furthermore increasingly explored for modulating the interface between implanted (bio)medical materials and host tissue. Herein, an overview of the biological barriers and host-material interfaces encountered by DDS upon oral, intravenous, and local administration is provided, and material engineering advances at different time and space scales to exemplify how current and future DDS can contribute to improved disease treatment are highlighted.

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.