|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
THE STEM CELL NICHE |
Departments of aCell and Molecular Physiology and
bBiomedical Engineering,
cProgram in Molecular Biology and Biotechnology, and
dCancer Center and Center for Gastrointestinal and Biliary Disease Biology, University of North Carolina School of Medicine, Chapel Hill, North Carolina, USA;
eDepartment of Medicinal Chemistry, Center for Therapeutic Biomaterials, Salt Lake City, Utah, USA
Key Words. Nuclear magnetic resonance • Metabolomics • Footprinting • Human • Hepatic stem cells • Hepatoblasts • Extracellular matrix • Hyaluronans • Collagens
Correspondence: Correspondence: William S. Turner, Ph.D., University of North Carolina School of Medicine, Campus Box 7038, Glaxo Building, Rooms 32-35, Chapel Hill, North Carolina 27599, USA. Telephone: 919-966-0347; Fax: 919-6112; e-mail: wsturner{at}email.unc.edu; or Lola M. Reid, Ph.D., University of North Carolina School of Medicine, Campus Box 7038, Glaxo Building, Rooms 32–35, Chapel Hill, North Carolina 27599, USA. Telephone: 919-966-0347; Fax: 919-6112; e-mail: lmreid{at}email.unc.edu
Received on October 19, 2007;
accepted for publication on January 30, 2008.
First published online in STEM CELLS EXPRESS March 6, 2008.
| ABSTRACT |
|---|
|
|
|---|
-fetoprotein (AFP); those in hydrogels were primarily hHBs, expressing AFP, albumin, and urea. Both hHpSCs and hHBs used energy provided by anaerobic metabolism. Variations in hyaluronan-matrix chemistry resulted in distinct profiles correlating with growth or with differentiative responses. Metabolomic footprinting offers noninvasive and nondestructive assessment of physiological states of stem/progenitor cells ex vivo. Disclosure of potential conflicts of interest is found at the end of this article.
| INTRODUCTION |
|---|
|
|
|---|
Human hepatic stem cells (hHpSCs) and human hepatoblasts (hHBs) are small cells (both less than 12 µm) with overlapping but distinct antigenic profiles; they are located in vivo within or next to the ductal plates of fetal livers and the canals of Hering in adult livers [5, 6]. The methods for identification, isolation, and culture of cells in monolayer cultures or in purified hyaluronan hydrogels have been described in prior reports [5, 7–12] Cells were cultured in Kubota's medium (KM), a serum-free medium designed for hepatic progenitors, and either on culture plastic or in hyaluronan hydrogels mixed with specific extracellular matrix components. Conditioned medium from cultures was collected daily and subjected to metabolomic profiling. Cell biological characterizations of cells in these conditions are described elsewhere (Turner et al., manuscript in preparation). This study focuses upon use of NMR to evaluate the cellular metabolome of the cells under these conditions.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Hyaluronan-Matrix Complexes
The complexes comprised two forms of hyaluronan (HA) [13]. Hyaluronan modified with a 3,3'-dithiobis(propanoic dihydrazide) thiol group (HA-DTPH) is a chemically modified HA possessing multiple thiols for cross-linking with disulfide bridges [14] or with a bivalent electrophile [15]. Carbylan-S (CMHA-S) is a carboxymethylated HA derivative further modified with multiple thiols for cross-linking [16]. The cross-linking was achieved by treatment with polyethylene glycol diacrylate. We prepared complexes of HA-DTPH or CMHA-S with type I collagen (Sigma-Aldrich, St. Louis, http://www.sigmaaldrich.com) and/or laminin (Sigma-Aldrich). Seven variations of substratum were tested: (a) culture plastic (PL), (b) HA-DTPH + collagen in a gelatin form (H + G), (c) HA-DTPH + gelatin + laminin (H + G + L), (d) CMHA-S (C), (e) CMHA-S + laminin (C + L), (f) CMHA-S + gelatin (C + G), and (g) CMHA-S + gelatin + laminin (C + G + L). Additional details on the preparation of the hydrogels are given in the supplemental online data.
Hyaluronan Matrices as Substrata
The mixtures of extracellular matrix components (sECMs) were placed into six-well polystyrene culture dishes (Falcon; Becton, Dickinson and Company, Franklin Lakes, NJ, http://www.bd.com) just after addition of the cross-linker and prior to gelation and incubated for 16 hours incubation at 37°C in a CO2 incubator (Forma Scientific, Thermo Scientific, Waltham, MA, http://www.thermo.com) to allow gelation to occur. One milliliter of KM was added to each hydrogel followed by hHpSCs and hHBs suspended in KM at seeding densities of 1 million cells per hydrogel (volume, 0.5 ml). If embedded into the hydrogels, cells were seeded prior to cross-linking. Cells embedded in the hydrogels remained as aggregates of cells throughout the hydrogel. Cells were cultured for 4 weeks, with media changes every day.
NMR Spectroscopy
Sample Collection. Samples were taken in triplicate from each experimental condition and each control on days 2, 8, and 12. Samples were stored at –80°C immediately following collection.
Sample Preparation. After thawing at room temperature, an aliquot of 540 µl of each sample was added to 5-mm magnetic resonance tubes (Wilmad Lab-Glass Inc., Buena, NJ, http://www.wilmad-labglass.com) along with 60 µl of a D2O (Cambridge Isotope Laboratories, Andover, MA, http://www.isotope.com) solution containing 81.84 mM formate (Alfa Aesar, Ward Hill, MA, http://www.alfa.com) for ratiometric peak area quantification, 8.94 mM 3-(trimethylsily)propionic-2,2,3,3-day4 acid sodium salt (TSP) (Sigma-Aldrich) for a chemical shift reference, and 0.2% NaN3 (Sigma-Aldrich) to inhibit bacterial growth.
1H NMR Spectroscopy. All 1H NMR spectra were obtained using a Varian INOVA 400 MHz NMR spectrometer (Varian Inc., Palo Alto, CA, http://www.varianinc.com). A standard one-pulse sequence with 60° pulse angle and a 1.5-second H2O presaturation was used, with an interpulse delay of 6 seconds, an acquisition time of 2.558 seconds, a sweep width of 4,650.08 Hz, and the transmitter offset centered on the H2O signal. The spectra were digitized with 24,000 points.
Spectral Processing. Processing of all NMR spectra was done in ACD/1D NMR Manager, version 8.0 (Advanced Chemistry Development, Inc., Toronto, ON, Canada, http://www.acdlabs.com). Exponential line broadening of 0.1 Hz and zero filling to 32,000 points were applied to each spectrum. After Fourier transformation, spectra were autophased and baseline corrected using uniform settings applied to all spectra. All spectra were referenced to the TSP peak at 0.00 ppm. Regions of the spectra associated with TSP (0.25 ppm and upfield), formate (8.23–8.51 ppm), and water (4.55–5.03 ppm) were excluded from the spectra for the purposes of spectral binning. The spectra were integrated using the Intelligent Bucketing method in the ACD software, with bin sizes of 0.02–0.06 ppm. This resulted in 250 bins. Bin area values for each bin were then used as inputs into principal component analysis (PCA).
Principal Component Analysis. PCA was used to identify the effects of various substrates and times of incubation on the 1H NMR spectra of the media. PCA was performed using standard procedures. Each spectral bin was mean-centered and autoscaled by subtracting its mean and dividing by its SD. The covariance matrix was computed using all spectrum samples from all three livers, seven substrates, and three days. The 250 eigenvectors of this covariance matrix were sorted in descending order by their eigenvalues. These eigenvectors define the new dimensions, or principal components, of the 250-dimensional data space originally defined by the 250 spectral bins. Unlike the coordinate system defined by the values of individual bins, the new coordinate system—defined by the principal components—maximizes the data variance along the first principal components. For analyses, the data samples (represented in the data space by single points) are projected onto a two-dimensional plane defined by two among the first several principal components so as to reveal any groupings of samples with respect to substrate of incubation time.
Post-PCA Spectral Binning Modifications. Full batch manual integrations were applied to the PCA-identified regions of interest to allow for sufficient bin broadening to fully capture integral values that may have been diminished due to the upper bin limit of 0.06 ppm. Care was taken to ensure that resonances adjacent to bins of interest were not captured in this process.
Spectral Assignments. Spectral assignments were made by reference to literature values of chemical shifts in various media and biological fluids.
Statistical Analyses
For each metabolite, a mixed model was fit with metabolite as the response. Times of 2, 8, and 12 days were used, and the following scaffolds were used: PL, C, C + L, C + G, C + G + L, H, and H + L. Fixed time points were used for assessing substrata effects, so that data for different matrices were compared at the same time point. Major variables that were impossible to control in the samples were the distinctions in the human livers from one donor to another. Those variations were managed in the statistical analyses as a random effect. Five hundred sixty pairwise comparisons were made between scaffolds all falling in the classes of C versus CL, C versus GC, C versus plastic, CL versus GCL, GC versus H, GCL versus HL, GCL versus plastic, H versus HL, H versus plastic, and HL versus plastic for each of the samples for three times, as well as for overall comparisons for all times combined. Since 10 comparisons were being made for each time point, a Bonferroni correction was done to account for multiple comparisons, and consequently only those p values less than .05/10= .005 were considered statistically significant.
Model
The statistical model can be expressed as the equation Yijkl=
+ βi +
j +
k +
i·j +
ijkl, where
is the overall mean response of the metabolite; β is a fixed substrate effect; subscript i indicates the seven substrates, numbered 1–7;
is a fixed effect of time and is denoted by subscript j, which is numbered 1–3 and indicates the day, where j= 1 is day 2, j= 2 is day 8, and j= 3 is day 12;
is the random effect of liver, which is interpreted as the assumption that all livers in the study are considered to be random samples of all possible livers; subscript k, which is labeled 1–3, indicates which of the three livers in the study is being used;
is the interaction effect of substrate and day and consequently has the subscripts i and j;
is the random error associated with the lth repetition of any set of substrate, day, and liver; and subscript l also registers 1–3 to indicate the repetition.
| RESULTS |
|---|
|
|
|---|
|
|
|
135°] shift of points; Figs. 2, 3) are the bins that in the Figure 4 loading plots are shifted in the same (
135°) or opposite (
215°) direction from the plot origin. In other words, the spectral bins sensitive to incubation time are located in the top left and bottom right corners of the loading plot. Among these, bins were found containing the following metabolites: glucose, alanine, lactate, arginine, glutamine, and glutamate. These bins are labeled accordingly in the Figure 4 plot. Thus, PCA suggests that concentrations of the metabolites in the supernatant change systematically from day 2 to day 12.
|
|
Figure 5B is the histogram of alanine. Again the x-axis represents the substratum, which is the type of matrix-hyaluronan scaffold, and the y-axis is the percentage formate region under the curve of the spectra. Alanine is produced in high quantity at day 2 of the culturing process and subsequently declines in production by day 12. There is a shift between the first two graphed sECMs, CMHA-S and CMHA-S-laminin, in which production falls and consumption begins by day 12.
Figure 5C represents the histogram for lactate. The level of lactate within KM is at a low level (0.0157% formate peak). However, cells by day 2 show a minimum of an eightfold increase in lactate. Tissue culture plastic is the substratum with the highest induction of lactate metabolites in the cells. Over the course of the experiment, lactate production subsequently falls to close to the amount found in KM. This suggests that the cells are recovering from the processing procedures that have injured them, relying less on anaerobic glycolysis for energy, and have subsequently settled down to a more normal metabolism indicated by the restoration of lactate to low levels.
Figure 5D represents arginine. The amino acid metabolite is present in KM at low levels compared with that found by day 2. Cells cultured on all forms of scaffolds/substrata demonstrate a high level of arginine, which subsequently falls in concentration quickly. By day 12, arginine levels almost reached a concentration equivalent to the starting amount introduced in culture media.
Figure 5E represents glutamine, which has an overall consumption rate on all substrates, as recognized between the media and the cultures on days 2, 8, and 12. However, net consumption of the glutamine goes down over the course of time and appears to stabilize within a narrow range by day 8. A common set of metabolites for this experiment is listed in supplemental online Table 1, giving a basis of reference for phenotyping the footprint of the human hepatic progenitor metabolites.
The statistical model used (described in Materials and Methods) contains the predicted trends for each treatment over time with consideration of substratum. The model is constructed by using the actual values to predict what occurs after all factors are taken into account. As with any of the metabolites, the model is considered individualized. The overall p values for this model are significant for time but not for substratum or the day/substratum interactions, which consider all the factors at once. The statistical model shows that the overall mean response of any metabolite is the combination of a fixed substratum effect, a fixed time effect, a random liver effect, and the interaction of the substratum and time variables. These effects are documented in pairwise comparisons for significance, found in Table 1.
|
|
-Ketoglutarate in the Krebs cycle can be affected by transamination reactions with alanine and aspartate to yield glutamate, which in turn is an intermediate toward the production of glutamine and/or albumin. Oxaloacetate from the tricarboxylic acid cycle is the byproduct of the interaction of the
-ketoglutamate with the aspartate, as well as the conversion from malate, and in similar fashion, the combination of oxaloacetate and glutamate forms aspartate and
-ketoglutarate. Aspartate is funneled into the urea cycle, leading to the production of arginine and, subsequently, urea. Pyruvate can be metabolized into alanine or lactate, with reversal of the process being plausible. Measured end products, urea and albumin, demonstrate that the cells are parenchyma, as only these cells synthesize these compounds. The relationship between arginine and urea production is indicated by the data in Figure 6B. The x-axis represents the normalized arginine concentration compared with the net urea analysis, where the y-axis represents the day. Across all substrata tested, the trend from day 2 to day 12 is a decrease in production of both arginine and urea. The base amount of arginine in KM is 200 mg/l. From day 2 through day 12, the amount of arginine production fell from 178.7516 to 30.7449 mg/dl arginine. The paths of the arginine (Fig. 6B, dotted line) and the urea (Fig. 6B, unbroken line) are similar, because they are directly tied to a coupled reaction in which arginine is used to make urea. As the system slows and less arginine is produced, less is available for the production of urea.
Similarly, albumin production in Figure 6C can be tied to metabolites detected in the phenotype through glutamate. Glutamate production is highest early in the culturing period (Fig. 5F) on all substrata and decreases over the time course. However, net production is shown over all substrata when removing the amount introduced by KM. The net effect when paths are skewed according to culturing techniques shows that as time progresses throughout the experiment, the amount of the cells and the conversion of glutamate to albumin or glutamine are opposite. The albumin production of the system decreases as the glutamine production increases.
| DISCUSSION |
|---|
|
|
|---|
7–10 days [9]. Daily metabolomic phenotyping by NMR assays of the hHBs cultured in the hydrogels provided a global survey of metabolic pathways in the cells. This is the first report, to our knowledge, of metabolomic phenotyping of a stem/progenitor cell population via unbiased high-throughput analytical methods using computational mathematics to interpret the metabolic footprint. The novelty of the approach is furthered indicated by the fact that there have been relatively few studies on metabolomics, even of mature hepatocytes, but none at all of hepatic progenitors. Most of the published reports are with respect to drug metabolism [2, 18, 19] or on a specific metabolic pathway, such as carbohydrates [20] or lipids [21]. Those studies characterize transcription factors or signaling pathways but do not analyze the complement of small molecules produced by the cells. It is possible that metabolomics were ignored because of the need for improved methods for doing metabolomic analyses on cells in culture. In this study, we carefully modified the extracellular matrix composition of the hydrogels to identify their influence on metabolomic differences indicative of a stem cell versus differentiative phenotype.
The extracellular matrix in combination with soluble signals (autocrine, paracrine, and endocrine) synergistically regulates all aspects of cells [22, 23]. Matrix effects on mature hepatocytes have long been known to regulate specific gene expression and, by inference, metabolic profiles [18, 23–27]. Changes in the chemistry of the matrix and soluble signals parallel the maturational stages of the cells [28–31]. The matrix chemistry near or in the liver's stem cell niches is composed of hyaluronans and their receptors, CD44H, laminins and certain integrins (
4β6), type III and IV collagen, and forms of chondroitin sulfate proteoglycans (PGs) with few sulfate residues [5, 11, 12, 29, 30]. The matrix chemistry transitions to that associated with fully mature cells and composed of collagen type I, some collagen type III, fibronectin, and forms of highly sulfated proteoglycans (e.g. heparin-PGs) but not laminin, hyaluronans, or type IV collagen [3–8, 23, 29, 30, 32, 33]. Because of the microenvironment [11, 12], the hHpSCs and hHBs are less likely to come into contact with the shear stresses from blood flow or with the hormonal factors contained within the blood of the sinusoids than their progeny adjacent to fenestrated endothelial cells of the sinusoids [34]. The fenestrated endothelia offer cellular boundaries in some regions and only matrix chemistry at the sites of the fenestrae. The knowledge of the matrix chemistry was used to define the variations in hyaluronan hydrogel matrix tested.
Applying metabolomic analytical methods, extensive multivariate statistical analyses were performed on the NMR data; these analyses further illustrated that the changes in the matrix chemistry can give rise to significant changes in the production of certain metabolites, including consumption of specific compounds. Implications of the changes in the matrix chemistry and comparison with what occurs with the control cultures (plastic) are indicative of what is going on within the cells on matrix scaffolds, especially those in 3D cultures. Table 1 is indicative of the comparison of cells on plastic and matrix-hyaluronan mixes. Of the 560 paired comparisons, those between the CMHA-S and plastic show the most significance in metabolic divergence. This is not to be confused as an identification of the best culturing conditions. The choice of culturing conditions is dictated, in part, by the desired biological responses of the cells. However, using this method of analysis, CMHA-S does exhibit properties that make it the preferred system, of those assayed, for the maintenance of hepatoblasts with minimal differentiation.
The relevance of being three-dimensional is also strongly implicated by our findings. The hHpSCs and hHBs in serum-free KM in monolayer cultures showed a strong selection for hHpSCs within 1 or 2 weeks [7, 9, 11, 12]. In the hyaluronan hydrogels, the selection proved true for hepatoblasts under all variations of hyaluronan and matrix components (e.g., HA-DTPH alone, HA-DTPH with laminin, and HA-DTPH with collagen). The variations in matrix components added to the hyaluronans affected expansion rates or differentiated functions in subtle ways, all the conditions assayed were equivalent in the ability to sustain human hepatic progenitors at the hepatoblast lineage stage. Three-dimensionality of cells when embedded in type I collagen resulted in more rapid differentiation, whereas in hyaluronans, it resulted in more rapid expansion. Therefore, the response of the cells is defined by the chemistry of the matrix components and by the medium conditions, as well as by being in a three-dimensional state. Alternatively, it may be a factor of whether a matrix component(s) is presented on one or multiple cell surfaces.
Serum was avoided, given that it is not defined, results in differentiation of mesenchymal and parenchymal progenitors, and results in selective overgrowth of mesenchymal cell populations [35]. The serum-free conditions contributed to greater stability and uniformity of gene expression, corroborating similar findings in prior studies [18].
The form of cross-linking of the hydrogels proved important in their influence on the cells. Aldehyde-cross-linked hydrogels were found to be supportive, albeit resulting in slow growth and with technical difficulties in retrieval of the cells [36]. The hydrogels used here are cross-linked by disulfide bridges and were found to be permissive for rapid growth and for ease of cell recovery. This permissiveness is assumed to be due to the cells' ability to enzymatically modify the scaffolding, a requirement for cells to go through cell division.
Metabolomic profiles (Fig. 5; Table 1; also described in supplemental online data) readily indicate cells in aerobic versus anaerobic metabolism. Our data indicate that the hHpSCs on culture plastic and the hHBs in the hydrogels rely primarily on anaerobic glycolysis, especially initially in the incubation. Alanine and lactate are indicative of anaerobic glycolysis [37]. Others have reported that low oxygen concentrations in mammalian embryonic tissues are normal. Thus, maintenance of cultures at low oxygen tensions should be selective conditions for stem cells and early progenitors [38–40]. Cells engaged primarily in anaerobic, not aerobic, metabolism are highly indicative of stem cells [41]. This correlates with the known oxygen conditions in vivo, in which the site of the liver's stem cell niche has significantly lower levels of oxygen than that next to the more mature parenchymal cells [40].
Also, hHpSCs and hHBs are remarkably tolerant of ischemia, surviving for hours in the livers of asystolic donors, whereas mature parenchymal cells typically die within an hour following cardiac arrest [9, 42]. Mammalian cells act in response to a hypoxic environment by increasing their carbohydrate consumption and switching to anaerobic respiration, a pattern also found in stem cells [40, 43]. Glycolytic enzymes are present in the cytosol and do not require mitochondria. Glycolytic metabolic pathways, the "Warburg effect," are known to be associated with cancer cells, now recognized to be transformed stem cells. Lactate is an end product of anaerobic glycolysis and is indicative of low oxygen availability, a relatively undifferentiated state [37], or consumption of it by specific cell types. Anaerobic use of glucose with low production of lactate has been reported within the hepatic culturing system. Specific stem cell populations, such as hemopoietic progenitors, are able to survive in this state by using glycolytic pathways [40, 44]. Grayson et al. [40] and Cipolleschi et al. [44] report that human mesenchymal stem cells have increased glucose consumption and lactate production under such conditions. Certain lineage stages of cells, especially stem cells, are able to survive under hypoxic conditions over extended periods of time [40], tolerance for ischemia [42], and relatively undifferentiated state.
The data here show that glucose consumption is high and that lactate production is high (Fig. 4). Each of the metabolites is reduced over time, which can be explained as an initial inflammatory process and acute injury responses when the cells are first plated and then adaptation to culture with a gradual shift to lineage stages optimally supported by the culture conditions used. Transamination of pyruvate, which leads to increased amounts of alanine, occurs to help decrease systems of lactate.
Conversion from anaerobic to aerobic metabolic pathways is indicative of differentiation [41] Differentiation is also indicated by metabolomic profiles implicating gluconeogenesis and/or the urea cycle pathway. Gluconeogenic amino acids are alanine, threonine, tyrosine, and phenylalanine [2]. Profiles for these can be constructed using the NMR phenotyping, where reduction could be indicative of that metabolic pathway. The urea cycle pathway includes biomarkers of arginine, glutamine, and urea [45]. Arginine is taken up by the cells and is used to produce urea [2]. As a direct player in production (Fig. 6), the log-log graph of arginine compared with production of urea follows the same path. Furthermore, the impact of this pathway and the albumin production pathway is complicated where the production of glutamine from glutamate is concerned. Glutamate can be shuttled toward one pathway or another, and in the production of glutamine, it is used as a cationic scavenger [2]. This is plausible for the graphic illustrated in Figure 6, in which glutamine and albumin follow opposing trends in production and consumption.
| CONCLUSION |
|---|
|
|
|---|
This application is directly applicable to cells in 3-D cultures or in bioreactors in which one can better mimic the multifaceted complexity of in vivo situations (2, 11, 41, 42). The database being established should prove ideal for rapid daily assessment of cell populations to see whether they are surviving, expanding, undergoing lineage restriction and/or differentiating to specific adult fates.
Bioreactors are being created with useful features (e.g. dissolvable boundaries, see through portals) and are complemented by platforms of analyses to overcome the limitations of access to cells within the devices. Metabolomic profiles offer non-invasive and non-destructive assays to assess the cells in bioreactors or other forms of 3-D culture conditions.
| DISCLOSURE OF POTENTIAL CONFLICTS OF INTEREST |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| STEM CELLS | THE ONCOLOGIST | CME | ALPHAMED PRESS JOURNALS |