Abstract
Organoids, the condensed three-dimensional (3D) tissues emerged at the early stage of organogenesis, are a promising approach to regenerate functional and vascularized organ mimics. While incorporation of heterotypic cell types, such as human mesenchymal stem cells (hMSCs) and human induced pluripotent stem cells (hiPSCs)-derived neural progenitors aid neural organ development, the interactions of secreted factors during neurogenesis have not been well understood. The objective of this study is to investigate the impact of the composition and structure of 3D hybrid spheroids of hiPSCs and hMSCs on dorsal cortical differentiation and the secretion of extracellular matrices and trophic factors in vitro. The hybrid spheroids were formed at different hiPSC:hMSC ratios (100:0, 75:25, 50:50, 25:75, 0:100) using direct mixing or pre-hiPSC aggregation method, which generated dynamic spheroid structure. The cellular organization, proliferation, neural marker expression, and the secretion of extracellular matrix proteins and the cytokines were characterized. The incorporation of MSCs upregulated Nestin and β-tubulin III expression (the dorsal cortical identity was shown by Pax6 and TBR1 expression), matrix remodeling proteins, and the secretion of transforming growth factor-β1 and prostaglandin E2. This study indicates that the appropriate composition and structure of hiPSC-MSC spheroids promote neural differentiation and trophic factor and matrix secretion due to the heterotypic cell–cell interactions.
Introduction
C
Introducing heterotypic cell–cell interactions in the transplanted cells has shown promising outcomes. For example, transplantation of hiPSC-derived cardiomyocytes, endothelial cells, and smooth muscle cells showed much better cell engraftment than cardiomyocytes alone in a large animal model.10,11 Similarly, cotransplantation of activated Schwann cells and MSCs improved functional recovery of hindlimbs and remyelination of the injured axons with reduced formation of glial scar. 12 However, cotransplantations of a cell mixture in suspension lack the tissue architecture that has been shown to significantly influence secretory properties and cell survival. 13 Presenting the cells in three-dimensional (3D) tissue format may overcome the drawback of single-cell suspension.
Attempts have been made to transplant microtissues or spheroids to enhance cellular functional properties compared with single-cell suspension. For example, the formation of spheroids upregulates the secretion of trophic factors and cytokines for NPCs. 14 For MSCs, the formation of spheroids has been found to enhance the expression of antiapoptotic and anti-inflammatory proteins such as tumor necrosis factor-inducible gene 6 protein (TSG-6), STC-1, IL-6, and prostaglandin E2 (PGE-2).15–19 MSC spheroids also enhance extracellular matrix (ECM) secretion, such as fibronectin and laminin,20–22 which promotes NPC differentiation, prolongs cell survival, and provides better protection from nutrient/growth factor deprivation. 23 However, modulation of cellular microenvironment indicated by the secretion of trophic factors and ECM proteins has not been well understood during spheroid formation of multiple cell types.
Recently, Takebe et al. reported a generalized method to form organ buds of different types of tissues (e.g., liver, kidney, brain, etc.) using MSCs, endothelial cells, and hiPSC-derived progenitor cells. 24 The transplanted organ buds were shown to efficiently self-organize into multiple functional and vascularized organs, such as liver or kidney.24,25 Moreover, complex 3D organoids have been created recently from hiPSCs, such as cortical or cerebral organoids,26–28 to predict neurotoxicity of different drugs, 29 and to improve transplantation outcome.25,30 By promoting the secretion of endogenous growth factors and ECM proteins, the formation of heterogeneous 3D tissue architecture recapitulates essential cellular microenvironment and cell–matrix interactions of functional tissues.14,31 However, the precise role of 3D cellular organization and structural environment on organoids' biological function have been underexplored.
Given these recent advancements, this study constructs spheroids of human MSCs (hMSCs) and hiPSC-derived NPCs and investigates the impact of cell organization on the proliferation, secretion of trophic factors and ECMs, and neural differentiation in vitro. It is hypothesized that the heterotypic cell–cell interactions and the secretion of endogenous ECM/trophic factors in the 3D hybrid spheroids of MSCs and hiPSC-derived NPCs promote neural differentiation. Two hybrid spheroid structures fabricated by direct mixing or pre-hiPSC aggregation were assessed to understand the influence of tissue structure and the secreted trophic factors and ECMs on neural differentiation. Direct mixing enables maximum heterotypic cell–cell interactions, whereas pre-hiPSC aggregation promotes the formation of different spatial cellular organization, such as the core–shell structure. The outcome of this study advances our understanding of the interactions of multiple cell types and the modulation of extracellular microenvironment during neurogenesis toward the application of tissue organoids for in vitro disease modeling and cell therapy.
Materials and Methods
Undifferentiated human iPSC culture
Human iPSK3 cells were derived from human foreskin fibroblasts transfected with plasmid DNA encoding reprogramming factors OCT4, NANOG, SOX2, and LIN28 (kindly provided by Dr. Stephen Duncan, Medical College of Wisconsin, and Dr. David Gilbert, Department of Biological Sciences of Florida State University).32,33 Human iPSK3 cells were maintained in mTeSR serum-free medium (StemCell Technologies, Inc., Vancouver, Canada) on six-well plates coated with growth factor-reduced Geltrex (Life Technologies). The cells were passaged by Accutase dissociation every 5–6 days and seeded at 1 × 106 cells per well of a six-well plate in the presence of 10 μM Y27632 (Sigma) for the first 24 h.34–36
hMSC culture
Standardized frozen hMSCs from multiple donors were obtained from the Tulane Center for Gene Therapy and cultured as previously described.37,38 The hMSCs were isolated from the bone marrow of healthy donors ranging in age from 19 to 49 years based on plastic adherence, negative for CD34, CD45, CD117 (all <2%) and positive for CD29, CD44, CD49c, CD90, CD105, and CD147 markers (all >95%), and possess trilineage differentiation potential upon induction in vitro.39,40 Briefly, hMSCs were expanded at a density of 1.7 × 103 cells/cm2 using αMEM (Invitrogen) medium supplemented with 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin in a standard 5% CO2 incubator. At ∼80% confluence, adherent cells were harvested with 0.25% trypsin-EDTA (Sigma-Aldrich) and further propagated.
Hybrid hiPSC-hMSC spheroid formation and neural differentiation
Before forming hybrid hiPSC-hMSC spheroids, undifferentiated iPSK3 cells (or preformed hiPSC aggregates) were labeled with CellTracker Green (2.5 μM, Life Technologies) and hMSCs were labeled with CellTracker Red (2.5 μM) for 30 min. The labeling allows for imaging and tracking of the self-organization and location of two cell types after the mixing.
Method A—direct mixing
A total number of 2 × 104 cells were seeded into U-bottom Ultra-low attachment (ULA) 96-well plates (Corning Incorporated, Corning, NY) at different iPSC to MSC ratios (100:0, 75:25, 50:50, 25:75, and 0:100, n = 3–12) in neural differentiation medium composed of Dulbecco's modified Eagle's medium/Nutrient Mixture F-12 (DMEM/F-12) plus 2% B27 serum-free supplement (Life Technologies). Y27632 (10 μM) was added during the seeding and removed after 24 h. At day 1, the cell aggregates were treated with dual SMAD signaling inhibitors 10 μM SB431542 (Sigma) and 100 nM LDN193189 (Sigma) (Supplementary Fig. S1; Supplementary Data are available online at www.liebertpub.com/tea).35,41,42 After 7 days, cells were grown in neural medium without growth factors. After 14–20 days in suspension, the 3D NPC spheroids were replated onto Geltrex-coated surface and characterized for various neural markers (β-tubulin III, Nestin, and TBR1, etc.) and ECM remodeling proteins.
Method B—pre-hiPSC aggregation
Human iPSCs (0.5–2 × 104 cells) were seeded into U-bottom ULA 96-well plates in neural differentiation medium composed of DMEM/F-12 plus 2% B27 to construct hiPSC aggregates first (Supplementary Fig. S1). Y27632 (10 μM) was added during the seeding and removed after 24 h. The aggregates were treated with 10 μM SB431542 and 100 nM LDN193189. At day 5, hMSCs were added in the wells and cultured to form core–shell aggregates. For each group, a total number of 2 × 104 cells were seeded at five different iPSC to MSC ratios (100:0, 75:25, 50:50, 25:75, and 0:100, n = 3–12). After another 7 days, cells were grown in neural medium without growth factors. At 14–20 days after the addition of MSCs, similar characterizations were performed as those for the direct mixing method.
Aggregate size distribution
The images of various spheroids were captured over the culture time (7 days) by a phase-contrast microscopy. The captured images were converted to binary images using ImageJ software (http://rsb.info.nih.gov/ij) and analyzed with the “particle analysis tool.” Through particle analysis in ImageJ software, the Feret's diameter of each aggregate in the images can be calculated, which provides the size distribution of the aggregates. 43
Biochemical assays
MTT assay
The spheroids under different conditions were incubated with 5 mg/mL 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT; Sigma) solution at day 7 after being incorporated with MSCs. The absorbance of the samples was measured at 490 nm using a microplate reader (Bio-Rad, Richmond, CA).
DNA assay
The DNA content of the individual spheroids at day 7 after the addition of MSCs was determined. DNA standard was prepared by dissolving salmon testes DNA in TEX (10 mM Tris, 1 mM EDTA, 0.1% Triton X-100 at pH 8) and a standard curve was constructed for each assay. The aggregates were lysed with 0.1 mg/mL proteinase K (Fisher Scientific, Pittsburgh, PA) at 50°C overnight. The lysates (100 μL) were mixed with 100 μL of PicoGreen (Molecular Probes, Eugene, OR) in a 96-well plate. The plate was incubated for 5 min in the dark and then read on a fluorescent plate reader (FLX800; Bioinstrument, Inc., Winooski, VT).
LIVE/DEAD assay
The spheroids were evaluated for viability using the Live/Dead® Staining Kit (Molecular Probes). After 7 days, the cells were incubated in DMEM-F12 containing 1 μM calcein AM (green) and 2 μM ethidium homodimer I (red) for 30 min. The samples were imaged under a fluorescent microscope.
Enzyme-linked immunosorbent assay
To quantify the growth factors secreted by different spheroids, culture supernatants were collected at day 4 after adding MSCs. Concentrations of prostaglandin E2 (PGE2) and transforming growth factor (TGF)-β1 were measured by enzyme-linked immunosorbent assay (ELISA) according to the manufacturer's instructions (R&D Systems, Minneapolis, MN for PGE2; Life Technologies for TGF-β1). Briefly, the samples were added into 96-well microplates and incubated with primary/secondary antibody solution conjugated with horseradish peroxidase (HRP) for 2–3 h. After washing, 3,3′,5,5′-tetramethylbenzidine (TMB) substrate solution was added for 30 min and stopped by the stop solution. The absorbance units (AU) were measured using a microplate reader (Bio-Rad) at a wavelength of 450 nm with background subtraction at 540–570 nm. The AU values were corrected by subtracting the values of negative control stained with HRP-IgG only. The concentrations of PGE2 or TGF-β1 were normalized by the ratios of MSC numbers to total cell numbers.
Histology
For histology, the hybrid spheroids of hiPSCs (prelabeled with CellTracker Green) and hMSCs (prelabeled with CellTracker Red) were fixed in 10% formalin, dehydrated, and embedded in paraffin wax. The sections of 10 μm were cut and fluorescent images were captured with an Olympus IX70 microscope with MagnaFire SP 2.1B software to show the localization of hiPSCs and hMSCs. 19
Immunocytochemistry
Briefly, the samples were fixed with 4% paraformaldehyde (PFA) and permeabilized with 0.2–0.5% Triton X-100. The samples were then blocked for 30 min and incubated with various mouse or rabbit primary antibodies (Supplementary Table S1) for 4 h. After washing, the cells were incubated with the corresponding secondary antibody: Alexa Fluor® 488 goat anti-Mouse IgG1 or Alexa Fluor 488 or 594 goat anti-Rabbit IgG (Life Technologies) for 1 h. The samples were counterstained with Hoechst 33342 and visualized using a fluorescence microscope (Olympus IX70, Melville, NY) or a confocal microscope (Zeiss LSM 880). For ECM expression of laminin and Collagen IV, replated aggregates onto 0.1% gelatin-coated surface (to reduce the interference from Geltrex) were fixed and stained. MSCs were labeled with CellTracker Red for visualizing the colocalization of MSCs and ECM proteins. The images from five independent fields (800–1000 cells) were analyzed using ImageJ software. The intensity was calculated based on the area of marker of interest normalized to the nuclei, indicating the relative expression among different conditions.
Flow cytometry
To quantify the levels of various neural marker expression, the cells were harvested by trypsinization and analyzed by flow cytometry. 44 Briefly, 1 × 106 cells per sample were fixed with 4% PFA and washed with staining buffer (2% FBS in phosphate buffered saline). The cells were permeabilized with 100% cold methanol, blocked, and then incubated with primary antibodies against Nestin and β-tubulin III followed by the corresponding secondary antibody Alexa Fluor 488 goat anti-Mouse IgG1 (for β-tubulin III) and Alexa Fluor 488 goat anti-Rabbit IgG (for Nestin). The cells were acquired with BD FACSCanto™ II flow cytometer (Becton Dickinson) and analyzed against isotype controls using FlowJo software.
Reverse transcription–polymerase chain reaction analysis
Total RNA was isolated from different cell samples using the RNeasy Mini Kit (Qiagen, Valencia, CA) according to the manufacturer's protocol followed by the treatment of the DNA-Free RNA Kit (Zymo, Irvine, CA).
45
Reverse transcription was carried out using 2 μg of total RNA, anchored oligo-dT primers (Operon, Huntsville, AL), and Superscript III (Invitrogen, Carlsbad, CA) (according to the protocol of the manufacturer). Primers specific for target genes (Table 1) were designed using the software Oligo Explorer 1.2 (GeneLink, Hawthorne, NY; Table 1). The gene β-actin was used as an endogenous control for normalization of expression levels. Real-time reverse transcription–polymerase chain reaction (RT-PCR) reactions were performed on an ABI7500 instrument (Applied Biosystems, Foster City, CA), using SYBR1 Green PCR Master Mix (Applied Biosystems). The amplification reactions were performed as follows: 2 min at 50°C, 10 min at 95°C, and 40 cycles of 95°C for 15 s and 55°C for 30 s, and 68°C for 30 s. Fold variation in gene expression was quantified by means of the comparative Ct method:
Statistical analysis
Each experiment was carried out at least three times. The representative experiments were presented and the results were expressed as mean ± standard deviation. To assess the statistical significance, one-way analysis of variance followed by Fisher's least significant difference post hoc tests were performed. A p-value <0.05 was considered statistically significant.
Results
Aggregation and localization of hiPSC-hMSC hybrid spheroids
Method A—direct mixing
For the four hiPSC to hMSC ratios (100:0, 75:25, 50:50, and 25:75) formed by direct mixing, the cells formed structured aggregates and showed the proliferation by aggregate diameter increase (e.g., for 50:50, from 687 ± 25 to 903 ± 19 μm) during day 1–7 (Fig. 1Ai, ii). Most cells were viable based on LIVE/DEAD viability assay (Supplementary Fig. S2). By contrast, the size of hMSC (0:100) spheroids decreased from 753 ± 22 to 446 ± 22 μm during day 1–7 (Fig. 1Aii). Notably, hMSC condensation in the center of spheroids was observed after day 2 (Fig. 1Ai). To confirm the observation by phase-contrast images, hMSCs were labeled with CellTracker Red and hiPSCs were labeled with CellTracker Green to visualize spatial cell distribution of two different cell types (Fig. 2A, Supplementary Fig. S3 and Supplementary Fig. S4). At day 1, hMSCs and hiPSCs were randomly distributed in the spheroids. After day 2, red hMSCs were gradually condensed in the core, whereas hiPSCs were mainly found around hMSC core. In addition, the relative area occupied by the condensed hMSCs decreased gradually from day 3 to 7.

