Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Neural Cell 3D Microtissue Formation Is Marked by Cytokines' Up-Regulation

  • Yinzhi Lai ,

    Contributed equally to this work with: Yinzhi Lai, Amish Asthana

    Affiliation Cellular Bioengineering Laboratory, Department of Biological and Agricultural Engineering, Faculty of Engineering, Driftmier Engineering Center, University of Georgia, Athens, Georgia, United States of America

  • Amish Asthana ,

    Contributed equally to this work with: Yinzhi Lai, Amish Asthana

    Affiliation Cellular Bioengineering Laboratory, Department of Biological and Agricultural Engineering, Faculty of Engineering, Driftmier Engineering Center, University of Georgia, Athens, Georgia, United States of America

  • Ke Cheng,

    Current address: Heart Institute, Cedars-Sinai Medical Center, Los Angeles, California, United States of America

    Affiliation Cellular Bioengineering Laboratory, Department of Biological and Agricultural Engineering, Faculty of Engineering, Driftmier Engineering Center, University of Georgia, Athens, Georgia, United States of America

  • William S. Kisaalita

    williamk@engr.uga.edu

    Affiliation Cellular Bioengineering Laboratory, Department of Biological and Agricultural Engineering, Faculty of Engineering, Driftmier Engineering Center, University of Georgia, Athens, Georgia, United States of America

Abstract

Cells cultured in three dimensional (3D) scaffolds as opposed to traditional two-dimensional (2D) substrates have been considered more physiologically relevant based on their superior ability to emulate the in vivo environment. Combined with stem cell technology, 3D cell cultures can provide a promising alternative for use in cell-based assays or biosensors in non-clinical drug discovery studies. To advance 3D culture technology, a case has been made for identifying and validating three-dimensionality biomarkers. With this goal in mind, we conducted a transcriptomic expression comparison among neural progenitor cells cultured on 2D substrates, 3D porous polystyrene scaffolds, and as 3D neurospheres (in vivo surrogate). Up-regulation of cytokines as a group in 3D and neurospheres was observed. A group of 13 cytokines were commonly up-regulated in cells cultured in polystyrene scaffolds and neurospheres, suggesting potential for any or a combination from this list to serve as three-dimensionality biomarkers. These results are supportive of further cytokine identification and validation studies with cells from non-neural tissue.

Introduction

Providing a 3D spatial microenvironment for cells to grow in, is the sole criterion that has traditionally been associated with three-dimensional cell culture. However, with recent advances in the field in the past decade, the meaning of 3D cell culture has been extended to providing the “total microenvironment” that supports the formation of microtissue that exhibits “complex” physiological relevance (CPR) or better emulation of the in vivo functionality in a manner not possible in 2D cultures [1]. Three main categories or microenvironment factors (MEFs) or “three-dimensions” from the literature include: 1) chemical or biochemical composition, 2) spatial (geometric 3D) and temporal dimensions, and 3) force and substrate physical properties [1][4]. However, there is still a lack of a quantifiable entity which can establish if the cellular response in a 3D culture is actually physiologically relevant and in vivo-like or just different from 2D. The identification and validation for this entity or a potential three-dimensionality biomarker is necessary due to three compelling reasons. First, apart from the concept of “three-dimensional matrix adhesion” originally proposed by Cukierman et al. [5] as a possible indication or “diagnosis” or marker for a culture state of three-dimensionality, the fields of tissue engineering and/or cell-based biosensors have not provided knowledge on the basis of which a consensus for three-dimensionality and the associated complex physiological relevance could be established. Because of this, claims of “physiologically more relevant” are readily made for cells cultured on any surface or scaffold that provides loosely defined 3D geometry, either at the nano- or micro- structure levels or their combinations, as long as the resulting cell phenotypes are different between the 2D and 3D geometries. Second, the concept of using combinatorial approaches to fabricate libraries of polymers or other material scaffolds [6], [7] for tissue engineering or cell-based drug discovery call for high throughput assay by which “hit materials” can be quickly identified for further development. Cell-material interaction outcome can potentially guide the development of such assays or biosensors [8]. An interaction with a material which yields cells that emulate in vivo conditions would be most desirable. Three-dimensionality biomarkers would provide the intellectual basis for material discovery platform development. Third, in order to lower the costs associated with 3D platforms and make them more accessible for high throughput screening (HTS) applications, simplification of the platform without giving up the physiologically relevant behavior of the cells is necessary, as discussed in detail by Lai et al. [4].

