Stem Cells 2002;20:119-138
www.StemCells.com
© 2002 AlphaMed Press
Ligand/Receptor Signaling Threshold (LIST) Model Accounts for gp130-Mediated Embryonic Stem Cell Self-Renewal Responses to LIF and HIL-6
Sowmya Viswanathana,
Tania Benatara,
Stefan Rose-Johnb,
Doug A. Lauffenburgerc,
Peter W. Zandstraa
a Institute of Biomaterials and Biomedical Engineering, University of Toronto, Toronto, Ontario, Canada;
b Department of Biochemistry, Christian-Albrechts-Universität zu Kiel, Kiel, Germany;
c Biotechnology Process Engineering Center, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
Key Words. Embryonic stem cells • gp130 • Ligand-receptor signaling complexes • Thresholds • Self-renewal • LIF
Peter W. Zandstra, Ph.D., Institute of Biomaterials and Biomedical Engineering, Room 407, Roseburgh Building, 4 Taddle Creek Road, Toronto, Ontario, M5S 3G9, Canada. Telephone: 416-978-8888; Fax: 416-978-4317; e-mail: peter.zandstra{at}utoronto.ca; website: http://stemcell.ibme.utoronto.ca
 |
ABSTRACT
|
|---|
We previously demonstrated that embryonic stem (ES) cell self-renewal required sustained signaling by leukemia inhibitory factor (LIF) in a concentration-dependent manner, allowing us to hypothesize that thresholds in ligand-receptor signaling modulate stem cell differentiation control. To test this hypothesis, we have experimentally and computationally compared the abilities of two gp130-signaling cytokines (LIF and Hyper-interleukin-6 [HIL-6]) to sustain ES cell self-renewal. Quantitative measurements of ES cell phenotypic markers (stage-specific embryonic antigen-1 and E-cadherin), functional assays (alkaline phosphatase activity and embryoid body formation efficiency), and transcription factor (Oct-4) expression over a range of LIF and HIL-6 concentrations demonstrated a superior ability of LIF to maintain ES cell pluripotentiality at higher concentrations (
500 pM). Additionally, we observed distinct qualitative differences in the ES cell self-renewal dose response profiles between the two cytokines. A computational model permitted calculation of the number of signaling complexes as a function of receptor expression, ligand concentration, and ligand/receptor-binding properties, generating predictions for the degree of self-renewal as a function of cytokine concentration by comparison of these calculated complex numbers to experimentally determined threshold cytokine concentrations. Model predictions, consistent with experimental data, indicated that differences in the potencies of these two cytokines were based primarily on differences in receptor-binding stoichiometries and properties. These results support a ligand/receptor signaling threshold model of ES cell fate modulation through appropriate types and levels of cytokine stimulation. Insights from these results may be more generally applicable to tissue-specific stem cells and could aid in the development of stem cell-based technologies.
 |
INTRODUCTION
|
|---|
Human embryonic stem (ES) cells [14] have the potential to serve as renewable sources of functional tissue-specific cells in clinical tissue-replacement therapies. A prerequisite for the use of these and other stem cells in a clinical context is a thorough understanding of the mechanisms that regulate their decision to either self-renew and remain in the stem cell compartment or differentiate and exit the stem cell compartment. Recent evidence from hematopoietic stem cells (HSC) [58] and murine ES cells [911] suggests that stem cell fate may be amenable to exogenous control by appropriate types and levels of ligand-mediated receptor stimulation.
We recently proposed a ligand/receptor signaling threshold (LIST) model for stem cell fate determination where parametric values of the appropriate cell-surface signaling cytokine-receptor complexes determined whether a gene program associated with a particular cell fate was activated or not [12]. According to this conceptual model, a threshold number of cell-surface cytokine-receptor complexes determines whether or not a stem cell undergoes a self-renewal or differentiation response. Although this model is consistent with a large body of data, its applicability to stem cell differentiation control has not been rigorously examined.
Murine ES cells [13, 14] respond to cytokines from the interleukin-6 (IL-6) family, which maintain them in an undifferentiated state by signaling through the gp130 signal-transducing receptor component (a molecule that is also important in regulating other systems such as HSC [5, 8, 1517] and primordial germ cells [1821]). Members of the IL-6 family of cytokines, including leukemia inhibitory factor (LIF) [9, 10], IL-6 in combination with the soluble IL-6 receptor (sIL-6R) [17, 22] (ES cells lack membrane-bound IL-6 receptor [22]), cardiotrophin-1 (CT-1) [23], ciliary neurotrophic factor (CNTF) [24] and oncostatin M (OSM) [25, 26], can all maintain ES cells in an undifferentiated state for prolonged periods in vitro, albeit their efficacy varies [23, 25, 26]. Members of the IL-6 family of cytokines signal via multimeric receptor complexes usually comprised of a ligand-specific
-receptor subunit (IL-6R, IL-11R, and CNTFR) and gp130, the common signal-transducing subunit [27, 28]. IL-6 induces gp130 homodimerization [29], whereas CNTF [30] and LIF [3032] signal via heterodimers of gp130 and LIF receptors (LIFR). Membrane-bound
-receptors can be functionally replaced by soluble forms that lack the cytoplasmic domain, but may still function agonistically, as in the case of sIL-6R [27, 33, 34].
We had previously developed quantitative phenotypic (stage-specific embryonic antigen-1 [SSEA-1] [35, 36]) and functional (alkaline phosphatase [ALP] activity [3739] and the capacity of ES cells to form complex cell aggregates called embryoid bodies (EBs) [4043]) assays to identify ES cells and their immediately differentiated progeny [11]. These assays, along with measurements of E-cadherin (a cell adhesion molecule important in early development [44, 45]) and Oct-4 (a transcription factor critical for maintenance of ES cell pluripotentiality [46, 47]) expression described herein, provide the quantitative tools necessary to conduct mechanistic investigations of factors that govern early ES cell fate decisions.
We have taken advantage of differences in the signaling mechanisms of the two IL-6-type cytokines, LIF and Hyper-IL-6 (HIL-6), to examine the validity of our LIST model in regulating stem cell fate. HIL-6, a fusion protein formed by covalently linking human IL-6 to a truncated human sIL-6R by a flexible glycine linker [48], is active at 100-1,000-fold lower concentrations than unlinked IL-6 and IL-6R on responsive cells because it has increased stability and may be more resistant to internalization [5, 4850]. Our results suggest that cell-surface interactions between signaling cytokines and their corresponding receptors provide a mechanistic explanation of the dose dependency of ES cell self-renewal responses to exogenous cytokine stimuli. Differences in the abilities of LIF and HIL-6 to sustain ES cell pluripotentiality can be attributed to the distinct stoichiometries involved in forming functional cell-surface signaling ligand-receptor complexes. Our LIST model explains the ability of LIF to maintain ES cells in an undifferentiated state over a relatively wide range of concentrations, while HIL-6 is able to act equivalently over a relatively narrow range of doses, a result also recently shown to be true for HSC [5]. Our results also confirm that the actions of these cytokines cannot be attributed to differential mitogenic requirements of these cells [11], and that these cytokines may be acting, at least in part, by modulating the differential survival of particular subset(s) of cells. Our ability to correlate differences in cell-specific and measurable parameters (for example, receptor occupancy) to ES cell responses has provided us a mechanistic insight into the regulation of stem cell fate, a critical first step in being able to exploit these cells for clinical gains.
 |
