Modeling of replicative senescence in hematopoietic development

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.

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.

Wnt-5A/B Signaling in Hematopoiesis throughout Life

Wnt signaling is well-known to play major roles in the hematopoietic system, from embryogenesis to aging and disease. In addition to the main β-catenin-dependent pathway, it is now clear that Wnt5a and the structurally related Wnt5b are essential for hematopoiesis, bone marrow colonization and the final steps of hematopoietic stem cell (HSC) maturation via β-catenin-independent signaling. Wnt5a and Wnt5b ligands prevent hematopoietic exhaustion (by maintaining quiescent, long-term HSCs), induce the proliferation of progenitors, and guide myeloid development, in addition to being involved in the development of aging-related alterations. The aim of this review is to summarize the current knowledge on these roles of Wnt5a and Wn5b signaling in the hematopoietic field.

Personal response to immune checkpoint inhibitors of patients with advanced melanoma explained by a computational model of cellular immunity, tumor growth, and drug

Immune checkpoint inhibitors, such as pembrolizumab, are transforming clinical oncology. Yet, insufficient overall response rate, and accelerated tumor growth rate in some patients, highlight the need for identifying potential responders. To construct a computational model, identifying response predictors, and enabling immunotherapy personalization. The combined dynamics of cellular immunity, pembrolizumab, and the melanoma cancer were modeled by a set of ordinary differential equations. The model relies on a scheme of T memory stem cells, progressively differentiating into effector CD8⁺ T cells, and additionally includes T cell exhaustion, reinvigoration and senescence. Clinical data of a pembrolizumab-treated patient with advanced melanoma (Patient O’) were used for model calibration and simulations. Virtual patient populations, varying in one parameter or more, were generated for retrieving clinical studies. Simulations captured the major features of Patient O’s disease, displaying a good fit to her clinical data. A temporary increase in tumor burden, as implied by the clinical data, was obtained only when assuming aberrant self-renewal rates. Variation in effector T cell cytotoxicity was sufficient for simulating dynamics that vary from rapid progression to complete cure, while variation in tumor immunogenicity has a delayed and limited effect on response. Simulations of a-specific clinical trial were in good agreement with the clinical results, demonstrating positive correlations between response to pembrolizumab and the ratio of reinvigoration to baseline tumor load. These results were obtained by assuming inter-patient variation in the toxicity of effector CD8+ T cells, and in their intrinsic division rate, as well as by assuming that the intrinsic division rate of cancer cells is correlated with the baseline tumor burden. In conclusion, hyperprogression can result from lower patient-specific effector cytotoxicity, a temporary increase in tumor load is unlikely to result from real tumor growth, and the ratio of reinvigoration to tumor load can predict personal response to pembrolizumab. Upon further validation, the model can serve for immunotherapy personalization.

Replicative cellular age distributions in compartmentalized tissues

The cellular age distribution of hierarchically organized tissues can reveal important insights into the dynamics of cell differentiation and self-renewal and associated cancer risks. Here, we examine the effect of progenitor compartments with varying differentiation and self-renewal capacities on the resulting observable distributions of replicative cellular ages. We find that strongly amplifying progenitor compartments, i.e. compartments with high self-renewal capacities, substantially broaden the age distributions which become skewed towards younger cells with a long tail of few old cells. For several of these strongly amplifying compartments, the age distribution becomes virtually independent of the influx from the stem cell compartment. By contrast, if tissues are organized into many downstream compartments with low self-renewal capacity, the shape of the replicative cell distribution in more differentiated compartments is dominated by stem cell dynamics with little added variation. In the limiting case of a strict binary differentiation tree without self-renewal, the shape of the output distribution becomes indistinguishable from that of the input distribution. Our results suggest that a comparison of cellular age distributions between healthy and cancerous tissues may inform about dynamical changes within the hierarchical tissue structure, i.e. an acquired increased self-renewal capacity in certain tumours. Furthermore, we compare our theoretical results to telomere length distributions in granulocyte populations of 10 healthy individuals across different ages, highlighting that our theoretical expectations agree with experimental observations.

