Hematopoiesis in numbers

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.

Mathematical modeling of endogenous and exogenously administered T cell recirculation in mouse and its application to pharmacokinetic studies of cell therapies

Introduction

In vivo T cell migration has been of interest to scientists for the past 60 years. T cell kinetics are important in the understanding of the immune response to infectious agents. More recently, adoptive T cell therapies have proven to be a most promising approach to treating a wide range of diseases, including autoimmune and cancer diseases, whereby the characterization of cellular kinetics represents an important step towards the prediction of therapeutic efficacy.

Methods

Here, we developed a physiologically-based pharmacokinetic (PBPK) model that describes endogenous T cell homeostasis and the kinetics of exogenously administered T cells in mouse. Parameter calibration was performed using a nonlinear fixed-effects modeling approach based on published data on T cell kinetics and steady-state levels in different tissues of mice. The Partial Rank Correlation Coefficient (PRCC) method was used to perform a global sensitivity assessment. To estimate the impact of kinetic parameters on exogenously administered T cell dynamics, a local sensitivity analysis was conducted.

Results

We simulated the model to analyze cellular kinetics following various T cell doses and frequencies of CCR7+ T cells in the population of infused lymphocytes. The model predicted the effects of T cell numbers and of population composition of infused T cells on the resultant concentration of T cells in various organs. For example, a higher percentage of CCR7+ T cells among exogenously administered T lymphocytes led to an augmented accumulation of T cells in the spleen. The model predicted a linear dependence of T cell dynamics on the dose of adoptively transferred T cells.

Discussion

The mathematical model of T cell migration presented here can be integrated into a multi-scale model of the immune system and be used in a preclinical setting for predicting the distribution of genetically modified T lymphocytes in various organs, following adoptive T cell therapies.

Dynamic Tracking of Native Polyclonal Hematopoiesis in Adult Mice

Hematopoietic dysfunction has been associated with a reduction in the number of active precursors. However, precursor quantification at homeostasis and under diseased conditions is constrained by the scarcity of available methods. To address this issue, we optimized a method for quantifying a wide range of hematopoietic precursors. Assuming the random induction of a stable label in precursors following a binomial distribution, the estimation depends on the inverse correlation between precursor numbers and the variance of precursor labeling among independent samples. Experimentally validated to cover the full dynamic range of hematopoietic precursors in mice (1 to 105), we utilized this approach to demonstrate that thousands of precursors, which emerge after modest expansion during fetal-to-adult transition, contribute to native and perturbed hematopoiesis. We further estimated the number of precursors in a mouse model of Fanconi Anemia, showcasing how repopulation deficits can be segregated into autologous (cell proliferation) and non-autologous causes (lack of precursor). Our results support an accessible and reliable approach for precursor quantification, emphasizing the contemporary perspective that native hematopoiesis is highly polyclonal.

Structural characterization of Russula griseocarnosa polysaccharide and its improvement on hematopoietic function

The hematopoietic function of a polysaccharide derived from Russula griseocarnosa was demonstrated in K562 cells, and subsequently purified through chromatography to obtain RGP1. RGP1 is a galactan composed of 1,6-α-D-Galp as the main chain, with partial substitutions. A -CH3 substitution was detected at O-3 of 1,6-α-D-Galp. The possible branches at O-2 of 1,6-α-D-Galp was α-L-Fucp. In mice with cyclophosphamide (CTX)-induced hematopoietic dysfunction, RGP1 alleviated bone marrow damage and multinucleated giant cell infiltration of the spleen, increased the number of long-term hematopoietic stem cells, and regulated the levels of myeloid cells in the peripheral blood. Furthermore, RGP1 promoted the differentiation of activated T cells and CD4+ T cells without affecting natural killer cells and B cells. Proteomic analysis, detection of cytokines, and western blotting revealed that RGP1 could alleviate hematopoietic dysfunction by promoting the activation of CD4+ T cells and the Janus kinase/ signal transducer and activator of transcription 3 pathway. The present study provides experimental evidence to support the application of RGP1 in CTX-induced hematopoietic dysfunction.

Bone marrow inflammation in haematological malignancies