Kinetics of hybrid hiPSC-hMSC aggregate size for direct mixing

Spatial cell distribution of hiPSC-hMSC spheroids formed by direct mixing (method A) and pre-hiPSC aggregation (method B) at different culture times.
Method B—pre-hiPSC aggregation
Undifferentiated hiPSCs (0.5–2 × 104 cells) were seeded for 5 days to form hiPSC aggregates. At day 5, hMSCs were added onto preformed hiPSC aggregates and the core–shell (hiPSC core–hMSC shell) structures were observed under different iPSC:MSC ratios (75:25, 50:50, and 25:75) (Fig. 1Bi, ii). Most cells were viable based on LIVE/DEAD viability assay (Supplementary Fig. S2). Slight increase in aggregate size was observed during day 2–7 (e.g., for 75:25 ratio, from 843 ± 47 to 914 ± 11 μm). Adversely, the size of hMSC (0:100) spheroids decreased in diameter from 894 ± 22 to 558 ± 3 μm during day 1–7. In addition, iPSC-MSC hybrid aggregates (e.g., 75:25) displayed core–shell spatial distribution as demonstrated by CellTracker labeling (Fig. 2A, Supplementary Fig. S3 and Supplementary Fig. S4), with hMSCs as shell and preformed hiPSC aggregates as core. Interestingly, hMSCs showed the tendency to diffuse into the preformed hiPSC aggregates and move toward the spheroid center during the culture period (e.g., day 2 after merging, 50:50 ratio).
Cell proliferation and metabolic activity of hiPSC-hMSC hybrid spheroids
MTT activity showed cell viability and metabolic activity of the cells in the spheroids (Fig. 3A). Higher MTT activity was observed for the hiPSC spheroids compared with the other four groups, either using the direct mixing or the pre-hiPSC aggregation methods. Spheroids of 0:100 hMSC and 25:75 iPSC-MSC generated through direct mixing method showed the relatively lower MTT activity compared with the other groups, whereas in the pre-hiPSC aggregation method, hybrid spheroids showed comparable MTT activity, which is much higher than hMSC spheroids.