A dynamical systems model for the measurement of cellular senescence

Senescent cells provide a good in vitro model to study ageing. However, cultures of ‘senescent’ cells consist of a mix of cell subtypes (proliferative, senescent, growth-arrested and apoptotic). Determining the proportion of senescent cells is crucial for studying ageing and developing new anti-degenerative therapies. Commonly used markers such as doubling population, senescence-associated β-galactosidase, Ki-67, γH2AX and TUNEL assays capture diverse and overlapping cellular populations and are not purely specific to senescence. A newly developed dynamical systems model follows the transition of an initial culture to senescence tracking population doubling, and the proportion of cells in proliferating, growth-arrested, apoptotic and senescent states. Our model provides a parsimonious description of transitions between these states accruing towards a predominantly senescent population. Using a genetic algorithm, these model parameters are well constrained by an in vitro human primary fibroblast dataset recording five markers at 16 time points. The computational model accurately fits to the data and translates these joint markers into the first complete description of the proportion of cells in different states over the lifetime. The high temporal resolution of the dataset demonstrates the efficacy of strategies for reconstructing the trajectory towards replicative senescence with a minimal number of experimental recordings.

Mathematical modeling of bone marrow - peripheral blood dynamics in the disease state based on current emerging paradigms, part II

The cancer stem cell hypothesis has gained currency in recent times but concerns remain about its scientific foundations because of significant gaps that exist between research findings and comprehensive knowledge about cancer stem cells (CSCs). In this light, a mathematical model that considers hematopoietic dynamics in the diseased state of the bone marrow and peripheral blood is proposed and used to address findings about CSCs. The ensuing model, resulting from a modification and refinement of a recent model, develops out of the position that mathematical models of CSC development, that are few at this time, are needed to provide insightful underpinnings for biomedical findings about CSCs as the CSC idea gains traction. Accordingly, the mathematical challenges brought on by the model that mirror general challenges in dealing with nonlinear phenomena are discussed and placed in context. The proposed model describes the logical occurrence of discrete time delays, that by themselves present mathematical challenges, in the evolving cell populations under consideration. Under the challenging circumstances, the steady state properties of the model system of delay differential equations are obtained, analyzed, and the resulting mathematical predictions arising therefrom are interpreted and placed within the framework of findings regarding CSCs. Simulations of the model are carried out by considering various parameter scenarios that reflect different experimental situations involving disease evolution in human hosts. Model analyses and simulations suggest that the emergence of the cancer stem cell population alongside other malignant cells engenders higher dimensions of complexity in the evolution of malignancy in the bone marrow and peripheral blood at the expense of healthy hematopoietic development. The model predicts the evolution of an aberrant environment in which the malignant population particularly in the bone marrow shows tendencies of reaching an uncontrollable equilibrium state. Essentially, the model shows that a structural relationship exists between CSCs and non-stem malignant cells that confers on CSCs the role of temporally enhancing and stimulating the expansion of non-stem malignant cells while also benefitting from increases in their own population and these CSCs may be the main protagonists that drive the ultimate evolution of the uncontrollable equilibrium state of such malignant cells and these may have implications for treatment.

Modeling large fluctuations of thousands of clones during hematopoiesis: The role of stem cell self-renewal and bursty progenitor dynamics in rhesus macaque

