- Research article
- Open Access
TGF-β stimulation in human and murine cells reveals commonly affected biological processes and pathways at transcription level
BMC Systems Biology volume 8, Article number: 55 (2014)
The TGF-β signaling pathway is a fundamental pathway in the living cell, which plays a key role in many central cellular processes. The complex and sometimes contradicting mechanisms by which TGF-β yields phenotypic effects are not yet completely understood. In this study we investigated and compared the transcriptional response profile of TGF-β1 stimulation in different cell types. For this purpose, extensive experiments are performed and time-course microarray data are generated in human and mouse parenchymal liver cells, human mesenchymal stromal cells and mouse hematopoietic progenitor cells at different time points. We applied a panel of bioinformatics methods on our data to uncover common patterns in the dynamic gene expression response in respective cells.
Our analysis revealed a quite variable and multifaceted transcriptional response profile of TGF-β1 stimulation, which goes far beyond the well-characterized classical TGF-β1 signaling pathway. Nonetheless, we could identify several commonly affected processes and signaling pathways across cell types and species. In addition our analysis suggested an important role of the transcription factor EGR1, which appeared to have a conserved influence across cell-types and species. Validation via an independent dataset on A549 lung adenocarcinoma cells largely confirmed our findings. Network analysis suggested explanations, how TGF-β1 stimulation could lead to the observed effects.
The analysis of dynamical transcriptional response to TGF-β treatment experiments in different human and murine cell systems revealed commonly affected biological processes and pathways, which could be linked to TGF-β1 via network analysis. This helps to gain insights about TGF-β pathway activities in these cell systems and its conserved interactions between the species and tissue types.
The transforming growth factor-beta1 (TGF-β1) signaling pathway is a fundamental pathway in the living cell, which plays a role in many central cellular processes. The TGF-β superfamily contains over 30 different proteins, such as BMPs, Activins, Inhibins, and the TGF-β1 isoforms [1–3]. The pathway contributes to regulation of various cellular processes, such as apoptosis, cell differentiation, cell growth as well as tumor suppression and immune regulation processes .
There are three TGF-β isoforms (TGF-β1, TGF-β2, TGF-β3) which have different physiological and pathological effects on epithelial, endothelial, lymphatic, myeloid and mesenchymal tissues . The TGF-β pathway is one of the most studied pathways [6–10]. However, the complex and sometimes contradicting mechanisms by which TGF-β yields phenotypic effects is not yet completely understood . The classical TGF-β1 pathway is already well established since several years . However, the identification of alternative signaling pathways that contain different receptors and Smad proteins has increased the overall complexity of the TGF-β1 signaling pathway . Additional file 1: Figure S1 shows a simplified cartoon sketch comprising mainly Smads cascades in the TGF-β1 signaling pathway.
In this study we investigated and compared downstream effects of TGF-β1 stimulation on the dynamical response of gene expression in mouse and human in different cell and tissue types. Two types of mouse hematopoietic progenitor cells were used: multipotent progenitor (MPP) and dendritic cell (DC) committed progenitors, referred to as common dendritic progenitor (CDP) cells. CDP differentiate from MPP and give rise to two types of DC: plasmacytoid DC (pDC) and conventional DC (cDC). MPP and CDP were obtained from bone marrow by in vitro culture with a specific cytokine cocktail and FACS sorting [12, 13]. Furthermore, we employed human mesenchymal stromal cells (MSC), which differentiate into osteocytes, chondrocytes or adipocytes [14–16]. Finally, primary murine hepatocytes (HPC) and immortalized human hepatocytes (human HPC, HepG2) cells were used. We have taken these different cell types for three reasons: (i) All these cells are highly responsive to TGF-β. (ii) The different cell types reflect different degrees of differentiation. (iii) The different cells show a variable response to TGF-β. While in hepatocytes TGF-β induces apoptosis, multipotent progenitors initiate a differentiation programme in response to TGF-β.
Very little and vague information is known about the detailed influence of TGF-β1 in these different cell systems. For example, TGF-β1 is known to be necessary for MSC proliferation. It is essential for chondrogenic differentiation. On the other hand, TGF-β1 participates in inhibition of adipogenic and osteogenic differentiation. Furthermore, there are evidences, that TGF-β1 contributes to supporting myogenic differentiation of MSC [17–19]. There are also evidences that the TGF-β pathway play a role in the induction of cellular senescence in MSC . Although TGF-β1 triggers primary early responses (e.g. Smad activation) and EMT in human HPC (HepG2) cells, cell cycle arrest and apoptosis are generally not promoted by TGF-β1 [21, 22]. Furthermore, TGF-β1 is known to be crucial for development of Langerhans cells, the cutaneous contingent of migratory dendritic cells, both in vivo and in vitro and it evidently contributes in accelerating their differentiation and directing their subsets specification toward cDCs [12, 23–25].
We used a panel of bioinformatics methods, ranging from statistical testing over functional and promoter sequence analysis to clustering for pattern discovery in our gene expression time series data. Only one gene, the SKI-like oncogene (Skil), was commonly found to be differentially expressed (DE) in all cell types. Skil is a component of the SMAD-pathway, which regulates cell growth and differentiation. Moreover, Smad7 that blocks TGF-β receptor activity seems to play a major common role, because it was identified as DE in most cell types. Despite of the differences on the level of individual genes we observed a conserved effect of TGF-β1 stimulation on a number of biological processes and pathways. Moreover, we could identify a few overrepresented transcription factor binding sites, which were commonly found in several cell types. Specifically EGR1 seems to have major relevance for the transcriptional stimulation response in mouse and human.
By analysis of an independent dataset on human A549 lung adenocarcinoma cells (CRL) from GEO (access No. GSE17708)  we were able to reproduce a highly significant proportion of the commonly identified biological processes, pathways and transcriptional factors in our datasets. Network analysis suggests explanations, how TGF-β1 stimulation could lead to the observed effects.
Results and discussion
Time series transcriptome measurements
All cell types were treated with TGF-β in three biological replicates. TGF-β treatment concentrations were optimized in each cell type to show a maximal effect. Extracted RNA samples were hybridized to microarrays (Affymetrix Gene 1.0 ST) for genome-wide transcriptome analysis. Mouse progenitor cells and HepG2 cells were measured at 6 successive time points, mouse primary HPC cells at 5, and human MSCs at 4 different time points. Additional file 2: Table S1 gives an overview of our experiments and the measured time-points, the “Methods” section gives details about cell cultures, stimulation, RNA-isolation and array hybridization in our experiments.
Differential gene expression
Transcriptional response is highly tissue specific on gene level
We employed the “betr” method  to quantify the probability of differential expression of genes in whole time-courses (see Methods). Using this approach we were able to assess differential gene expression for each gene in each cell type in a comparable manner. We considered a gene to have differential time-course expression (DE), if it had a probability of >99% and was at least two-fold up- or down-regulated at one time point minimum (Additional file 1: Figures S2 a & b, Additional file 2: Tables S2 & S8).
The strongest stimulatory effect of TGF-β1 was observed in CDP cells (614 genes). Eight out of these genes in CDP are already known to play a role in the TGF-β pathway (Tgfb3, Smad7, Thbs1, Tgfbr1, Smurf1, Smad3, Smad6, Tgfbr2). In mouse HPC a significantly lower number of DE genes were found compared to other cell types.
We conducted set comparisons of DE genes across cell types. It is worth mentioning in this context that comparisons between mouse and human genes were done on the basis of homologous genes (see Methods). Not surprisingly, the found overlap was particularly high among mouse hematopoietic progenitor cells (MPP and CDP). These were 173 genes, which equals a harmonic mean of above 41% of DE genes in both cell types (Additional file 1: Figure S2 a). Only two of these genes, namely Smad7 and Tgfbr1 are known to play a role in the TGF-β pathway. Three genes (Lox, Pmepa1, Skil) are found to be DE in all mouse cells (CDP, MPP and HPC). Pmepa1 (Prostate Transmembrane Protein) is known to interact with Smad and suppress the TGF-β pathway [28, 29]. Only the protein-coding gene Skil (Ski-like-oncogene) that encodes a protein in the SMAD-pathway [30, 31] was found to have a DE time-course in all cell types. In addition, the gene Smad7 was commonly found in all cell types except mouse HPC cells. 18 genes including ROR1, C10orf10, SMAD7, FSTL3, GADD45B, JUNB, ZFP36, OLFM2, SPTLC3, ID1, LMCD1, SLC38A3, GXYLT2, SKIL, HES1, RASGEF1B, CITED2 and PDGFA were DE in all human cells (MSC, HepG2). The heatmaps in Additional file 1: Figure S3 visualize patterns of temporal behavior for particular groups of genes. Here again we see similarity in gene expressions between mouse dendritic cells.
These findings on one hand stress the similarity of the transcriptional response in MPP and CDP, which is not very surprising given the fact that these cells were both derived from bone marrow. On the other hand they highlight that TGF-β1 treatment affects by far not only genes within the canonical TGF-β1 pathway, but leads to a large number of diverse secondary downstream effects, which are only partially overlapping across different cell types. In other words there is a high tissue specificity of the transcriptional TGF-β stimulation response on the level of individual genes.
TGF-β1 pathway genes react time-dependant and tissue-specific
We had a closer look at genes, which are known to play a role in the TGF-β1 pathway, such as Bmp(s), Smad(s) and Id(s). In Figure 1 the log2 fold changes of 17 genes involved in the TGF-β1 pathway, which are DE in at least one cell type, are depicted. It can be noticed that almost all genes show time-dependant transcriptional effects. These effects are distinct between early and later time points, with moderate activities until 4 h and mostly higher activities at late times. It can also be noticed that cells of similar origin are more alike. For example, Bmp2, Bmp4, Bmp6, Cdkn2b and Comp are dys-regulated (i.e. significantly differ from 0 level according to “betr”) only in human and not in mouse tissues. Fs1 is similar to these genes, but also shows activity in mouse HPC. Id1 in human cells is up-regulated at earlier time points and a down-regulated after 4 h. Inhba shows activity only in MSC cells where its expression after 1 hour constantly increases. Smad3, Smad6 and Smad7 reveal similar time courses in mouse MPP and CDP cells and in human MSCs. Smad3 is increasingly down-regulated over time and the other two genes are always up-regulated. Smurf1 is always over-expressed and shows a curve that is opposite to Smad3, Smad6 and Smad7. Tgfb3 is over-expressed at later time points in MPP and CDP cells and shows almost no activity in the other cell types. Thbs1 is highly active in all cell types. However, while it is underrepresented in MPP and CDP, it shows elevated expression in mouse and human HPC. Tgfbr1 and Tgfbr2 behave similar, in particular in mouse progenitor cells, where Tgfbr2 is less up-regulated than Tgfbr1.
Time point specific analysis confirms highly tissue specific expression changes on gene expression level
In order to cross-validate our previous analysis, which considers time series as a whole, we conducted also a time-point specific analysis of differential gene expression using linear models for microarray data (Limma). For this purpose we compared the gene expression at 4 hours after stimulation to the initial expression at time point 0. The time period of 4 hours was chosen because at least short-time relevant effects are expected in all cell types after this period.
In the context of this time point specific analysis of transcriptional effects we considered a gene to be differentially expressed (DE) in a given cell type, if FDRBH < 0.01 and the absolute fold changes was (logFC) > 1. The overlap analysis of DE genes at 4 h agrees with the time-course analysis. There are no or very few genes in common between the different cell types except in the case of mouse dendritic cells (Additional file 1: Figure S5 A). Moreover, the direction of regulation (up or down) differs between cell types (Additional file 2: Table S9).
The heatmap in Additional file 1: Figure S4 depicts the log fold changes of all genes, which are DE in at least one cell type. The plot indicates two gene sets, which clearly show a similar behavior in mouse MPP and CDP cell types. The first set contains 36 genes that are over-expressed. The other set (42 genes) is under-expressed. Interestingly, the 36 genes being up-regulated in MPP and CDP cells are not regulated by TGF-β1 in other cell types. Although not DE genes in every cell type, the genes Smad7, Pmepa1 (beside the gene Skil) seem to be up-regulated in all the cells. The rest of the genes are regulated in a rather cell-type specific manner.
Cluster analysis reveals functionally similar gene groups in different cell types
We conducted a time series cluster analyses in order to find groups of DE genes showing similar expression changes over time (observed as within cell-type temporal behavior shown in Additional file 1: Figures S3, between cell-types similarities shown in Additional file 1: Figure S4). The cluster analyses yielded 12 clusters in MPP and mouse HPC, 20 in CDP and 11 in human MSC (Additional file 2: Table S3). Genes contained in individual clusters can be found in Additional file 2: Table S14. Figure 2 depicts the mean curves for each of these clusters in each cell type. We investigated the functional similarity of genes across different clusters. For this purpose the R-Package “GOSemSim”  utilizing the semantic similarity measure proposed by Wang et al.  was employed. Semantic similarities are a means to compare GO annotations of gene pairs in a quantitative manner, for example on the basis of the information content of GO terms. We refer the reader to  for an excellent overview.
A heatmap depicting these GO semantic similarities suggested a high functional similarity of genes in several clusters from different cell types (Additional file 1: Figure S6 Additional file 2: Tables S15, S16). In particular cluster B (MPP), and cluster B (CDP) are highly similar to each other (semantic similarity > 0.7). Time-course log2 fold changes of the corresponding genes are shown in (Figure 3 top). As can be noticed the clusters are of different size, but have several genes in common (13 genes). Functional analysis revealed that genes in these clusters are enriched for cell adhesion molecules (CAMs), valine, leucine and isoleucine biosynthesis, Pantothenate and CoA biosynthesis and regulation of cellular extravasation. Enrichment analysis was conducted here via the R-package GOstats , which employs a hypergeometric test taking into account the dependency structure among GO terms.
The second group of functionally similar clusters (Figure 3, bottom) contains cluster K (human HPC) and cluster B (MSCs). Genes in these clusters play (among others) a role in TGF-β and Notch signaling pathways (Additional file 2: Tables S15, S16).
Taken together our cluster analyses showed that despite evident differences on the level of individual genes, functionally similar clusters of genes can be identified across cell types.
Enrichment analysis reveals commonly affected biological processes, pathways and transcription factors in all cell types
Motivated by our previous findings we asked, whether there were common functional patterns across all cell types. For this purpose we scanned GO terms and KEGG pathways for significant association with differential time course gene expression in each cell type (Additional file 2: Tables S4, S10, S11).
Our analysis brought up 6 KEGG pathways and 11 GO terms, which were significantly associated to all cell types (FDR < 5%, Figure 4). The 6 KEGG pathways associated to all cell types were: Metabolic pathways, Glutathione metabolism, Lysosome, Purine metabolism, Peroxisome and PPAR signaling pathway. The 11 GO terms associated to all cells were: oxidation-reduction process, innate immune response, positive regulation of transcription from RNA polymerase II promoter, negative regulation of apoptotic process, angiogenesis, lipid metabolic process, positive regulation of cell proliferation, positive regulation of cell migration, proteolysis, positive regulation of transcription DNA-dependent and response to drug. The role of TGF-β in apoptosis, cell proliferation as well as immune response is well known. Moreover, an effect of TGF-β perturbation on PPAR signaling has been described in skin fibroblasts . In  the authors describe TGF-β mediated oxidative stress and decreased glutathione concentration in fibrosis models. Finally, there is evidence that TGF-β has an effect on angiogensis and cell migration . Hence, our findings largely fit to the current biological knowledge about TGF-β.
Conserved role of EGR1/2 transcription factors
We analyzed DE genes with respect to overrepresented sequence motifs in their promoter regions with the XXmotif tool . Significant motifs were then compared to known position weight matrices (TRANSFAC) of transcription factors (TFs) via STAMP .
The analysis in each cell type predicted between 11 and 21 regulating TFBS in the time-course analysis (Additional file 2: Tables S6, S17, S18), except for mouse HPC, where no overrepresented TFBS could be detected. This may be attributed to the small number of 16 DE genes in this cell type. Overlaps were particularly high within mouse MPP and CDP and within human cells.FOXP1, KROX, TEF, POU6F1, FOX and PITX binding sites were commonly identified in mouse MPP and CDP. KROX, HFH4 and PAX4 were found in all human cells. FOX, FOXP1, KROX and TEF were found to be themselves representatives of DE genes. Figure 5 shows a network representation of all eight TFBS together with the set of DE genes containing respective binding sites. The plot reveals a relative clear difference between mouse and human cells with the exception of the KROX TFBS, which appears in all four cell types. KROX represents EGR1 and EGR2.
Network analysis suggests possible signal transduction pathways in mouse and human
In order to better understand, how TGF-β may influence the commonly identified transcription factor, biological processes and the PPAR-pathway we conducted a network analysis. Using protein-protein interaction information from the BioGRID database  we reconstructed a mouse and a human specific network depicting dys-regulated paths from TGF-β to SKIL, SMAD7, EGR1 as well as genes involved into glutathione metabolism, purine metabolism, PPAR signaling, oxidation-reduction process, innate immune response, negative regulation of apoptotic process, angiogenesis, positive regulation of cell proliferation and positive regulation of cell migration (Figures 6 and 7; see further details in Methods part).
Our network analysis suggests pathways, by which TGF-β stimulation is possibly propagated via protein-protein interactions to our commonly identified biological processes. Due to the organism specificity of interactome information these pathways show certain differences: Far less protein-protein interactions are known in mouse than in human. In human, for example, negative regulation of apoptosis might be mediated via SMAD3 and SOX9 . In contrast, the GO and network analysis in mouse suggests a direct role of TGFBR1.
Enrichment of biological processes, pathways and tfbss is reproducible on independent dataset
To validate the central finding from our data, namely the existence of commonly affected biological processes, pathways and transcription factors in all cell types, we downloaded an independent dataset comprising gene expression data measured at 9 time points (0, 0.5, 1, 2, 4, 8, 16, 24, 72 h) after TGF-β stimulation in human A549 lung adenocarcinoma cell-lines (CRL, GSE17708). The dataset was analyzed in the same manner as described for our data before. High fractions of the 11 GO terms and 6 KEGG pathways commonly identified in all of our cell types were also found in GSE17708 (Figure 8, Additional file 2: Table S10). Out of the KEGG pathways and GO terms associated to all of our human cells 70% and 74%, respectively could be reproduced on the independent dataset (Figure 8, Additional file 2: Table S11). Notably, 11 (61%) out of the 18 genes which exhibiting differential time courses in both our human MSC and HPC cells were found also to have differential time-courses in GSE17708 cells, these were ROR1, SMAD7, FSTL3, GADD45B, JUNB, ZFP36, ID1, LMCD1, GXYLT2, SKIL and HES1. This corresponding fraction is significantly larger than expected by chance (p < 1E-9, hyper-geometric test).
The KROX TFBS (corresponding to TFs EGR1 and EGR2), which was enriched in all of our cell types, was also found in GSE17708. Moreover, the other two TFBS that we identified in our human cells (HFH4, PAX4) were also enriched in the A549 lung cancer cell line (Additional file 2: Tables S17).
Taken together this analysis reveals a high reproducibility of our commonly identified biological processes, pathways as well as TFs.
We have conducted an in-depth comparison of the dynamical TGF-β1 response profile on gene expression level across several cell types. Despite of a generally high degree of cell type specificity, there appears to be a common functional response, which is conserved across cell types and species (i.e. mouse and human). Our analysis suggests a common effect of TGF-β1 stimulation on apoptosis, cell proliferation, immune response, angiogenesis, cell migration, PPAR signaling, oxidative stress as well as purine and glutathione metabolism. Network analysis gives hints to possible pathways, by which these effects could be mediated.
On the level of individual genes the SKI-like oncogene and Smad7 were differentially expressed in most (Smad7) or all (SKI-like oncogene) cell types and thus appear to play a major role. Smad7 is involved into the canonical TGF-β pathway . It is a general antagonist of the TGF-β family (for review see ). The SKI-like oncogene is a direct target gene of Smad2, which regulates its transcription . It plays a role in cell growth and differentiation.
Notably, a high fraction of the biological processes, pathways and TFBS that we identified to be enriched in all our cell types was found also in an independent dataset from a lung cancer cell line. This strengthens the confidence into our results.
In summary our findings indicate that despite a high variability of transcriptional response across cell types and organisms there appears to be a set of commonly affected processes and pathways. In addition, the TFBS analysis suggested a major role of the transcription factor EGR1 in the TGF-β response in human and mouse. Indeed the induction of EGR1 via TGF-β stimulation has been already reported earlier  and thus fits to the existing knowledge about TGF-β1 − induced transcriptional response in other cell systems.
Previous studies of TGF-β stimulation were mainly limited to one specific cell type (e.g. fibroblasts). In this paper we went beyond this point and conducted experiments in different cell types under as much as possible comparable conditions. In consequence we were able to compare transcriptional responses across cell types and organisms, which revealed common patterns. The identification of common and specific signal transduction pathways that are affected by TGF-β in human and mice will allow us to define potential therapeutic targets and will further enable us to characterize gene expression patterns and complex regulatory networks. In addition, future work using our and other transcriptome data can, for example, address the identification of TGF-β dependent mesenchymal or epithelial gene signatures or the definition of cell specific cancer signatures.
Animal experiments required for obtaining murine MPP, CDP and HPC were approved by local authorities (Landesamt für Natur, Umwelt und Verbraucherschutz Nordrhein-Westfalen - LANUV NRW) in compliance with the German animal protection law. Experiments were performed at the Institute for Animal Research of the RWTH Aachen University Hospital, under projects 10534A4 and 10718A4, entitled “Untersuchungen zur Hämatopoese aus adulten Blutstammzellen”. Animal maintenance, handling, and anesthesia were performed according to the Federation for Laboratory Animal Science Associations FELASA recommendations.
All samples of human MSC cells were used after patient’s written consent using guidelines approved by the Ethic Committee of the RWTH University of Aachen (Permit number: EK128/09).
Mouse hematopoietic progenitor cells (MPP & CDP)
After sorting, MPP and CDP were treated with 10 ng/mL recombinant human TGF-β1 (R&D Systems, Minneapolis, USA) for 2, 4, 8, 12 and 24 h as described in  or left untreated. Cells were lysed in 350 μl TRI-Reagent and stored at −80°C.
RNA was isolated using the MagMAX-96 for Microarrays kit (Life Technologies, Darmstadt, Germany) according to manufacturer’s protocol.
Primary mouse hepatocytes and human HepG2 cells (HPC)
Hepatocytes (HPC) represent the most prominent cell population in the liver. Primary HPC are sensitive to TGF-β1, and express the corresponding type I (ALK5), type II (TβRII), and type III (betaglycan) receptors. TGF-β1 promotes cell cycle arrest and apoptosis of primary HPC. In addition, in vitro TGF-β1 provokes epithelial-to-mesenchymal transition (EMT)-like processes in this hepatic cell subpopulation, which most likely do not occur in vivo .
Primary murine HPC were isolated from male C57BL/6 mice according to the collagenase method of Seglen . Cells were plated in collagen coated 6-well dishes at a density of 1.2 × 106 cells using HepatoZYME-SFM (Gibco, Life Technology, Darmstadt, Germany). Four hours after seeding the medium was renewed and cells were grown for a further 24 hrs culture period.
HepG2 (DSMZ: DSM ACC180) were cultured in RPMI (PAA, Pasching, Austria) containing 10% fetal calf serum (PAA), 1 × Penicillin/Streptomycin (Lonza, Cologne, Germany). Medium was renewed every second day. For the experiment, cells were passaged and plated in 6-well dishes using accutase (PAA) at a density of 4 × 105 cells. One day before the experiment, cells were washed with PBS (1×), medium changed to HepatoZYME-SFM (Gibco) and cultured for further 24 hrs.
One hour before the experiment, the medium was exchanged and cells stimulated with 1 ng/mL recombinant human TGF-β1 (R&D Systems, Minneapolis, USA) for indicated time intervals (HepG2: 0 min, 20 min, 1 h, 2 h, 4 h, 24 h; murine HPC: 0 min, 1 h, 2 h, 4 h). The cells were harvested using Qiazol for cell lysis (Qiagen, Hilden, Germany), directly frozen and stored at −80°C.
RNA was isolated using the RNeasy Kit system (Qiagen), performing a DNAse digestion according to the manufacturer’s protocol.
Human Mesenchymal Stromal Cells (MSC)
Mesenchymal stromal cells (MSC) are found in all supportive tissue as in fat tissue, bone marrow and cord. MSC are characterized by their plastic adherence and their differentiation potential towards adipogenic, osteogenic and chondrogenic lineages. All MSC express the surface markers CD29, CD73, CD90 and CD105 and they lack the expression of CD14, CD31, CD34 and CD45 [14, 51].
Isolation and expansion
MSCs were isolated from mononuclear cells (MNCs) by plastic adherence as described before [16, 52, 53]. In brief, bone fragments from the caput femoris of patients undergoing femoral head prosthesis were flushed with phosphate- buffered saline (PBS) and washed twice with PBS. MNC were then resuspended in culture medium and seeded into tissue culture flasks. Cells were cultured at 37°C in a humidified atmosphere with 5% CO2. The first medium exchange was performed after 48 h to remove nonadherent cells. Thereafter, media changes were performed twice per week and MSCs were passaged when reaching 80-90% of confluence.
MSC from three different donors were used in an early passage (p3-5) for stimulation with TGF-b1. 1x106 MSC were seeded into 6-well culture plates. When the cells were attached after 24 h 1 ng/mL recombinant TGF-β1 (R&D Systems) was added to the culture media at different time points. The cells were harvested at the same time point with Qiazol (Qiagen) and directly frozen and stored at −80°C.
RNA was isolated via phenol/chloroform extraction using the miRNeasy Kit (Qiagen), performing a DNAse digestion.
Human samples were assayed using Affymetrix GeneChip® type “Human Gene 1.0 ST Array” with 34,760 probe-sets and mouse samples were assayed in Affymetrix GeneChip type “Mouse Gene 1.0 ST Array” with 32,321 probe-sets. Hybridization, wash and staining were done according to manufacturer’s recommended standard techniques.
Normalization and preprocessing
Raw probe intensities were normalized and summarized to expression levels using the FARMS algorithm which utilizes a factor analysis approach . A rigorous quality assessment confirmed a fairly good quality of the chips with exception of mouse HPC chips where Initial chip quality assessment revealed a strong batch effect and one bad chip (replicate no. 1 at time 1 hour). The bad chip was excluded and batch adjustment was performed to alleviate that effect on those chips via the “ComBat” method .
Affymetrix probe IDs were mapped to Entrez gene IDs using the Bioconductor annotation packages “mogene10sttranscriptcluster.db” in mouse chips and “hugene10sttranscriptcluster.db” in human chips [56, 57].
Differential gene expression
Time point specific analysis
Differential gene expression analyses via “Limma” Linear Models for Microarray Data  using empirical Bayes method  was performed by comparing samples at each time point after TGF-β stimulation to the unstimulated cells at time point 0. Statistical dependencies of samples between time points and replicates were considered via a factorial design matrix in “Limma” using a “time” and a “replicate” factor, and contrasts are considered for interaction effects. Corrections for multiple testing was done using the Benjamini & Hochberg’s method . Significant differentially expressed genes are considered those with FDRBH < 0.01 and absolute log2 fold change value (logFC) ≥ 1.
Analysis of whole time courses
The small number of replicates in our experiments limits the power of statistical testing procedures for assessing differential gene expression at individual time points. Furthermore, the number of measured time points is not the same for each cell type, which complicates any further meta-analysis. Therefore, we employed the “betr” method to analyze whole time series at once . The algorithm of this method uses a random-effects model together with the empirical Bayes method to estimate probabilities for differential expression of whole time courses. Genes were considered to be significant at a probability cutoff of > 99% for the whole time-course analysis and absolute log2 fold change (logFC) ≥ 1. Since “betr” requires the same number of replicates per time point and one chip in mouse HPC had to be omitted due to low quality (see Normalization & Preprocessing) unfortunately in this particular cell line we had to exclude time point 1 h completely from the time-course analysis.
Clustering of gene expression time series was done via the MFDA method proposed in . The method assumes gene expression time series within each given cluster to follow a mean curve plus some Gaussian noise. It decides cluster allocation via a Gaussian Mixture Model utilizing an EM algorithm for parameter estimation and decides the optimal number of clusters via the Bayesian Information Criterion (BIC).
It is worth mentioning that we applied MFDA not on raw gene expression data here, but on log fold-changes relative to time point 0. The reason was that we wanted to group genes not on the basis of their absolute expression values, but on the basis of similar response to the stimulus.
KEGG and go analysis
Analyses of pathways in KEGG  and biological processes in Gene Ontology project (GO)  were performed as follows: The –log(p-value) of all genes in the individual time point analysis and –log(1 – probabilities) of all genes in the whole time-course analysis, respectively, were taken as a ranking score for each transcript. Gene sets of KEGG pathways and GO terms were then tested for their association with these ranking scores via a univariate logistic regression based test [63, 64]. Resulting p-values of KEGG pathways and GO terms were adjusted according to Benjamini & Yekutieli’s false discovery rate control under dependency , and significant KEGG pathways and GO terms reported at a FDRBY cutoff value of 5%.
Transcription factor binding sites analysis
Analysis of transcription factors binding sites (TFBSs) was performed using the de novo sequence motif detection method XXmotifs . Identified sequence motifs were then aligned to known TRANSFAC TFBS via STAMP  and the top match is considered. The XXmotif method uses BLAST  all-against-all comparisons to mask regions of local homology in order to avoid false positives. The method then performs an enrichment analysis after transforming the found patterns to position weight matrices (PWMs). The STAMP method utilizes a global or un-gapped local alignment to detect DNA motifs similarities to defined PWMs. Furthermore, it considers familial binding profiles, thus improving transcription factors (TF) classification accuracy. TFBSs analysis was done using these methods in each cell type for those genes, which according to the time-course analysis showed a probability of > 0.99% for differential expression. Promoter sequences of the genes under consideration (2Kbp upstream of transcription start site) were obtained from the Ensembl database via “biomaRt” [69, 70]. Only the top matching motif for each TRANSFAC TFBS was considered and significant TFBSs were reported at an E-Value threshold of 1e-3.
Mapping of TFBS to individual transcription factors was performed via manual inspection of TRANSFAC PWMs. We obtained all proteins which had been used to construct each of the PWMs. With the help of the commercial software GeneGo Metacore®  we then mapped protein names to Entrez gene IDs. As a consequence we found for the TFBS FOXP1 the differentially expressed genes Foxp2 and FoxP1 (MPP, CDP mouse cells). For the TFBS FOX the human gene FAU was identified (human HPC). For KROX we found Egr1/EGR1 (mouse HPC, human MSC) and EGR2 (human MSC). For TEF we identified Klf3/KLF3 (MPPs, CDPs, and human HPC), TRIM37 and USP7 (human HPC, MSC).
Identification of homologous genes
Human homologs of mouse genes were identified via the KEGG Sequence Similarity Data Base (SSDB), which contains local alignments of amino acid sequences for protein coding genes from different species. We here considered two genes to be homologs, if the alignment E-value was below 1e-30 (bit-score > 112). In case of more than one homologous gene, all are considered.
Information about protein-protein interactions was collected separately for human and mouse from the BioGRID database version 3.2.109 . Correspondingly, a network comprising 16,011 nodes and 140,471 physical interactions was constructed for human. For mouse the network consisted of 6,233 nodes and 16,100 physical interactions. Nodes in these networks were weighted by the average probability (mean over all cell types from the same organism) for differential time course expression according to the betr model (see above). A “distance” for each edge was then calculated as 2 minus the sum of its incident nodes’ weights. Hence: the smaller the distance the higher the weight of its incident nodes. We used Dijkstra’s algorithm to search for minimum distance (i.e. maximum node weight) path connecting TGFB1 with each of SKIL, SMAD7, EGR1, PPARG and all genes annotated to glutathione metabolism, purine metabolism, oxidation-reduction process, innate immune response, negative regulation of apoptotic process, angiogenesis, positive regulation of cell proliferation and positive regulation of cell migration. For each of the last mentioned terms we only kept those genes as representatives which showed the minimum distance to TGFB1. If there were several paths of the same minimum distance, all of them were considered.
In the network for mouse Tgfb1 was not identified and hence we started our analysis with Tgfbr1 instead.
Availability of supporting data
The datasets supporting the results of this article are included within the article and its supplemental material files.
In addition, our microarray datasets for the different cell types analyzed here are deposited in the Gene Expression Omnibus database  under accession numbers GSE45942 (murine HPC), GSE45945 (human HPC), GSE46019 (human MSC) and GSE46109 (mouse MPP and CDP).
Transforming growth factor-beta
Common dendritic progenitor(s)
Linear models for microarray data
Mesenchymal stromal cell(s)
Plasmacytoid dendritic cell
Conventional dendritic cell
Lung adenocarcinoma cell-line
Bayesian estimation of temporal regulation
Linear models for microarray data
Transcription factor binding site
Kyoto encyclopedia of genes and genomes.
Renner U, Paez-Pereda M, Arzt E, Stalla GK: Growth factors and cytokines: function and molecular regulation in pituitary adenomas. Front Horm Res. 2004, 32: 96-109.
Giacomini D, Páez-Pereda M, Theodoropoulou M, Labeur M, Refojo D, Gerez J, Chervin A, Berner S, Losa M, Buchfelder M, Renner U, Stalla GK, Arzt E: Bone morphogenetic protein-4 inhibits corticotroph tumor cells: involvement in the retinoic acid inhibitory action. Endocrinol. 2006, 147: 247-256. 10.1210/en.2005-0958.
Hinck AP: Structural studies of the TGF-βs and their receptors - insights into evolution of the TGF-β superfamily. FEBS Lett. 2012, 586: 1860-1870. 10.1016/j.febslet.2012.05.028.
Oomizu S, Chaturvedi K, Sarkar DK: Folliculostellate cells determine the susceptibility of lactotropes to estradiol’s mitogenic action. Endocrinol. 2004, 145: 1473-1480. 10.1210/en.2003-0965.
Hentges S, Sarkar DK: Transforming growth factor-beta regulation of estradiol-induced prolactinomas. Front Neuroendocrinol. 2001, 22: 340-363. 10.1006/frne.2001.0220.
Alexandrow MG, Moses HL: Transforming growth factor beta 1 inhibits mouse keratinocytes late in G1 independent of effects on gene transcription. Cancer Res. 1995, 55: 3928-3932.
Massagué J: TGF-β signal transduction. Annu Rev Biochem. 1998, 67: 753-791. 10.1146/annurev.biochem.67.1.753.
Massagué J: The transforming growth factor-beta family. Annu Rev Biochem. 1990, 6: 597-641. 10.1146/annurev.cb.06.110190.003121.
Roberts ABSM: Physiological actions and clinical applications of transforming growth factor-beta (TGF-beta). Growth Factors. 1993, 8: 1-9. 10.3109/08977199309029129.
Roberts AB, Heine UI, Flanders KC SM: Transforming growth factor-beta. Major role in regulation of extracellular matrix. Ann N Y Acad Sci. 1990, 580: 225-232. 10.1111/j.1749-6632.1990.tb17931.x.
Orlova VV, Liu Z, Goumans M-J, Ten Dijke P: Controlling angiogenesis by two unique TGF-β type I receptor signaling pathways. Histol Histopathol. 2011, 26: 1219-1230.
Felker P, Seré K, Lin Q, Becker C, Hristov M, Hieronymus T, Zenke M: TGF-beta1 accelerates dendritic cell differentiation from common dendritic cell progenitors and directs subset specification toward conventional dendritic cells. J Immunol. 2010, 185: 5326-5335. 10.4049/jimmunol.0903950.
Seré K, Baek J-H, Ober-Blöbaum J, Müller-Newen G, Tacke F, Yokota Y, Zenke M, Hieronymus T: Two distinct types of Langerhans cells populate the skin during steady state and inflammation. Immunity. 2012, 37: 905-916. 10.1016/j.immuni.2012.07.019.
Dominici M, Le Blanc K, Mueller I, Slaper-Cortenbach I, Marini F, Krause D, Deans R, Keating A, Prockop D, Horwitz E: Minimal criteria for defining multipotent mesenchymal stromal cells. The Int Society for Cellular Therapy position statement. Cytotherapy. 2006, 8: 315-317. 10.1080/14653240600855905.
Wagner W, Ho AD: Mesenchymal stem cell preparations–comparing apples and oranges. Stem Cell Rev. 2007, 3: 239-248. 10.1007/s12015-007-9001-1.
Walenda G, Hemeda H, Schneider RK, Merkel R, Hoffmann B, Wagner W: Human Platelet Lysate Gel Provides a Novel 3D-Matrix for Enhanced Culture Expansion of Mesenchymal Stromal Cells. Tissue Eng Part C Methods. 2012, 1: 1-32.
Post S, Abdallah BM, Bentzon JF, Kassem M: Latent TGF-beta binding proteins (LTBPs)-1 and −3 coordinate proliferation and osteogenic differentiation of human mesenchymal stem cells. Bone. 2008, 43: 679-688. 10.1016/j.bone.2008.06.016.
Wang D, Park JS, Chu JSF, Krakowski A, Luo K, Chen DJ, Li S: Proteomic profiling of bone marrow mesenchymal stem cells upon transforming growth factor beta1 stimulation. J Biol Chem. 2002, 277: 36045-36051. 10.1074/jbc.M200575200.
Gao L, McBeath R, Chen CS: Stem cell shape regulates a chondrogenic versus myogenic fate through Rac1 and N-cadherin. Stem Cells. 2010, 28: 564-572.
Ito T, Sawada R, Fujiwara Y, Seyama Y, Tsuchiya T: FGF-2 suppresses cellular senescence of human mesenchymal stem cells by down-regulation of TGF-beta2. Biochem Biophys Res Commun. 2007, 359: 108-114. 10.1016/j.bbrc.2007.05.067.
Buenemann CL, Willy C, Buchmann A, Schmiechen A, Schwarz M: Transforming growth factor-beta1-induced Smad signaling, cell-cycle arrest and apoptosis in hepatoma cells. Carcinog. 2001, 22: 447-452. 10.1093/carcin/22.3.447.
Xu XM, Yuan GJ, Li QW, Shan SL JS: Hyperthermia Inhibits Transforming Growth Factor Beta-Induced Epithelial-Mesenchymal Transition (EMT) in HepG2 Hepatocellular Carcinoma Cells. Hepatogastroenterology. 2012, 59 (119): 2059-63.
Borkowski BTA, Letterio JJ, Farr AG, Udey MC: A Role for Endogenous Transforming Growth Factor 1 in Langerhans Cell Biology. Skin Trans. 1996, 184: 1-6.
Strobl H, Riedl E, Scheinecker C, Bello-Fernandez C, Pickl WF, Rappersberger K, Majdic O, Knapp W: TGF-beta 1 promotes in vitro development of dendritic cells from CD34+ hemopoietic progenitors. J Immunol. 1996, 157: 1499-1507.
Strobl H, Bello-Fernandez C, Riedl E, Pickl WF, Majdic O, Lyman SD, Knapp W: flt3 ligand in cooperation with transforming growth factor-beta1 potentiates in vitro development of Langerhans-type dendritic cells and allows single-cell dendritic cell cluster formation under serum-free conditions. Blood. 1997, 90: 1425-1434.
Sartor MA, Mahavisno V, Keshamouni VG, Cavalcoli J, Wright Z, Karnovsky A, Kuick R, Jagadish HV, Mirel B, Weymouth T, Athey B, Omenn GS: ConceptGen: a gene set enrichment and gene set relation mapping tool. Bioinformatics. 2010, 26: 456-463. 10.1093/bioinformatics/btp683.
Aryee MJ, Gutiérrez-Pabello JA, Kramnik I, Maiti T, Quackenbush J: An improved empirical bayes approach to estimating differential gene expression in microarray time-course data: BETR (Bayesian Estimation of Temporal Regulation). BMC Bioinform. 2009, 10: 409-10.1186/1471-2105-10-409.
Xu LL, Shanmugam N, Segawa T, Sesterhenn IA, McLeod DG, Moul JW, Srivastava S: A novel androgen-regulated gene, PMEPA1, located on chromosome 20q13 exhibits high level expression in prostate. Genomics. 2000, 66: 257-263. 10.1006/geno.2000.6214.
Xu LL, Shi Y, Petrovics G, Sun C, Makarem M, Zhang W, Sesterhenn IA, McLeod DG, Sun L, Moul JW, Srivastava S: PMEPA1, an androgen-regulated NEDD4-binding protein, exhibits cell growth inhibitory function and decreased expression during prostate cancer progression. Cancer Res. 2003, 63: 4299-4304.
Nomura N, Sasamoto S, Ishii S, Date T, Matsui M, Ishizaki R: Isolation of human cDNA clones of ski and the ski-related gene, sno. Nucleic Acids Res. 1989, 17: 5489-5500. 10.1093/nar/17.14.5489.
Cohen SB, Zheng G, Heyman HC, Stavnezer E: Heterodimers of the SnoN and Ski oncoproteins form preferentially over homodimers and are more potent transforming agents. Nucleic Acids Res. 1999, 27: 1006-1014. 10.1093/nar/27.4.1006.
Yu G, Li F, Qin Y, Bo X, Wu Y, Wang S: GOSemSim: an R package for measuring semantic similarity among GO terms and gene products. Bioinform. 2010, 26: 976-978. 10.1093/bioinformatics/btq064.
Wang JZ, Du Z, Payattakool R, Yu PS, Chen C-F: A new method to measure the semantic similarity of GO terms. Bioinform. 2007, 23: 1274-1281. 10.1093/bioinformatics/btm087.
Pesquita C, Faria D, Falcão AO, Lord P, Couto FM: Semantic similarity in biomedical ontologies. PLoS Comput Biol. 2009, 5: e1000443-10.1371/journal.pcbi.1000443.
Falcon S, Gentleman R: Using GOstats to test gene lists for GO term association. Bioinform. 2007, 23: 257-258. 10.1093/bioinformatics/btl567.
Ghosh AK, Bhattacharyya S, Lakos G, Chen S-J, Mori Y, Varga J: Disruption of transforming growth factor beta signaling and profibrotic responses in normal skin fibroblasts by peroxisome proliferator-activated receptor gamma. Arthritis Rheum. 2004, 50: 1305-1318. 10.1002/art.20104.
Liu R-M, Gaston Pravia KA: Oxidative stress and glutathione in TGF-beta-mediated fibrogenesis. Free Radic Biol Med. 2010, 48: 1-15. 10.1016/j.freeradbiomed.2009.09.026.
Yang EY, Moses HL: Transforming growth factor beta 1-induced changes in cell migration, proliferation, and angiogenesis in the chicken chorioallantoic membrane. J Cell Biol. 1990, 111: 731-741. 10.1083/jcb.111.2.731.
Luehr S, Hartmann H, Söding J: The XXmotif web server for eXhaustive, weight matriX-based motif discovery in nucleotide sequences. Nucleic Acids Res. 2012, 40 (Web Server issue): W104-W109.
Mahony S, Benos PV: STAMP: a web tool for exploring DNA-binding motif similarities. Nucleic Acids Res. 2007, 35 (Web Serv Issue): W253-W258.
Stark C, Breitkreutz B-J, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006, 34 (Database issue): D535-D539.
Yanagisawa K, Osada H, Masuda A, Kondo M, Saito T, Yatabe Y, Takagi K, Takahashi T: Induction of apoptosis by Smad3 and down-regulation of Smad3 expression in response to TGF-beta in human normal lung epithelial cells. Oncogene. 1998, 17: 1743-1747. 10.1038/sj.onc.1202052.
Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.
Yan X, Liu Z, Chen Y: Regulation of TGF-b signaling by Smad7 Overview of TGF-b Signaling Pathways. Acta Biochim Biophys Hung. 2009, 41: 263-272.
Lee KL, Lim SK, Orlov YL, Yit LY, Yang H, Ang LT, Poellinger L, Lim B: Graded Nodal/Activin Signaling Titrates Conversion of Quantitative Phospho-Smad2 Levels into Qualitative Embryonic Stem Cell Fate Decisions. PLoS Genet. 2011, 7: 19-
Chen S-J, Ning H, Ishida W, Sodin-Semrl S, Takagawa S, Mori Y, Varga J: The early-immediate gene EGR-1 is induced by transforming growth factor-beta and mediates stimulation of collagen gene expression. J Biol Chem. 2006, 281: 21183-21197. 10.1074/jbc.M603270200.
Chu AS, Diaz R, Hui J-J, Yanger K, Zong Y, Alpini G, Stanger BZ, Wells RG: Lineage tracing demonstrates no evidence of cholangiocyte epithelial-to-mesenchymal transition in murine models of hepatic fibrosis. Hepatol. 2011, 53: 1685-1695. 10.1002/hep.24206.
Aden DP, Fogel A, Plotkin S, Damjanov I, Knowles BB: Controlled synthesis of HBsAg in a differentiated human liver carcinoma-derived cell line. Nature. 1979, 282: 615-616. 10.1038/282615a0.
Knowles BB, Howe CC, Aden DP: Human hepatocellular carcinoma cell lines secrete the major plasma proteins and hepatitis B surface antigen. Science. 1980, 209: 497-499. 10.1126/science.6248960.
Seglen P: Preparation of isolated rat liver cells. Methods Cell Biol. 1976, 13: 29-83.
Horwitz EM, Le Blanc K, Dominici M, Mueller I, Slaper-Cortenbach I, Marini FC, Deans RJ, Krause DS, Keating A: Clarification of the nomenclature for MSC: The International Society for Cellular Therapy position statement. Cytotherapy. 2005, 7: 393-395. 10.1080/14653240500319234.
Lohmann M, Walenda G, Hemeda H, Joussen S, Drescher W, Jockenhoevel S, Hutschenreuter G, Zenke M, Wagner W: Donor age of human platelet lysate affects proliferation and differentiation of mesenchymal stem cells. PLoS One. 2012, 7: e37839-10.1371/journal.pone.0037839.
Koch CM, Reck K, Shao K, Lin Q, Joussen S, Ziegler P, Walenda G, Drescher W, Opalka B, May T, Brummendorf TH, Zenke M, Saric T, Wagner W: Pluripotent stem cells escape from senescence-associated DNA methylation changes. Genome Res. 2013, 23: 248-259. 10.1101/gr.141945.112.
Hochreiter S, Clevert D-A, Obermayer K: A new summarization method for Affymetrix probe level data. Bioinform. 2006, 22: 943-949. 10.1093/bioinformatics/btl033.
Johnson WE, Li C, Rabinovic A: Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostat Oxford England. 2007, 8: 118-127. 10.1093/biostatistics/kxj037.
Arthur L: mogene10sttranscriptcluster.db: Affymetrix Mouse Gene 1.0-ST Array Transcriptcluster Revision 8 annotation data (chip mogene10sttranscriptcluster). R package version 8.0.1. 2010, http://www.bioconductor.org,
Arthur L: hugene10sttranscriptcluster.db: Affymetrix Mouse Gene 1.0-ST Array Transcriptcluster Revision 8 annotation data (chip hugene10sttranscriptcluster). R package version 8.0.1. 2010, http://www.bioconductor.org,
Smyth GK: Limma: Linear Models for Microarray Data. Bioinform. 2005, pag: pages:397-pages:420.
Casella G: An Introduction to Empirical Bayes Data Analysis. Am Stat. 1985, 39: 83-87. 10.2307/2682801.
Benjamini Y, Hochberg Y: Controling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J Royal Stat Society Series B. 1995, 57: 289-300.
Ma P, Castillo-Davis CI, Zhong W, Liu JS: A data-driven clustering method for time course gene expression data. Nucleic Acids Res. 2006, 34: 1261-1269. 10.1093/nar/gkl013.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Sartor MA, Leikauf GD, Medvedovic M: LRpath: a logistic regression approach for identifying enriched biological groups in gene expression data. Bioinform. 2009, 25: 211-217. 10.1093/bioinformatics/btn592.
Montaner D, Dopazo J: Multidimensional gene set analysis of genomic data. PLoS ONE. 2010, 5: e10348-10.1371/journal.pone.0010348.
Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple hypothesis testing under dependency. Ann Stat. 2001, 29: 1165-1188. 10.1214/aos/1013699998.
Mahony S, Auron PE, Benos PV: DNA Familial Binding Profiles Made Easy: Comparison of Various Motif Alignment and Clustering Strategies. PLoS Comput Biol. 2007, 3: 14-10.1371/journal.pcbi.0030014.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mole Biol. 1990, 215: 403-410. 10.1016/S0022-2836(05)80360-2.
Flicek P, Amode MR, Barrell D, Beal K, Brent S, Chen Y, Clapham P, Coates G, Fairley S, Fitzgerald S, Gordon L, Hendrix M, Hourlier T, Johnson N, Kähäri A, Keefe D, Keenan S, Kinsella R, Kokocinski F, Kulesha E, Larsson P, Longden I, McLaren W, Overduin B, Pritchard B, Riat HS, Rios D, Ritchie GRS, Ruffier M, Schuster M, et al: Ensembl 2011. Nucleic Acids Res. 2011, 39 (Database issue): D800-D806.
Durinck S, Moreau Y, Kasprzyk A, Davis S, De Moor B, Brazma A, Huber W: BioMart and Bioconductor: a powerful link between biological databases and microarray data analysis. Bioinform. 2005, 21: 3439-3440. 10.1093/bioinformatics/bti525.
Durinck S, Spellman PT, Birney E, Huber W: Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protocol. 2009, 4: 1184-1191. 10.1038/nprot.2009.97.
GeneGo Metacore Software. [http://www.genego.com]
Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo/]
The authors would like to acknowledge the funding institutions that supported the different research projects of this study and promoted the co-operations under the umbrella of the interregional project “Alma-in-Silico”. These include: The Interreg IV European Fund and the committed regions and projects partners. The state of North Rhine-Westphalia represented in the German Stem Cells Network (http://www.stemcells.nrw.de), the German Ministry of Education and Research (BMBF) through researching grants of the German Research Foundation DFG : WA 1706/2-1, SFB/TRR57 P13, SFB542/B10, and SSP1356 Ze432/5-2.
HF is a member of the DFG funded excellence cluster ImmunoSensation.
The sponsors had no role in the study design, analysis, preparing the manuscripts or the decision to publish the work.
The authors declare that they have no competing interests.
HF, MZ, RW, KH, BS and KA conceived, designed and coordinated the study. GW, SM, KS, QL, BS and KF prepared the samples and carried out microarray experiments. KA, NM and HF performed data analysis. KA, HF, MZ, RW helped to interpret the results and to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.