Metabolic activity and DNA content of hybrid hiPSC-hMSC neural spheroids formed by different aggregation methods.
Similarly, the DNA content of the spheroids correlated with relative abundance of hiPSCs for the direct mixing method (Fig. 3Bi-a). hiPSC spheroids (100:0) showed the highest DNA content, whereas hMSC spheroids (0:100) showed the lowest DNA content. For pre-hiPSC aggregation method, DNA content showed similar trend to that observed in direct mixing method, which is dependent of the relative ratio of iPSC in spheroids, except 25:75 iPSC-MSC, which contained comparable DNA content to 50:50 iPSC-MSC (Fig. 3Bii-a). Again, the hMSC group showed the lowest DNA content, owing to the increased apoptosis due to upregulated caspase3/7 expression and altered mitochondria bioenergetics on 3D aggregation. 46 If the DNA content was converted into total cell number, the fold increase of 50:50 and 25:75 iPSC-MSC was higher than other groups for direct mixing method (Fig. 3Bi-b). For pre-iPSC aggregation, 25:75 iPSC-MSC group had the highest fold increase, indicating that this condition may promote hiPSC proliferation (Fig. 3Bii-b).
Expression of ECM and remodeling proteins in hybrid hiPSC-hMSC spheroids
To interrogate the influence of the coculture on ECM distribution in the hybrid spheroids, the expression of collagen type IV (Col IV) and Laminin (LN) was investigated by immunocytochemistry (Fig. 4). Both Collagen IV and laminin were homogeneously expressed in hiPSC aggregates and hMSC aggregates (Supplementary Fig. S5). In hybrid iPSC-MSC spheroids, MSCs were labeled with CellTracker Red to show colocalization with ECM proteins. Using the direct mixing method, the hybrid aggregates showed Col IV and LN expression (Fig. 4A). The group of 25:75 iPSC-MSC had higher Col IV expression, whereas LN expression seemed to increase with the ratio of hMSCs (Fig. 4C). For pre-hiPSC aggregation method, a higher level of Col IV expression was observed for the group of 25:75 iPSC-MSC. Higher LN expression was observed for the groups with higher MSC portion (Fig. 4B, D). Extensive expression of Col IV and LN was found around the hMSCs, indicating that differential hybrid spheroid structure may affect endogenous ECM dynamics.

