Cindy
Chung
ab,
Beth L.
Pruitt
b and
Sarah C.
Heilshorn
*a
aMaterials Science and Engineering, McCullough Building, 476 Lomita Mall, Stanford, CA, USA. E-mail: heilshorn@stanford.edu; Fax: +650 498 5596; Tel: +650 723 3763
bMechanical Engineering, Durand Building, 496 Lomita Mall, Stanford, CA, USA
First published on 15th July 2013
Cellular therapies have great potential to provide alternative treatment options for those suffering from heart disease. In order to optimize cell delivery for therapeutic efficacy, a greater understanding of parameters that impact stem cell differentiation, survival, growth, and development are needed. In this study, we examine the role of hydrogel crosslink density on spontaneous cardiomyocyte (CM) differentiation of murine embryoid bodies (EBs). CM differentiation was accelerated in hydrogels of low crosslink density, where 100% of the hydrogels were positive for CM differentiation compared to only 53% in the high crosslink density group after 8 days of culture. DNA microarray data suggests that enhanced CM differentiation in the low crosslink density hydrogels was not tissue specific but rather a result of favoured EB development and cell proliferation. Additionally, enhanced EB growth and differentiation in low crosslink density hydrogels was independent of RGD ligand density and not a consequence of enhanced diffusion. We also demonstrate that matrix metalloproteinase activity is required for spontaneous CM differentiation in 3D hydrogels. Low hydrogel crosslink density regulates spontaneous EB differentiation by promoting EB growth and development. Elucidating the effects of microenvironmental cues on cell differentiation can aid in the optimization of stem cell-based therapies for tissue regeneration.
With the ability for indefinite self-renewal and differentiation into multiple cell lineages, embryonic stem cells (ESCs) have the potential to be used as a viable cell source for clinical applications. However, in order for ESCs to be useful for medical therapies, a greater understanding of factors that impact their differentiation is essential. Increasing evidence points to the importance of the extracellular matrix (ECM) in mediating stem cell behavior and fate decisions.7,8In vivo, the ECM provides a complex network of physical, chemical, and mechanical cues. However, with the development of tunable engineered matrices, we are now able to systematically investigate cellular response to specific matrix properties.
In particular, matrix crosslink density has been shown to play a major role in regulation of cell morphology, growth, and function.9–13 Highly crosslinked poly(ethylene glycol) hydrogels physically inhibit cell growth in the absence of proteolytic sites.14 In addition, increased crosslink density has also been shown to inhibit endogenous ECM production and distribution by mesenchymal stem cells undergoing chondrogenic differentiation.15
Changes in hydrogel crosslink density can often alter several other hydrogel properties, e.g., modulus, diffusivity, and mesh size. Assuming a homogenous and amorphous hydrogel structure, increases in hydrogel crosslink density would result in increased modulus, decreased diffusivity, and decreased mesh size. Increasing hydrogel crosslink density may also limit the accessibility of the incorporated cell-adhesive ligands, e.g., RGD peptide. RGD peptide is often incorporated into hydrogel networks to promote integrin-mediated adhesion of cells to their surrounding extracellular matrix and subsequent signal transduction from the matrix to cells.16–18 All of these factors must be considered when analyzing cell response to changes in hydrogel crosslink density.
Previously, we used recombinant elastin-like protein (ELP) hydrogels19,20 to demonstrate that matrix crosslink density can be used to transiently inhibit the spontaneous contraction of pre-differentiated embryoid bodies (EBs) encapsulated in 3D matrices. When contracting EBs were encapsulated in ELP hydrogels of high crosslink density, contractile behaviour was transiently suspended. However, embryonic stem cell-derived cardiomyocytes were capable of adapting to the higher crosslinked network with recovered contractility after a period of culture time, suggesting an adaptation of the cells to their surrounding ECM. Here, we examine the role of matrix crosslink density and RGD-mediated integrin adhesion on CM differentiation with interest in how crosslink density affects cell gene expression, growth, and matrix remodeling with a 3D hydrogel environment.
![[R with combining low line]](https://www.rsc.org/images/entities/char_0052_0332.gif)
![[G with combining low line]](https://www.rsc.org/images/entities/char_0047_0332.gif)
SPASSAA). An ELP variant with a scrambled RDG sequence (labeled as 0 mM RGD) was used for comparison in experiments isolating RGD cell-adhesion effects.
![]() | ||
Fig. 1 Schematic of ELP hydrogel formation. ELP consists of 4 repeats of a cassette containing 1 bioactive domain and 3 identical elastin-like domains, which include lysine residues within the elastin-like domain to provide amine reactive crosslinking sites that form covalent bonds with THPC (A). Elastic moduli of 5 wt% ELP hydrogels (n = 4–5) with varying THPC : ELP reactive group stoichiometry (B). Effective diffusion coefficient for 70 kDa FITC-dextran in 5 wt% ELP hydrogels of varying crosslink density (C). All values are reported as mean ± SEM, where * denotes significant difference and n.s. denotes no significant difference, p < 0.05. | ||
:
1, 0.75
:
1, and 1
:
1). Acellular 5 wt% ELP hydrogels were tested in unconfined compression on an ARG2 rheometer (TA Instruments). Briefly, ELP hydrogels (50 μl) were crosslinked in cylindrical molds (8 mm diameter) for 10 min at room temperature followed by 10 min at 37 °C. Molds were removed, and hydrogels were compressed at 2 μm s−1 at 37 °C in a custom-made phosphate buffered saline (PBS) bath. Elastic modulus in compression was determined by taking the slope of the stress-strain curve for small strains (<15%). Effective diffusion coefficients were determined from fluorescent recovery after photobleaching (FRAP) experiments. Fluorescein isothiocyanate (FITC)-labeled dextran (70 kDa, Sigma) was encapsulated at 4 mg ml−1 in 5 wt% ELP hydrogels of different crosslink densities in custom molds (3 mm in diameter and 0.5 mm in height). FRAP experiments were performed on a Leica SPE confocal microscope. A 100 micron by 100 micron region of the image was photobleached with a 488 nm laser at 100% intensity for 15 seconds. A time series was captured every 1.6 seconds for 2 minutes to capture fluorescent recovery. The time series was analyzed using open source matlab code “frap_analysis” to determine the effective diffusion coefficient within the hydrogels.22 Values are reported as mean ± SEM.
:
1 and 1
:
1 crosslink densities in 5 wt% ELP hydrogels with 5.3 mM RGD ligands was performed. Using a false discovery rate (FDR) <0.05, a gene list of differentially regulated genes was created and examined using Ingenuity Pathway Analysis to determine the top up-regulated and down-regulated biological functions.
:
1 and 1
:
1 crosslinking groups, respectively (Fig. 1B). As expected, elastic moduli increases as THPC to ELP reactive group ratio approaches 1
:
1. To address the potential effects of limitations in nutrient and waste diffusion within hydrogels of different crosslink densities, FRAP experiments were conducted. For all hydrogels, similar effective diffusion coefficients for 70 kDa FITC-labeled dextran was observed, demonstrating that observed cellular effects among hydrogels of different crosslink densities would not be a consequence of diffusional limitations (Fig. 1C).
:
1 THPC to ELP reactive group stoichiometry, exhibited accelerated EGFP expression when compared to higher crosslink density groups. Significant differences in EGFP expression were observed at days 4, 6, and 8 (Fig. 2B). By day 8, 100% of the 0.5
:
1 hydrogels were positive for EGFP expression compared to only 53% in the 1
:
1 group.
![]() | ||
Fig. 2 Representative fluorescent images of EGFP expression within hydrogels after 2, 4, and 8 days of culture for crosslink densities of 0.5 : 1, 0.75 : 1, and 1 : 1 THPC to ELP reactive groups (A). Yield of EGFP positive hydrogels (n = 15–18), where * denotes significant difference among groups (p < 0.05) by χ2 test for binomial distributions (B). | ||
Contractility rate was used as an assessment of CM function. Despite the delayed CM differentiation in the higher crosslinking groups, no significant differences were observed in beat rate among all groups when corrected for the first day of spontaneous contraction (Fig. 3B). Even though no differences were observed among the groups, a general increase in beat rate was observed with time for all conditions. To delve into the underlying causes for the differences in CM differentiation rate, we looked at murine EB gene expression through DNA microarrays, murine EB outgrowth, and the ability of the cells to remodel their surrounding matrix.
![]() | ||
| Fig. 3 Yield of contracting hydrogels (A) and contractility rate (B) plotted for each hydrogel-encapsulated embryoid body exhibiting contractile behavior for up to 8 days. Data are displayed as the mean and SEM for each group (n = 8–14). Significant differences among groups (p < 0.05) using χ2 test for binomial distributions is denoted by *. No significant differences in contractility rates were observed among all groups. | ||
:
1 and 1
:
1 crosslink densities. Genes exhibiting significantly different gene regulation (p < 0.05) were analyzed for functional pathway enrichment. In comparing the 0.5
:
1 to 1
:
1 crosslink densities, genes associated with organism, embryonic, and tissue development, along with cellular assembly and organization were significantly up-regulated, while genetic programs related to cell death and disease were significantly down-regulated (Fig. 4). No significant differences were observed between 0.5
:
1 and 1
:
1 hydrogels crosslink density groups when looking for specific cardiac markers: GATA-4 (p = 0.4), α-myosin heavy chain (p = 0.15), cardiac troponin I (p = 0.93), connexin 43 (p = 0.42). Though differential expression for transcription factor NKX2-5 (p = 0.04) was observed, relative gene expression was not substantially up-regulated (0.5
:
1/1
:
1 = 1.08), where a 0.5
:
1/1
:
1 = 1 would represent no fold change.
![]() | ||
Fig. 4 Functional pathway enrichment analysis of genes differentially regulated (p < 0.05) between 0.5 : 1 and 1 : 1 crosslink density groups. Top five up and down regulated categories are presented, where the magnitude of functional increase (green) or decrease (red) of 0.5 : 1 group over 1 : 1 group for respective processes are indicated by activation z-score. | ||
:
1 group compared to the 1
:
1 crosslink density hydrogels with significant differences noted at later time points (Fig. 5A and B). Further support for increased murine EB outgrowth in hydrogels of lower crosslink density was observed when cell nuclei were stained with DAPI to visualize cell coverage within the hydrogels, where greater cell coverage was found for the 0.5
:
1 group compared to the 1
:
1 group (Fig. 5C).
![]() | ||
Fig. 5 Metabolic activity (n = 5–6) (A) and DNA content (n = 11–12) (B) for murine EBs encapsulated in 5 wt% ELP hydrogels of 0.5 : 1 and 1 : 1 crosslink density, with 5.3 or 0 mM RGD ligand density, at days 1, 4, and 8. DAPI coverage (n = 7–11) of murine EBs encapsulated after 8 days of culture (C). Significant difference (p < 0.05) compared to 1 : 1 counterparts determined by Tukey's posthoc test is denoted by *. No significant differences were observed between 0 (using scrambled RDG ELP variant) and 5.3 mM RGD hydrogels at the same crosslink density. | ||
By using a negative control ELP containing a non-adhesive, scrambled RDG sequence for comparison, control experiments were performed to isolate the effects of the RGD ligand. Interestingly, these studies clearly demonstrate that murine EB survival, proliferation, and growth were independent of RGD ligand presence (Fig. 5). No significant differences were found among ELP hydrogels with the same crosslink density that contained the RGD adhesive sequence or the inactive RDG sequence.
:
1 crosslink density group also favors greater MMP activity per cell.
![]() | ||
Fig. 6 MMP activity within 5 wt% ELP hydrogels with 5.3 or 0 mM RGD ligand density at 0.5 : 1 and 1 : 1 crosslink density (n = 3). | ||
To investigate the role of MMPs in spontaneous CM differentiation, MMP inhibitors were added to the culture media of 0.5
:
1 hydrogels. Compared to the control group without the addition of MMP inhibitors, the addition of soluble MMP inhibitors (PD166793 and doxycycline) impeded cell outgrowth and resulted in delayed CM differentiation (Fig. 7). Interestingly, this delay was not as pronounced as that observed when crosslink density of the hydrogel is increased (Fig. 2B). MMP-inhibited 0.5
:
1 hydrogels continued to permit CM differentiation with greater yields than that of hydrogels with higher crosslink densities.
Engineered matrices have been developed as tools to investigate extracellular matrix (ECM) effects on cell behavior. Using modular designs, matrix parameters can be independently tuned to isolate specific matrix effects in hopes of deconstructing complex micro-environmental cues. In particular, recombinant ELPs, attractive for their ability to mimic the modulus of native cardiac tissue and their resilience to withstand repeated strains during cell contraction, can be crosslinked into hydrogels of varying crosslink density while maintaining constant RGD cell-adhesive ligand density and protein weight percent.20 Previously, we have shown that ELP hydrogel crosslink density could modulate transient contractile inhibition of human embryoid bodies that had been allowed to spontaneously differentiate into cardiomyocytes prior to encapsulation.19 In this study, we investigated the spontaneous CM differentiation of mouse EBs within hydrogels of varying crosslink density, focusing on how crosslink density affects gene expression, cell growth, and matrix remodeling.
CM differentiation was temporally accelerated in hydrogels with lower crosslink density (0.5
:
1) (Fig. 2). DNA microarray data indicated significant differences in gene expression levels between murine EBs encapsulated in ELP hydrogels of 0.5
:
1 and 1
:
1 crosslink densities. Surprisingly, differentially regulated genes were not cell-lineage specific, but rather more indiscriminate, demonstrating that hydrogels of 0.5
:
1 crosslink density provided a more supportive environment for overall EB development (Fig. 4). Genes associated with growth, assembly, and organization (from the cell to organism scale) were significantly up-regulated, while death and disease were significantly down-regulated. Thus, the observed increases in CM differentiation in hydrogels with lower crosslink densities (0.5
:
1), could be a reflection of general EB growth and development.
Despite good viability of mouse EBs in all groups as previously reported,21 large differences were observed in how the cells grew and spread (Fig. 5) within hydrogels of low and high crosslink density. While hydrogels with high crosslink densities of other materials have been reported to experience diffusional limitations, we found that the diffusion of 70 kDa dextran was not restricted within the 1
:
1 hydrogels (Fig. 1C). As most paracrine signals secreted by cells are smaller than 70 kDa, these signals are likely smaller than the hydrogel mesh size and are capable of diffusing within hydrogels of low and high crosslink densities at similar rates. Thus, the effects of diffusional limitations in cell signaling and nutrient and waste transport are expected to be negligible in dictating cellular response.
Increasing hydrogel crosslink density may also limit the accessibility of the incorporated RGD ligands. To address this concern, ELPs including a scrambled, non-cell adhesive RDG domain were used for comparison. In some poly (ethylene glycol) systems, the addition of RGD ligands reduced EB aggregation and promoted endothelial differentiation; however, crosslink density was not specifically controlled.34 Our results (Fig. 5) show that the presence of RGD ligand was not essential for mouse EB metabolic activity, growth, and proliferation in the ELP hydrogels, demonstrating that crosslink density was the key parameter in EB growth regulation. Increases in crosslink density, or decreases in mesh size, have been shown to inhibit cell growth in other cell-laden 3D hydrogel systems.14,35 Accordingly, hydrogels of lower crosslink density contain greater voids and spaces, permitting cells to grow, spread, and subsequently differentiate. Crosslink density can also affect the local deposition of EB secreted proteins, which may drive EB differentiation patterns.36 It is important to note the differences among hydrogel culture systems. In this study, we focused on the differentiation of individual EBs, and we were able to control crosslink density while varying RGD ligand density independently to determine cell adhesive ligand effects.
To determine if the delayed CM differentiation was the result of the restricted ability for cells to remodel the surrounding matrix, mouse EB-laden hydrogels were assayed for MMP activity. MMPs represent the main group of regulating proteases in the ECM. MMPs induce ECM remodeling, which is essential for cell migration, proliferation, and differentiation during normal development.37,38 Hydrogels of lower crosslink density exhibited greater MMP activity per DNA content (Fig. 6). Increased crosslink density can result in restricted matrix distribution within hydrogels.15 Chen et al. has shown that ESCs are dependent on endogenous ECM for survival and differentiation.39 Thus, higher crosslink density in the 1
:
1 hydrogels could limit endogenous matrix deposition, and subsequently lower MMP remodeling activity. Additionally, during CM differentiation, cardiomyocyte progenitor cells (CPCs) have shown increased expression of genes associated with extracellular matrix synthesis (e.g. collagen type I, III, and IV, fibronectin, and elastin) and matrix remodeling (e.g. MMP −2 and −9, and tissue inhibitor of metalloproteinase (TIMP) −1, −2, and −4) compared to non-differentiating CPCs.40 Thus, the observed increase in MMP activity within the 0.5
:
1 hydrogels may be correlated to CM differentiation.
To investigate the role of matrix remodeling, encapsulated mouse EBs in 0.5
:
1 hydrogels were treated with exogenous, broad-spectrum MMP inhibitors, PD166793 and doxycycline. The addition of soluble MMP inhibitors delayed and decreased the yield of CM differentiation (Fig. 7). This effect was most notable early in the culture. As the cells continued to grow and proliferate in the 0.5
:
1 hydrogels with culture time, EB development could trigger increased matrix remodeling activity that may overcome the effects of the MMP inhibitors. Overall, the delayed differentiation suggests that cell-dictated matrix remodeling does play a role in spontaneous cardiac differentiation. This matrix remodeling is likely local to the EBs, as no visible bulk gel degradation was observed after two weeks of culture in previous studies with ELP hydrogels.19
![]() | ||
Fig. 7 Cardiomyocyte differentiation yield within 5 wt% ELP hydrogels with 0.5 : 1 crosslink density and 5.3 mM RGD ligand density quantified by EGFP expression (A) and contractility (B) with or without the addition of MMP inhibitors, PD166793 and doxycycline (n = 11–14). Significant differences among groups (p < 0.05) using χ2 test for binomial distributions is denoted by *. | ||
With changes in crosslink density, ELP hydrogel modulus is also altered (Fig. 1B). Studies in 2D have shown that substrate modulus can influence pluripotency and differentiation of embryonic stem cells.41,42 Substrates with elasticity matching native tissue of the desired cell type are thought to favor cell fate decisions towards the desired cell-lineage and cell maturation.43–45
In a 3D matrix metalloproteinase-sensitive poly (ethylene glycol) hydrogel system, compliant hydrogels ∼300 Pa, mimicking the elasticity of embryonic cardiac tissue, demonstrated increased expression of cardiac transcription factor, NKX2-5, while stiffer hydrogels ∼4000 Pa showed decreased expression.10 Here we demonstrate similar results, where our compliant hydrogel promoted cardiomyocyte differentiation. Given the compliant nature of the ELP hydrogel system, we would expect greater differentiation towards neural, cardiac, and muscle lineages as opposed to cartilage or bone.
| This journal is © The Royal Society of Chemistry 2013 |