Embryonic stem cells are conventionally differentiated by modulating specific growth factors in the cell culture media. Recently the effect of cellular mechanical microenvironment in inducing phenotype specific differentiation has attracted considerable attention. We have shown the possibility of inducing endoderm differentiation by culturing the stem cells on fibrin substrates of specific stiffness [1]. Here, we analyze the regulatory network involved in such mechanically induced endoderm differentiation under two different experimental configurations of 2-dimensional and 3-dimensional culture, respectively. Mouse embryonic stem cells are differentiated on an array of substrates of varying mechanical properties and analyzed for relevant endoderm markers. The experimental data set is further analyzed for identification of co-regulated transcription factors across different substrate conditions using the technique of bi-clustering. Overlapped bi-clusters are identified following an optimization formulation, which is solved using an evolutionary algorithm. While typically such analysis is performed at the mean value of expression data across experimental repeats, the variability of stem cell systems reduces the confidence on such analysis of mean data. Bootstrapping technique is thus integrated with the bi-clustering algorithm to determine sets of robust bi-clusters, which is found to differ significantly from corresponding bi-clusters at the mean data value. Analysis of robust bi-clusters reveals an overall similar network interaction as has been reported for chemically induced endoderm or endodermal organs but with differences in patterning between 2-dimensional and 3-dimensional culture. Such analysis sheds light on the pathway of stem cell differentiation indicating the prospect of the two culture configurations for further maturation.
References
[1]
Jaramillo M, Singh S, Velankar S, Kumta PN, Banerjee I (2011) Inducing endoderm differentiation by modulating mechanical properties of soft substrates. Pittsburgh: University of Pittsburgh.
[2]
Murry CE, Keller G (2008) Differentiation of Embryonic Stem Cells to？Clinically Relevant Populations: Lessons from Embryonic Development. Cell 132: 661–680.
[3]
Engler AJ, Sen S, Sweeney HL, Discher DE (2006) Matrix Elasticity Directs Stem Cell Lineage Specification. Cell 126: 677–689.
[4]
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al. (2005) Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America 102: 15545–15550.
[5]
Sill M, Kaiser S, Benner A, Kopp-Schneider A (2011) Robust biclustering by sparse singular value decomposition incorporating stability selection. Bioinformatics 27: 2089–2097.
[6]
Denton A, Wu J, Townsend M, Sule P, PruSZ B (2008) Relating gene expression data on two-component systems to functional annotations in Escherichia coli. BMC Bioinformatics 9: 294.
[7]
Hartigan JA (1972) Direct Clustering of a Data Matrix. Journal of the American Statistical Association 67: 123–129.
[8]
Cheng Y, Church GM (2000) Biclustering of expression data. 8th International Conference on Intelligent Systems for Molecular Biology. pp. 93–103.
[9]
Getz G, Levine E, Domany E (2000) Coupled two-way clustering analysis of gene microarray data. Proceedings of the National Academy of Sciences 97: 12079–12084.
[10]
Wang H, Wang W, Yang J, Yu P (2002) Clustering by pattern similarity in large data sets. ACM-SIGMOD international conference on management of data. pp. 418–427.
[11]
Abdullah A, Hussain A (2006) A new biclustering technique based on crossing minimization. Neurocomputing 69: 1882–1896.
[12]
Reiss D, Baliga N, Bonneau R (2006) Integrated biclustering of heterogeneous genome-wide datasets for the inference of global regulatory networks. BMC Bioinformatics 7: 280.
[13]
Hochreiter S, Bodenhofer U, Heusel M, Mayr A, Mitterecker A, et al. (2010) FABIA: factor analysis for bicluster acquisition. Bioinformatics 26: 1520–1527.
[14]
Madeira SC, Oliveira AL (2004) Biclustering algorithms for biological data analysis: a survey. IEEE/ACM Trans Comput Biol Bioinform 1: 24–45.
[15]
Orling J (1977) Containment in Graph Theory: Covering Graphs with Cliques Nederl Akad Wetensch Indag Math 39: 8.
[16]
Divina F, Aguilar-Ruiz JS (2006) Biclustering of expression data with evolutionary computation. Knowledge and Data Engineering, IEEE Transactions on 18: 590–602.
[17]
Yang E, Foteinou PT, King KR, Yarmush ML, Androulakis IP (2007) A novel non-overlapping bi-clustering algorithm for network generation using living cell array data. Bioinformatics 23: 2306–2313.
[18]
Willems E, Leyns L, Vandesompele J (2008) Standardization of real-time PCR gene expression data from independent biological replicates. Analytical Biochemistry 379: 127–129.
[19]
Politis DN, Romano JP (1994) The Stationary Bootstrap. Journal of the American Statistical Association 89: 1303–1313.
[20]
Stine R (1989) An Introduction to Bootstrap Methods. Sociological Methods & Research 18: 243–291.
[21]
Efron B, Tibshirani RJ (1993) An Introduction to the Bootstrap: Chapman \& Hall.
[22]
Huebsch N, Arany PR, Mao AS, Shvartsman D, Ali OA, et al. (2010) Harnessing traction-mediated manipulation of the cell/matrix interface to control stem-cell fate. Nat Mater 9: 518–526.
[23]
Cao L, Gibson JD, Miyamoto S, Sail V, Verma R, et al. (2011) Intestinal lineage commitment of embryonic stem cells. Differentiation 81: 1–10.
[24]
Anderson WJ, Zhou Q, Alcalde V, Kaneko OF, Blank LJ, et al. (2008) Genetic targeting of the endoderm with claudin-6CreER. Developmental Dynamics 237: 504–512.
[25]
Costa RH, Kalinichenko VV, Holterman A-XL, Wang X (2003) Transcription factors in liver development, differentiation, and regeneration. Hepatology 38: 1331–1347.
[26]
Herbst R, Nielsch U, Sladek F, Lai E, Babiss LE, et al. (1991) Differential regulation of hepatocyte-enriched transcription factors explains changes in albumin and transthyretin gene expression among hepatoma cells. New Biol 3: 289–296.
[27]
Pani L, Overdier DG, Porcella A, Qian X, Lai E, et al. (1992) Hepatocyte nuclear factor 3 beta contains two transcriptional activation domains, one of which is novel and conserved with the Drosophila fork head protein. Mol Cell Biol 12: 3723–3732.
[28]
Cockell M, Stolarczyk D, Frutiger S, Hughes G, Hagenbuchle O, et al. (1995) Binding sites for hepatocyte nuclear factor 3 beta or 3 gamma and pancreas transcription factor 1 are required for efficient expression of the gene encoding pancreatic alpha-amylase. Mol Cell Biol 15: 1933–1941.
[29]
Martinez Barbera JP, Clements M, Thomas P, Rodriguez T, Meloy D, et al. (2000) The homeobox gene Hex is required in definitive endodermal tissues for normal forebrain, liver and thyroid formation. Development 127: 2433–2445.
[30]
Denson LA, McClure MH, Bogue CW, Karpen SJ, Jacobs HC (2000) HNF3[beta] and GATA-4 transactivate the liver-enriched homeobox gene, Hex. Gene 246: 311–320.
[31]
Lehner F, Kulik U, Klempnauer J, Borlak J (2010) Mapping of liver-enriched transcription factors in the human intestine. World J Gastroenterol 16: 3919–3927.
[32]
Rojas A, Schachterle W, Xu S-M, Martín F, Black BL (2010) Direct transcriptional regulation of Gata4 during early endoderm specification is controlled by FoxA2 binding to an intronic enhancer. Developmental Biology 346: 346–355.
[33]
Lopez-Sanchez C, Garcia-Masa N, Ga？an CM, Garcia-Martinez V (2009) Movement and commitment of primitive streak precardiac cells during cardiogenesis. Int J Dev Biol 53: 1445–1455.
[34]
Nagy P, Bisgaard HC, Thorgeirsson SS (1994) Expression of hepatic transcription factors during liver development and oval cell differentiation. The Journal of Cell Biology 126: 223–233.
[35]
Morrison GM, Oikonomopoulou I, Migueles RP, Soneji S, Livigni A, et al. (2008) Anterior Definitive Endoderm from ESCs Reveals a Role for FGF Signaling. Cell Stem Cell 3: 402–415.
[36]
Banerjee I, Ierapetritou MG (2003) Development of an adaptive chemistry model considering micromixing effects. Chemical Engineering Science 58: 4537–4555.
[37]
Banerjee I, Maiti S, Parashurama N, Yarmush M (2010) An integer programming formulation to identify the sparse network architecture governing differentiation of embryonic stem cells. Bioinformatics 26: 1332–1339.
[38]
Banerjee I, Ierapetritou MG (2006) An Adaptive Reduction Scheme to Model Reactive Flow. Combust & Flame 144:
[39]
Goldberg D (1989) Genetic Algorithms in Search, Optimization, and Machine Learning: Addison-Wesley Professional.
[40]
Bieniasz LK, Rabitz H (2006) Extraction of Parameters and Their Error Distributions from Cyclic Voltammograms Using Bootstrap Resampling Enhanced by Solution Maps:？Computational Study. Analytical Chemistry 78: 8430–8437.
