Optimality conditions for cell-fate heterogeneity that maximize the effects of growth factors in PC12 cells

Recently, the heterogeneity that arises from stochastic fate decisions has been reported for several types of cancer-derived cell lines and several types of clonal cells grown under constant environmental conditions. However, the relation between this stochasticity and the responsiveness to extracellular stimuli remains largely unknown. Here we focused on the fate decisions of the PC12 cell line, which was derived from rat pheochromocytoma, and is a model system to study differentiation into sympathetic neurons. Whereas epidermal growth factor (EGF) stimulates the proliferation of populations of PC12 cells, nerve growth factor (NGF) promotes the differentiation of neurites to neuron-like cells. We found that phenotypic heterogeneity increased with time at several surrounding serum concentrations, suggesting stochastic cell-fate decisions in single cells. We made a simple mathematical model assuming Markovian transitions of the cell fates, and estimated the transition rates based on Bayes' theorem. The model suggests that depending on the serum concentration, EGF (NGF) even directs differentiation (proliferation) at the single-cell level. The maximum effects of the growth factors were ensured when the transition rates were appropriately controlled by the serum concentration to produce a nonextremal, moderate amount of cell-fate heterogeneity. Our model was validated by the experimental finding that the means and variances of the local cell densities obey a power-law relationship. These results suggest that even when efficient responses to growth factors are observed at the population level, the growth factors stochastically direct the cell-fate decisions in different directions at the single-cell level.