Effects of different cell aggregation methods on collagen IV and laminin colocalization with MSCs.
Matrix metalloproteinases (MMPs) play a critical role in neural cell proliferation, migration, and differentiation. 47 In this study, MMP9 (a gelatinase) expression was observed for spheroids of different iPSC-MSC ratios (Fig. 5A). Moreover, the gene expression of MMP2 (another gelatinase) and MMP3 (a broader ECM cleavage enzyme) were analyzed by RT-PCR (Fig. 5). For direct mixing method, significantly higher expression of MMP2 (i.e., 39.7 ± 1.2, 45.7 ± 1.6, and 63.3 ± 0.5-fold for 75:25, 50:50, and 25:75, respectively) was observed for hybrid spheroids compared with the hiPSC spheroids (Fig. 5Bi). For MMP3, higher expression for hybrid spheroids (i.e., 7.5 ± 0.1, 2.4 ± 0.1, and 9.4 ± 2.7-fold for 75:25, 50:50, and 25:75, respectively) than 100:0 hiPSC group was observed. For pre-hiPSC aggregation method, higher expression of MMP2 (i.e., 6.7 ± 0.2, 9.4 ± 0.9, fold for 75:25 and 50:50, respectively) was observed for hybrid spheroids compared with hiPSC spheroids, and MMP3 expression was much higher (i.e., 23.7 ± 0.9, 41.3 ± 0.8, fold for 75:25 and 50:50, respectively) (Fig. 5Bii). These results indicate that incorporating MSCs in the hybrid spheroids upregulated MMP2 and MMP3 gene expression.