In a recent clone-tracking experiment, millions of uniquely tagged hematopoietic stem cells (HSCs) and progenitor cells were autologously transplanted into rhesus macaques and peripheral blood containing thousands of tags were sampled and sequenced over 14 years to quantify the abundance of hundreds to thousands of tags or “clones.” Two major puzzles of the data have been observed: consistent differences and massive temporal fluctuations of clone populations. The large sample-to-sample variability can lead clones to occasionally go “extinct” but “resurrect” themselves in subsequent samples. Although heterogeneity in HSC differentiation rates, potentially due to tagging, and random sampling of the animals’ blood and cellular demographic stochasticity might be invoked to explain these features, we show that random sampling cannot explain the magnitude of the temporal fluctuations. Moreover, we show through simpler neutral mechanistic and statistical models of hematopoiesis of tagged cells that a broad distribution in clone sizes can arise from stochastic HSC self-renewal instead of tag-induced heterogeneity. The very large clone population fluctuations that often lead to extinctions and resurrections can be naturally explained by a generation-limited proliferation constraint on the progenitor cells. This constraint leads to bursty cell population dynamics underlying the large temporal fluctuations. We analyzed experimental clone abundance data using a new statistic that counts clonal disappearances and provided least-squares estimates of two key model parameters in our model, the total HSC differentiation rate and the maximum number of progenitor-cell divisions.

Hematopoiesis of virally tagged cells in rhesus macaques is analyzed in the context of a mechanistic and statistical model. We find that the clone size distribution and the temporal variability in the abundance of each clone (viral tag) in peripheral blood are consistent with (i) stochastic HSC self-renewal during bone marrow repair, (ii) clonal aging that restricts the number of generations of progenitor cells, and (iii) infrequent and small-size samples. By fitting data, we infer two key parameters that control the level of fluctuations of clone sizes in our model: the total HSC differentiation rate and the maximum proliferation capacity of progenitor cells. Our analysis provides insight into the mechanisms of hematopoiesis and a framework to guide future multiclone barcoding/lineage tracking measurements.

Understanding Hematological Toxicities Using Mathematical Modeling

Balancing antitumor efficacy with toxicity is a significant challenge, and drug-induced myelosuppression is a common dose-limiting toxicity of cancer treatments. Mathematical modeling has proven to be a powerful ally in this field, scaling results from animal models to humans, and designing optimized treatment regimens. Here we outline existing mathematical approaches for studying bone marrow toxicity, identify gaps in current understanding, and make future recommendations to advance this vital field of safety research further.

Cell Density-Dependent Cytological Stage Profile and Its Application for a Screen of Cytostatic Agents Active Toward Leukemic Stem Cells

Proliferation and expansion of leukemia is driven by leukemic stem cells (LSCs). Multidrug resistance (MDR) of LSCs is one of the main reasons of failure and relapses in acute myeloid leukemia (AML) treatment. In this study, we show that maintaining HL-60 at low cell culture density or applying a 240-day treatment with anthrapyridazone (BS-121) increased the percentage of primitive cells, which include LSCs determining the overall stage profile. This change manifested in morphology, expression of both cell surface markers and redox-state proteins, as well as mitochondrial potential. Moreover, four sublines were generated, each with unique and characteristic stage profile and cytostatic sensitivity. Cell density-induced culture alterations (affecting stage profiles) were exploited in a screen of anthrapyridazones. Among the compound tested, C-123 was the most potent against primitive cell stages while generating relatively low amounts of reactive oxygen species (ROS). Furthermore, it had low toxicity in vivo and weakly affected blood morphology of healthy mice. The cell density-dependent stage profiles could be utilized in preliminary drug screens for activity against LSCs or in construction of patient-specific platforms to find drugs effective in case of AML relapse (drug extrapolation). The correlation between ROS generation in differentiated cells and toxic effect observed in HL-60 has a potential application in myelotoxicity predictions. The discovered properties of C-123 indicate its potential application in AML treatment, specifically in conditioned myeloablation preceding allogeneic transplantation and/or ex vivo treatment preceding autologous transplantation.

Quantifying replicative senescence as a tumor suppressor pathway and a target for cancer therapy

To study quantitatively replicative senescence as a tumor suppressor mechanism, we investigate the distribution of a growing clonal cell population restricted by Hayflick’s limit. We find that in the biologically relevant range of parameters, if the imbalance between cell division and death is moderate or low (high death-to-birth ratio), senescence offers significant protection against cancer by halting abnormal cell proliferation at early pre-diagnostic stages of tumor development. We also find that by the time tumors are typically detected, there is a high probability that telomerase is activated, even if the cell of origin was telomerase negative. Hence, the fact that most cancers are positive for telomerase is not necessarily an indication that cancer originated in a telomerase positive cell. Finally, we discuss how the population dynamics of cells can determine the outcomes of anti-telomerase cancer therapies, and provide guidelines on how the model could potentially be applied to develop clinically useful tools to predict the response to treatment by telomerase inhibitors in individual patients.

