Abstract
Human induced pluripotent stem cells (hiPSCs) emerge as a promising source to construct human brain-like tissues, spheroids, or organoids in vitro for disease modeling and drug screening. A suspension bioreactor can be used to generate large size of brain organoids from hiPSCs through enhanced diffusion, but the influence of a dynamic bioreactor culture environment on neural tissue patterning from hiPSCs has not been well understood. The objective of this study is to assess the influence of a suspension bioreactor culture on cortical spheroid (i.e., forebrain-like aggregates) formation from hiPSCs. Single undifferentiated hiPSK3 cells or preformed embryoid bodies were inoculated into the bioreactor. Aggregate size distribution, neural marker expression (e.g., Nestin, PAX6, β-tubulin III, and MAP-2), and cortical tissue patterning markers (e.g., TBR1, BRN2, SATB2, and vGlut1) were evaluated with static control. Bioreactor culture was found to promote the expression of TBR1, a deep cortical layer VI marker, and temporally affect SATB2, a superficial cortical layer II–IV marker that appears later according to inside-out cortical tissue development. Prolonged culture after 70 days showed layer-specific cortical structure in the spheroids. Differential expression of matrix metalloproteinase-2 and −3 was also observed for bioreactor and static culture. The altered expression of cortical markers by a suspension bioreactor indicates the importance of culture environment on cortical tissue development from hiPSCs.
Introduction
R
A dynamic spinner bioreactor culture enhances diffusion in neural progenitor cell (NPC) spheres derived from hiPSCs and allows further growth and maturation of large size aggregates (up to 4 mm in diameter).10,18 While the ultimate solution to generate large size of forebrain tissues needs the vascularization of the constructs,4,19 bioreactor culture provides an immediate and simple approach for the human PSC (hPSC)-derived aggregate development into mini-tissues through enhanced diffusion. 20 The suspension-based stirred bioreactor is simple, scalable, and has the potential to generate a large quantity of cortical organoids for downstream applications in drug screening and disease modeling. Generation of stem cell spheroids or organoids in a suspension bioreactor has been demonstrated for mesenchymal stem cells, hPSC-derived cardiac cells, and hPSC-derived NPCs.20–24 The size of stem cell aggregates can be regulated by agitation speed/shear stress, which indirectly affects the diffusion of biomolecules.25–27 Collectively, the simple setup of suspension bioreactor allows us to delineate the complex biological problems involved in cortical tissue development from hPSCs. However, the influence of suspension bioreactors on brain tissue structure has not been well explored.
The objective of this study is to investigate the effect of dynamic bioreactor culture on the formation of cortical spheroids from hiPSCs. Current systems have derived cortical cells or spheroids from hPSCs in static cultures.2,28 How dynamic bioreactor culture impacts the development of cortical tissues from hPSCs remains unclear. In addition, current methods to derive cortical organoids from hPSCs involve long-term culture (up to 140 days). It would be advantageous if the developmental time is shortened by modulating the culture system. In this study, it is hypothesized that cortical tissue development from hiPSCs may be accelerated by the dynamic culture due to enhanced diffusion of patterning factors. So this study compared the aggregate formation from hiPSCs, neural patterning marker expression, and the expression of cortical layer-specific genes in a spinner bioreactor and a static culture. To evaluate the temporal effect of bioreactor culture, the inoculation with single hiPSC suspension and preformed embryoid bodies (EBs) was investigated. The outcome of this study should provide the knowledge on cortical tissue development from hPSCs and a transformative platform to generate cortical spheroids or organoids for disease modeling, drug screening, and possibly 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).29,30 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 six-well plate in the presence of 10 μM Y27632 (Sigma) for the first 24 h. 31
Neural differentiation and cortical spheroid formation from hiPSCs in static culture
For neural differentiation in static culture, human iPSK3 cells were seeded into ultra-low attachment 24-well plates (Corning Incorporated, Corning, NY) at 4.0–4.5 × 105 cells per well in 1 mL of differentiation medium composed of Dulbecco's Modified Eagle Medium/Nutrient Mixture F-12 (DMEM/F-12) plus 2% B27 serum-free supplement (Life Technologies). 32 Y27632 (10 μM) was added during the seeding and removed after 24 h. On day 1, the cells formed EBs and were treated with no growth factor or with dual SMAD signaling inhibitors: 10 μM SB431542 (Sigma) and 100 nM LDN193189 (Sigma), hereafter referred to as LDN/SB. After 7 days, the cells were treated with cyclopamine (1 μM) and fibroblast growth factor-2 (10 ng/mL) until day 14. The suspension cultures were maintained for 7 days in short-term study or for 21–71 days in long-term study. For immunocytochemistry characterizations, the aggregates were replated onto Geltrex-coated surface for 5–7 days.
Preparation of spinner bioreactors
A 50 mL glass spinner bioreactor (#356875; Wheaton) was used for studying NPC differentiation and neural tissue patterning from hiPSCs. Before the inoculation, the vessel was immersed with 2 mL Sigmacote (#SL2; Sigma) overnight to avoid cell attachment on the glass surface. Then, the vessel was dried at room temperature. The spinner bioreactor was autoclaved at 121°C for 20 min before using. Wheaton Micro-Stir platform (#W900701-A; Wheaton) was used for bioreactor culture. The Micro-Stir platform allows for the control of agitation speed and the pattern (e.g., intermittent agitation with different cycles). 33 The Micro-Stir platform was put into a 37°C, 5% CO2 incubator before setting up bioreactor culture (Fig. 1A).

Short-term bioreactor culture for no growth factor induction.
Neural induction and cortical spheroid formation from hiPSCs in bioreactor culture
The bioreactor was inoculated with day 0 undifferentiated hiPSK3 single cell suspension in the presence of Y27632 or day 8 EBs, in parallel to static control. For undifferentiated cells, the seeding density was 4.0–4.5 × 105 cells/mL. Alternatively, the EBs formed with the equivalent cell density to static culture were put in the bioreactor. The working volume of the bioreactor was 15 mL. After single-cell inoculation, the bioreactor was agitated at an intermittent mode (15-min on/15-min off cycle for 12 h) at 80 rpm. After that, a constant agitation (100 rpm) was maintained until the harvest time (day 7–67 in bioreactor). For EB inoculation, the bioreactor was agitated at a constant speed (100 rpm) until the harvest time. Complete medium change was performed every 3–4 days. During the culture time, samples were taken to monitor aggregate size distribution, glucose consumption, and lactate production. The cell numbers were determined at the harvest time by counting a 2 mL aggregate sample. The harvested aggregates were either replated to Geltrex-coated surface for another 4–7 days for immunocytochemistry, flow cytometry, or reverse transcription–polymerase chain reaction (RT-PCR) analysis, or kept in low-attachment plates (up to a total of 71-day culture) for confocal microcopy analysis.
Aggregate size distribution
The images of aggregates (both the bioreactor and static culture in the side-by-side comparison of the same experiment) were captured over the culture time by 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 provided the size distribution of the aggregates.
Biochemical assays: cell number and glucose/lactate assay
Cell numbers at the harvest time of bioreactor were determined using a hemocytometer after cell trypsinization of the aggregates. The supernatants were analyzed for glucose and lactate concentrations with a YSI 2700 Select™ Bioanalyzer (Yellow Spring Instruments, Yellow Spring, OH). The lactate yield from glucose (mol/mol) was calculated based on the ratio of lactate production to glucose consumption.
Bromodeoxyuridine (BrdU) assay
Briefly, the aggregates were incubated in a medium containing 10 μM BrdU (Sigma) for four hours. 34 The cells were then fixed with 70% cold ethanol, followed by a denaturation step using 2 N HCl/0.5% Triton X-100 for 30 min in the dark. The samples were reduced with 1 mg/mL sodium borohydride for 5 min and incubated with mouse anti-BrdU (1:100; Life Technologies) in a blocking buffer (0.5% Tween 20/1% bovine serum albumin in phosphate-buffered saline [PBS]), followed by Alexa Fluor® 488 goat anti-Mouse IgG1 (Molecular Probes). The cells were counterstained with Hoechst 33342 and viewed under a confocal microscope (Zeiss LSM 880).
Immunocytochemistry
Briefly, the samples were fixed with 4% paraformaldehyde (PFA) and permeabilized with 0.2–0.5% Triton X-100 for intracellular markers. The samples were then blocked and incubated with various mouse or rabbit primary antibodies (Supplementary Table S1; Supplementary Data are available online at www.liebertpub.com/tea). After washing, the cells were incubated with the corresponding secondary antibody: Alexa Fluor 488 goat anti-Mouse IgG or Alexa Fluor 594 goat anti-Rabbit IgG (Life Technologies). The samples were stained with Hoechst 33342 and visualized using a fluorescent microscope (Olympus IX70, Melville, NY) or a confocal microscope (Zeiss LSM 880). The images from five independent fields (800–1000 cells) were analyzed using ImageJ software. The proportion of positive cells was calculated based on the area of marker of interest normalized to the nuclei using ImageJ analysis, indicating the relative expression among different conditions.
Flow cytometry
To quantify the levels of various marker expressions, the cells from the bioreactor and static control were harvested by trypsinization and analyzed by flow cytometry. 35 Briefly, 1 × 106 cells per sample were fixed with 4% PFA and washed with a staining buffer (2% fetal bovine serum in PBS). The cells were permeabilized with 100% cold methanol for intracellular markers, blocked, and then incubated with primary antibodies against Nestin, PAX6, β-tubulin III, or SATB2 followed by the corresponding secondary antibody (Supplementary Table S1). The cells were acquired with BD FACSCanto™ II flow cytometer (Becton Dickinson) and analyzed against isotype controls using FlowJo software. At least 2 × 105 cells were analyzed for each sample.
RT-PCR
Total RNA was isolated from neural cell samples using the RNeasy Mini Kit (Qiagen, Valencia, CA) according to the manufacturer's protocol followed by the treatment of DNA-Free RNA Kit (Zymo, Irvine, CA). 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 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 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
To assess the statistical significance of the experimental results, Student's t-test was performed on the measurements from three replicates of bioreactor and static samples. For better comparison, the data from the side-by-side differentiation for bioreactor and static culture in the same experiment were presented. The results were expressed as [mean ± standard deviation]. A p-value <0.05 was considered statistically significant.
Results
Short-term bioreactor culture for no growth factor induction and LDN/SB induction
Neural differentiation of hiPSK3 cells in the bioreactor was first evaluated for 7 days using a neural medium without any growth factor and inoculated with single hiPSK3 cell suspension (Fig. 1). On day 1, bioreactor condition had larger EB-like aggregates than static condition (both seeded at 0.42 × 106 cells/mL). With the culture time, the aggregate size in both conditions increased and the difference was not significant (Fig. 1B). The static condition had larger aggregates by day 7 than the bioreactor condition, probably due to the merge of aggregates since the bioreactor reduced the frequency of EB merging due to agitation. The size distribution of aggregates was analyzed through image analysis (Fig. 1C). By day 7, the average diameter was 301 ± 67 μm for bioreactor condition and 394 ± 100 μm for static condition (Supplementary Fig. S1A). Lower lactate (mol)/glucose (mol) ratio was observed for bioreactor culture than static culture (1.70 ± 0.01 vs. 1.94 ± 0.08) (e.g., Run 1 in Table 2 and Supplementary Fig. S1B). However, the range of lactate/glucose ratio showed large variations in different runs (e.g., Run 3 in Table 2, 2.13 ± 0.08 vs. 2.73 ± 0.04 for bioreactor vs. static). The harvested cell density (day 7) was comparable for the two conditions (e.g., 1.13–1.26 × 106 cells/mL in Run 1 of Table 3). Variations from different runs were observed (e.g., 0.79–0.83 × 106 cells/mL for Run 2). If cell source was not good, no good aggregates were formed (data not shown). The aggregates were harvested and replated for another 7 days. The expression of Nestin and PAX6 was similar for both conditions (Fig. 1D, E). The cells also expressed patterning marker TBR1 (a forebrain layer VI cortical marker) and HOXB4 (a hindbrain/spinal cord marker) (Fig. 1E). RT-PCR analysis for HOXB4, HB9 (motor neuron marker), TBR1, and vGluT1 (a glutamatergic neuron marker) showed higher TBR1 expression and lower vGluT1 expression for bioreactor culture than static culture (Fig. 1F). The expression of HOXB4 and HB9 was comparable.
Statistical difference (p-value <0.05) between bioreactor and static cultures for no growth factor method.
Statistical difference (p-value <0.05) between bioreactor and static cultures for LDN/SB induction method.
Then, the neural differentiation was induced using LDN193189 and SB431542 for 7 days and inoculated with single cell suspension (Fig. 2). Similarly, on day 1, bioreactor condition had larger size of aggregates than static condition, but the difference was not significant at later time points. With time, the aggregate size increased and the static condition had larger aggregates by day 7 than the bioreactor condition (274 ± 122 μm vs. 416 ± 164 μm) (Fig. 2A and Supplementary Fig. S2). The lactate (mol)/glucose (mol) ratio for the bioreactor culture was lower than the static culture (1.58 ± 0.21 vs. 2.10 ± 0.17) (e.g., Run 1 in Table 2, similar observation was found for additional runs) and the harvested cell density (day 7) was comparable (0.65–0.89 × 106 cells/mL) (Table 3). The harvested aggregates were replated for neural marker analysis (Fig. 2B). The cells expressed patterning markers of Nestin, PAX6, TBR1m, and HOXB4, as well as the more mature neural marker β-tubulin III (Fig. 2C, D). The quantification results showed higher PAX6 and TBR1 expression for bioreactor culture than static culture, but no difference in Nestin and HOXB4 expression was observed at this time point (Fig. 2B, C), showing the influence of culture system on some neuronal markers. The expression of β-tubulin III was comparable (41–46%), evaluated at 10 days after replating of day 7 aggregates (Fig. 2D). These results indicate that bioreactor culture may alter the expression of forebrain patterning markers.

Short-term bioreactor culture for LDN/SB induction.
Long-term bioreactor culture for single cell inoculation
Prolonged neural differentiation to generate cortical spheroids was then performed for 21 days in the bioreactor in the presence of LDN193189 and SB431542 (Fig. 3A). With the increased culture time, the aggregate size increased, and the bioreactor condition had larger aggregates than the static condition after day 10 (Fig. 3B, C). By day 21, the average diameter was 1047 ± 174 μm for bioreactor condition and 803 ± 167 μm for static condition (Supplementary Fig. S3). The harvested cell density (day 21) was 0.93–1.38 × 106 cells/mL for bioreactor culture and static culture (Table 3). The structure of cortical spheroids (day 21) was evaluated by confocal microscopy (Fig. 4Ai). The expression of TBR1 was observed in the spheroids, but the expression of SATB2 (a marker for superficial cortical II–IV layer) was minimal. The aggregates also expressed HOXB4 and β-tubulin III. The expression of presynaptic marker synapsin I and postsynaptic marker PSD95 was observed (day 28) (Fig. 4Aii). RT-PCR analysis showed higher HOXB4 and TBR1 expression, but lower vGluT1 expression for bioreactor culture than static culture (Fig. 4B). Gene expression of matrix metalloproteinase (MMP)-2 and −3 was evaluated to show the influence of culture environment on matrix remodeling (Fig. 4C and Supplementary Fig. S4). MMP-2 expression was lower, but MMP-3 expression was higher for bioreactor culture than static culture. These results demonstrate the development of early stage cortical spheroids from hiPSCs, and bioreactor culture increased the gene expression of HOXB4 and TBR1 and altered MMP-2/-3 expression.

Long-term bioreactor culture for single cell inoculation—aggregate characterization.

Long-term bioreactor culture for single cell inoculation—neural differentiation characterization. The day 21 aggregates from bioreactor culture and static culture were harvested and characterized.
Long-term bioreactor culture for EB inoculation
To evaluate if the exposure to bioreactor culture has a temporal effect, the bioreactor was inoculated with day 8 EBs after the ectoderm induction by dual SMAD inhibitors (Fig. 5A). With the increased culture time, the aggregate size increased. The bioreactor condition had larger aggregates than static condition, probably due to better nutrient diffusion (Fig. 5B, C). By day 27, the average diameter was 1121 ± 91 μm for bioreactor condition and 997 ± 152 μm for static condition (Supplementary Fig. S3). The replated spheroids expressed cortical markers of TBR1 and glutamate, as well as synaptic markers synapsin I and PSD95 (Fig. 6A). HOXB4 expression was also visible. RT-PCR analysis showed higher HOXB4 and TBR1 expression, but slightly lower vGluT1 expression for bioreactor culture than static culture (Fig. 6B). The expression of SATB2 (day 32) showed 37% for bioreactor culture and 14% for static culture (Fig. 6C. In Supplementary Fig. S5, SATB2 is 32 ± 15% for bioreactor and 13 ± 2% for static culture). High expression of β-tubulin III (78–92%) was observed for both cultures. These data indicate that bioreactor culture increases the expression of TBR1, SATB2, and HOXB4 at the neural patterning stage after day 8 of ectoderm induction.

Long-term bioreactor culture for EB inoculation—aggregate characterizations.

Long-term bioreactor culture for EB inoculation—neural differentiation characterization. On day 27, aggregates from bioreactor and static cultures were replated and characterized.
Characterization of long-term cortical spheroids derived from hiPSCs
Long-term culture of cortical spheroids in bioreactor inoculated with the preformed EBs was performed (Fig. 7). The structure of cortical spheroids was evaluated by confocal microscopy for the expression of TBR1 and SATB2 (Fig. 7A). While SATB2 expression appeared on day 43, no clear layer structure was observed. The expression of β-tubulin III was around 88–89% for both cultures (Fig. 7B), distributed throughout the spheroids (Supplementary Fig. S6). The expression of SATB2 was 40.1% for bioreactor and 50.6% for static culture on day 44, indicating that static culture reached the similar level of SATB2 expression after a longer culture period (Fig. 7B). Another superficial layer II–IV marker BRN2 (layer III) was observed on day 48 (Fig. 7C). Prolonged culture was performed up to day 71 and the layered cortical structure was more visible (Fig. 7D), with a specific superficial layer of SATB2+ cells formed next to the deep layer of TBR1+ cells. The expression of TBR1 and SATB2/BRN2 at different time points demonstrates that the cortical spheroids followed “inside-out” developmental pattern of human cortex (Supplementary Fig. S7). Moreover, co-localization of BRN2 and BrdU was observed, indicating the proliferation capability of superficial layer cells (Fig. 7E). By contrast, TBR1 and BrdU expression were not co-localized, indicating the formation of deep layer cells at the early stage of development before the appearance of superficial layer cells. 2 These results demonstrate the formation of layer-specific structure of cortical spheroids from hiPSCs in a long-term culture.

Layer-specific structural characterizations of long-term cortical spheroids derived from human induced pluripotent stem cells.
Additional characterizations confirmed the expression of cortical markers BRN2 and TBR1 (Supplementary Fig. S8A). GABAergic interneurons were also observed based on the expression of GABA and vGAT (about 15–20%) (Supplementary Fig. S8B). In addition, mature neuron markers MAP-2 and tau, the microtubule-associated proteins involved in pathological development of various neurological diseases,36,37 were expressed in the replated spheroids (Supplementary Fig. S8C). The cell populations also have some astrocytes indicated by the expression of GFAP (Supplementary Fig. S8D). These data further indicate the cortical tissue properties of the derived spheroids.
Discussion
Recent developments in recapitulating neurodegenerative diseases and the findings of differential susceptibility of different neuronal subtypes demonstrate the critical needs for novel in vitro brain-like tissue models for disease modeling and drug screening.13,14,38–40 Forebrain cortical neurons, spheroids, and organoids have wide applications in neurological disorders, including Alzheimer's disease, microcephaly, brain injury, and stroke.2,9,18,41 Human cortex has six layers with the deep layer VI neurons (TBR1+ early-born neurons) being generated first and the superficial layer II–IV neurons (SATB2+ late-born neurons) appearing later, 42 following an “inside-out” developmental pattern. 43 However, the brain-like tissue models not only need to contain relevant cell populations but also the 3D organization of appropriate structure. The forebrain-like 3D cortical tissue development from hPSCs has been investigated in vitro in static culture environment.28,42,44 Comparison of 3D aggregate-based neural differentiation with 2D monolayer differentiation showed that aggregate-based neural differentiation resulted in higher MAP-2+ cells (93% vs. 45% at day 40). 45
Despite rapid development of in vitro generation of 3D mini-brain or organoid models from hPSCs, the size and the complexity of the neural structures are limited by the culture systems and inadequate supply of nutrients and oxygen. 46 Bioreactors can enhance the diffusion of biomolecules and patterning factors and were used for generating large size of organoids in a couple of studies,10,18 but the exact role of bioreactor system in brain organoid development remains unknown. Our study indicates that bioreactor culture promotes the expression of TBR1 (a forebrain cortical layer marker) and HOXB4 (a hindbrain/spinal cord marker), affecting neural tissue patterning in brain spheroids/organoids derived from hiPSCs.
The bioreactor culture has some possible effects on cortical tissue development based on our results. On one hand, the dynamic culture environment exposes shear stress to the surface of the aggregates, which causes signal transduction of mechanical forces, leading to the activation of intracellular pathways such as Wnt (Supplementary Fig. S9),26,27 which may contribute to the upregulation of HOXB4. 47 On the other hand, the dynamic culture environment enhances the diffusion of neural patterning factors (e.g., cyclopamine), nutrients, and oxygen and promotes cortical tissue development,48,49 which may contribute to the upregulation of TBR1.
Based on the geometry of the bioreactor used in this study, the estimated maximum shear stress is 2.3–3.4 dyne/cm2, 50 although most spheroids would not experience the impeller maximum shear stress. Shear stress in bioreactors alters the kinetics of ligand-receptor binding and changes the frequency of cell collision, which may also lead to altered cell phenotype.25,51 The fluidic pattern in a spinner bioreactor shows the maximum flow velocity at the edge of the impeller, and the lowest flow velocity in the center under the impeller based on our computational fluidic dynamics (CFD) analysis (Supplementary Fig. S10), 52 which contributes to the heterogeneity of aggregate size. The detailed analysis of predicting the actual stresses experienced by the spheroids needs more complex CFD analysis to be performed in future.
In addition, bioreactor culture alters the concentration gradients inside the aggregates in a more controllable way48,49 and thus may exert differential impact on neural patterning from hPSCs compared to static culture. For example, the sonic hedgehog (Shh) antagonist cyclopamine would show a deeper penetration into the aggregates for bioreactor culture than static culture, influencing the ventralization process; so the inhibition of Shh may enhance forebrain cortical marker TBR1 expression. Moreover, bioreactor culture enhances oxygen diffusion, which plays a role in cell metabolism. Further studies are required to establish the direct correlation of oxygen and/or other metabolic state factors with cortical tissue development. These differences in culture environment may collectively contribute to the difference of cortical tissue development from hiPSCs between the bioreactor culture and static culture.
The bioreactor culture of day 8 preformed EBs also showed the elevated gene expression of HOXB4 and TBR1, indicating that the temporal exposure to dynamic culture environment is important during neural patterning stage (day 8–15) rather than neural ectoderm induction stage (day 0–7). Neural patterning stage involves multiple signaling pathways to modulate the regional specificity of brain tissue along anterior-posterior (A-P) axis and dorsal-ventral (D-V) axis.53,54 It was thought that the gradients of signaling molecules during neural patterning stage play a critical role in spatially defining the neuronal subtypes in the 3D neural spheroids/organoids. 54
Besides the bulk fluid dynamic environment, it was thought that the interstitial flow (0.1–1.0 dyne/cm2) inside the aggregates may play an important role in tissue morphogenesis and development. 55 The parenchymal cells of brain are surrounded by interstitial fluid. 56 Regulation of the fluid composition and volume is achieved by the barriers between central nervous system and blood. The interstitial space is usually occupied by ECMs, which determine the resistance to fluid flow. Therefore, interstitial flow is very slow and generates low levels of shear stress. Besides enhancing nutrient diffusion, the interstitial flow affects ECM deposition and cytokine production from the cells. Shear stress at interstitial levels (0.135 dyne/cm2) may affect the nuclear localization of transcriptional coactivator with PDZ binding motif (TAZ), 57 an effector protein and signaling mediator of mechanotransduction along with Yes-associated protein,58,59 which could impact cellular differentiation. 57
Based on the knowledge of interstitial flow, it is speculated that stirred agitation in a spinner bioreactor may provide interstitial flow inside the aggregates that is felt by the cells and promotes nutrient diffusion in our study. The embryonic stem cell (ESC) aggregate porosity (or void space) is estimated to be about 0.26–0.48.48,60 If the aggregate surface experiences 3–15 dyne/cm2 (the range of blood flow shear), the flow penetration into the aggregates (tissue-like structure) would be at the interstitial level (0.1–0.2 dyne/cm2) due to the resistance from ECMs among the cells in the aggregates. 55 Such interstitial flow would be minimal in static culture. It is possible that the interstitial flow inside the aggregates of bioreactor culture contributes to the alter expression of neural patterning makers (TBR1 and HOXB4) observed in this study.
The results from this study indicate that MMP-2 and MMP-3 expression can be altered by the culture system. MMP-2 expression was slightly decreased, while MMP-3 expression was increased in bioreactor culture compared to static culture. MMP-2, also known as gelatinase A, mainly degrades type IV collagen, the most abundant basement membrane protein. 61 MMP-3, also known as stromelysin-1, degrades collagen types II, III, IV, IX, and X, proteoglycans, fibronectin, laminin, and elastin. MMP-3 is also an activator for pro-MMP-9, which regulates MMP-9 activity. 62 MMP-3 deficiency has been shown to result in the delayed granule cell radial migration and disturbance in synaptogenesis in Purkinje cells of cerebellar cortex.63,64 Our study demonstrates the effect of dynamic culture on forebrain-like cortical tissue development from hiPSCs by temporally modulating ECM remodeling proteins. The dynamic culture of ESC aggregates has been shown to modulate the ECM (e.g., laminin) and growth factor (e.g., bone morphogenetic protein 4) expression pattern compared to static culture, during the development of mesodermal lineages. 65 The ECMs in brain tissues consist of basement membrane proteins (e.g., laminin, fibronectin, collagen IV, and heparan sulfate proteoglycans), perineuronal nets (e.g., hyaluronic acid and chondroitin sulfate proteoglycans), and interstitial matrix (e.g., laminin and tenascin R).66,67 Investigation of ECM remodeling and composition affected by the culture systems may reveal the mechanisms of brain tissue memory and neurological disease progression, which needs further study.
It shall be noted that the reported system in the current form cannot fully address the robustness in the production of spheroids. Large variations exist in metabolic activity and aggregate size, where cell source may also play a role. In the future, several strategies may be explored to address the issues related to the variability in spheroid/organoid production. For examples, controlled input cell population and defined matrix components may be used together with optimized feeding strategy. They may help reduce the spheroid-spheroid or organoid-organoid variability and improve the reproducibility of the experimental systems.
Conclusions
This study indicates that bioreactor culture alters the expression of neural patterning markers during the cortical spheroid formation from hiPSCs compared to static culture. The results demonstrate the importance of culture environment in cortical tissue development from hiPSCs during the establishment of in vitro 3D brain-like tissue models. This study has significance in deriving brain organoids from hPSCs for disease modeling and drug discovery.
Footnotes
Acknowledgments
The authors would like to thank Ms. Ruth Didier in 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, and Dr. Ravindran Chella of FSU Department of Chemical and Biomedical Engineering for his help in CFD analysis. This work is supported by FSU start up fund, FSU Bridge Fund, and National Science Foundation (Grant No. 1652992).
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.