Effects of different cell aggregation methods on the expression of ECM remodeling proteins.
Cytokine secretion by hybrid hiPSC-hMSC spheroids
MSCs are significant source of TGF-β1,48–50 which exerts antiapoptotic and growth/migration-enhancing activities on NPCs. 51 PGE2 is a factor released by MSCs and has been shown to regulate multiple functions of different immune cells. 52 Both TGF-β1 and PGE2 were measured in this study. The results indicate that the 75:25 iPSC-MSC and 0:100 hMSC groups formed by direct mixing had higher normalized TGF-β1 concentration than the other groups (Fig. 6Ai). Using pre-hiPSC aggregation method, the normalized TGF-β1 concentration increases with the relative abundance of hMSCs (Fig. 6Aii). For PGE2, higher secretion was observed for 75:25 and 25:75 iPSC-MSC groups than the other groups for direct mixing method (Fig. 6Bi). For pre-hiPSC aggregation method, the normalized PGE2 concentration increases with the relative ratio of hMSCs in hybrid spheroids (Fig. 6Bii). These results indicate that incorporating MSCs at an appropriate ratio in the hybrid spheroids upregulates TGF-β1 and PGE2 expression compared with hiPSC spheroids.

TGF-β1 and PGE2 secretion by hybrid hiPSC-hMSC neural spheroids formed by different aggregation methods.
Neural differentiation of hybrid hiPSC-hMSC spheroids
Neural marker expression was assessed by immunocytochemistry after 14 days of neural differentiation of hybrid spheroids. Upregulated Pax6 (forebrain neural progenitors) and Nestin expression was observed for hybrid spheroids compared with hiPSC aggregates (Fig. 7A–D). Forebrain cortical neuron marker, TBR1, was also expressed for the hybrid spheroids at this early stage. The expression of Pax6 and TBR1 defines dorsal cortical identity of neural cells in forebrain. 53 Extensive expression of β-tubulin III, a marker for mature neurons, was observed for the aggregates of different iPSC-MSC ratios. Importantly, the neurite density of hybrid spheroids was increased with the ratio of hMSCs for pre-hiPSC aggregation method. However, little expression of Nestin, Pax6, β-tubulin III, and TBR1 was observed for the hMSC group. The expression of GFAP, a marker for astrocytes, was expressed for all five groups.