Tissue inflammation is a hallmark of tumour microenvironments. In the bone marrow, tumour-associated inflammation impacts normal niches for haematopoietic progenitor cells and mature immune cells and supports the outgrowth and survival of malignant cells residing in these niche compartments. This Review provides an overview of our current understanding of inflammatory changes in the bone marrow microenvironment of myeloid and lymphoid malignancies, using acute myeloid leukaemia and multiple myeloma as examples and highlights unique and shared features of inflammation in niches for progenitor cells and plasma cells. Importantly, inflammation exerts profoundly different effects on normal bone marrow niches in these malignancies, and we provide context for possible drivers of these divergent effects. We explore the role of tumour cells in inflammatory changes, as well as the role of cellular constituents of normal bone marrow niches, including myeloid cells and stromal cells. Integrating knowledge of disease-specific dynamics of malignancy-associated bone marrow inflammation will provide a necessary framework for future targeting of these processes to improve patient outcome.

An integrative mechanistic model of thymocyte dynamics

Background

The thymus plays a central role in shaping human immune function. A mechanistic, quantitative description of immune cell dynamics and thymic output under homeostatic conditions and various patho-physiological scenarios are of particular interest in drug development applications, e.g., in the identification of potential therapeutic targets and selection of lead drug candidates against infectious diseases.

Methods

We here developed an integrative mathematical model of thymocyte dynamics in human. It incorporates mechanistic features of thymocyte homeostasis as well as spatial constraints of the thymus and considerations of age-dependent involution. All model parameter estimates were obtained based on published physiological data of thymocyte dynamics and thymus properties in mouse and human. We performed model sensitivity analyses to reveal potential therapeutic targets through an identification of processes critically affecting thymic function; we further explored differences in thymic function across healthy subjects, multiple sclerosis patients, and patients on fingolimod treatment.

Results

We found thymic function to be most impacted by the egress, proliferation, differentiation and death rates of those thymocytes which are most differentiated. Model predictions also showed that the clinically observed decrease in relapse risk with age, in multiple sclerosis patients who would have discontinued fingolimod therapy, can be explained mechanistically by decreased thymic output with age. Moreover, we quantified the effects of fingolimod treatment duration on thymic output.

Conclusions

In summary, the proposed model accurately describes, in mechanistic terms, thymic output as a function of age. It may be further used to perform predictive simulations of clinically relevant scenarios which combine specific patho-physiological conditions and pharmacological interventions of interest.

Impact of age and sex on myelopoiesis and inflammation during myocardial infarction

Of all the different risk factors known to cause cardiovascular disease (CVD), age and sex are considered to play a crucial role. Aging follows a continuum from birth to death, and therefore it inevitably acts as a risk for CVD. Along with age, sex differences have also been shown to demonstrate variations in immune system responses to pathological insults. It has been widely perceived that females are protected against myocardial infarction (MI) and the protection is quite apparent in young vs. old women. Acute MI leads to changes in the population of myeloid and lymphoid cells at the injury site with myeloid bias being observed in the initial inflammation and the lymphoid in the late-resolution phases of the pathology. Multiple evidence demonstrates that aging enhances damage to various cellular processes through inflamm-aging, an inflammatory process identified to increase pro-inflammatory markers in circulation and tissues. Following MI, marked changes were observed in different sub-sets of major myeloid cell types viz., neutrophils, monocytes, and macrophages. There is a paucity of information regarding the tissue and site-specific functions of these sub-sets. In this review, we highlight the importance of age and sex as crucial risk factors by discussing their role during MI-induced myelopoiesis while emphasizing the current status of myeloid cell sub-sets. We further put forth the need for designing and executing age and sex interaction studies aimed to determine the appropriate age and sex to develop personalized therapeutic strategies post-MI.

Modeling early treatment response in AML from cell-free tumor DNA

Monitoring disease response after intensive chemotherapy for acute myeloid leukemia (AML) currently requires invasive bone marrow biopsies, imposing a significant burden on patients. In contrast, cell-free tumor DNA (ctDNA) in peripheral blood, carrying tumor-specific mutations, offers a less-invasive assessment of residual disease. However, the relationship between ctDNA levels and bone marrow blast kinetics remains unclear. We explored this in 10 AML patients with NPM1 and IDH2 mutations undergoing initial chemotherapy. Comparison of mathematical mixed-effect models showed that (1) inclusion of blast cell death in the bone marrow, (2) transition of ctDNA to peripheral blood, and (3) ctDNA decay in peripheral blood describes kinetics of blast cells and ctDNA best. The fitted model allows prediction of residual bone marrow blast content from ctDNA, and its scaling factor, representing clonal heterogeneity, correlates with relapse risk. Our study provides precise insights into blast and ctDNA kinetics, offering novel avenues for AML disease monitoring.