MicroRNAs and Metabolites in Serum Change after Chemotherapy: Impact on Hematopoietic Stem and Progenitor Cells

Hematopoietic regeneration after high dose chemotherapy necessitates activation of the stem cell pool. There is evidence that serum taken after chemotherapy comprises factors stimulating proliferation and self-renewal of CD34⁺ hematopoietic stem and progenitor cells (HSPCs) – however, the nature of these feedback signals is yet unclear. Here, we addressed the question if specific microRNAs (miRNAs) or metabolites are affected after high dose chemotherapy. Serum taken from the same patients before and after chemotherapy was supplemented for in vitro cultivation of HSPCs. Serum taken after chemotherapy significantly enhanced HSPC proliferation, better maintained a CD34⁺ immunophenotype, and stimulated colony forming units. Microarray analysis revealed that 23 miRNAs changed in serum after chemotherapy – particularly, miRNA-320c and miRNA-1275 were down-regulated whereas miRNA-3663-3p was up-regulated. miRNA-320c was exemplarily inhibited by an antagomiR, which seemed to increase proliferation. Metabolomic profiling demonstrated that 44 metabolites were less abundant, whereas three (including 2-hydroxybutyrate and taurocholenate sulphate) increased in serum upon chemotherapy. Nine of these metabolites were subsequently tested for effects on HSPCs in vitro, but none of them exerted a clear concentration dependent effect on proliferation, immunophenotype and colony forming unit formation. Taken together, serum profiles of miRNAs and metabolites changed after chemotherapy. Rather than individually, these factors may act in concert to recruit HSPCs into action for hematopoietic regeneration.

Mathematical Modelling as a Tool to Understand Cell Self-renewal and Differentiation

Mathematical modeling is a powerful technique to address key questions and paradigms in a variety of complex biological systems and can provide quantitative insights into cell kinetics, fate determination and development of cell populations. The chapter is devoted to a review of modeling of the dynamics of stem cell-initiated systems using mathematical methods of ordinary differential equations. Some basic concepts and tools for cell population dynamics are summarized and presented as a gentle introduction to non-mathematicians. The models take into account different plausible mechanisms regulating homeostasis. Two mathematical frameworks are proposed reflecting, respectively, a discrete (punctuated by division events) and a continuous character of transitions between differentiation stages. Advantages and constraints of the mathematical approaches are presented on examples of models of blood systems and compared to patients data on healthy hematopoiesis.

A multicompartment mathematical model of cancer stem cell-driven tumor growth dynamics

Tumors are appreciated to be an intrinsically heterogeneous population of cells with varying proliferation capacities and tumorigenic potentials. As a central tenet of the so-called cancer stem cell hypothesis, most cancer cells have only a limited lifespan, and thus cannot initiate or reinitiate tumors. Longevity and clonogenicity are properties unique to the subpopulation of cancer stem cells. To understand the implications of the population structure suggested by this hypothesis--a hierarchy consisting of cancer stem cells and progeny non-stem cancer cells which experience a reduction in their remaining proliferation capacity per division--we set out to develop a mathematical model for the development of the aggregate population. We show that overall tumor progression rate during the exponential growth phase is identical to the growth rate of the cancer stem cell compartment. Tumors with identical stem cell proportions, however, can have different growth rates, dependent on the proliferation kinetics of all participating cell populations. Analysis of the model revealed that the proliferation potential of non-stem cancer cells is likely to be small to reproduce biologic observations. Furthermore, a single compartment of non-stem cancer cell population may adequately represent population growth dynamics only when the compartment proliferation rate is scaled with the generational hierarchy depth.

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.