Taken together, the subfield or field of 3D culture needs ubiquitous validated biomarkers. As a first step, in search for three-dimensionality biomarkers, we initiated a cytokine expression comparative transcriptomic study with neural progenitor (NP) cells grown on 2D flat surfaces, 3D polymeric scaffolds and as neurospheres (NS). NS were used as the in vivo surrogate, since they have been shown to emulate many in vivo functions that have not been possible in 2D cultures [9], [10]. Cytokines are involved in many crucial cell functions like innate and adaptive inflammatory host defenses, cell growth, differentiation, cell death, angiogenesis, and development and repair processes [11]. Based on the structural homologies of their receptors they can be broadly classified into families like Colony Stimulating Factors, Interleukins, Interferons, TGF (transforming growth factor) family, TNF (tumor necrosis factor) superfamily, PDGF (platelet-derived growth factor) family and Chemokines [11]. Although cytokines have been extensively studied in the field of immunology and oncology, tissue or cell-based biosensor engineers have paid little attention to these small proteins that have potential to revolutionize the field. The evidence for their existence in 3D cultures is compelling but they have not yet been looked at as candidates for potential 3D biomarkers. However, they were an ideal family to explore for the search. The rationale behind their choice was based on the fact that, in a 3D microenvironment cells are surrounded by homotypic neighbors forming a loosely bound disorganized aggregate. When compared to in vivo, such a scenario exists only during avascular tumorogenesis or early stages of inflammatory wound healing and both these phenomenon are regulated by the same molecules – Cytokines [12]. So in vitro, the cells growing in 3D relate to any of those two models depending upon their type – malignant or primary, respectively, and therefore upregulation of their cytokine levels was physiologically relevant.

Materials and Methods

2.1 Materials and reagents

Neural progenitor cells were obtained from Regenerative Bioscience Center at University of Georgia as well as Millipore (ENStem-A™, Billerica, MA). Polystyrene, chloroform, ammonium bicarbonate were obtained from Sigma (St. Louis, MO). Neural basal media, penicillin/streptomycin, L-glutamine, recombinant human leukemia inhibitory factor (hLIF), basic fibroblast growth factor (bFGF), B-27 supplement and phosphate buffered saline (PBS) were obtained from Gibco (Gaithersburg, MD). RNeasy mini kit was obtained from Qiagen (Valencia, CA). Human Genome U133 plus 2.0 microarray chips were obtained from Affymetrix (Santa Clara, CA).

2.2 Scaffold fabrication and neural progenitor cell culture

The scaffold fabrication process in Cheng et al. [13] was followed. Briefly, a viscous polymer solution was prepared by dissolving polystyrene (PS) in chloroform. Sieved ammonium bicarbonate particles in the range of 40–60 µm were added to the polymer solution and mixed thoroughly. These particles generated larger scaffold pore sizes from 60–100 µm, probably due to particle agglomeration. The paste mixture of polymer/salt/solvent was cast into the wells of a standard glass cell culture vessel. The optimal polymer/salt/solvent ratio of 1∶20∶5 (w/w/v) as determined in Cheng et al. [13] was used to generate the paste mixture. Casting was followed by chloroform evaporation, vacuum drying, and baking at 80°C overnight. Before use for cell culture, the scaffolds were sterilized by immersing in 70% alcohol overnight. PS scaffolds offer advantages of low cost, transparency for optical detection, and/or compatibility with existing instrumentation platforms in High Throughput Screening (HTS) applications.

Human neural stem cells or neural progenitors (NP) were maintained in neural basal media (Invitrogen, PA) supplemented with penicillin/streptomycin, L-glutamine, recombinant human leukemia inhibitory factor (hLIF), bFGF and B-27, at 37°C in a 5% CO2 humidified incubator. Medium was half changed every 48 hours. For differentiation, cells were allowed to grow until 90% confluent, after which they were exposed to differentiation medium which was half changed every 24 hours. The composition of differentiation medium was similar to the subculture medium described above with the exception of bFGF. For subculturing, 90% confluent cells were aspirated by pipetting and the subculture ratio was typically 1∶2 to 1∶3. Before cell seeding, both the 3D scaffolds and 2D substrates were coated with poly-ornithine and laminin to rule out any differences that may be caused by the polymer material itself. Neurospheres (NS) were formed by plating cells into non-coated dishes with shaking. Dishes without coating and shaking prevented NP cells from attaching to the surface and encourage them to attach to each other and form spheroid structures. The cell seeding density was 50,000 cells/cm2 for both 2D and 3D cultures.

2.3 Microarray gene expression analysis

Total RNA was isolated from all samples using Qiagen RNeasy Kits (Qiagen, Valencia, CA) according to the manufacturer's standard protocol. The quantity of mRNA isolated from each sample was determined using absorption at 260 and 280 nm. The purity of each sample was monitored using the A260/A280 ratio as well as housekeeping genes. A ratio of 1.8–2.1 was considered a “clean” sample and could be used in microarray experiments. Samples were sent to Affymetrix Core facility at Medical College of Georgia (MCG) for Human Whole Genome U133 plus 2.0 GeneChip Expression Analysis (Affymetrix, Santa Clara, CA). The expression data is publicly available on the GEO site as Series GSE13715.

The expression value of each gene was obtained by Expression Console (Affymetrix) with the RMA (Robust Multiple-array Average) algorithm. RMA is an algorithm widely applied to create expression values for Affymetrix data. The raw intensity values from Affymetrix genechips are background corrected, log2 transformed and then quartile normalized before a linear model was fit to obtain an expression measure for each probe set.

2.4 Statistical analysis

Linear discriminant analysis (LDA) with cross-validation was fulfilled by SAS DISCRIM procedure. Student t-test was used to compare the mean expression from two conditions (2D vs. 3D and 2D vs. NS) with significance level of 0.05. To determine whether cytokines as a group were significantly regulated by culture condition (i.e., 2D, 3D and NS), the permutation t-test method [14] was utilized.