MATERIALS AND METHODS
|
|---|
Cells and Cytokines
Murine ES cells (CCE) [42] were maintained at 37°C in humidified air with 5% CO2 as previously described [11]. All cells were used within 15 passages of the initial thawing. HIL-6 [48] was a generous gift from Stefan Rose-John (Christian-Albrechts-Universität; Kiel, Germany; http://www.uni-kiel.de/index-e.html). HIL-6 concentration was independently established by a human IL-6 enzyme-linked immunosorbent assay kit (CLB; Amsterdam, The Netherlands; http://www.clb.de). STO embryonic fibroblast cells were maintained in Dulbecco's modified Eagle's medium supplemented with 50 U/ml penicillin and 50 µg/ml streptomycin (Invitrogen Life Technologies; Rockville, MD; http://www.lifetech.com), 10% cosmic calf serum (HyClone; Logan, UT; http://www.hyclone.com), and 2 mM L-glutamine (Life Technologies).
Analysis of ES Cell Viability and Proliferation
Single cell suspensions of ES cells were initiated in triplicate in three gelatin-coated 24-well tissue culture plates at 7,500 cells/well and supplemented with indicated LIF and HIL-6 concentrations. On day 1, adherent cells were washed in phosphate-buffered saline (PBS), lysed in 0.1% Triton 100x, and incubated with 0.5 µM of SYTOX® Green (Molecular Probes; Eugene, OR; http://www.probes.com) for 10 minutes at room temperature. SYTOX® Green exhibited strong fluorescent enhancement when bound to cellular nucleic acid, a standard curve related cell number to fluorescent intensity (data not shown). Fluorescence emission was measured on a SPECTRAmax® spectrofluorometer (490 nm excitation, 524 nm emission with a cut-off filter at 515 nm) with SOFTmax® PRO software. A blank well containing all reagents except the cells served as a background control. Cell numbers were measured on days 2 and 3 in a similar manner (Fig. 1A
). ES cells were also cultured in 75 cm2 flasks for 5 days at indicated doses of LIF and HIL-6 and then replated in 24-well tissue culture plates (10,000 cells/well) at the corresponding ligand concentrations on day 5 in order to measure adherent cell numbers on days 6-8 as before (Fig. 1B
).

View larger version (24K):
[in this window]
[in a new window]
|
Figure 1. Experimental protocol. Single ES cell suspensions were either (A) directly seeded in 24-well gelatin-coated plates and cell numbers were determined consecutively on d 1-3 or (B) cultured in 75 cm2 dishes for d 1-5 and then seeded in 24-well gelatin-coated plates on d 5 and analyzed subsequently on d 6-8. (C) Cells were also seeded at equivalent densities on d 1, 3, 5, and 7 in 25 cm2 or 75 cm2 dishes at different doses of LIF and HIL-6 and harvested to assess differentiation status.
|
|
Flow Cytometric Analysis of SSEA-1, E-Cadherin, and Intracellular Oct-4 Expression
ES cells were seeded at 2 x 106 cells/25 cm2 flask or at 6 x 106 cells/75 cm2 flask at indicated LIF and HIL-6 concentrations and subcultured every other day for 7 days using 0.25% trypsin-EDTA (Sigma; St. Louis, MO; http://www.sigmaaldrich.com) to generate single-cell suspensions. Cultures were reinitiated on alternate days with 1 x 106 cells/25 cm2 flask or 3 x 106 cells/75 cm2 flask at different doses of LIF and HIL-6 to control for cell density (Fig. 1C
). An aliquot of 1 x 106 cells was washed in ice cold Hanks' HEPES-buffered saline solution (Life Technologies) containing 2% fetal calf serum (HF) and incubated with
-mouse CD16/CD32 receptor monoclonal antibody at 1 µg/100 µl (Pharmigen; San Diego, CA; http://www.pharmingen.com) to block nonspecific binding of immunoglobulins to the mouse Fc
III/II receptors. Cells to be stained for Oct-4 were fixed with 100 µl of IntraPrepTM Permeabilization Reagent 1 (Immunotech; Westbrook, ME) for 15 minutes at room temperature, washed 1x with HF, permeabilized with IntraPrepTM Permeabilization Reagent 2 and incubated with a 1:100 dilution mouse
-mouse Oct3/4 monoclonal antibody (Transduction Laboratories; Lexington, KY) for 15 minutes at room temperature. Cells were washed 2x with 1 ml of HF, followed by staining with 1:100 dilution of fluorescein isothiocyanate-conjugated (FITC-conjugated) goat anti-mouse IgG antibody (Sigma). Finally, the cells were washed 2x with HF, and resuspended in HF for analysis on a flow cytometer (XL; Beckman-Coulter; Miami, FL; http://www.coulter.com) using the ExpoADCXL4 software (Beckman-Coulter). Cells analyzed for SSEA-1 and E-cadherin were incubated first with a 1:10 dilution of monoclonal anti-SSEA-1 (Developmental Studies Hybridoma Bank; Iowa City, IA; http://www.uiowa.edu/~dshbwww) for 40 minutes at 4°C, washed 2x with ice-cold HF and incubated with a 1:500 dilution of phycoerythrin-conjugated goat anti-mouse IgM antibody (Immunotech). Cells were once again washed 2x with HF, resuspended in 100 ml of HF and subsequently stained with a 1:10 monoclonal anti-uvomorulin (rat IgG1 isotype) at 3.75 µg/100 µl (Sigma) for 40 minutes at 4°C, washed 2x with ice-cold HF, and incubated with a 1:100 dilution of FITC-conjugated goat anti-rat IgG antibody (Sigma) for another 40 minutes at 4°C. Finally, cells were washed 2x with HF, stained with 1 mg/ml 7-AAD (Molecular Probes) and resuspended in HF for analysis on a flow cytometer. Positive staining was defined as the emission of a level of fluorescence that exceeded levels obtained by >99.5% of cells from the control population (cultured at 500 pM of LIF) stained with only the secondary antibodies.
SSEA-1 Selection
Petri dishes (Optilux 1001; Becton Dickinson; Franklin Lakes, NJ; http://www.bd.com) were coated with 50 µg of goat anti-mouse IgM (Jackson ImmunoResearch Labs; Westgrove, PA; http://www.jacksonimmuno.com) diluted in panning buffer (0.05 M Tris, 0.15 M NaCl to pH = 9.5) for either 40 minutes at room temperature, or overnight at 4°C. Petri dishes were washed 4x with 3 ml of Hanks' balanced salt solution with phenol red (GIBCO; Rockville, MD), followed by 2x washes with HF, leaving the last wash on for 5 minutes at room temperature. Dishes were coated with monoclonal anti-SSEA-1 hybridoma supernatant (Developmental Studies Hybridoma Bank) diluted 1:2 with HF for 40 minutes at room temperature, and then washed. Between 1 and 2 x 107 cells (ES or STOs) were plated onto SSEA-1-coated dishes in HF for 1 hour at 4°C with a gentle swirl at 30 minutes, and an aliquot of cells was removed for flow cytometric analysis before the panning. Nonadherent cells were removed from the first two washes with HF and counted. Adherent cells were removed after vigorous pipetting of the remaining cells and counted. Aliquots from the two fractions were analyzed for their SSEA-1 and Oct-4 content in parallel by flow cytometry as above. It was not possible to analyze SSEA-1 and Oct-4 expression simultaneously as SSEA-1 expression was reduced when ES cells were permeabilized to stain for intracellular Oct-4. Oct-4 protein expression of the two fractions was also independently analyzed by Western analysis.
Functional Assays
The EB formation capacity assay and the ALP activity assay were performed as previously described [11]. We had previously demonstrated that EB formation capacity of ES cells correlated with other differentiation assays including SSEA-1 expression and ALP activity [11].
Oct-4 Protein Detection
ES cells were cultured as described above. Aliquots of 2 x 106 cells were washed 3x in PBS, followed by lysis in Solution A (10 mM HEPES, 10 mM KCl, 0.1 mM EDTA, 0.1 mM EGTA, 1 mM dithiothreitol, 0.5 mM phenylmethylsulfonyl fluoride (PMSF) plus protease (aprotinin and leupeptin [BioShop; Toronto, ON; http://www.bioshopcanada.com]) and phosphatase inhibitors (sodium fluoride, sodium orthovandate, sodium pyrophosphate [Sigma]) for 10 minutes at 4°C. Cytoplasmic proteins were removed after centrifugation at 13,000 g for 10 seconds. Nuclear proteins were eluted by resuspending in Solution B (20 mM HEPES, 0.4 M NaCl, 1 mM EGTA, 1 mM EDTA, 0.1 mM PMSF, with protease and phosphatase inhibitors) for 30 minutes at 4°C, followed by centrifugation at 13,000 g for 30 seconds. Quantification of protein levels was performed using Plus Protein Assay Reagent (Pierce; Rockford, IL; http://www.piercenet.com). Equal amounts of protein (28 µg) were separated on SDS-PAGE, and transferred onto polyvinylidene difluoride membranes. Determination of Oct-3/4 levels was performed by incubation of membranes with mouse anti-mouse Oct-3 antibody (Transduction Laboratories) in buffer containing 5% milk in tris buffered saline solution containing 0.1% Tween 20 (TBST), followed by incubation with goat anti-mouse IgG-HRP (Sigma). Proteins were visualized by enhanced chemiluminescence (ECL) using a FluorchemTM Imaging System (Alpha Innotech; San Leandro, CA; http://www.alphainnotech.com).
Statistical Analysis
All data are reported as mean ± standard deviation. Experiments measuring changes in ES cell marker expressions were repeated two to five times, and statistical significance was assessed using paired Student's t-test. Experiments measuring ES cell survival and proliferation were conducted over five times; a representative experiment using SYTOX® Green is shown.
 |
RESULTS
|
|---|
ES Cell Survival and Proliferation
There were no significant differences in ES cell growth rates when cultured at different concentrations of LIF and HIL-6 over a period of 3 days (Fig. 2A
). Survival efficiencies (percentage of cells relative to input) were comparable at all doses of LIF and HIL-6 after 1 day in culture, when most cells were predominantly undifferentiated (cells had high SSEA-1 and E-cadherin expression, data not shown) (Fig. 2B
). Cell density limitations prevented us from extending this assay to 7 days when downregulation of ES cell markers is known to occur and cells could not be subcultured as we were interested in measuring growth rates. To circumvent this problem, we cultured ES cells at different doses of LIF and HIL-6 in 75 cm2 tissue-culture flasks for 5 days and replated cells at low densities in 24-well plates (Fig. 1B
) to measure their growth and survival from day 6 to 8. Even after 6 to 8 days of differentiation, growth rates at different ligand concentrations were not statistically different (Fig. 2C
), however, there were significant survival differences (p < 0.05) between replated cells exposed to low (
5 pM) and high (
500 pM) ligand concentrations (Fig. 2D
). The first day after replating, only 13% of the cells survived in the absence of LIF or HIL-6, whereas survival rates were up to 65% and 40% at 500 pM LIF and HIL-6 respectively (p < 0.05). Cell survival efficiencies were greater in the presence of higher LIF concentrations (
50 pM) than at equivalent HIL-6 concentrations (p < 0.05) suggesting that LIF was more potent than HIL-6 in favoring survival of replated cells.
ES Cell Differentiation
There were dramatic changes in ES cell morphology at different doses of LIF and HIL-6 over a 7-day period (Fig. 3
) when subcultured on alternate days (to ensure that seeding cell density numbers were equivalent at different ligand concentrations [Fig. 1C
]). Undifferentiated ES cells formed homogeneous compact, spherical colonies, while differentiated cells were flattened and expressed a heterogeneous morphology. Differentiation in response to different concentrations of LIF and HIL-6 was quantitatively measured by several, independent assays of ES cell differentiation and, as expected, undifferentiated ES cells maintained above 50 pM LIF concentrations had significantly higher SSEA-1 and E-cadherin expression, ALP activity, and EB formation capacity than those maintained at lower LIF concentrations (
5 pM) (p < 0.05) (Figs. 4A-4D
). Surprisingly, ES cells did not show a similar monophasic dose dependency on different concentrations of HIL-6 (Figs. 4A-4D
). ES cells maintained at 50 pM HIL-6 had significantly higher ALP activity and EB formation capacity than those maintained in the absence of any ligand or in the presence of 500 pM of HIL-6 (p < 0.10). Further, there were no significant differences (p > 0.10) in the SSEA-1 and E-cadherin expression, ALP activity, and EB formation capacity between cells maintained in the presence of 500 pM HIL-6 and in the absence of any ligand, suggesting that HIL-6 exhibited a relatively narrow range of concentration at which it was most efficient at maintaining ES cell pluripotentiality. There were, however, significant differences between the abilities of LIF and HIL-6, at higher ligand concentrations (50-500 pM), to maintain ES cell pluripotentiality (p < 0.10) with LIF being more potent than HIL-6 (Figs. 4A-4D
). Taken together, these results suggested that ES cell self-renewal varied in a dose-dependent manner for both LIF and HIL-6, a monophasic dose dependency for LIF and a bell-shaped, biphasic dose dependency for HIL-6 (though due to variability in the assays it was not possible to definitively conclude this) with LIF being more potent than HIL-6 at higher concentrations (
500 pM). It was noteworthy that there was a substantial background level of self-renewal obtained in the absence of exogenous cytokine supplementation indicative of either the selective survival of cells with a particular phenotype (i.e., SSEA-1 and E-cadherin positive) or perhaps, of endogenous production of cytokines capable of signaling through gp130, as we [11] and others [51] have previously shown. To minimize the effects of autocrine LIF production (which has been shown to be lower than the levels required for ES cell self-renewal [52]), cell density was carefully monitored and cultures were reinitiated on alternate days.
To confirm and extend the above analysis, ES cell differentiation was also measured through the analysis of Oct-4 protein expression at different doses of LIF and HIL-6 over a 7-day period. To verify that flow cytometric measurements of intranuclear Oct-4 protein expression correlated with other established ES cell markers (for example, SSEA-1 and Western analysis of Oct-4), we conducted an SSEA-1 selection experiment using petri dishes coated with an anti-SSEA-1 antibody. Fluorescent-activated sorting, based on differences in mean SSEA-1 fluorescence, was initially explored as the primary means of correlating Oct-4 intranuclear expression with SSEA-1 expression, however, poor ES cell viability after cell sorting prevented us from forming reliable conclusions. Using an alternative approach, ES cells or a control population of embryonic fibroblasts (STO cells, which do not express SSEA-1 or Oct-4) were incubated in the anti-SSEA-1 antibody-coated plates and the adherent (i.e., SSEA-1hi) and nonadherent (i.e., SSEA-1low) fractions were separated as described for further analysis. Both the adherent (data not shown) and nonadherent fractions (
75% of the population) of STO cells were negative for both SSEA-1 and Oct-4 (Figs. 5A, 5B
). As expected, the adherent ES cell fractions (
52% of the total population) had higher levels of SSEA-1 expression (65.9% of ES cells were SSEA-1hi) than the nonadherent fractions (
20% of the population), which had lower levels of SSEA-1 expression (23.5% were SSEA-1hi) (Fig. 5A
). In fact, after the selection, the SSEA-1hi population was almost completely removed from the nonadherent populations (Fig. 5A
). Importantly, the ES adherent fractions also had higher levels of Oct-4 expression (72.2% of ES cells were Oct-4hi relative to STO cells) than the nonadherent fractions (22.9% were Oct-4hi) (Fig. 5B
). A separate experiment confirmed that the adherent cells had higher Oct-4 protein expression by Western analysis than the nonadherent fraction (Fig. 5C
). This suggested that flow cytometric measurements of intracellular Oct-4 protein levels correlated with other established assays of ES cell differentiation status.
Using this novel, flow cytometry-based assay, our results showed that nuclear Oct-4 protein expression was downregulated after 7 days in culture in a dose-dependent manner for both LIF and HIL-6 (Figs. 6A-6C
). Cells cultured with LIF had increased Oct-4 protein expression (as analyzed both by Western blots and flow cytometry) at higher concentrations (50-500 pM) than at lower concentrations (0-5 pM) (p < 0.05). As before, at 500 pM, LIF was more potent than HIL-6 at sustaining Oct-4 protein levels after 7 days in culture (p < 0.05). Oct-4 expression exhibited a moderately biphasic dose dependency on HIL-6 concentrations, although, in this series of experiments, no significant differences between cells cultured at 50 pM and 500 pM of HIL-6 (p > 0.10) were observed.
Our experiments have repeatedly shown, albeit to different extents, that the maintenance of high levels of SSEA-1, E-cadherin, and Oct-4 expression, the capacity to form EBs, and high ALP activity (all associated with undifferentiated ES cells) were more efficient when cultures were supplemented with high doses of LIF than with equivalent doses of HIL-6 (Figs. 4 and 6
) (p < 0.10), and this was not due to differential depletion of ligands, as experiments were conducted with frequent medium exchanges. Additionally, while there were significant differences on day 7 between 50 and 500 pM HIL-6 in supporting EB formation efficiency and ALP activity (p < 0.10) suggestive of a biphasic dose dependency, smaller, nonsignificant differences were observed between 50 and 500 pM HIL-6-dependent SSEA-1/E-cadherin (p > 0.10) and Oct-4 expression (p > 0.10). These differences may have been masked by differential kinetics exhibited by ES cells in the loss of SSEA-1, E-cadherin (data not shown), and Oct-4 expression (Fig. 6D
) when exposed to a range of HIL-6 and LIF doses. Due to variability within and between assays, it was difficult to definitively conclude that ES cell self-renewal responses had a biphasic dose dependency in response to HIL-6 concentrations. This led us to examine whether further insight into the mechanisms by which HIL-6 and LIF were acting could not be better elucidated using a computational modeling approach.
Ligand/Receptor Signaling Threshold Model for ES Cell Fate
A mechanistic model, based on signaling complex formation (LIF-LIFR-gp130 for LIF and gp130-HIL-6-gp130 [trimer] or gp130-HIL-6-gp130-HIL-6 [tetramer] for HIL-6), parametric values of which may be critical in determining ES cell fate, was developed. Mutagenesis studies have revealed that IL-6 and LIF are trivalent ligands [5357]. LIF is thought to bind two LIFR molecules and one gp130 molecule, but a functional LIF complex is thought to be trimeric consisting of LIF, LIFR, and gp130 in a stoichiometric ratio of 1:1:1 [30] (Fig. 7A
). There is some debate as to the composition of a functional IL-6 signaling unit; it may be a tetramer composed of IL-6, IL-6R, and gp130 in a stoichiometric ratio of 1:1:2 [5860] or, alternatively, it may be a hexamer composed of two molecules of each component in a stoichiometric ratio of 2:2:2 [58, 61, 62]. Since the IL-6R binding site is unavailable in HIL-6 [48], the functional unit may be a trimer composed of HIL-6 and gp130 in a 1:2 stoichiometric ratio (Fig. 7B
) or a tetramer consisting of HIL-6 and gp130 in a stoichiometric ratio of 2:2 (Fig. 7C
). A combination model for the competitive formation of both the HIL-6 tetramers and HIL-6 trimers was also proposed (Fig. 7D
).

View larger version (20K):
[in this window]
[in a new window]
|
Figure 7. Proposed cell surface stoichiometry of LIF and HIL-6 (see text). (A) LIF-LIFR-gp130 signaling complex. (B) HIL-6-gp130-gp130 trimer complex. (C) HIL-6-gp130-HIL-6-gp130 tetramer complex. (D) A combination model involving competition between HIL-6-gp130-gp130 trimers and HIL-6-gp130-HIL-6-130 tetramers.
|
|
As an initial approach in the development of our LIST model, we calculated cell-surface signaling complex numbers at steady state without taking into account variations in signaling complex numbers due to ligand/receptor trafficking. This assumed that initial binding and trafficking events, which may occur on the order of seconds to hours, did not impact cell fate decisions (typically occurring on the order of hours to days). The numbers of LIF (CLIF) and HIL-6 signaling complexes formed in the trimer (CHIL-6trimer), tetramer (CHIL-6tetramer), or combination (Ccombination) modes were calculated on a per cell basis using mass action kinetics equations and were a function of gp130 and LIFR numbers (R1 and R2 respectively), ligand concentration (L), binding affinity constants (kf and kr), and cross-linking affinity constants (kc and ku) (Appendix 1 and Table 1
[50, 56, 6365]). Some values in our calculations could not be found for this system in the published literature and, thus, were estimated based on published ranges for similar processes (Table 1
). Signaling complex numbers for individual cells were summed over the entire population to determine bulk concentration of signaling complexes and plotted at different ligand concentrations (Fig. 8A
); the dose dependency of this simulation closely correlated with our experimental data of day 7 ES cell responses to log-fold changes in exogenous cytokine concentrations. Model-predicted bulk concentrations of LIF signaling complexes followed a sigmoidal dose-response curve; model-predicted bulk concentrations of HIL-6-trimer-type signaling complexes followed a bell-shaped, biphasic dose-response curve (Fig. 8A
). Maximal HIL-6-trimer-type complex formation (25% receptor occupancy) occurred around 1,000 pM of HIL-6 (equivalent to the low affinity binding constant estimated for HIL-6; similar to the published value of the low affinity equilibrium binding constant for LIF [64]) and the curve was symmetric about this ligand concentration as predicted previously for interactions of bivalent ligands with receptors [66]. Model-predicted bulk concentrations of HIL-6-tetramer-type signaling complexes exhibited a monophasic dependency on HIL-6 concentration, similar to LIF (Fig. 8A
), however, maximal HIL-6-tetramer-type complex formation (
25% receptor occupancy) occurred only at very high ligand concentrations (
1 µM). A combination model of HIL-6 signaling competitively between the trimer and tetramer modes followed the trimer curve until 1,000 pM (equivalent to the low affinity binding constant estimated for HIL-6), and then switched to the tetramer mode, which out-competed the trimer mode of signaling at higher HIL-6 concentrations (Fig. 8A
). Our model correlated with experimentally observed data that showed that 500 pM LIF was significantly more potent than 500 pM HIL-6 in forming signaling complex numbers, and hence in sustaining ES cell pluripotentiality. Further, as we have shown experimentally, the model predicted a monophasic dependency of signal complex number formation on LIF concentrations, and a bi/multiphasic dependency on HIL-6 depending on which mode of signaling was predominant (trimer, tetramer, or a combination of the two).
An advantage of our computational approach was that it allowed us (by lumping and varying parameters over a range of values in model simulations [66]) to understand how variations in these parameters could effect cell-surface signaling complex formation (and potentially ES cell responsessee below). For example, we showed that by varying cell-surface receptor numbers (such as may occur due to changes in receptor internalization or synthesis rates), we affected the maximal number of signaling complexes formed (at steady-state) at a given ligand concentration (Fig. 8A
) (consequently affecting the percentage of cells in a given population above or below a signaling threshold). Varying binding affinity (KD) constants caused horizontal shifts in the curve (Fig. 8A
) suggesting that greater ligand concentrations were needed to achieve maximal receptor occupancy. A normalized curve of HIL-6-trimer complexes/gp130 receptor numbers (reflecting receptor occupancy) attained a maximal value when the ligand concentration (L) was equivalent to the binding affinity constant (KD2), and was a function of the product of ligand cross-linking constant (KC2) and gp130 receptor numbers (RT2) [66] (Fig. 8B
). Similar normalization for LIF complexes, HIL-6-tetramer, and HIL-6 combination complexes revealed that maximum complex numbers formed, and the ligand concentration at which these occurred depended on the lumped parameters (Appendix 1).
Calculated signaling complex numbers (Table 1
) were also used to predict whether or not an individual cell had exceeded a specified threshold for self-renewal. Although our model thresholds were based upon experimental observations that ES cells self-renewed at LIF concentrations greater than 100 pM while cells differentiated at LIF concentrations less than 0.5 pM (4-10 pM LIF was previously shown to be the half-maximally effective concentration to prevent differentiation [26]), our computational approach allowed us to effectively investigate this choice of threshold and other ligand specific parameters. The mean LIF-equivalent complexes formed at the selected concentrations were designated as the "high signal threshold" and the "low signal threshold," respectively. If the cell-surface complex number exceeded the "high signal threshold," it was given an individual probability of self-renewal (Pi) value of 1. If the cell-surface complex number on an individual cell fell below the "low signal threshold," the cell was assigned a Pi value of 0, i.e., it had zero probability of self-renewing. Cells were assigned a Pi value of 0.5 when their cell-surface complex numbers were between the high and low signal thresholds, signifying that they had a 50% chance of self-renewing and a 50% chance of differentiating upon cell division (alternatively, they could divide asymmetrically). Model results were computed with and without the Pi = 0.5 option (Fig. 9A
). The individual cell probability values were summed over a theoretical population of 106 cells to yield a population probability of ES cell self-renewal and plotted as a function of bulk ligand-receptor signaling complexes for LIF and HIL-6 (Fig. 9B
). The most noteworthy feature of the computational plots was the illustration that thresholds in signaling complex numbers were the important determinants of ES cell self-renewal responses; the model-computed threshold for 50% or more of the population to self-renew (calculated using values from Table 1
) was between 0.5 and 26 pM of LIF-equivalent signaling complex numbers. On replotting our experimental data from multiple assays as a function of signaling complex numbers (only data for SSEA-1/E-cadherin-positive cells are shown), we observed a similar threshold-dependent quantal effect (Fig. 9C
): the threshold for LIF-mediated self-renewal of 50% of the population was between 0.4 and 4 pM of signaling complexes; the threshold for HIL-6-trimer mediated self-renewal of 50% of the population was between 3 and 22 pM signaling complexes. The threshold-dependence of ES cell markers on cell-surface signaling complex numbers calculated by our computational model, at steady-state, supported our LIST model premise that signaling complex numbers are important determinants of stem cell fate decisions.
 |
DISCUSSION
|
|---|
Parametric values of ligand-receptor interactions have previously been shown to be relevant in determining cell fate for hematopoietic stem cells [7, 8], neural progenitors [67], T lymphocytes [68, 69], and even in determining cell fate during development in nonvertebrates [70]. While different stem cell systems may utilize different combinations of signaling ligands and receptors, critical levels of their interactions may ultimately determine the nature of the response that is elicited. A generalized model that illustrates how quantitative variations in ligand-receptor interactions, arising from interactions of the cell with its microenvironment, can result in alteration in cell fate choices would thus be extremely useful (Fig. 10
). Importantly, the number of signaling ligand-receptor complexes is not a static property of the cell; it is dynamically modulated by extrinsic factors including soluble factors and cell-cell and cell-matrix interactions, which can influence receptor up- or downregulation, autocrine secretion, ligand depletion, negative feedback signaling, etc.

View larger version (28K):
[in this window]
[in a new window]
|
Figure 10. Ligand-receptor signaling threshold model for ES cell and tissue-specific stem cell fate determination. The parametric value of signaling complex numbers on individual cells determines which side of the threshold, i.e., self-renewal versus differentiation, they fall on, and thereby determines their cell fate. In this example, cells with lower-than-the-threshold number of signaling complexes differentiate, while cells with higher-than-the-threshold number of signaling complexes self-renew.
|
|
Herein, we have obtained results indicating that thresholds in cell-surface ligand-receptor complexes can affect ES cell fate by examining ES cell responses to different doses of IL-6-type cytokines. ES cells are an ideal model system to study stem cell fate as their self-renewal status in vitro correlates with numerous quantitative assays, and decisions made by ES cells to self-renew or differentiate may be similar to those made by other stem cell systems such as HSC [5, 7, 8]. We found correlation between levels of molecular markers (Oct-4 protein expression), cell-surface phenotype markers (SSEA-1 and E-cadherin expression), and cytokine stimulation, implying that thresholds in levels of ligand-receptor activation at the cell surface were directly translated to levels of downstream transcription factors that were activated and were responsible for triggering different gene programs. This correlation might prove useful in designing protocols aimed at enriching stem cell population(s) based on specific cell-surface antigen(s) expression.
We used two cytokines from the IL-6 family of cytokines, LIF, which signals through gp130-LIFR heterodimers, and HIL-6, which signals through gp130-gp130 homodimers, to demonstrate how variations in receptor-binding properties of the signaling cytokines can influence cell fate as predicted by our LIST model. Importantly, differences in LIF and HIL-6 signaling may be additionally influenced by the differential signaling potencies of heterodimers and homodimers of gp130, which, though not shown conclusively, may explain some of the discrepancies between model predictions (Fig. 9B
) and experimental data (Fig. 9C
). There is some evidence to suggest that heterodimer signaling may be more potent in activating and amplifying janus kinase/signal transducers and activators of transcription (Jak/STAT) and Ras/MAPK (mitogen-activated protein kinase) pathways than homodimer signaling [26], and this may be linked to their differential sorting on internalization, as has been shown for ErbB receptor-tyrosine kinases [71]. However, LIFR is considered to be functionally similar to the gp130 receptor (although not identical [7274]) suggesting that the signaling potencies of gp130-gp130 homodimers and gp130-LIFR heterodimers may be equivalent, and other explanations, such as we have provided herein, likely play a role in the control of stem cell responses.
We demonstrated that attenuation of cell-surface complex numbers can be achieved by signaling through HIL-6, which required a higher stoichiometric ratio of signaling receptors to ligand than the simple 1:1 ratio required for LIF, and, accordingly, resulted in the formation of lower signaling complex numbers at ligand concentrations that far exceeded available receptor numbers. That is, for a hypothetical cell with "X" number of gp130 receptors, stimulation by LIF resulted in the formation of "X" number of signaling complexes at saturation (providing that LIFR numbers were not limiting), while stimulation by equimolar concentrations of HIL-6 resulted in the formation of X/2 signaling complexes at saturation. If the postulated threshold for self-renewal lies between X/2 and X signaling complexes, LIF would have been more efficient than HIL-6 in maintaining ES cell pluripotentiality. In fact, equimolar concentrations of LIF and HIL-6, although similar, were not equivalent in maintaining marker expression associated with undifferentiated cells in the various assays used to gauge ES cell differentiation status. This correlated with our LIST model predictions, which were based on simple receptor-ligand binding and dimerization mass action equations. The correlation between the LIST model predictions and experimental data was all the more remarkable since the model simulations were confined to cell-surface ligand-receptor interactions and had to contend with estimates of receptor distributions and binding/dimerization constants, and with steady-state approximations of dynamically regulated processes. This suggested that the number of cellsurface signaling complexes at steady-state may be critical determinants of downstream processes that were responsible for turning on different genetic programs associated with ES cell self-renewal or differentiation. Discrepancies between the model and experimental data also arose because we had to contend with correlating discrete probability values of self-renewal (0, 0.5, 1) assigned to individual cells and iterated for a given population, to experimental data that represented a continuous range of possible values of cell-specific phenotypic and molecular marker expression. Inherent variability within the assays diminished the correlation; however, when data from multiple assays were taken together, the quantal relationship between cells positive for ES cell self-renewal markers and the number of cell-surface signaling complex numbers was apparent, thus validating our core hypothesis that this parameter had an important (governing) effect on the decision of cells to self-renew. Our computational model, which explained cell fate specifications at the population level by taking individual cell properties (i.e., signaling ligand-receptor complex numbers) into account, is a first step toward developing a more comprehensive deterministic model that includes multiple interactions among individual cells within a population, as well as the dynamics and interplay between cell-signaling networks.
Our approach was distinct from stochastic models of stem cell differentiation control that typically assume that cell fate processes are random and are best described by statistical probability distributions [7581]. These models argue that cytokines support an overall probability of progenitor cell survival and/or proliferation by selecting for cells that have already made the decision to self-renew or differentiate by some unknown random process. Our deterministic model, on the other hand, proposed a specific mechanistic basis (i.e., signaling through the ligand-receptor complexes) to account for the self-renewal of individual mouse ES cells in vitro based on cell-surface-receptor occupancy. Given our data that LIF and HIL-6 also mediate the survival of a subset of cells independent of any effects on ES cell proliferation [7, 8, 11] (possibly through the differential engagement of STAT3, known to be important in both cell survival [8284] and maintenance of ES cell pluripotentiality [85, 86]), we cannot eliminate the possibility that IL-6 type cytokines also influence ES cell fate by enhancing cell survival in addition to promoting their self-renewal. Ongoing studies should reveal the relative contribution of these two mechanisms on ES cell fate control.
Our model predictions correlated with ES cell marker's dependence on LIF concentration. For HIL-6, on the other hand, certain ES cell markers (EB formation efficiency and ALP activity) had a biphasic dependency on HIL-6 doses, while other markers (SSEA-1/E-cadherin expression and Oct-4 protein expression) did not exhibit a strong biphasic dose dependency. These results suggested that the kinetics of these endpoint assays used showed variability in response to LIF versus HIL-6 stimulation. The observed biphasic dependency of EB formation efficiency and ALP activity was consistent with model predictions for HIL-6 signaling in the trimer mode and with other experimental data [5]. The self-antagonism of HIL-6 at higher concentrations, in agreement with data obtained for other bivalent cross-linking ligands, such as growth hormone [87, 88] and IL-6 [89] and with model predictions for other bivalent ligands [66], resulted from competition for gp130 receptors. The dose dependence of SSEA-1/E-cadherin and Oct-4 protein expression on HIL-6 concentrations could not be used to distinguish among the proposed trimer, tetramer, or combination modes of HIL-6 signaling. Signaling in the trimer mode may not have been apparent at the HIL-6 concentrations tested, which were not sufficiently high to observe significant antagonistic behavior. Alternatively, HIL-6 may have been signaling in the tetramer mode but, according to our model, signaling would only have been effective at very high ligand concentrations (in the µM range). It may also have been possible that HIL-6 signaled competitively in both the trimer and tetramer mode, oscillating from one mode to another with time and ligand concentration. ES cell responses would then demonstrate a partially biphasic (following the trimer mode of signaling up to 1,000 pM and then following the tetramer mode) dependency on HIL-6 concentration according to the combination model we developed, with the axis of symmetry lying at KD2, the equilibrium binding constant associated with HIL-6 binding to gp130. A different kind of combination model was previously suggested, which required the preformation of gp130 homodimers [59], however, there is recent evidence for the kind of sequential binding of HIL-6 to two gp130 receptors [60, 90] that can be interpreted in favor of all three models suggested here.
Our LIST model can be adapted to describe critical ligand-receptor interactions in different cell systems; for example, to explain the observed increase in osteoclast formation upon tetracycline-regulated upregulation of gp130 levels in a stromal/osteoblastic cell line upon stimulation with IL-6/sIL-6R, but not with OSM [91].
 |
SUMMARY
|
|---|
We have proposed a mechanistically based model for stem cell fate control where the level of cytokine signaling results in the activation of gene programs responsible for determining cell fate. We have shown that gp130-mediated signaling influences self-renewal and possible survival outcomes in a threshold-dependent manner that was independent of any effects on proliferation. We demonstrated that there were differences at higher ligand concentrations in the abilities of LIF and HIL-6 to sustain ES cell pluripotentiality, likely due to different receptor-ligand binding stoichiometries. An important message from this study is that outcomes of stem cell fate choices are amenable to exogenous control by appropriate types and levels of cytokine stimulation. Insights from our ES cell model should provide deterministic clues to the numerous, seemingly stochastic parameters that affect stem cell fate and, thus, help in developing clinically relevant protocols to either expand stem cell populations or direct their differentiation into appropriate tissues in vitro.
 |
APPENDIX
|
|---|
 |
ACKNOWLEDGMENT
|
|---|
SSEA-1 was developed by D. Solter and obtained from the Developmental Studies Hybridoma Bank maintained by The University of Iowa, Department of Biological Sciences, Iowa City, IA. CCE cells were obtained from George Daley (Whitehead Institute, Cambridge, MA). Thanks are also due to Ting Yin for expert technical assistance.
This study was possible through research support from a National Science Foundation (NSF) Engineering Research Center (ERC) grant to the Biotechnology Process Engineering Center (BPEC) at the Massachusetts Institute of Technology, and Natural Sciences and Engineering Research Council (NSERC) of Canada post-graduate award to S. Viswanathan. P. Zandstra is a Canada Research Chair in Stem Cell Bioengineering.
 |
REFERENCES
|
|---|
-
Thomson JA, Itskovitz-Eldor J, Shapiro SS et al. Embryonic stem cell lines derived from human blastocysts. Science
1998;282:11451147.[Abstract/Free Full Text]
-
Shamblott MJ, Axelman J, Wang S et al. Derivation of pluripotent stem cells from cultured human primordial germ cells. Proc Natl Acad Sci USA
1998;95:1372613731.[Abstract/Free Full Text]
-
Reubinoff BE, Pera MF, Fong CY et al. Embryonic stem cell lines from human blastocysts: somatic differentiation in vitro. Nat Biotechnol
2000;18:399404.[CrossRef][Medline]
-
Schuldiner M, Yanuka O, Itskovitz-Eldor J et al. From the cover: effects of eight growth factors on the differentiation of cells derived from human embryonic stem cells. Proc Natl Acad Sci USA
2000;97:1130711312.[Abstract/Free Full Text]
-
Audet J, Miller CL, Rose-John S et al. Distinct role of gp130 activation in promoting self-renewal divisions by mitogenically stimulated murine hematopoietic stem cells. Proc Natl Acad Sci USA
2001;98:17571762.[Abstract/Free Full Text]
-
Pawliuk R, Eaves C, Humphries RK. Evidence of both ontogeny and transplant dose-regulated expansion of hematopoietic stem cells in vivo. Blood
1996;88:28522858.[Abstract/Free Full Text]
-
Ramsfjell V, Bryder D, Bjorgvinsdottir H et al. Distinct requirements for optimal growth and In vitro expansion of human CD34(+)CD38(-) bone marrow long-term culture-initiating cells (LTC-IC), extended LTC-IC, and murine in vivo long-term reconstituting stem cells. Blood
1999;94:40934102.[Abstract/Free Full Text]
-
Zandstra PW, Conneally E, Petzer AL et al. Cytokine manipulation of primitive human hematopoietic cell self-renewal. Proc Natl Acad Sci USA
1997;94:46984703.[Abstract/Free Full Text]
-
Williams RL, Hilton DJ, Pease S et al. Myeloid leukaemia inhibitory factor maintains the developmental potential of embryonic stem cells. Nature
1988;336:684687.[CrossRef][Medline]
-
Smith AG, Heath JK, Donaldson DD et al. Inhibition of pluripotential embryonic stem cell differentiation by purified polypeptides. Nature
1988;336:688690.[CrossRef][Medline]
-
Zandstra PW, Le HV, Daley GQ et al. Leukemia inhibitory factor (LIF) concentration modulates embryonic stem cell self-renewal and differentiation independently of proliferation. Biotechnol Bioeng
2000;69:607617.[CrossRef][Medline]
-
Zandstra PW, Lauffenburger DA, Eaves CJ. A ligand-receptor signaling threshold model of stem cell differentiation control: a biologically conserved mechanism applicable to hematopoiesis. Blood
2000;96:12151222.[Abstract/Free Full Text]
-
Evans MJ, Kaufman MH. Establishment in culture of pluripotential cells from mouse embryos. Nature
1981;292:154156.[CrossRef][Medline]
-
Martin GR. Isolation of a pluripotent cell line from early mouse embryos cultured in medium conditioned by teratocarcinoma stem cells. Proc Natl Acad Sci USA
1981;78:76347638.[Abstract/Free Full Text]
-
Peters M, Muller AM, Rose-John S. Interleukin-6 and soluble interleukin-6 receptor: direct stimulation of gp130 and hematopoiesis. Blood
1998;92:34953504.[Free Full Text]
-
Sui X, Tsuji K, Tanaka R et al. gp130 and c-Kit signalings synergize for ex vivo expansion of human primitive hemopoietic progenitor cells. Proc Natl Acad Sci USA
1995;92:28592863.[Abstract/Free Full Text]
-
Yoshida K, Chambers I, Nichols J et al. Maintenance of the pluripotential phenotype of embryonic stem cells through direct activation of gp130 signalling pathways. Mech Dev
1994;45:163171.[CrossRef][Medline]
-
De Felici M, Dolci S. Leukemia inhibitory factor sustains the survival of mouse primordial germ cells cultured on TM4 feeder layers. Dev Biol
1991;147:281284.[CrossRef][Medline]
-
Matsui Y, Toksoz D, Nishikawa S et al. Effect of Steel factor and leukaemia inhibitory factor on murine primordial germ cells in culture. Nature
1991;353:750752.[CrossRef][Medline]
-
Matsui Y, Zsebo K, Hogan BL. Derivation of pluripotential embryonic stem cells from murine primordial germ cells in culture. Cell
1992;70:841847.[CrossRef][Medline]
-
Koshimizu U, Taga T, Watanabe M et al. Functional requirement of gp130-mediated signaling for growth and survival of mouse primordial germ cells in vitro and derivation of embryonic germ (EG) cells. Development
1996;122:12351242.[Abstract]
-
Saito M, Yoshida K, Hibi M et al. Molecular cloning of a murine IL-6 receptor-associated signal transducer, gp130, and its regulated expression in vivo. J Immunol
1992;148:40664071.[Abstract]
-
Piquet-Pellorce C, Grey L, Mereau A et al. Are LIF and related cytokines functionally equivalent? Exp Cell Res
1994;213:340347.[CrossRef][Medline]
-
Conover JC, Ip NY, Poueymirou WT et al. Ciliary neurotrophic factor maintains the pluripotentiality of embryonic stem cells. Development
1993;119:559565.[Abstract]
-
Rose TM, Weiford DM, Gunderson NL et al. Oncostatin M (OSM) inhibits the differentiation of pluripotent embryonic stem cells in vitro. Cytokine
1994;6:4854.[CrossRef][Medline]
-
Ernst M, Oates A, Dunn AR. Gp130-mediated signal transduction in embryonic stem cells involves activation of Jak and Ras/mitogen-activated protein kinase pathways. J Biol Chem
1996;271:3013630143.[Abstract/Free Full Text]