Myelodysplastic syndromes are diseases which are characterized by ineffective blood formation. There is accumulating evidence that they are caused by an aberrant hematopoietic stem cell. However, it is yet unclear how this malignant clone suppresses normal hematopoiesis. To this end, we generated mathematical models under the assumption that feedback signals regulate self-renewal and proliferation of normal and diseased stem cells. The simulations demonstrate that the malignant cells must have particularly higher self-renewal rates than normal stem cells – rather than higher proliferation rates. On the other hand, down-regulation of self-renewal by the increasing number of malignant cells in the bone marrow niche can explain impairment of normal blood formation. In fact, we show that serum of patients with myelodysplastic syndrome, as compared to serum of healthy donors, stimulates proliferation and moderately impacts on maintenance of hematopoietic stem and progenitor cells in vitro. Thus, aberrant high self-renewal rates of the malignant clone seem to initiate disease development; suppression of normal blood formation is then caused by a rebound effect of feedback signals which down-regulate self-renewal of normal stem and progenitor cells as well.

Biphasic modulation of cancer stem cell-driven solid tumour dynamics in response to reactivated replicative senescence

OBJECTIVES

Cell senescence is a physiological programme of irreversible mitotic arrest that is triggered after a variety of intracellular and extracellular events. Its purpose is to protect tissue integrity by disabling mitosis in stressed or damaged cells. The senescence program serves as a tumour suppressor, and cancer cells are believed to bypass senescence to advance to malignancy. Recent studies have shown that senescence can be reactivated in cancer cells through a number of external perturbations, including oncogene activation, tumour suppressor gene withdrawal and irradiation.

MATERIALS AND METHODS

We have developed an agent-based model of solid tumour growth whose input population composition is based on the cancer stem-cell hypothesis. It is used to show how cancer stem cells can drive tumour progression, while non-stem cancer cells (CCs) interfere with this by impeding cancer stem-cell dynamics.

RESULTS

Here we show that intratumoural competition between the two cell types may arise to modulate tumour progression and ultimately cancer presentation risk. Model simulations reveal that reactivation of the replicative senescence programme in CCs initially increases total tumour burden, as attrition from cell death is partially averted, but evolves to provide tumour control in the long-term through increasing constraints on stem-cell compartment kinetics.

CONCLUSIONS

Reactivation of replicative senescence can prolong CC competition with cancer stem cells, thereby ultimately inhibiting malignant progression regardless of tumour size.

Clonal selection and therapy resistance in acute leukaemias: mathematical modelling explains different proliferation patterns at diagnosis and relapse

Recent experimental evidence suggests that acute myeloid leukaemias may originate from multiple clones of malignant cells. Nevertheless, it is not known how the observed clones may differ with respect to cell properties, such as proliferation and self-renewal. There are scarcely any data on how these cell properties change due to chemotherapy and relapse. We propose a new mathematical model to investigate the impact of cell properties on the multi-clonal composition of leukaemias. Model results imply that enhanced self-renewal may be a key mechanism in the clonal selection process. Simulations suggest that fast proliferating and highly self-renewing cells dominate at primary diagnosis, while relapse following therapy-induced remission is triggered mostly by highly self-renewing but slowly proliferating cells. Comparison of simulation results to patient data demonstrates that the proposed model is consistent with clinically observed dynamics based on a clonal selection process.

Stochasticity and determinism in models of hematopoiesis

This chapter represents a novel view of modeling in hematopoiesis, synthesizing both deterministic and stochastic approaches. Whereas the stochastic models work in situations where chance dominates, for example when the number of cells is small, or under random mutations, the deterministic models are more important for large-scale, normal hematopoiesis. New types of models are on the horizon. These models attempt to account for distributed environments such as hematopoietic niches and their impact on dynamics. Mixed effects of such structures and chance events are largely unknown and constitute both a challenge and promise for modeling. Our discussion is presented under the separate headings of deterministic and stochastic modeling; however, the connections between both are frequently mentioned. Four case studies are included to elucidate important examples. We also include a primer of deterministic and stochastic dynamics for the reader's use.

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.