Results

The development of stem cell-based, HTS compatible, in vivo tissue-emulating biosensor platform requires a substrate or scaffold that is not only easy to fabricate but also provides stem cells with a 3D micro environment. For this purpose, 3D porous scaffolds were prepared from polystyrene (PS), a polymer material that has been widely applied in fabricating cell culture vessels. The physical and chemical properties of the scaffolds used in this study have been described by Cheng et al. [13]. Mouse superior cervical ganglion (SCG) and NP cells cultured in these scaffolds, and NP cells cultured as neurospheres (spheroids) have exhibited similar responses as freshly dissected SCG tissue in terms of voltage gated calcium channel and resting membrane potential, while 2D cultured cells exhibited significantly higher responses [13]. Based on Tuj staining (neuronal marker) the neural progenitor cells in scaffolds and neurospheres differentiated into nerve cells [15]. Additionally, spheroids have successfully emulated in vivo drug resistance [9], providing rationale for use of the neurospheres as an in vivo surrogate.

3.1 Quality of transcript expressions

Five experimental groups, with four biological replicates each, were chosen for this study. The five experimental groups were: 1) NP cells before differentiation induction cultured on 2D culture vessel (B2D), 2) NP cells before differentiation in 3D scaffolds (B3D), 3) one week culture of NP cells after differentiation induction on 2D surfaces (A2D), 4) one week culture of NP cells after differentiation in 3D scaffolds (A3D), and 5) neurospheres (NS). In order to minimize the cell passage and other environmental effects on cells, the samples were generated in 2D–3D paired manner, which means at least one 2D and one 3D sample were generated at the same time and from the same batch of cells. Intensity results of the transcripts were examined and normalized to exclude background signals and reported in arbitrary units using RMA (Robust Multiple-array Average) algorithm. RMA is an algorithm widely applied to create expression values for Affymetrix genechip data. The raw intensity values from Affymetrix genechips were first background corrected, log2 transformed and then quartile normalized. The final expression measure was obtained by fitting a linear model for each probe set on each array.

Before we analyze the expression data, quality control tests were performed. No array results were accepted if the correlation coefficient (R2 value) between replicates was less than 0.95. Overall, samples exhibited good correlations with each other. Differentiation status and culture vessel difference didn't affect the overall gene expression pattern (Figure 1, top). We also examined the overall pattern of gene expression based on signal intensity box-plot (Figure 1, bottom). As is the case in box-plot data presentations, the lowest bar represents the observed sample minimum, the base of the box represents the lower quartile while the top represents the upper quartile, the bar inside the box represents the median, and the top bar represents the highest sample observation. The box-plot displays differences between populations without making any assumptions of the underlying statistical distribution. With the exception of a few samples, most of the five values represented in all the data sets (Figure 1, bottom), are in agreement. The third sample from B3D group (B3D_3) exhibited a very distinctive distribution of the expression value in which it had the lowest sample maximum and the smallest “distance” between upper quartile and lower quartile and had a very low correlation with other samples grown under similar conditions, as shown in Figure 1 (top). In order for array data to be accepted for analysis, the overall gene expression pattern has to be similar with no wide divergence between samples. Large variations indicate either hybridization errors or problems with the quality of the RNA used. For this reason B3D_3 was excluded from further analysis. Also, the third sample from A3D group had overall low expression values and a noticeable shift in intensity histogram (data not shown). It had a low correlation coefficient with sample grown under similar condition (Figure 1, top) and was also excluded from further analysis.

thumbnail
Figure 1. Microarray data quality.

The quality of the microarray data was accessed by Pearson's correlation of the samples and the overall distribution of the mRNA expression. The upper panel shows the Pearson's correlation for each pair of the samples. B2D indicate the samples from “Before differentiation, 2D culture condition”, A2D “After differentiation, 2D culture”, B3D “Before differentiation, 3D culture”, A3D “After differentiation, 3D culture”, and NS “Neurospheres”. Pearson's correlation coefficient is a widely used similarity measure for gene expression data. It measures the similarity between two profiles by calculating the linear relationship of the distributions of the two corresponding random variables. The Pearson's correlation value is indicated the color bar color scale (values close to 1 indicate good correlation). The lower panel shows the box-plot of each sample. The box shows the range of the middle 50% with a line in the center for the median value. Additional lines indicate the overall range of the data.

https://doi.org/10.1371/journal.pone.0026821.g001

3.2 Cytokines' transcript and culture classification