Neural marker expression of hybrid hiPSC-hMSC neural spheroids formed by different aggregation methods. Day 14 spheroids were replated for 3 days and immunocytochemistry was performed.
Neural markers, Nestin and β-tubulin III, were further quantified by flow cytometry (Fig. 8). For direct mixing method, Nestin expression was found higher for the hybrid spheroids (i.e., 47.4% ± 2.0%, 45.3% ± 0.1%, and 49.3% ± 5.3% for 75:25, 50:50, and 25:75, respectively) compared with the iPSC only group (38.7% ± 2.2%) (Fig. 8Ai, C). Similarly, for pre-hiPSC aggregation method, higher Nestin expression (i.e., 47.4% ± 2.0%, 40.5% ± 1.9%, and 43.0% ± 4.2% for 75:25, 50:50, and 25:75, respectively) was observed for the hybrid spheroids than the iPSC only group (35.7% ± 0.8%) (Fig. 8Aii, C). For β-tubulin III expression, higher levels were observed for hybrid spheroids (i.e., 59.0% ± 4.0%, 56.3% ± 9.8%, and 48.6% ± 2.1% for 75:25, 50:50, and 25:75, respectively) than the iPSC only group (42.4% ± 1.4%) using direct mixing method (Fig. 8Bi, D). Similarly, for pre-hiPSC aggregation method, β-tubulin III expression was higher for hybrid spheroids (i.e., 52.5% ± 0.8%, 73.2% ± 6.6%, and 61.4% ± 20.8% for 75:25, 50:50, and 25:75, respectively) than the iPSC only group (44.3% ± 5.7%) (Fig. 8Bii, D).