Minimizing the risk of cancer: tissue architecture and cellular replication limits

Normal somatic cells are capable of only a limited number of divisions, which prevents unlimited cell proliferation and the onset of tumours. Cancer cells find ways to circumvent this obstacle, typically by expressing the enzyme telomerase and less often by alternative recombination strategies. Given this fundamental link between cellular replication limits and cancer, it is important to understand how a tissue's architecture affects the replicative capacity of a cell population. We define this as the average number of remaining divisions at equilibrium. The lower the replication capacity, the lower the chances to escape the replication limit during abnormal growth when a tumour develops. In this paper, we examine how the replication capacity is influenced by defining characteristics of cell lineages, such as the number of intermediate cell compartments, self-renewal capability of cells and division rates. We describe an optimal tissue architecture that minimizes the replication capacity of dividing cells and thus the risk of cancer. Interestingly, some of the features that define an optimal tissue architecture have been documented in a variety of tissues, suggesting that they may have evolved as a cancer-protecting strategy in multicellular organisms.

Aging induced decline in T-lymphopoiesis is primarily dependent on status of progenitor niches in the bone marrow and thymus

Age-related decline in the generation of T cells is associated with two primary lymphoid organs, the bone marrow (BM) and thymus. Both organs contain lympho-hematopoietic progenitor/stem cells (LPCs) and non-hematopoietic stromal/niche cells. Murine model showed this decline is not due to reduced quantities of LPCs, nor autonomous defects in LPCs, but rather defects in their niche cells. However, this viewpoint is challenged by the fact that aged BM progenitors have a myeloid skew. By grafting young wild-type (WT) BM progenitors into aged IL-7R^−/− hosts, which possess WT-equivalent niches although LPCs are defect, we demonstrated that these young BM progenitors also exhibited a myeloid skew. We, further, demonstrated that aged BM progenitors, recruited by a grafted fetal thymus in the in vivo microenvironment, were able to compete with their young counterparts, although the in vitro manipulated old BM cells were not able to do so in conventional BM transplantation. Both LPCs and their niche cells inevitably get old with increasing organismal age, but aging in niche cells occurred much earlier than in LPCs by an observation in thymic T-lymphopoiesis. Therefore, the aging induced decline in competence to generate T cells is primarily dependent on status of the progenitor niche cells in the BM and thymus.

Epigenetic-aging-signature to determine age in different tissues

All tissues of the organism are affected by aging. This process is associated with epigenetic modifications such as methylation changes at specific cytosine residues in the DNA (CpG sites). Here, we have identified an Epigenetic-Aging-Signature which is applicable for many tissues to predict donor age. DNA-methylation profiles of various cell types were retrieved from public data depositories - all using the HumanMethylation27 BeadChip platform which represents 27,578 CpG sites. Five datasets from dermis, epidermis, cervical smear, T-cells and monocytes were used for Pavlidis Template Matching to identify 19 CpG sites that are continuously hypermethylated upon aging (R > 0.6; p-value <10⁻¹³). Four of these CpG sites (associated with the genes NPTX2, TRIM58, GRIA2 and KCNQ1DN) and an additional hypomethylated CpG site (BIRC4BP) were implemented in a model to predict donor age. This Epigenetic-Aging-Signature was tested on a validation group of eight independent datasets corresponding to several cell types from different tissues. Overall, the five CpG sites revealed age-associated DNA-methylation changes in all tissues. The average absolute difference between predicted and real chronological age was about 11 years. This method can be used to predict donor age in various cell preparations - for example in forensic analysis.

Replicative senescence of mesenchymal stem cells causes DNA-methylation changes which correlate with repressive histone marks