The focus on cytokines in this study was based on previous proteomic and transcriptomic evidences from other 3D culture studies indicating that cytokines' expression levels are altered by culture conditions [16][19]. The list of cytokines analyzed was obtained from the Immunology Database and Analysis Portal (ImmPort) (https://www.immport.org/immportWeb/queryref/geneListSummary.do). There were a total of 681 probesets found in our microarray corresponding to the list of human cytokines. The permutation method was used to determine whether cytokines as a group was significantly regulated by culture conditions (i.e. 2D and 3D). Although we reduced the number of probesets to be analyzed from the genechip total of 54000 to 681, by focusing only on cytokines, the number of variables (transcript expression level in our case) still greatly exceeded the number of samples. Permutation testing is superior to standard tools for estimating statistical significance for multiple hypotheses testing due to its ability to empirically determine whether the observed class distinction could be obtained by chance. Briefly, the assay samples were re-iteratively reassigned randomly to the classification label (2D, 3D or NS). This process was done by randomly sampling function from R console's base package. After each reiteration, a list of p values were calculated measuring the significance of newly designated culture condition label for each probesets (2D vs. 3D or 2D vs. NS). After 1000 permutations, summary statistics (we used the median p value) with the permuted class distinction were generated and compared to those obtained in the experimental data [14], [20]. Results showed that cytokines as a group was significant when comparing 2D with 3D or NS (at a significance level of 0.05).

Next, we used linear discriminant analysis (LDA) to test if cytokine transcriptomic expressions can form a good criterion to separate 2D samples from 3D samples. LDA is a multivariate statistical technique commonly used to build a linear prescriptive or descriptive model which can characterize or classify each observation into two or more groups. LDA is also commonly used as a dimension reduction method for further classification analysis. Compared to cluster analysis such as Principle Component Analysis (PCA), LDA requires prior knowledge of the classes, which fitted our purpose of study.

The LDA results are presented in Table 1. As shown cytokines successfully classified 2D from 3D samples with acceptable low error counts. After-differentiation had a lower error rate in comparison to before-differentiation condition. This indicated that if 3D biomarker cytokines do exist, they may be culture time dependent. As pointed out by Lai et al. [4], under such circumstances biomarking with a profile as opposed to an end- point measurement would be the most practical and would increase the robustness of the biomarker. NS exhibited 0 error classification count, suggesting that it was more different from all other culture conditions.

thumbnail
Table 1. Classification of samples by linear discrimination analysis after cross-validation with 681cutokine probesets.

https://doi.org/10.1371/journal.pone.0026821.t001

In Table 2 we list 16 cytokines selected by stepwise discriminant analysis with a cut off p-value of 0.01. With these 16 cytokines, the classification of 2D, 3D and NS is achieved with 0 error rate after cross-validation. Ten of these cytokines are commonly up-regulated in both 3D and NS culture conditions (however not statistically significant). Two are down-regulated both in 3D and NS. Although statistically these cytokines can discriminate 2D from 3D effectively, they are not as a group necessarily the best candidates for 3D biomarkers. This is because, the multivariate nature of the data makes the transcript expression levels correlated with each other, which means each of the variables (cytokine transcript expression levels in our case) can be represented by a group of other variables. Other cytokines may have the same or similar statistical power of these 16 cytokines to discriminate 2D from 3D, with better biological meanings. For this reason, cytokines significantly up-regulated in 3D and NS culture conditions were examined in detail.

thumbnail
Table 2. Most influential genes in LDA Classification by stepwise selection (p<0.01).

https://doi.org/10.1371/journal.pone.0026821.t002

3.3 Cytokines up-regulated in 3D and NS

In Tables S1 and S2, we present cytokine significantly up-regulated (p≤0.05) in 3D and NS cultures. Forty and ninety-one probesets were up-regulated in 3D and NS culture conditions, respectively. The higher number in NS cultures is possibly due to inability to control the size of neurospheres. Many NS were observed to be larger than the 3D structure pore sizes (e.g., three times the maximum pore size of 100 µm) [13]. With large NS, the core of the cellular aggregate may experience hypoxia to the extent that genes not observed in 3D are up-regulated. Evidence in support of this comes from MIP-2 gene (inflammatory protein-2) induced by hypoxia [21] that was also found to be up-regulated in NS but not in 3D conditions in our study. Table 3 lists 13 genes (ANGTL7/CDT6, ARMET/MANF, BMP8B/OP2, CCL13/MCP-4, FGF5, GHRL, IL-11, IL-1B/IL-1F2, NOV/IBP-9, PDGFB, STC1, TGFA, and VEGF-A), whose transcripts were up-regulated in both 3D and NS culture conditions, and examples of their functions. We particularly focused on these genes in the following physiologic function discussion because they were less likely to be up-regulated because of conditions like hypoxia that may be present in NS but absent in 3D conditions.

thumbnail
Table 3. Cytokines up-regulated in both 3D and neurospheres with role examples in cells of nerve tissue origin.

https://doi.org/10.1371/journal.pone.0026821.t003

Discussion