Nestin and β-tubulin III expression quantified by flow cytometry for hybrid hiPSC-hMSC spheroids formed by different aggregation methods. Flow cytometry histograms for spheroids formed at different iPSC to MSC ratios (100:0, 75:25, 50:50, and 25:75) at day 14 of neural differentiation.
Discussion
The microenvironment of neuropermissive matrix, trophic factors, and adherent molecules plays critical roles in neural tissue regeneration by forming guiding strands that direct axon regeneration.4,54 The cocultivation of MSCs with PSCs in 2-D promotes neural differentiation through synergistic paracrine interactions. 55 However, it remains unclear if direct cell–cell interactions take part in the neuropromoting effect of the coculture. On the other hand, the formation of 3D spheroids significantly enriched the secreted ECM proteins and growth factors with potent neurogenic functions (e.g., laminin and vascular endothelial growth factor [VEGF]).14,21,56 Moreover, the arrangement and the interactions between stromal cells and PSCs is changed in 3D aggregates, leading to the modulation of cell–cell contacts' downstream signaling (e.g., Notch). 57 However, their influences on neural differentiation are largely unknown. Therefore, this study investigates the influence of heterotypic cell–cell interactions on neural differentiation and the secretory functions of both cell populations, which have been poorly understood to date.
Dynamic regulation of heterotypic cell–cell interactions
In vivo, the complex tissues are integrated in a spatially restricted manner and provide dictatorial signals that organize sequential developmental events with tissue-specific function. 58 The in vitro culture aims to mimic the similar organization events, in which cells reorganize to minimize the surface energy of the 3D structure according to differential adhesion hypothesis (DAH) and the differential surface tension hypothesis. 17 In this study, hMSC condensation in the core of hybrid aggregates was observed during culture period. In direct mixing, different cell populations self-organize based on their relative degree of adhesiveness, and hMSCs may have stronger homotypic cell–cell interactions than hiPSC-NPCs. The observed core–shell spheroids are consistent with DAH predictions, which suggest that hMSCs with greater cohesiveness formed a core. In contrast, when mixed with preaggregated hiPSCs, hMSCs were found at the outer layer of the spheroids initially (with iPSC-NPC core and hMSCs as shell), which suggests that differential cell cortical tension between hiPSC-NPCs and hMSCs may dominate over homotypic cadherin–cadherin interactions. After the addition of MSCs, the formation of a cohesive hybrid spheroid requires additional energy input to destabilize the hiPSC-NPC aggregate structure, which results in prolonged kinetics to reach energy equilibrium. Therefore, mixing monodispersed hMSCs and precultured hiPSC spheroids suggested a potential method to control relative cell position in complex heterotypic 3D spheroids.
The formation of heterotypic spheroids with different ratios was reported to not influence the cell sorting process and the final state of self-assembly. 59 However, the kinetics to reach the surface energy balance is altered and correlated with the relative abundance of the two cell types. 59 Consistently, it was found that the 75:25 group required more time to reach the equilibrium than the 25:75 group for pre-hiPSC aggregation method. In addition, differential proliferation rate between the hiPSCs and hMSCs changed the ratio of mixed cell population over time and consequently the surface energy balance. 60 In direct mixing, cell proliferation had minimal influence on the final state, since surface energy-minimized 3D structure was already reached by day 2. However, in the preaggregation methods, the sustained imbalance in the cell ratio induced by the proliferation of hiPSCs in the 75:25 and 25:75 groups inversely correlated with the time required to reach the surface energy equilibrium. Taken together, the results suggest that both cell mixing ratios and aggregation methods regulate the kinetics of heterospheroid formation. The direct mixing method resulted in spheroids with an outer hiPSC layer and hMSC core. By contrast, the preaggregation method enabled a metastable surface energy status, where a ring of hMSCs covers hiPSCs.
The role of ECM secretion in heterotypic cell–cell interactions
HMSCs are an important source of ECM. The culture of hMSCs as spheroids was reported to enhance the secretion of ECM proteins such as fibronectin and laminin.20–22 Similarly, PSC-derived NPC aggregates produce fibronection, Col IV, and vitronectin. 44 Using the direct mixing method, Col IV and LN are distributed evenly within the heterotypic spheroids, suggesting that hMSCs and hiPSC-NPCs contribute cooperatively to the structure of heterotypic spheroids. In the pre-hiPSC aggregation method, a higher level of Col IV was observed in the aggregates formed with the higher hMSC percentage (25:75), indicating that the invasion of hMSC into the preformed hiPSC spheroids requires significant ECM remodeling.
MMPs are a group of proteolytic enzymes, which mediate ECM degradation and remodeling. 47 The cocultivation of hiPSCs and hMSCs in heterotypic spheroids revealed differential regulation of MMP-2 and MMP-3, and that MSCs are the major contributor to the MMP expression (Fig. 5). In the direct mixing method, MMP-2 expression increased with the portion of hMSCs in hybrid spheroids, which may contribute to the MSC condensation in the core of the hybrid spheroids. On the other hand, for the preaggregation method, a correlation between MMP-3 expression and the relative abundance of hMSCs in the heterotypic spheroids were observed, suggesting an important contribution of MMP-3 to facilitate hMSC migration and invasion in preaggregated hiPSC-NPCs.
The role of cytokine secretion in heterotypic cell–cell interactions
HMSCs and NPCs are important sources of soluble molecules, which play significant roles in autocrine and paracrine effects in 3D cultures. For instance, MSC spheroids formed through hanging drop or forced aggregation have been found to enhance the expression of antiapoptotic and anti-inflammatory proteins such as TSG-6 and STC-1, anticancer molecules such as IL-24 and CD82,15–18 as well as the angiogenic factors such as VEGF.61,62 Similarly, NPCs secrete platelet-derived growth factor, bone morphogenic protein, epidermal growth factor, basic fibroblast growth factor, and Wnt ligands, which regulate proliferation and promote early differentiation into neuronal and glial lineages. 63
In this study, the presence of MSCs (e.g., 75:25 iPSC-MSC) in hybrid spheroids showed significantly higher secretion of PGE-2, either using direct mixing or the preaggregation method compared with the hiPSC only group. In addition, the secretion of TGF-β1 increased with the relative abundance of hMSCs. Mild hypoxia or autophagy may be involved in the enhanced secretion of PGE-2 and TGF- β1 by hMSCs.64,65 The specific location of hMSCs (i.e., the core) in the hybrid aggregate may be linked with the increased soluble factor secretion. Moreover, nitric oxide is a potent inducer of PGE-2 secretion by NPCs, which can also be activated by the hypoxia induced in large aggregates. 66 In addition, PGE-2 induces profound ECM remodeling mediated by the induction of MMPs, 67 which play a role in the growth factor secretion by hMSCs. 68 Thus, it is likely that PGE-2 secretion precedes ECM remodeling in the hybrid spheroids.
The molecular mechanisms of the increased expression of PGE-2 in hMSC aggregates may involve the regulation of NF-κB upon aggregation (e.g., through caspase and IL-1α signaling). 69 Indeed, the inhibition of NF-κB signaling by small molecules (e.g., QNZ, a NF-κB inhibitor) was found to abrogate PGE-2 secretion in hMSC aggregates. 69 Of interest, it was found that the secretion of TGF-β1 and PGE-2 in neural progenitors is mediated by inflammatory molecules (i.e., TNF-α and/or IFN-γ), which triggers the NF-κB signaling. 70 Consequently, the increased secretion of PGE-2 suggested cooperative role between hMSCs and hiPSC-NPCs in the hybrid spheroids, which may be promoted by the activation of NF-κB signaling in the two cell types.
The role of MSC in neural differentiation
MSC-derived ECM proteins can enhance neural differentiation, prolong cell survival, and provide better protection from nutrient and growth factor deprivation. 23 However, modulation of cellular microenvironment indicated by the secretion of trophic factors and ECM proteins has not been well understood during hybrid spheroid formation. In this study, it was found that hiPSCs are the main cell type capable of cortical neuronal differentiation in the hybrid spheroids, while the expression of neural markers was low in hMSCs. Using direct mixing and preaggregation methods, the presence of hMSCs in the hybrid spheroids promoted neuronal differentiation (i.e., Nestin and β-tubulin III expression) compared with the hiPSC only group. Our differentiation procedure should promote neurons of dorsal cortical identity of forebrain (expressing markers Pax6 and TBR1).42,53 These results suggest an important paracrine role of hMSCs for the regulation of hiPSC fates, which are mediated by direct interactions between the two cell types.
ECM remodeling may also be an important contributing factor that enhances neural differentiation. Hybrid spheroids composed of a higher ratio of hMSC-promoted MMP expression, which is characterized by a balance between degradation and neosynthesis of structural ECM proteins. Indeed, the formation of hMSC aggregates upregulates the expression of several types of MMPs (MMP-2, -9, and -1/13). 71 Moreover, activation of MMPs was reported to enhance stem cell differentiation toward neuronal lineage, through the NF-κB signaling. 72 On the other hand, the formation of hMSC spheroids was also reported to enhance ECM secretion.20–22 The secretion of ECM proteins by hMSCs such as collagen IV and laminin in the heterotypic spheroids may affect hiPSC fate and promote neuronal differentiation.
Previous studies demonstrated that the secretome of hMSCs also promoted neural differentiation of NPCs.73,74 In particular, TGF-β1 exerts antiapoptotic and growth/migration-enhancing activities on NPCs through SMAD2 (for growth and migration) and ERK1/2 (for antiapoptotic behavior) signaling, 51 which can promote β-tubulin III expression. 75 In our study, TGF- β1 may partially contributes to the enhanced neuronal differentiation of hiPSCs. Alternatively, PGE-2 or proinflammatory environment was demonstrated to promote neuronal differentiation, through the regulation of the NF-κB signaling.76,77 These observations are consistent with results of the present study, which demonstrate that the presence of hMSC in the hybrid spheroid promotes the levels of Nestin and β-tubulin III.
Conclusions
The modulation of heterotypic interactions between hMSCs and hiPSCs, imposed by different mixing methods and relative population abundance, dynamically regulates cellular self-organization in the hybrid dorsal cortical spheroids. The dynamics of cellular spatial organization in the hybrid spheroids promote the gene expression of ECM remodeling proteins and the secretion of PGE-2 and TGF-β1. The microenvironment changes induced by ECM and soluble factor secretion (higher at 25:75 iPSC to MSC ratio in general) due to MSC-hiPSC interactions cooperatively promote cortical neural differentiation of hiPSCs.
Footnotes
Acknowledgments
The authors thank Ms. Ruth Didier of FSU Department of Biomedical Sciences for her help with flow cytometry analysis, Dr. Brian K. Washburn and Kristina Poduch in FSU Department of Biological Sciences for their help with RT-PCR analysis, Dr. Stephen Duncan at Medical College of Wisconsin and Dr. David Gilbert in FSU Department of Biological Sciences for human iPSK3 cells. This work is supported by FSU start up fund, FSU Bridge Grant, NSF Career Award (No. 1652992 to Y.L.), and by NIH R03EB020770 (Y.L.). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