Cells in culture undergo replicative senescence. In this study, we analyzed functional, genetic and epigenetic sequels of long-term culture in human mesenchymal stem cells (MSC). Already within early passages the fibroblastoid colonyforming unit (CFU-f) frequency and the differentiation potential of MSC declined significantly. Relevant chromosomal aberrations were not detected by karyotyping and SNP-microarrays. Subsequently, we have compared DNA-methylation profiles with the Infinium HumanMethylation27 Bead Array and the profiles differed markedly in MSC derived from adipose tissue and bone marrow. Notably, all MSC revealed highly consistent senescence-associated modifications at specific CpG sites. These DNA-methylation changes correlated with histone marks of previously published data sets, such as trimethylation of H3K9, H3K27 and EZH2 targets. Taken together, culture expansion of MSC has profound functional implications - these are hardly reflected by genomic instability but they are associated with highly reproducible DNA-methylation changes which correlate with repressive histone marks. Therefore replicative senescence seems to be epigenetically controlled.

Cooperative effects of Akt-1 and Raf-1 on the induction of cellular senescence in doxorubicin or tamoxifen treated breast cancer cells

Escape from cellular senescence induction is a potent mechanism for chemoresistance. Cellular senescence can be induced in breast cancer cell lines by the removal of estrogen signaling with tamoxifen or by the accumulation of DNA damage induced by the chemotherapeutic drug doxorubicin. Long term culturing of the hormone-sensitive breast cancer cell line MCF-7 in doxorubicin (MCF-7/DoxR) reduced the ability of doxorubicin, but not tamoxifen, to induce senescence. Two pathways that are often upregulated in chemo- and hormonal-resistance are the PI3K/PTEN/Akt/mTOR and Ras/Raf/MEK/ERK pathways. To determine if active Akt-1 and Raf-1 can influence drug-induced senescence, we stably introduced activated ΔAkt-1(CA) and ΔRaf-1(CA) into drug-sensitive and doxorubicin-resistant cells. Expression of a constitutively-active Raf-1 construct resulted in higher baseline senescence, indicating these cells possessed the ability to undergo oncogene-induced-senescence. Constitutive activation of the Akt pathway significantly decreased drug-induced senescence in response to doxorubicin but not tamoxifen in MCF-7 cells. However, constitutive Akt-1 activation in drug-resistant cells containing high levels of active ERK completely escaped cellular senescence induced by doxorubicin and tamoxifen. These results indicate that up regulation of the Ras/PI3K/PTEN/Akt/mTOR pathway in the presence of elevated Ras/Raf/MEK/ERK signaling together can contribute to drug-resistance by diminishing cell senescence in response to chemotherapy. Understanding how breast cancers containing certain oncogenic mutations escape cell senescence in response to chemotherapy and hormonal based therapies may provide insights into the design of more effective drug combinations for the treatment of breast cancer.

Stability analysis of multi-compartment models for cell production systems

We study two- and three-compartment models of a hierarchical cell production system with cell division regulated by the level of mature cells. We investigate the structure of equilibria with respect to parameters as well as local stability properties for the equilibria. To interpret the results we adapt the concept of reproduction numbers, which is well known in ecology, to stem cell population dynamics. In the two-compartment model, the positive equilibrium is stable wherever it exists. In the three-compartment model, we find that the intermediate stage of differentiation is responsible for the emergence of an instability region in the parameter plane. Moreover, we prove that this region shrinks as the mortality rate for mature cells increases and discuss this result.

Expansion of adipose mesenchymal stromal cells is affected by human platelet lysate and plating density

The composition of mesenchymal stromal cells (MSCs) changes in the course of in vitro culture expansion. Little is known how these cell preparations are influenced by culture media, plating density, or passaging. In this study, we have isolated MSCs from human adipose tissue in culture medium supplemented with either fetal calf serum (FCS) or human platelet lysate (HPL). In addition, culture expansion was simultaneously performed at plating densities of 10 or 10,000 cells/cm(2). The use of FCS resulted in larger cells, whereas HPL significantly enhanced proliferation. Notably, HPL also facilitated expansion for more population doublings than FCS (43 ± 3 vs. 22 ± 4 population doubling; p < 0.001), while plating density did not have a significant effect on long-term growth curves. To gain further insight into population dynamics, we conceived a cellular automaton model to simulate expansion of MSCS. It is based on the assumptions that the number of cell divisions is limited and that due to contact inhibition proliferation occurs only at the rim of colonies. The model predicts that low plating densities result in more heterogeneity with regard to cell division history, and favor subpopulations of higher migratory activity. In summary, HPL is a suitable serum supplement for isolation of MSC from adipose tissue and facilitates more population doublings than FCS. Cellular automaton computer simulations provided additional insights into how complex population dynamics during long-term expansion are affected by plating density and migration.