Consistent with our results, up-regulation of cytokines in 3D cultures compared to 2D has been reported by several transcriptomic studies using cells from the four main tissue types (nerve, muscle, connective, and epithelial) cultured in a wide variety of platforms. For example, Klapperich and Bertozzi [19] showed that seven cytokines (IL-8, CXCL1, CXCL2, CXCL3, CXCL5, VEGF, LIF) were up-regulated in human fetal lung fibroblasts (IMR-90) cultured in a collagen–glycosaminoglycan (collagen/GAG) 3D mesh when compared to 2D surfaces. Also, up-regulation of six cytokines (CXCL1- 3, IL-8, MIP-3a, Angiopoetin like4) by a melanoma cell line (NA8) cultured on poly-2-hydroxyethyl methacrylate (polyHEMA) plates was reported by Ghosh at al. [18]. Transcriptomic findings such as those in the above examples have been further substantiated by studies at the protein level. For example, Enzerink et al. [16] have reported induction of chemokine (CCL2-5, CXCL1-3, CXCL8) secretion due to clustering of cells in five different fibroblast cell lines cultured in agarose. Also, Fischbach et al. [17] cultured tumor cells (oral squamous cell carsinoma) in a 2D and 3D RGD-alginate system and reported a dramatic enhancement of IL-8 levels in 3D. Another study by the same group showed that when the same cells were grown in Matrigel (lrBM) there was up-regulation of cytokines when compared to 2D. This observation is of particular importance as cells grown on Matrigel have already been shown to produce an outcome similar to in vivo, like the formation of mammary gland acinus and milk-like secretion into lumen [22] proving that Matrigel can provide all the relevant microenvironmental factors. This suggests that the up-regulation of cytokines in 3D compared to 2D is not a random differential response but is pertinent as a similar response is elicited when a proven physiologically relevant microenvironmental platform is provided.

We have used the themes of tumorogenesis, inflammation and development as shown in Table 3 to closely examine the13 up-regulated transcripts in both 3D and NS culture conditions; cells in a 3D culture in vitro relate to in vivo phenomenon like avascular tumor progression, early stages of inflammatory wound healing or embryonic development depending upon their type- malignant, primary or stem [4]. As these conditions are regulated by autocrine and paracrine cytokine signaling in vivo, the up-regulation of cytokines in 3D culture is physiologically relevant.

Angiogenesis is recurrent function among the genes listed in Table 3, which is not surprising. As already mentioned, as a microtissue grows beyond a certain size, nutrient and oxygen depletion become limiting factors leading to the inhibition of cell proliferation and initiation of angiogenic signaling. Oxygen concentration in 3D tissues depends on the balance between oxygen delivery and consumption. In vivo, this balance is tightly regulated by evenly distributed capillary networks but in vitro homotypic 3D microtissues lack vasculature and therefore develop a hypoxic core as their size increases. This event leads to the cells producing chemical signals (cytokines) for angiogenesis and is quite similar to the response occurring in normal hypoxic tissues where balanced signaling cascades lead to vascular remodeling and angioadaptation until the tissue oxygen concentration is back within its normal range [23]. The up-regulation of VEGF-A and other genes with angiogenic function like ANGPTL7/CDT6, FGF5, NOV/IBP-9, and TGFA is pointing to a functional class of cytokines with great potential as three-dimensionality biomarkers.

The above genes may play other roles besides being factors involved with hypoxia induced angiogenesis. For example, VEGF-A has been shown to regulate neuronal development, survival, neurite growth and it also possesses gliotrophic properties. Also, it's up-regulation is not limited to neuronal cells; it has been shown to be up-regulated in 3D cultures of a variety of cell and tissue types including human fetal lung fibroblasts [19], oral squamous cell carcinoma, glioblastoma, breast cancer [17], neonatal rat cardiomyocytes and neonatal mouse cardiomyocytes [24]. This up-regulation in a non-cell type specific manner in 3D cultures of VEGF-A, and probably other members of the class, lends credibility to the notion of ubiquity of these potential biomarkers in different cells derived from different tissue types.

Another recurrent gene function theme in Table 3 is terminal differentiation, a process by which cells commit to being part of a particular tissue or organ and perform a particular function. Terminal differentiation is preceded by inhibition of proliferation and cell cycle arrest. This exerts endoplasmic reticulum stress on the cell, but anti-apoptotic factors like Bcl-2 protect the cell from apoptosis and commit it to differentiation. In 3D, much like in organs in vivo, there is a spatial constraint on the microtissue, exerted by the defined pore size of the scaffold, which prevents them from proliferating freely, maintaining them in a quiescent state. Evidence in support of this comes from the fact that cells grown in 3D have shown lower proliferation rates than their 2D counterparts and higher expression of cyclin dependent kinase inhibitors (CDI) p21 [25]. Also integrin mediated adhesion to the ECM leads to activation of Bcl-2 family of genes and makes the 3D cells more resistant to apoptosis [26]. Such conditions are favorable for the cell to undergo terminal differentiation and this can be clearly seen by up-regulation in a number of differentiation and survival factor in 3D and NS compared to 2D, like PDGFB and STC1. Neural progenitor cells cultured in 3D have been shown to differentiate and produce a heterogeneous population consisting of neurons, astrocytes and oligodendrocytes [27]. Such cultures essentially behave as a coculture where cells are in symbiotic relationship with each other and might produce factors that can modulate the functions of the other cell types. This has been confirmed by intermediate filament glial fibrillary acidic protein (GFAP), a marker for astrocytic differentiation that was up-regulated in 3D and NS compared to 2D cultures (data not shown). It is a well known fact that astrocytes play a trophic role in supporting neurons. Factors secreted by astrocytes generally belong to the FGF, TGF and EGF families that play an important role in early neurogenesis [28]. Members of these families act as potent mitogens for multipotential neural progenitors and have been implicated in the regulation of several aspects of neurogenesis. Therefore up-regulation of members of these cytokine families in 3D and NS in this study is not surprising. The identification and validation of a few cytokines as three-dimensionality biomarkers that are ubiquitous among cells from different tissue types needs to be done.