Breast cancer remotely imposes a myeloid bias on haematopoietic stem cells by reprogramming the bone marrow niche

Myeloid cell infiltration of solid tumours generally associates with poor patient prognosis and disease severity1-13. Therefore, understanding the regulation of myeloid cell differentiation during cancer is crucial to counteract their pro-tumourigenic role. Bone marrow (BM) haematopoiesis is a tightly regulated process for the production of all immune cells in accordance to tissue needs14. Myeloid cells differentiate during haematopoiesis from multipotent haematopoietic stem and progenitor cells (HSPCs)15-17. HSPCs can sense inflammatory signals from the periphery during infections18-21 or inflammatory disorders22-27. In these settings, HSPC expansion is associated with increased myeloid differentiation28,29. During carcinogenesis, the elevation of haematopoietic growth factors supports the expansion and differentiation of committed myeloid progenitors5,30. However, it is unclear whether cancer-related inflammation also triggers demand-adapted haematopoiesis at the level of multipotent HSPCs. In the BM, HSPCs reside within the haematopoietic niche which delivers HSC maintenance and differentiation cues31-35. Mesenchymal stem cells (MSCs) are a major cellular component of the BM niche and contribute to HSC homeostasis36-41. Modifications of MSCs in systemic disorders have been associated with HSC differentiation towards myeloid cells22,42. It is unknown if MSCs are regulated in the context of solid tumours and if their myeloid supportive activity is impacted by cancer-induced systemic changes. Here, using unbiased transcriptomic analysis and in situ imaging of HSCs and the BM niche during breast cancer, we show that both HSCs and MSCs are transcriptionally and spatially modified. We demonstrate that breast tumour can distantly remodel the cellular cross-talks in the BM niche leading to increased myelopoiesis.

Dynamics of Human Hematopoietic Stem and Progenitor Cell Differentiation to the Erythroid Lineage

Erythropoiesis, the development of erythrocytes from hematopoietic stem cells, occurs through four phases: erythroid progenitor development, early erythropoiesis, terminal erythroid differentiation (TED), and maturation. According to the classical model that is based on immunophenotypic profiles of cell populations, each of these phases comprises multiple differentiation states that arise in a hierarchical manner. After segregation of lymphoid potential, erythroid priming begins during progenitor development and progresses through progenitor cell types that have multilineage potential. Complete separation of the erythroid lineage is achieved during early erythropoiesis with the formation of unipotent erythroid progenitors: burst forming unit-erythroid and colony forming unit-erythroid. These erythroid committed progenitors undergo TED and maturation, which involves expulsion of the nucleus and remodeling to form functional biconcave, hemoglobin-filled erythrocytes. In the last decade or so, many studies employing advanced techniques such as single cell RNA-sequencing (scRNA-seq) as well as the conventional methods, including colony forming cell assays and immunophenotyping, have revealed heterogeneity within the stem, progenitor, and erythroblast stages, and uncovered alternate paths for segregation of erythroid lineage potential. In this review, we provide an in-depth account of immunophenotypic profiles of all cell types within erythropoiesis, highlight studies that demonstrate heterogeneous erythroid stages, and describe deviations to the classical model of erythropoiesis. Overall, although scRNA-seq approaches have provided new insights, flow cytometry remains relevant and is the primary method for validation of novel immunophenotypes.

DRAG in situ barcoding reveals an increased number of HSPCs contributing to myelopoiesis with age

Ageing is associated with changes in the cellular composition of the immune system. During ageing, hematopoietic stem and progenitor cells (HSPCs) that produce immune cells are thought to decline in their regenerative capacity. However, HSPC function has been mostly assessed using transplantation assays, and it remains unclear how HSPCs age in the native bone marrow niche. To address this issue, we present an in situ single cell lineage tracing technology to quantify the clonal composition and cell production of single cells in their native niche. Our results demonstrate that a pool of HSPCs with unequal output maintains myelopoiesis through overlapping waves of cell production throughout adult life. During ageing, the increased frequency of myeloid cells is explained by greater numbers of HSPCs contributing to myelopoiesis rather than the increased myeloid output of individual HSPCs. Strikingly, the myeloid output of HSPCs remains constant over time despite accumulating significant transcriptomic changes throughout adulthood. Together, these results show that, unlike emergency myelopoiesis post-transplantation, aged HSPCs in their native microenvironment do not functionally decline in their regenerative capacity.

Histone H2A deubiquitinases in the transcriptional programs of development and hematopoiesis: a consolidated analysis