Haematopoietic stem and progenitor cells in rheumatoid arthritis

RA is the prototypic chronic inflammatory disease, characterized by progressive articular cartilage and bone destruction. The systemic nature of RA is evidenced by the increased risk of atherosclerosis and lymphoproliferative disorders. Components of both the innate and adaptive immune system are implicated in the pathophysiology of the articular and extra-articular manifestations of the disease. A fundamental process in the onset of RA is the breakdown in self-tolerance. Accelerated ageing of immune cells (immunosenescence) appears to be a major mechanism favouring the disruption of tolerance. Telomere erosion, a hallmark of immunosenescence, is present in lymphoid (naïve and memory T cells) and myeloid (granulocytes) cells in RA. The premature ageing process also involves the haematopoietic stem and progenitor cells (CD34(+) HSPC), thus extending the RA immunopathogenesis to include early events in the shaping of the immune system. This review summarizes current concepts of HSPC ageing and its impact on immune regeneration, highlighting the phenotypic and functional similarities between elderly and RA HSPC.

DNA methylation pattern changes upon long-term culture and aging of human mesenchymal stromal cells

Within 2–3 months of in vitro culture-expansion, mesenchymal stromal cells (MSC) undergo replicative senescence characterized by cell enlargement, loss of differentiation potential and ultimate growth arrest. In this study, we have analyzed DNA methylation changes upon long-term culture of MSC by using the HumanMethylation27 BeadChip microarray assessing 27 578 unique CpG sites. Furthermore, we have compared MSC from young and elderly donors. Overall, methylation patterns were maintained throughout both long-term culture and aging but highly significant differences were observed at specific CpG sites. Many of these differences were observed in homeobox genes and genes involved in cell differentiation. Methylation changes were verified by pyrosequencing after bisulfite conversion and compared to gene expression data. Notably, methylation changes in MSC were overlapping in long-term culture and aging in vivo. This supports the notion that replicative senescence and aging represent developmental processes that are regulated by specific epigenetic modifications.

Replicative senescence-associated gene expression changes in mesenchymal stromal cells are similar under different culture conditions

BACKGROUND

Research on mesenchymal stromal cells has created high expectations for a variety of therapeutic applications. Extensive propagation to yield enough mesenchymal stromal cells for therapy may result in replicative senescence and thus hamper long-term functionality in vivo. Highly variable proliferation rates of mesenchymal stromal cells in the course of long-term expansions under varying culture conditions may already indicate different propensity for cellular senescence. We hypothesized that senescence-associated regulated genes differ in mesenchymal stromal cells propagated under different culture conditions.

DESIGN AND METHODS

Human bone marrow-derived mesenchymal stromal cells were cultured either by serial passaging or by a two-step protocol in three different growth conditions. Culture media were supplemented with either fetal bovine serum in varying concentrations or pooled human platelet lysate.

RESULTS

All mesenchymal stromal cell preparations revealed significant gene expression changes upon long-term culture. Especially genes involved in cell differentiation, apoptosis and cell death were up-regulated, whereas genes involved in mitosis and proliferation were down-regulated. Furthermore, overlapping senescence-associated gene expression changes were found in all mesenchymal stromal cell preparations.

CONCLUSIONS

Long-term cell growth induced similar gene expression changes in mesenchymal stromal cells independently of isolation and expansion conditions. In advance of therapeutic application, this panel of genes might offer a feasible approach to assessing mesenchymal stromal cell quality with regard to the state of replicative senescence.