Overall, cytokines gene expression results in this study support the notion that 3D cultured cells in various formats are different from their 2D counterparts. Furthermore, up-regulated cytokines' transcripts, independent of culture format, have been identified; this group of 13 cytokines commonly up-regulated in cells cultured in polystyrene scaffolds and neurospheres are suggesting potential for any or a combination from this list to serve as three-dimensionality biomarkers. These results are supportive of further cytokine identification and in vitro/in vivo validation studies with cells from non-neural tissue.

Supporting Information

Table S1.

Significantly up-regulated genes in 3D.

https://doi.org/10.1371/journal.pone.0026821.s001

(DOCX)

Table S2.

Significantly up-regulated genes in neurospheres.

https://doi.org/10.1371/journal.pone.0026821.s002

(DOCX)

Acknowledgments

The authors thank Drs. Lina Wang, Andrew Sornborger and Paul Schliekelman for technical assistance.

Author Contributions

Conceived and designed the experiments: YL KC WSK. Performed the experiments: YL. Analyzed the data: YL AA WSK. Wrote the paper: YL AA WSK.

References

  1. 1. Kisaalita WS (2010) 3D cell-based biosensors in drug discovery programs: Microtissue engineering for high throughput screening. Boca Raton: CRC Press. 77 p.
  2. 2. Green JA, Yamada KM (2007) Three-dimensional microenvironments moderate fibroblast signaling responses. Adv Drug Deliv Rev 59: 1293–1298.
  3. 3. Griffith LG, Swartz MA (2006) Capturing complex 2D tissue physiology in vitro. Nat Rev Mol Biol 7: 211–224.
  4. 4. Lai Y, Asthana A, Kisaalita WS (2011) Biomarkers for simplifying HTS 3D cell culture platforms for drug discovery: the case for cytokines. Drug Discov Today 16(7/8): 293–297.
  5. 5. Cukierman E, Pankov R, Stevens DR, Yamada KM (2001) Taking cell-matrix adhesions to the third dimension. Science 294(5547): 1708–1712.
  6. 6. Simon CG Jr, Stephens JS, Dorsey SM, Becker ML (2007) Fabrication of combinatorial polymer scaffold libraries. Rev Sci Instrum 78(7): 2755–2761.
  7. 7. Yang Y, Bolikal D, Becker ML, Kohn J, Zeiger DN, et al. (2008) Combinatorial polymer scaffold libraries for screening cell-biomaterial interactions in 3D. Advanced Materials 20(11): 2037–2043.
  8. 8. Yang Y, Dorsey SM, Becker ML, Lin-Gibson S, Schumacher GE, et al. (2008) X-ray imaging optimization of 3D tissue engineering scaffolds via combinatorial fabrication methods. Biomaterials 29(12): 1901–1911.
  9. 9. Yamada KM, Cukierman E (2007) Modeling tissue morphogenesis and cancer in 3D. Cell 130(4): 601–610.
  10. 10. Kato-Negishi M, Tsuda Y, Onoe H, Takeuchi S (2010) A neurospheroid network-stamping method for neural transplantation to the brain. Biometarials 31: 8939–8945.
  11. 11. Oppenheim JJ (2001) Cytokines: past, present, and future. Int J Hematol 74(1): 3–8.
  12. 12. Coussens LM, Werb Z (2002) Inflammation and cancer. Nature 420(6917): 860–867.
  13. 13. Cheng K, Lai Y, Kisaalita WS (2008) Three-dimensional polymer scaffolds for high throughput cell-based assay systems. Biomaterials 29(18): 2802–2812.
  14. 14. Good PI (2005) Permutation, parametric and bootstrap tests of hypotheses, 3rd ed. New York: Springer. 128 p.
  15. 15. Wu ZZ, Kisaalita WS, Wang L, Zachman AL, Zhao Y, et al. (2010) Effects of topography on the functional development of human neural progenitor cells. Biotechnol Bioeng 106(4): 649–659.
  16. 16. Enzerink A, Salmenpera P, Kankuri E, Vaheri A (2009) Clustering of fibroblasts induces proinflammatory chemokine secretion promoting leukocyte migration. Mol Immunol 46(8–9): 1787–1795.
  17. 17. Fischbach C, Kong HJ, Hsiong SX, Evangelista MB, Yuen W, et al. (2009) Cancer cell angiogenic capability is regulated by 3D culture and integrin engagement. Proc Natl Acad Sci U S A 106(2): 399–404.
  18. 18. Ghosh S, Spagnoli GC, Martin I, Ploegert S, Demougin P, et al. (2005) Three-dimensional culture of melanoma cells profoundly affects gene expression profile: a high density oligonucleotide array study. J Cell Physiol 204(2): 522–531.
  19. 19. Klapperich CM, Bertozzi CR (2004) Global gene expression of cells attached to a tissue engineering scaffold. Biomaterials 25(25): 5631–5641.
  20. 20. Cheok MH, Yang W, Pui CH, Downing JR, Cheng C, et al. (2003) Treatment-specific changes in gene expression discriminate in vivo drug response in human leukemia cells. Nat Genet 34(1): 85–90.
  21. 21. Zampetaki A, Mitsialis SA, Pfeilschifter J, Kourembanas S (2004) Hypoxia induces macrophage inflammatory protein-2 (MIP-2) gene expression in murine macrophages via NF-kappaB: the prominent role of p42/p44 and PI3 kinase pathways. FASEB J 18(10): 1090–1092.
  22. 22. Weaver VM, Petersen OW, Wang F, Larabell CA, Briand P, et al. (1997) Reversion of the malignant phenotype of human breast cells in three-dimensional culture and in vivo by integrin blocking antibodies. J Cell Biol 137(1): 231–245.
  23. 23. Pries AR, Reglin B, Secomb TW (2001) Structural adaptation of microvascular networks: functional roles of adaptive responses. Am J Physiol Heart Circ Physiol 281(3): H1015–H1025.
  24. 24. Kelm JM, Ehler E, Nielsen LK, Schlatter S, Perriard JC, et al. (2004) Design of artificial myocardial microtissues. Tissue Eng 10(1–2): 201–214.
  25. 25. Birgersdotter A, Sandberg R, Ernberg I (2005) Gene expression perturbation in vitro–a growing case for three-dimensional (3D) culture systems. Semin Cancer Biol 15(5): 405–412.
  26. 26. Frisch SM, Ruoslahti E (1997) Integrins and anoikis. Curr Opin Cell Biol 9(5): 701–706.
  27. 27. Watanabe K, Nakamura M, Okano H, Toyama Y (2007) Establishment of three-dimensional culture of neural stem/progenitor cells in collagen Type-1 Gel. Restor Neurol Neurosci 25(2): 109–117.
  28. 28. Kuhn HG, Winkler J, Kempermann G, Thal LJ, Gage FH (1997) Epidermal growth factor and fibroblast growth factor-2 have different effects on neural progenitors in the adult rat brain. J Neurosci 17(15): 5820–5829.
  29. 29. Peek R, Kammerer RA, Frank S, Otte-Holler I, Westphal JR (2002) The angiopoietin-like factor cornea-derived transcript 6 is a putative morphogen for human cornea. J Biol Chem 277(1): 686–693.
  30. 30. Apostolou A, Shen Y, Liang Y, Luo J, Fang S (2008) Armet, a UPR-upregulated protein, inhibits cell proliferation and ER stress-induced cell death. Exp Cell Res 314(13): 2454–2467.
  31. 31. Yu YQ, Liu LC, Wang FC, Liang Y, Cha DQ, et al. (2010) Induction profile of MANF/ARMET by cerebral ischemia and its implication for neuron protection. J Cereb Blood Flow Metab 30(1): 79–91.
  32. 32. Lindholm P, Peranen J, Andressoo JO, Kalkkinen N, Kokaia Z, et al. (2008) MANF is widely expressed in mammalian tissues and differently regulated after ischemic and epileptic insults in rodent brain. Mol Cell Neurosci 39(3): 356–371.
  33. 33. Allerstorfer S, Sonvilla G, Fischer H, Spiegl-Kreinecker S, Gauglhofer C, et al. (2008) FGF5 as an oncogenic factor in human glioblastoma multiforme: autocrine and paracrine activities. Oncogene 27(30): 4180–4190.
  34. 34. Lindholm D, Harikka J, da Penha Berzaghi M, Castren E, Tzimagiorgis G, et al. (1994) Fibroblast growth factor-5 promotes differentiation of cultured rat septal cholinergic and raphe serotonergic neurons: comparison with the effects of neurotrophins. Eur J Neurosci 6: 244–252.
  35. 35. Reuss B, Hertel M, Werner S, Unsicker K (2000) Fibroblast growth factors-5 and -9 distinctly regulate expression and function of the gap junction protein connexin43 in cultured astroglial cells from different brain regions. Glia 30(3): 231–241.
  36. 36. Volante M, Rosas R, Ceppi P, Rapa I, Cassoni P, et al. (2009) Obestatin in human neuroendocrine tissues and tumours: expression and effect on tumour growth. J Pathol 218(4): 58–466.
  37. 37. Ito Y, Yamamoto M, Li M, Mitsuma N, Tanaka F, et al. (2000) Temporal expression of mRNAs for neuropoietic cytokines, interleukin-11 (IL-11), oncostatin M (OSM), cardiotrophin-1 (CT-1) and their receptors (IL-11Ralpha and OSMRbeta) in peripheral nerve injury. Neurochem Res 25(8): 1113–1118.
  38. 38. Yanagisawa M, Nakashima K, Arakawa H, Ikenaka K, Yoshida K, et al. (2000) Astrocyte differentiation of fetal neuroepithelial cells by interleukin-11 via activation of a common cytokine signal transducer, gp130, and a transcription factor, STAT3. J Neurochem 74(4): 1498–1504.
  39. 39. Mehler MF, Rozental R, Dougherty M, Spray DC, Kessler JA (1993) Cytokine regulation of neuronal differentiation of hippocampal progenitor cells. Nature 362(6415): 62–65.
  40. 40. Nagashima G, Suzuki R, Asai J, Fujimoto T (2002) Immunohistochemical analysis of reactive astrocytes around glioblastoma: an immunohistochemical study of postmortem glioblastoma cases. Clin Neurol Neurosurg 104(2): 125–131.
  41. 41. Sharief MK (1998) Cytokines in multiple sclerosis: pro-inflammation or pro-remyelination? Mult Scler 4(3): 169–173.
  42. 42. Merrill JE, Benveniste EN (1996) Cytokines in inflammatory brain lesions: helpful and harmful. Trends Neurosci 19(8): 331–338.
  43. 43. Planque N, Perbal B (2003) A structural approach to the role of CCN (CYR61/CTGF/NOV) proteins in tumourigenesis. Cancer Cell Int 3(1): 15.
  44. 44. Kubota S, Takigawa M (2007) CCN family proteins and angiogenesis: from embryo to adulthood. Angiogenesis 10(1): 1–11.
  45. 45. Lin CG, Chen CC, Leu SJ, Grzeszkiewicz TM, Lau LF (2005) Integrin-dependent functions of the angiogenic inducer NOV (CCN3): implication in wound healing. J Biol Chem 280(9): 8229–8237.
  46. 46. Perbal B (2001) NOV (nephroblastoma overexpressed) and the CCN family of genes: structural and functional issues. Mol Pathol 54(2): 57–79.
  47. 47. Chevalier G, Yeger H, Martinerie C, Laurent M, Alami J, et al. (1998) novH: differential expression in developing kidney and Wilm's tumors. Am J Pathol 152(6): 1563–1575.
  48. 48. Calzolari F, Appolloni I, Tutucci E, Caviglia S, Terrile M, et al. (2008) Tumor progression and oncogene addiction in a PDGF-B-induced model of gliomagenesis. Neoplasia 10(12): 1373–1382.
  49. 49. Sjoborg M, Pietz K, Ahgren A, Yamada N, Lindvall O, et al. (1998) Expression of platelet-derived growth factor after intrastriatal ibotenic acid injury. Exp Brain Res 119(2): 245–250.
  50. 50. Williams BP, Park JK, Alberta JA, Muhlebach SG, Hwang GY, et al. (1997) A PDGF-regulated immediate early gene response initiates neuronal differentiation in ventricular zone progenitor cells. Neuron 18(4): 553–562.
  51. 51. Nakao N, Brundin P, Funa K, Lindvall O, Odin P (1994) Platelet-derived growth factor exerts trophic effects on rat striatal DARPP-32-containing neurons in culture. Exp Brain Res 101(2): 291–296.
  52. 52. Koide Y, Sasaki T (2006) Stanniocalcin-1 (STC-1) as a molecular marker for human cancer. Rinsho Byori 54(3): 213–220.
  53. 53. Zhang K, Lindsberg PJ, Tatlisumak T, Kaste M, Olsen HS, et al. (2000) Stanniocalcin: A molecular guard of neurons during cerebral ischemia. Proc Natl Acad Sci U S A 97(7): 3637–3642.
  54. 54. Serlachius M, Zhang KZ, Andersson LC (2004) Stanniocalcin in terminally differentiated mammalian cells. Peptides 25(10): 1657–1662.
  55. 55. Dunn IF, Heese O, Black PM (2000) Growth factors in glioma angiogenesis: FGFs, PDGF, EGF, and TGFs. J Neurooncol 50(1–2): 121–137.
  56. 56. Junier MP (2000) What role(s) for TGFalpha in the central nervous system? Prog Neurobiol 62(5): 443–473.
  57. 57. Plate KH, Breier G, Weich HA, Mennel HD, Risau W (1994) Vascular endothelial growth factor and glioma angiogenesis: coordinate induction of VEGF receptors, distribution of VEGF protein and possible in vivo regulatory mechanisms. Int J Cancer 59(4): 520–529.
  58. 58. Zagzag D, Lukyanov Y, Lan L, Ali MA, Esencay M, et al. (2006) Hypoxia-inducible factor 1 and VEGF upregulate CXCR4 in glioblastoma: implications for angiogenesis and glioma cell invasion. Lab Invest 86(12): 1221–1232.
  59. 59. Nag S, Takahashi JL, Kilty DW (1997) Role of vascular endothelial growth factor in blood-brain barrier breakdown and angiogenesis in brain trauma. J Neuropath Exp Neurol 56(8): 912–921.
  60. 60. Krum JM, Rosenstein JM (1998) VEGF mRNA and its receptor flt-1 are expressed in reactive astrocytes following neural grafting and tumor cell implantation in the adult CNS. Exp Neurol 154(1): 57–65.
  61. 61. Rosenstein JM, Krum JM (2004) New roles for VEGF in nervous tissue–beyond blood vessels. Exp Neurol 187(2): 246–253.