Monoubiquitinated lysine 119 of histone H2A (H2AK119ub) is a highly abundant epigenetic mark, associated with gene repression and deposited on chromatin by the polycomb repressor complex 1 (PRC1), which is an essential regulator of diverse transcriptional programs in mammalian development and tissue homeostasis. While multiple deubiquitinases (DUBs) with catalytic activity for H2AK119ub (H2A-DUBs) have been identified, we lack systematic analyses of their roles and cross-talk in transcriptional regulation. Here, we address H2A-DUB functions in epigenetic regulation of mammalian development and tissue maintenance by conducting a meta-analysis of 248 genomics datasets from 32 independent studies, focusing on the mouse model and covering embryonic stem cells (ESCs), hematopoietic, and immune cell lineages. This covers all the publicly available datasets that map genomic H2A-DUB binding and H2AK119ub distributions (ChIP-Seq), and all datasets assessing dysregulation in gene expression in the relevant H2A-DUB knockout models (RNA-Seq). Many accessory datasets for PRC1-2 and DUB-interacting proteins are also analyzed and interpreted, as well as further data assessing chromatin accessibility (ATAC-Seq) and transcriptional activity (RNA-seq). We report co-localization in the binding of H2A-DUBs BAP1, USP16, and to a lesser extent others that is conserved across different cell-types, and also the enrichment of antagonistic PRC1-2 protein complexes at the same genomic locations. Such conserved sites enriched for the H2A-DUBs and PRC1-2 are proximal to transcriptionally active genes that engage in housekeeping cellular functions. Nevertheless, they exhibit H2AK119ub levels significantly above the genomic average that can undergo further increase with H2A-DUB knockout. This indicates a cooperation between H2A-DUBs and PRC1-2 in the modulation of housekeeping transcriptional programs, conserved across many cell types, likely operating through their antagonistic effects on H2AK119ub and the regulation of local H2AK119ub turnover. Our study further highlights existing knowledge gaps and discusses important directions for future work.

Hematologic dysfunction in cancer: Mechanisms, effects on antitumor immunity, and roles in disease progression

With the major advances in cancer immunology and immunotherapy, it is critical to consider that most immune cells are short-lived and need to be continuously replenished from hematopoietic stem and progenitor cells. Hematologic abnormalities are prevalent in cancer patients, and many ground-breaking studies over the past decade provide insights into their underlying cellular and molecular mechanisms. Such studies demonstrate that the dysfunction of hematopoiesis is more than a side-effect of cancer pathology, but an important systemic feature of cancer disease. Here we review these many advances, covering the cancer-associated phenotypes of hematopoietic stem and progenitor cells, the dysfunction of myelopoiesis and erythropoiesis, the importance of extramedullary hematopoiesis in cancer disease, and the developmental origins of tumor associated macrophages. We address the roles of many secreted mediators, signaling pathways, and transcriptional and epigenetic mechanisms that mediate such hematopoietic dysfunction. Furthermore, we discuss the important contribution of the hematopoietic dysfunction to cancer immunosuppression, the possible avenues for therapeutic intervention, and highlight the unanswered questions and directions for future work. Overall, hematopoietic dysfunction is established as an active component of the cancer disease mechanisms and an important target for therapeutic intervention.

Why are cell populations maintained via multiple compartments?

We consider the maintenance of 'product' cell populations from 'progenitor' cells via a sequence of one or more cell types, or compartments, where each cell's fate is chosen stochastically. If there is only one compartment then large amplification, that is, a large ratio of product cells to progenitors comes with disadvantages. The product cell population is dominated by large families (cells descended from the same progenitor) and many generations separate, on average, product cells from progenitors. These disadvantages are avoided using suitably constructed sequences of compartments: the amplification factor of a sequence is the product of the amplification factors of each compartment, while the average number of generations is a sum over contributions from each compartment. Passing through multiple compartments is, in fact, an efficient way to maintain a product cell population from a small flux of progenitors, avoiding excessive clonality and minimizing the number of rounds of division en route. We use division, exit and death rates, estimated from measurements of single-positive thymocytes, to choose illustrative parameter values in the single-compartment case. We also consider a five-compartment model of thymocyte differentiation, from double-negative precursors to single-positive product cells.

Plasma iron controls neutrophil production and function

