Individual fates of mesenchymal stem cells in vitro
© Krinner et al. 2010
Received: 27 October 2009
Accepted: 27 May 2010
Published: 27 May 2010
Skip to main content
© Krinner et al. 2010
Received: 27 October 2009
Accepted: 27 May 2010
Published: 27 May 2010
In vitro cultivated stem cell populations are in general heterogeneous with respect to their expression of differentiation markers. In hematopoietic progenitor populations, this heterogeneity has been shown to regenerate within days from isolated subpopulations defined by high or low marker expression. This kind of plasticity has been suggested to be a fundamental feature of mesenchymal stem cells (MSCs) as well. Here, we study MSC plasticity on the level of individual cells applying a multi-scale computer model that is based on the concept of noise-driven stem cell differentiation.
By simulation studies, we provide detailed insight into the kinetics of MSC organisation. Monitoring the fates of individual cells in high and low oxygen culture, we calculated the average transition times of individual cells into stem cell and differentiated states. We predict that at low oxygen the heterogeneity of a MSC population with respect to differentiation regenerates from any selected subpopulation in about two days. At high oxygen, regeneration becomes substantially slowed down. Simulation results on the composition of the functional stem cell pool of MSC populations suggest that most of the cells that constitute this pool originate from more differentiated cells.
Individual cell-based models are well-suited to provide quantitative predictions on essential features of the spatio-temporal organisation of MSC in vitro. Our predictions on MSC plasticity and its dependence on the environment motivate a number of in vitro experiments for validation. They may contribute to a better understanding of MSC organisation in vitro, including features of clonal expansion, environmental adaptation and stem cell ageing.
The generation and maintenance of replenishing tissues relies on an appropriately regulated balance between self-renewal and differentiation within a relatively small population of adult stem cells. According to the common stem cell paradigm this balance can be explained assuming a strict differentiation hierarchy and irreversible fate decisions [1, 2]. However, the organisation of stem cell populations is strongly influenced by environmental factors such as specific cell-cell interactions, growth factor and oxygen supply, as well as the geometry and mechanical properties of the local environment [3, 4]. Accordingly, it has been suggested that stemness represents a particular regulatory cell state rather than an entity and that this state may be approached in principle by any cell [5, 6]. Supporting these ideas, recent experimental results in hematopoietic systems demonstrated that stem cell populations can actually regenerate from more differentiated subpopulations [7, 8]. Currently, there is an ongoing debate on fundamental dynamics underlying this kind of cell plasticity. In particular, it remains open whether de-differentiation is prerequisite to lineage changes. A thorough understanding of this phenomenon is expected to make an important contribution to the development of novel therapeutic strategies for treating degenerative disease, injury and neoplasia.
Mesenchymal stem cells (MSCs) are multi-potent cells that persist in adult life in some tissue types, such as bone-marrow stroma, fat, skeletal muscle, and synovium without loosing their capacity to proliferate and differentiate [9, 10]. Under appropriate culture conditions, they can multiply and transform into specialized cell types in vitro. Plasticity of MSCs of the 3T3 T type linked to de-differentiation has already been demonstrated in the Eighties . More recently, also differentiation of adult human MSC was found to be at least partially reversible . In fact plasticity has been suggested to represent a fundamental feature of MSC .
Recently, we have introduced a multi-scale computer model of MSC expansion, lineage commitment and differentiation which consistently explains a panel of experimental results regarding the oxygen dependence of these processes and predicts optimal culture conditions . This model utilises the concept of noise-driven stem cell differentiation  which is based on the functional stem cell approach to tissue organisation by Roeder & Loeffler [5, 16]. According to this concept, MSC plasticity bases on permanent fluctuations of the differentiation state of each individual cell, which enables more differentiated cells to re-gain stem cell properties and subsequently to switch lineage (details see below).
Here we aim at quantitative predictions on MSC organisation in vitro based on our former results. For this purpose we performed "experiments in silico" using our novel multi-scale model. We monitored the fates of individual MSCs under different culture conditions. Linking intracellular regulation of the differentiation state to cell biomechanics our computer simulations provide insight into possible mechanisms of how cell-cell and cell-substrate interaction can affect stem cell functionality. Thereby, our computer simulations were designed as MSC protocols in silico such that they can be directly tested in vitro.
In the following we first give a brief description of the model of MSC organisation in vitro introduced by Krinner et al.  and provide the experimentally validated data set used throughout this study. Subsequently, we present our simulation results on MSC plasticity and discuss the potential and the limits of our approach.
In our model cell differentiation is defined as the loss of stem cell properties. Cell differentiation is quantified by a continuous state variable α that can adopt values between zero (full stem cell competency) and one (completely differentiated cell). Each value of α may represent a set of regulatory network activation patterns. From the molecular point of view, α may depend on the abundance and sub-cellular localization of proteins and RNAs, as well as other types of signalling and metabolic molecules . Cell differentiation is assumed to occur independently of cell proliferation .
MSC differentiation involves lineage-priming . This process implies particular cellular decisions, which can be modelled considering a second state variable . We here assume that differentiation and de-differentiation dynamics do not depend on these decisions. However, a switch from one into another specific lineage may require a defined degree of stemness as suggested for the chondrogenic lineage . In this case, differentiation stabilises lineages and the described capability of de-differentiation is synonymous with MSC plasticity in general. We here focus on that kind of MSC plasticity.
where σ0 denotes the fluctuation strength in stem cell states and f is a Hill function approaching 0 and 1 at low and high pO2, respectively.
Cell proliferation is assumed to depend on the differentiation state α of a cell. In our model, it is restricted to intermediate differentiation states αp with: 0 < αs < αp < αd < 1 (Figure 1b). These states are termed 'progenitor states' in the following. For these proliferative states we assume an identical doubling rate r = 1/τ and average growth time τ. 'Stem cells' (α < αs) and 'differentiated cells' (α > αd) do not proliferate. The state fluctuations cause the cells to switch frequently between proliferative and non-proliferative state, which results in an effective average growth time larger than τ.
In order to simulate the spatio-temporal dynamics of MSC populations we use an IBM where the cells are modelled as elastic adhesive spheres . We assume that the cell volume in suspension cannot be smaller than a minimum value V0. A cell can move actively by migration and passively by being pushed, it can deform, adhere to other cells or a substrate, and it can grow and divide. A proliferating cell divides if its volume has grown to twice the volume V0.
In the first term on the right hand side νi denotes the Poisson's ratio of the interaction partner i (i = 1,2), Ei its Young modulus, Ri its radius (substrate radius R = ∞) and δ the surface deformation. The second term models adhesion proportional to the Hertz contact area, where ε12 is the anchorage given as adhesion energy per unit area.
Cell proliferation is modelled assuming a two phase cell cycle: During the interphase, a cell doubles its volume by stochastic increments. During the mitotic phase, a cell divides into two daughter cells of equal volume. This growth process results in an approximately Γ-distributed growth time τ of the cells . A cell undergoes a growth arrest if the sum of the deformation forces on it exceeds a critical value Fc.
where the sums run over all neighbouring cells j in direct contact to cell i. F ij Hertz denotes the Hertz force between cell i and cell j and Fi stoch the stochastic Langevin force on cell i. The friction coefficients γis and γij describe friction between cell i and the substrate and between cell i and cell j, respectively. These coefficients are assumed to be proportional to the respective contact areas. Details can be found in .
with transition probability and constant randomization rate R. Transition times ϑ(α) from an initial α into the regimes of stem cells (α < αs) or differentiated cells (α > αd) were computed using an absorbing boundary approach .
Our model of MSC differentiation dynamics depends on parameters describing intracellular regulation; the randomization rate R, the stem cell state fluctuation strength σ0, the parameters of the Hill function (n and k) and those specifying the proliferation rate (r and αs with αd = 1-αs). The IBM of spatio-temporal organisation of growing MSC populations depends on parameters specifying cell-cell and cell-substrate interaction, as the Poisson's ratio, the Young modulus, and the friction constants. Combining these models in a particular application one has to adjust a large parameter set.
Parameter set used in the simulations
Randomization Rate R
2.5 × 104 s-1
Stem Cell State Fluctuation Strength σ0
Hill Coefficient n
Dissociation Constant k
Differentiation Threshold αs (αd = 1-αs)
Minimal Cell Radius R0
Minimal Cell Volume V0
Proliferation Rate r = 1/τ
Young Modulus E
Contact Inhibition Threshold Fmax
1 × 10-9 N
Poisson's Ratio ν
Friction Coefficients γij, γis
3 × 107 Ns/m3 , 1 × 1011 Ns/m3 *
Cellular Diffusion Coefficient DCell
4 × 10-12 cm2/s
Cell-Cell Anchorage ε
6 × 10-5 N/m
Cell-Plane Anchorage ε
6 × 10-5 N/m
Quiescence Threshold Fq
Using the IBM the fates of individual cells in growing populations can be monitored. We simulated individual α-trajectories and compared the cell differentiation dynamics at low (5%) and high (20%) oxygen concentrations. The genealogies of two selected clones in α space are shown in Figure 1c and 1d, for low and high oxygen, respectively.
Our results demonstrate that at low oxygen a frequent exchange between the subpopulations occurs on a time scale of about 2 days. At high oxygen the average transition time for stem cells into the pool of differentiated cells increases to about 4 days. Transition times for differentiated cells into the stem cell pool at high oxygen are much larger (>100 days), indicating quasi-deterministic cell differentiation behaviour. We confirmed our results using the master equation approach. In Figure 2c the fraction of cells having entered the stem cell pool at 20% pO2 is shown as a function of the initial α value and the simulation time. Only in this particular case, the fraction of absorbed cells grows too slowly to calculate the average transition times. In the three other cases, they were computed with high precision (less than 10-12 of all cells remain to be absorbed).
Since stem cell states are more easily accessible at low oxygen compared to high oxygen we predict MSC plasticity to be more pronounced under these conditions.
In vitro validation of the above results would require single cell tracking of MSCs and techniques to identify the differentiation state of the tracked cells. Currently, considerable effort is taken in order to establish tracking techniques for stem cell systems [25, 26]. Unfortunately, MSCs are particularly hard to track, because they tend to aggregate; a phenomenon known as mesenchymal condensation [27, 28]. Thus, in the following we present results on MSC plasticity as seen on the population level which can be validated in simpler experimental setups.
Chang et al.  studied how fast the distribution of differentiation marker expression within a cell population regenerates from subpopulations with defined expression level. They performed the following experiment: a population of precursor cells was generated under standard conditions and characterised by the expression level of a particular differentiation marker. Subpopulations of cells with defined expression levels of the differentiation marker were separated. These subpopulations were cultivated under standard conditions and regeneration of the distribution of expression levels in the population was monitored over time by FACS.
These simulation results implicate that if regeneration refers to the growth of a few large clones, as in the case of differentiated cells at high oxygen, the effect of contact inhibition becomes more relevant for population regeneration. The α-distribution in large clones significantly differs from that of a low-density culture. Moreover, due to the increased number of differentiated cells, these populations show a lower CFU capacity (compare ).
In general, 'self-renewal' of the stem cell population appears in our model as steady occupancy of stem cell states due to a particular population dynamics. Thus, additional information on MSC organisation in vitro can be obtained by performing the regeneration experiments described above in parallel for all subpopulations. Splitting the mother population into a number of subpopulations according to the expression of a differentiation marker, applying the 'regeneration protocol' suggested above to each of these subpopulations and quantifying the number of stem cells in each subpopulation after a fixed regeneration time would allow to quantify the fraction of stem cells in a MSC population descending from a particular subpopulation.
Recent experimental findings indicate that cells can regain stem cell properties under defined environmental conditions. These results challenge the commonly agreed stem cell paradigm. This paradigm treats 'stemness' as a fixed property intrinsic to stem cells and assumes a deterministic and irreversible differentiation scenario for each cell . As an alternative, novel concepts of functional stem cells have been developed that assign the interaction between cells and their growth environment a greater emphasis [5, 13, 16]. Treating stemness no longer as a fixed property, these concepts do not exclude certain preferred trends in the differentiation sequence, but allow reversible developments for individual cells.
We here provided the first quantitative predictions on the environmental dependent organisation of MSCs in vitro applying this novel concept. We predicted: i) the average transition times of individual cells into stem cell and differentiated states, ii) the time scales of the regeneration of the distribution of differentiation marker expression in a MSC population from subpopulations of stem and unspecific differentiated cells, and iii) the origin of the cells forming the in vitro stem cell pool of MSC. Moreover, we predicted that all these properties depend on the environment. Our results also provide estimates of the time scales of MSC adaptation to changed environmental conditions. They are in good agreement with experimental findings on MSC adaptation to low oxygen [30–32]. Particularly the work of Tang et al.  and Volkmer et al.  strengthens our modelling approach because the experimentally observed improvement of the functional competence of an entire MSC population within less than 24 hours can only hardly be explained by the expansion of residual stem cells as suggested by pedigree models.
In all our simulations, we considered an oxygen dependence of the state fluctuations.
In contrast, biophysical features, as cell-cell and cell-substrate interactions, were assumed to affect the regenerative potential of the MSC by interfering with their proliferation control mechanisms only. A direct feedback of these interactions on the noise amplitudes was not considered. However, recent results demonstrate that lineage specification and proliferation of MSC populations can be triggered by substrate elasticity  and substrate micro-structure . Thus, we here suggest performing the proposed experiments on MSC plasticity on substrates that vary with respect to their elasticity and microstructure. These experiments would provide information on whether mechano-signalling can affect the kinetics of state transitions in MSCs and thus, can be used to time regeneration processes in vitro.
Our results on the composition of the stem cells pool suggest that most of the stem cells in MSC populations expanding in vitro originate from progenitor states. Thus, their mother cells underwent differentiation and de-differentiation processes and were proliferative active. Recent experimental results suggest that these cellular activities result in changes in the cellular phenotype called stem cell ageing . A model that consistently describes these phenomena is currently lacking.
Most of our results could be validated by in vitro experiments on the population level. A number of suggestions were given in the text. However, more detailed studies would require tracking of individual cell fates in a single expanding MSC population. Such experiments would provide additional information on cell-cell communication in the expanding population, which was suggested to impact MSC expansion . As already mentioned above, the tracking of MSC involves particular problems. Long term monitoring of MSC fates will require therefore sophisticated marker systems for both the clonal origin and the differentiation state of the cells. A number of stem cell and differentiation markers of MSC have been suggested. Good candidates are early transcription factors [37, 38].
Long-term fluctuations in differentiation marker expression in single cells would directly proof our concept of noise-driven stem cell organisation. For the generality of our concept, we expect such fluctuations to underlie somatic stem cell organisation independent of tissue and species.
The impact of these fluctuations may vary between different stem cell systems according to functional requirements . Thus, individual stem cell systems may appear as more or less hierarchical organised. The MSC system may exhibit a pronounced flexibility, in order to be capable of instantaneous fate decisions in the course of development and in case of injury [39, 40].
Understanding single cell behaviour is prerequisite to unveil general principles of the organisation of stem cell populations. Stem cell maintenance, expansion and environmental adaptation may in particular rely on single cell plasticity. Currently only limited data on the in vitro plasticity of individual stem cells are available. We here presented for the first time quantitative simulation results on in vitro MSC plasticity applying our novel concept of noise driven stem cell differentiation. Thereby we demonstrate the suitability of the IBM approach for studying these phenomena. Challenging current views on stem cell organisation, our results predict a highly dynamic stem cell pool, whose maintenance involves permanent de-differentiation events.
The work was supported by the BMBF grant number 0313836 and 0313081F.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.