Low plasma iron (hypoferremia) induced by hepcidin is a conserved inflammatory response that protects against infections but inhibits erythropoiesis. How hypoferremia influences leukocytogenesis is unclear. Using proteomic data, we predicted that neutrophil production would be profoundly more iron-demanding than generation of other white blood cell types. Accordingly in mice, hepcidin-mediated hypoferremia substantially reduced numbers of granulocytes but not monocytes, lymphocytes, or dendritic cells. Neutrophil rebound after anti-Gr-1-induced neutropenia was blunted during hypoferremia but was rescued by supplemental iron. Similarly, hypoferremia markedly inhibited pharmacologically stimulated granulopoiesis mediated by granulocyte colony-stimulating factor and inflammation-induced accumulation of neutrophils in the spleen and peritoneal cavity. Furthermore, hypoferremia specifically altered neutrophil effector functions, suppressing antibacterial mechanisms but enhancing mitochondrial reactive oxygen species-dependent NETosis associated with chronic inflammation. Notably, antagonizing endogenous hepcidin during acute inflammation enhanced production of neutrophils. We propose plasma iron modulates the profile of innate immunity by controlling monocyte-to-neutrophil ratio and neutrophil activity in a therapeutically targetable system.

DNA methylation-mediated differential expression of DLX4 isoforms has opposing roles in leukemogenesis

Background

Previously, we reported the expression of DLX4 isoforms (BP1 and DLX7) in myeloid leukemia, but the functional role of DLX4 isoforms remains poorly understood. In the work described herein, we further determined the underlying role of DLX4 isoforms in chronic myeloid leukemia (CML) leukemogenesis.

Methods

The expression and methylation of DLX4 isoforms were detected by real-time quantitative PCR (RT-qPCR) and real-time quantitative methylation-specific PCR (RT-qMSP) in patients with CML. The functional role of DLX4 isoforms was determined in vitro and in vivo. The molecular mechanism of DLX4 isoforms in leukemogenesis was identified based on chromatin immunoprecipitation with high-throughput sequencing (ChIP-Seq)/assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-Seq) and RNA sequencing (RNA-Seq).

Results

BP1 expression was increased in patients with CML with unmethylated promoter, but DLX7 expression was decreased with hypermethylated promoter. Functionally, overexpression of BP1 increased the proliferation rate of K562 cells with S/G2 promotion, whereas DLX7 overexpression reduced the proliferation rate of K562 cells with G1 arrest. Moreover, K562 cells with BP1 overexpression increased the tumorigenicity in NCG mice, whereas K562 cells with DLX7 overexpression decreased the tumorigenicity. Mechanistically, a total of 91 genes including 79 messenger RNAs (mRNAs) and 12 long noncoding RNAs (lncRNAs) were discovered by ChIP-Seq and RNA-Seq as direct downstream targets of BP1. Among the downstream genes, knockdown of RREB1 and SGMS1-AS1 partially revived the proliferation caused by BP1 overexpression in K562 cells. Similarly, using ATAC-Seq and RNA-Seq, a total of 282 genes including 151 mRNA and 131 lncRNAs were identified as direct downstream targets of DLX7. Knockdown of downstream genes PTPRB and NEAT1 partially revived the proliferation caused by DLX7 overexpression in K562 cells. Finally, we also identified and validated a SGMS1-AS1/miR-181d-5p/SRPK2 competing endogenous RNA (ceRNA) network caused by BP1 overexpression in K562 cells.

Conclusions

The current findings reveal that DNA methylation-mediated differential expression of DLX4 isoforms BP1 and DLX7 plays opposite functions in leukemogenesis. BP1 plays an oncogenic role in leukemia development, whereas DLX7 acts as a tumor suppressor gene. These results suggest DLX4 as a therapeutic target for antileukemia therapy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s11658-022-00358-0.

Niche Regulation of Hematopoiesis: The Environment Is "Micro," but the Influence Is Large

Immune cell production is governed by a process known as hematopoiesis, where hematopoietic stem cells (HSCs) differentiate through progenitor cells and ultimately to the mature blood and immune cells found in circulation. While HSCs are capable of cell-autonomous regulation, they also rely on extrinsic factors to balance their state of quiescence and activation. These cues can, in part, be derived from the niche in which HSCs are found. Under steady-state conditions, HSCs are found in the bone marrow. This niche is designed to support HSCs but also to respond to external factors, which allows hematopoiesis to be a finely tuned and coordinated process. However, the niche, and its regulation, can become dysregulated to potentiate inflammation during disease. This review will highlight the architecture of the bone marrow and key regulators of hematopoiesis within this niche. Emphasis will be placed on how these mechanisms go awry to exacerbate hematopoietic contributions that drive cardiovascular disease.