- Research article
- Open Access
Analyzing M-CSF dependent monocyte/macrophage differentiation: Expression modes and meta-modes derived from an independent component analysis
BMC Bioinformatics volume 9, Article number: 100 (2008)
The analysis of high-throughput gene expression data sets derived from microarray experiments still is a field of extensive investigation. Although new approaches and algorithms are published continuously, mostly conventional methods like hierarchical clustering algorithms or variance analysis tools are used. Here we take a closer look at independent component analysis (ICA) which is already discussed widely as a new analysis approach. However, deep exploration of its applicability and relevance to concrete biological problems is still missing. In this study, we investigate the relevance of ICA in gaining new insights into well characterized regulatory mechanisms of M-CSF dependent macrophage differentiation.
Statistically independent gene expression modes (GEM) were extracted from observed gene expression signatures (GES) through ICA of different microarray experiments. From each GEM we deduced a group of genes, henceforth called sub-mode. These sub-modes were further analyzed with different database query and literature mining tools and then combined to form so called meta-modes. With them we performed a knowledge-based pathway analysis and reconstructed a well known signal cascade.
We show that ICA is an appropriate tool to uncover underlying biological mechanisms from microarray data. Most of the well known pathways of M-CSF dependent monocyte to macrophage differentiation can be identified by this unsupervised microarray data analysis. Moreover, recent research results like the involvement of proliferation associated cellular mechanisms during macrophage differentiation can be corroborated.
Since microarray technology has become one of the most popular approaches in the field of gene expression analysis, numerous statistical methods have been used to provide insights into the biological mechanisms of gene expression regulation. The high dimension of expression data and the complexity of the regulatory mechanisms leading to transcriptional networks still forces statisticians and bioinformaticians to examine available methods and to develop new sophisticated approaches. However, there are already appropriate methods using different approaches to examine the underlying biological mechanisms determining the gene expression signatures and profiles measured by microarray experiments. Supervised methods using prior knowledge like Gene Set Enrichment Analysis (GSEA) deliver useful results under certain conditions. But there is still a lack of reliable data needed for non-classical analysis. Widely used unsupervised approaches, like hierarchical clustering and k-means clustering, use correlations or other distance or similarity measures to identify genes with similar behavior under similar conditions. But these methods are not able to represent more complex structures and interdependencies in the regulatory machinery.
In contrast to the algorithms mentioned above, independent component analysis (ICA) explores higher-order statistics to decompose observed gene expression signatures (GES), which form the rows of the input data matrix, into statistically independent gene expression modes (GEM), which form the rows of matrix S according to the data model XT= AS. ICA solves blind source separation (BSS) problems, where it is known that the observed data set represents a linear superposition of underlying independent source signals. But it can more generally be considered a matrix decomposition technique which extracts informative features from multivariate data sets like, for example, biomedical signals like EEG (Electroencephalography) [1, 2], MEG (Magnetoencephalography)  and fMRI (functional magnetic resonance imaging) [4–6] recordings. ICA can also be considered a projective subspace technique appropriate for noise reduction [7, 8], or artifact removal [9, 10] if generated from independent sources.
In this work we will concentrate on the linear case, in which each single microarray GES is considered a linear superposition of unknown statistically independent GEM. To decompose these mixtures into statistically independent components, ICA algorithms like FastICA or JADE have been used. Typically, these GEMs can be interpreted as being characteristic of ongoing, largely independent biological regulatory processes. The philosophy behind can be expressed as: co-expression means co-regulation. But the complexity of gene regulation and the various interactions of cellular processes demands a new interpretation of our ICA-derived components. In the following we use these extracted GEMs to generate sub-modes, which may provide biological pathway information. The genes contained in these pathway-associated sub-modes can be regarded as more or less self-contained parts of larger regulatory networks, which can be represented by combining these sub-modes into meta-modes according to the functional role of the associated genes.
Here we used M-CSF dependent in vitro differentiation of human monocytes to macrophages as a model process to demonstrate that ICA is a useful tool to support and extend knowledge-based strategies and to identify complex regulatory networks or novel regulatory candidate genes.
The major known pathways associated to M-CSF receptor dependent signaling [11–13] include expansion of the role of the MAP-kinase pathway [14, 15] and Jun/Fos, Jak/Stat and PI-3 kinase [16–18] dependent signal transduction. Up-regulation of immune-regulatory components involved in innate immunity response (e.g. MHC), specific (e.g. Fc) [19–21] and nonspecific (CRP, complement, galectins) [22–26] opsonin receptors as well as charge and motif pattern recognition receptors (e.g. SR-family, LRP, Siglecs etc.) [27–30], is characteristic for monocyte/macrophage differentiation. Beyond this, an increase of membrane biogenesis, vesicular trafficking and metabolic pathways including amino acids, glucose, fatty acids and sterols, as well as increased activity of lysosomal hydrolases that enhance phagocytotic function [31, 32], autophagy  and recycling is triggered through M-CSF signaling as a hallmark of innate immunity . These mechanisms are tightly coupled to changes in cytokine/chemokine response  and red/ox signaling (NOS e.g. NADPH-Oxidase, Glutathione, Thioredoxin, Selenoproteins) that drive chemotaxis migration, inflammation (e.g. Nfκ B), apoptosis (eg. Caspases, TP53, Nfκ B, ceramide) and survival [36–42].
Results and Discussion
M-CSF dependent monocyte to macrophage differentiation involves the activation and regulation of many different cellular pathways. In this study we used several microarray experiments and combined them to a data set, which we analyzed using the JADE algorithm. The extracted GEMs were labeled from 1 to 14, according to decreasing energy. Note that the extracted GEMs show positive as well as negative components. They are partitioned into a sub-mode containing the negative signals only, denoted by i.1, and a corresponding sub-mode of the positive signals, denoted by i.2, respectively. These sub-modes were then combined into so-called meta-modes according to the following super categories deduced from the MeSH-filter used: Apoptosis, signal transduction, cell cycle and regulatory sequences, see table. Sub-classification and mapping to distinct pathways was then performed with the extracted sub-modes using the BiblioSphere MeSH- and GeneOntology-filter tools. Note that our method not only takes into account that one gene can be part of more than one pathway, but also that one pathway can be involved in more than one cellular event. This cannot be achieved with classical clustering tools.
Within the meta-mode Signal transduction four sub-modes, 3.2, 6.2, 12.2 and 13.2 were combined together. The MAP-kinase pathway (Figure 1) could be identified as the major signal transduction pathway in sub-modes 3.2 and 12.2. Sub-mode 6.2 encompassed the functions signal transduction and cell communication. The remaining sub-mode 13.2 could not be mapped to a defined pathway, but the majority of genes within this sub-mode are associated with innate immunity and defense functions. Among these we identified relevant genes, also related to signal transduction, like CD86, BLNK. The transcription factors LMO2 and FLI1 were unique in sub-mode 13.2 whereas MMP9, CD36, CTSK, C1QR1 and MYCL1 as a TF were also present in several other sub-modes.
The 12 and 18 respectively, identified MAPK-pathway genes were all unique within their sub-modes (Table 1), except IL8 and DUSP1, which were present in both sub-modes. IL8 is a member of the CXC chemokine family and thus one of the major mediators of the inflammatory response. It is also a potent angiogenic factor and has a signalling function in the FAS-pathway, whereas DUSP1 is assumed to play an important role in the human cellular response to environmental stress, as well as in the negative regulation of cellular proliferation. Another central gene of the MAPK-pathway is caspase-1 (CASP1), which was represented in sub-mode 12.2 (Figure 2). Caspase-1 is responsible for the maturation of the multi-functional cytokine interleukin-1β and as member of the FAS caspase cascade it is involved in FAS mediated cell death . Further remarkable genes associated with MAP-kinase in this sub-mode were S100A8, S100A9, GADD45B, CTSK, SOD2 and the transcription factors JUNB and ATF3, since they were all represented in other sub-modes or pathways, or play a central role in the MAPK-pathway.
Sub-mode 3.2 combined the MAPK-pathway with the thioredoxin (TRX) reductase/thioredoxin system. TRX is involved in a variety of oxidation reduction reactions that regulate cell growth and survival decisions . It reduces ligand binding and DNA interaction by oxidizing cysteine residues within the DNA binding domain of glucocorticoid hormone receptors. Furthermore, TXNDC14 and TXNRD1 were found in this sub-mode. TRX also seems to be up-regulated by NGF through MAPK1 . Other genes associated with the MAPK-pathway were: STK17A, SH3BP5, RPS6KA1, CD44, G6PD, IL1RN and the transcription factors EGR2.
In sub-mode 6.2 all of the 29 genes involved in signal transduction were also related to the MeSH-term cell communication. Five of those signalling genes CFLAR, TXNDC1, YWHAZ, NOTCH2 and PSEN1 were also involved in the negative regulation of cell death.
The MeSH-term regulatory sequences is described as nucleic acid sequences involved in gene expression regulation. This meta-mode combines genes mapped to the TP53-pathway (sub-mode 14.1) and genes related to the oncogenes JUN/FOS (sub-mode 4.1 and 10.1), which are members of a family of transcription factors containing the basic-region-leucine zipper or bZIP motif. The BiblioSphere software did not define a specific pathway for sub-mode 11.2, but there were a couple of peptidases and proteinases like LYZ, GGH and CPM as well as a remarkable number of classical targets for the SREBP transcription factors, regulating cholesterol and fatty acid metabolism: SQLE, CYP51A1, HMGCR, FDFT1, INSIG1, IDI1, SC5DL and LDLR.
Sub-mode 14.1 represented an intersection of genes involved in gene expression regulation and the TP53 pathway. Genes which fulfill both criteria were ADM, CCND2, CD59, CDC42, DUSP6, GADD45A, GCH1, IER3, NDUFV2, PIM1, SLC2A3 and UBE3A. Moreover, sub-mode 14.1 received high significance values (Z-Score) for the three other meta-mode categories and was also the sub-mode with the highest amount of genes represented in other sub-modes as well. This can be interpreted as an evidence for the complex and networked nature of gene expression regulation and the interactivity of cellular pathways.
The transcription factor JUN also known as c-Jun belongs to the family of c-Jun N-terminal kinases (JNKs) which are important for development and survival of macrophages . Sub-modes 4.1 and 10.1 combined twelve genes with a known relationship to the JUN/FOS pathway: CCND2, CREM, CXCL1, GADD45A, IL1RN, JUN, MAPK13, MARCKS, RALA, PLAU, S100A8 and SOD2.
Differentiation, Cell Cycle
The meta-mode cell cycle was completely governed by the TP53 pathway. Although all three sub-modes 5.2, 11.1 and 12.1 represented TP53 related genes, the intersection of genes was marginal. Only the genes DUSP6, PCNA and PRKCA were mapped to the TP53 pathway and were also present in the sub-modes 5.2 and 12.1. Sub-mode 11.1 represented genes specialized in cell cycle pathways regulating the interphase and in particular the G1 phase, since it contained the genes PPP1R15A, DUT, CD44, CDKN1A and SMC4L1. Sub-modes 5.2 and 12.1 mainly represented genes involved in cell growth and proliferation.
Sub-mode 5.2 was characterized by the TP53 related genes DHFR, VCAN, APP, EIF2AK2 and the transcription factor HMGB2 and HMGB3. Here, the latter has not been mapped to TP53 pathway but is mentioned here because of its strong relation to HMGB2.
The unique TP53 genes in sub-mode 12.1 were: CAMK1, CTSB, GSTN1, NME1, HMGCR, GSN, CYP51A1 and IL1RN.
Apoptosis related pathways play a major role during the differentiation of monocytes to macrophages. Here we introduce the term "survival/apoptosis" for the MeSH term apoptosis, because the identified apoptosis pathways here function as survival mechanisms for the differentiating cells. It has been shown, that an absence of M-CSF induces apoptosis in cultivated monocytes . Since apoptosis is regulated through many different pathways and regulatory mechanisms, we could identify seven sub-modes (2.1, 3.1, 6.1, 8.1, 9.2, 13.1, 4.2) related to apoptosis. These could be classified to four different pathways involved in the regulation of apoptosis: TP53 pathway, BAX pathway, FAS-pathway and calreticulin (CALR) regulated apoptosis. Three of these sub-modes represented only one pathway. Sub-modes 2.1, 6.1 were mapped to the TP53 pathway and sub-mode 4.2 is governed by CALR regulated apoptosis, whereas the others could be mapped to more than one pathway.
Due to the strongly networked nature of biological regulatory mechanisms, a lot of genes involved in more than one pathway can be regarded as connections between those. Toshiyuki and Reed  showed that the human BAX-gene is directly regulated by TP53 (TP53), whereas BAX is participating in the regulation of endoplasmatic reticulum Ca2+  as well. In this way it acts as a gateway for selected apoptotic signals. This was represented by the sub-modes 3.1, 8.1 and 13.1 which could comparably be mapped to the TP53 and BAX pathway. Sub-mode 8.1 here combined the most interesting combination of genes. The genes CCL3, CCND3, PAICS, FYB, AKAB1, IL1RN, CXCL1, MT1A and the TFs EGR2 and ATF3 could be implicated with BAX. These genes overlapped with five of the seven genes mapped to the TP53 pathway: ATF3, BAX, CSPG2, EIF5B and IL1RN. Furthermore, the metallo-thioneins which are suggested to regulate DNA binding activity of TP53, MT1A, MT1F, MT1B and MT1X were represented in this sub-mode .
The role of CALR as a major Ca2+-binding (storage) protein in the lumen of the endoplasmatic reticulum is well known . Consequently, one might imagine that CALR is involved in the regulation of apoptotic signals. The following genes of sub-mode 4.2 are related to CALR: SLC11A1, CD93, PROCR, NME1 and ATP2B1. All of these genes, except ATP2B1, passed the MeSH-filter apoptosis. The link to the TP53 pathway is the transcription factor FOXO1A (also found in sub-mode 6.1) and PRKCB1, which is also involved in various other cellular signaling pathways.
The member of the TNF-receptor superfamily FAS plays a central role in the regulation of programmed cell death. Sub-mode 9.2 contained eleven genes related to FAS: GSTM1, RALGDS, ALOX5, VCAN, S100A9, S100A8, VIL2, LY75, STAB1, HEBP2 and CD44.
Although not all sub-modes could be mapped to specific meta-modes, the remaining sub-modes still provide useful information. While the genes sorted to sub-modes 7.1 and 7.2 deliver no significant pathway information, they share common behavior. Genes of sub-mode 7.1 were all down-regulated in macrophages or up-regulated in monocytes, respectively, whereas genes of sub-mode 7.2 were up-regulated in macrophages. Among these, known marker-genes for the different cell types could be identified: MNDA, FCN1 and the S100 calcium binding proteins S100A8, S100A9 and S100A12 as monocyte and IGF2R, TSPAN4, MMP9, CTSK, MMD, TNS1 and CALR as macrophage genes.
Furthermore, the sub-modes 5.1, 4.1, 8.2 and 14.2 contained Major Histocompatibility Class (MHC) genes. Whereas the sub-mode 5.1 genes HLA-A and HLA-C belong to MHC class I, the MHC genes of the three other sub-modes belong to MHC class II which are: HLA-DQB1, HLA-DQA1, HLA-DPB1, HLA-DPA1 and HLA-DMB.
Additional meta-mode tables are shown in Additional file 1.
It has been stated [52, 53] that the use of ICA for the analysis of gene expression data is a promising tool, but there is still a lack of a careful discussion of the results. Here we emphasized the exploration of the biological relevance and obtained a detailed insight into the networked structure of the underlying regulatory mechanisms. Two MAP kinase related pathways could be identified as the main regulatory pathways during differentiation: the classical MAP kinase pathway and the JNK and p38 MAP kinase pathway, see figure 1. These results confirm expectations, according to which the MAP kinase pathway is activated by the M-CSF stimulus and functions as the main signal transduction pathway triggering macrophage differentiation and related pathways.
The conspicuous presence of TP53 associated pathways in M-CSF induced monocyte differentiation is associated with a dramatic regulation of cell-cycle and apoptosis related genes. This leads to the assumption that human mononuclear phagocytes, which are considered to be arrested to non-proliferating cells, still preserve proliferative potential .
Furthermore, we could show that ICA is able to distinguish between monocytes and macrophages concerning differential gene expression. This helpful attribute can be used to find specific marker genes not only for different cell types as it is shown here, but also for different tissues or normal and tumor cells.
Moreover, we were able to identify different regulatory mechanisms during M-CSF dependent differentiation. Although signal transduction pathways are mainly regulated by protein modifications like phosphorylation or acetylation, genes associated to specific pathways involved in macrophage differentiation could be separated into sub-modes only by analyzing gene expression signatures and their related gene expression modes. Furthermore, this analysis could be improved by combining gene expression sub-modes extracted from different microarray experiments into informative gene expression meta-modes. The results are in full agreement with the experimental literature on M-CSF dependent differentiation  and illustrate the potential power of such information-theory-based, unsupervised and data-driven analysis methods, see Figure 3.
To fully explore the potential of such information-theory-based unsupervised analysis tools and especially to determine the suitability and reliability of ICA for the analysis of microarray datasets, further investigations are needed. The algorithms still suffer from the fact, that the number of estimated independent components, i.e. the extracted gene expression modes, depends on the number of available gene expression signatures and the dimension of the related gene expression profiles. Therefor, the availability of greater datasets should lead to advancements, and as shown here, greater datasets can be obtained by the careful combination of smaller datasets.
For our analysis we combined the gene-chip results from three different experimental settings. In each experiment human peripheral blood monocytes were isolated from healthy donors (experiment 1 and 2) and from donors with Niemann-Pick type C disease (experiment 3). Monocytes were differentiated to macrophages for 4 days in the presence of M-CSF (50 ng/ml, R&D Systems). Differentiation was confirmed by phase contrast microscopy. Gene expression profiles were determined using Affymetrix HG-U133A (experiment 1 and 2) and HG-U133plus2.0 (experiment 3) GeneChips covering 22215 probe sets and about 18400 transcripts (HG-U133A). Probe sets only covered by HG-U133plus2.0 array were excluded from further analysis. In experiment one pooled RNA was used for hybridization, while in experiment two and tree RNA from single donors were used. The final data set consisted of seven monocyte and seven macrophage expression profiles and contained 22215 probe sets. After filtering out probe sets which had at least one absent call, 5969 probe sets remained for further analysis. The complete data set is publicly available in the NCBI Gene Expression Omnibus  through the accession number GSE9801.
The bulk of preprocessing has been done using the Affymetrix GeneChip Operating Software (GCOS), where default presets were used. Additionally, we applied a logarithmic correction to the data. This has been done because effects with multiplicative behavior, which may contain biological relevant information, become linear after logarithmic transformation. Another reason is that, untransformed microarray expression profiles have a strongly skewed, hence unbalanced distribution. This means, that there is a large amount of expression values near zero whereas only very few genes show high expression levels (Figure 4). To avoid adverse effects caused by such unbalanced distributions, we applied a logarithmic transformation. The final data are usually represented as a data matrix whose columns represent expression signatures of N genes while the rows represent M corresponding gene expression profiles.
JADE-based extraction of gene expression modes
The algorithm JADE is based on fourth-order cumulant tensors T z of pre-whitened input data z = Qx given by
with the kurtosis = Cum(z i z i z i z i ) being the corresponding autocumulant. Associated with these cumulants is a fourth-order signal space (FOSS) which defines the range of all mappings T z : M → T z (M)
The corresponding matrices [T z (M)] ij will be called cumulant matrices in the following. Note that the dimensionality m of the FOSS equals at most the number of sources.
A spectral representation of the cumulant matrices can be obtained using the column vectors of the whitened mixing matrix with the corresponding eigenvalues related to the kurtosis of the independent components. This spectral representation can be used to obtain an eigenmatrix decomposition of the cumulant tensor according toT z (E(q)) = μ q E(q)
with 0 ≤ q ≤ m2 symmetric eigenmatrices E(q)= u q and u q the q-th column of the mixing matrix U, and μ q being a scalar eigenvalue.
After whitening, a m × m – dimensional orthogonal matrix D = [d(0) ... d(m-1)], which jointly diagonalizes all eigenmatrices of T z , is found by maximizing the joint diagonality criterion
where Diag(·) denotes the vector of diagonal matrix elements. The joint diagonalizer D is then equivalent to the whitened mixing matrix U, hence the unknown independent component expression mode can be estimated easily.
Sub-modes and meta-modes
As result of an ICA analysis of a set of gene expression signatures representing the rows of the transpose data matrix XT, we obtain a matrix S of independent components (the rows of S) which represent independent gene expression modes (GEMs) as well as a matrix A of basis vectors of the new feature space. To deduce meaningful biological information from the GEMs, the discovery of specific biological processes, which determine the modes, is the goal of our expression mode analysis. After decomposing the data matrix with ICA, each GEM has been split into two sub-modes which can be considered to feature genes which are co-expressed, thus co-regulated by the underlying regulatory process. A GEM consists of scores of gene contributions to the sub-modes which account for the observation that excitatory as well as inhibitory regulations exist. In order to extract the most significant genes, various statistical tools can be applied which, however, often suffer from the small M large N case. Therefore, in most cases a threshold is simply applied, or, after ranking, a fixed number of top and bottom genes are chosen and further analyzed . The rational behind these methods is that each extracted gene expression sub-mode is best represented by its most active genes. However, the choice of threshold or number of active genes is non-trivial, and will influence the results considerably. In this study we assume instead that mapping to distinct pathways is most non-ambiguous by using a relatively small number of genes.
Here, we took a different approach by selecting genes that are extremal with respect to some probabilistic model. For each GEM y(i) ∈ ℝ, where i indexes the genes, we calculated the first four central moments corresponding to mean, standard deviation, skewness and kurtosis of the underlying data distribution. These shape parameters are then used to fit a density according to the Pearson family  using maximum-likelihood, see Figure 5. We chose a Pearson density as prior since it allows for flexible modeling with respect to these first four moments, which seemed crucial as for example skewness varies considerably between modes, see Figure 5, and high kurtotic as well as close-to-Gaussian modes were present.
We then used the estimated Pearson densities to determine the 1 - α and α percentiles for α = 1%. Samples that lie below the 1-percentile are denoted as significantly down-regulated genes, and genes above the 99-percentile as significantly up-regulated genes. The corresponding sub-modes were labeled as i.1 for down-regulation and i.2 for up-regulation. In Table 2 we list the number of significant genes in each sub-mode.
We analyzed the gene sub-modes with BiblioSphere . BiblioSphere is a data mining tool intended to provide gene relationships from literature databases and genome-wide promoter analysis. The probe sets were mapped to transcripts and to known genes with use of the Genomatix database. To uncover the biological meaning of the genes in the sub-mode, we applied the MeSH-Filter (Medical Subject Headings) to our data, which is the National Library of Medicine's controlled vocabulary thesaurus. We decided to use the category biological sciences as filter criterion. Co-citations between the genes of the sub-mode were taken into account by using the literature mining tool of the BiblioSphere software. Interesting terms were identified through Z-Scores which indicate over-representation of genes in the referring biological categories. Z-Scores are given by Z - Score = (n - )/σ n where n is the number of observed genes meeting any given criterion, is the corresponding expected number and σ n gives the standard deviation of n. All terms mentioned in this work are significant with respect to the Genomatix guidelines.
Depending on our filter analysis we defined several meta-modes, where we combined sub-modes with similar categories. In some cases we subclassified sub-modes within one meta-mode. In this way 4 meta-modes could be generated, whereas 17 of 28 sub-modes could be mapped to at least one meta-mode. For some meta-modes we displaced the MeSH-Term category with additional categories with respect to the underlying biology.
Additionally we used the KEGG pathway database for biochemical pathway analysis to more thoroughly characterize the biological relevance of a meta-mode. The genes corresponding to the meta-modes were mapped on database pathways using Pathway-Express which is part of the Onto-Tools provided by Intelligent Systems and Bioinformatics Laboratory .
Makeigh S, Bell A, Jung TP, Sejnowski T: Independent component analysis of electroencephalographic data. In Advances in Neural Information Processing Systems 8 (NIPS'1995). Cambridge, Masachusetts: MIT Press; 1995.
Habl M, Bauer C, Ziegaus C, Lang EW: Perspectives in Neuroscience: Artificial Neural Networks in Medicine and Biology. Springer Publishers Berlin 2000 chap. Analyzing Brain Tumor Related EEG Signals With ICA Algorithms; 131–136.
Vigario R, Jousmäki V, Hämäläinen M, Hari R, Oja E: Independent component analysis for identification of artifacts in Magnetoencephalographic recordings. In Advances in Neural Information Processing Systems 10 (NIPS'1997). Cambridge, Masachusetts: MIT Press; 1997:229–235.
Yang K, Rajapakse JC: ICA gives higher-order functional connectivity of brain. Neural Information Processing – Letters and Reviews 2004, 2: 27–32.
Keck IR, Theis FJ, Gruber P, Lang EW, Specht K, Puntonet CG: 3D spatial analysis of fMRI Data on a word perception task. In Lecture Notes in Computer Science, LNCS 3195. Edited by: Puntonet CG, Prieto A. Berlin: Springer Verlag; 2004:977–984.
Theis FJ, Meyer-Baese A, Lang EW: Second-order blind source separation based on multi-dimensional autoco-variances. In Lecture Notes in Computer Science (LNCS 3195), Proc. ICA'2004. Edited by: CG Puntonet AP. Berlin: Springer Verlag; 2004:726–733.
Tomé AM, Teixeira AR, Lang EW, Stadlthanner K, Rocha A: Blind Source Separation Using Time-delayed Signals. In Proceedings of the International Joint Conference on Neural Networks, IJCNN'2004, Volume CD. Budapest, Hungary; 2004.
Gruber P, Stadlthanner K, Böm M, Theis FJ, Lang EW, Tomé AM, Teixeira AR, Puntonet CG, Saéz JMG: Denoising using local projective subspace methods. Neurocomputing 2006.
Stadlthanner K, Theis F, Lang EW, Tomé AM, Gronwald W, Kalbitzer HR: A Matrix Pencil Approach to the Blind Source Separation of Artifacts in 2D NMR Spectra. Neural Information Processing – Letters and Reviews 2003, 1: 103–110.
Stadlthanner K, Theis F, Lang EW, Tomé AM, Gronwald W, Kalbitzer HR: Separation of water artefacts in 2D NOESY protein spectra using congruent matrix pencils. Neurocomputing 2005.
Shi C, Simon D: Integrin signals, transcription factors, and monocyte differentiation. Trends Cardiovasc Med 2006, 16: 146–152. 10.1016/j.tcm.2006.03.002
Pixley FJ, Stanley ER: CSF-1 regulation of the wandering macrophage: complexity in action. Trends in Cell Biology 2004, 14: 628–638. 10.1016/j.tcb.2004.09.016
Ross FP, Teitelbaum SL: α v β3and macrophage colony-stimulating factor: partners in osteoclast biology. Immunological Reviews 2005, 208: 88–105. 10.1111/j.0105-2896.2005.00331.x
Wada T, Penninger JM: Mitogen-activated protein kinases in apoptosis regulation. Oncogene 2004, 23: 2838–2849. 10.1038/sj.onc.1207556
Bogoyevitch M, Boehm I, Oakley A, Ketterman A, Barr R: Targeting the JNK MAPK cascade for inhibition: basic science and therapeutic potential. Biochim Biophys Acta 2004, 1697: 89–101.
BehreDagger G, Whitmarsh AJ, Matthew P, Coghlan MP, Hoang T, Carpenter CL, Zhang DE, Davis RJ, Tenen DG: c-Jun Is a JNK-independent Coactivator of the PU.1 Transcription Factor. J Biol Chem 1999, 274(8):4939–4946. 10.1074/jbc.274.8.4939
Fox S, Haque S, Lovibond A, Chambers T: The possible role of TGF-beta-induced suppressors of cytokine signaling expression in osteoclast/macrophage lineage commitment in vitro. Journal of Immunology 2003, 170: 3679–3687.
Stephens L, Ellson C, Hawkins P: Roles of PI3Ks in leukocyte chemotaxis and phagocytosis. Curr Opin Cell Biol 2002, 14: 203–213. 10.1016/S0955-0674(02)00311-3
Houde M, Bertholet S, Gagnon E, Brunet S, Goyette G, Laplante A, Princiotta MF, Thibault P, Sacks D, Desjardins M: Phagosomes are competent organelles for antigen cross-presentation. Nature 2003, 425: 402–406. 10.1038/nature01912
Vieira O, Botelho R, Grinstein S: Phagosome maturation: aging gracefully. Biochem J 2002, 366: 689–704.
Booth J, Trimble W, Grinstein S: Membrane dynamics in phagocytosis. Semin Immunol 2001, 13: 357–64. 10.1006/smim.2001.0332
Sobota A, Strzelecka-Kiliszek A, Gladkowska E, Yoshida K, Mrozinska K, Kwiatkowska K: Binding of IgG-Opsonized Particles to FcγR Is an Active Stage of Phagocytosis That Involves Receptor Clustering and Phosphorylation. Journal of Immunology 2005, 175: 4450–4457.
Swanson JA, Hoppe AD: The coordination of signaling during Fc receptor-mediated phagocytosis. Journal of Leukocyte Biology 2004, 76: 1093–1103. 10.1189/jlb.0804439
Mina-Osorio P, Ortega E: Signal regulators in FcR-mediated activation of leukocytes? Trends Immunol 2004, 25: 529–35. 10.1016/j.it.2004.08.008
Lau D, Mollnau H, Eiserich JP, Freeman BA, Daiber A, Gehling UM, Br"ummer J, Rudolph V, M"unzel T, Heitzer T, Meinertz T, Baldus S: Myeloperoxidase mediates neutrophil activation by association with CD11b/CD18 integrins. PNAS 2005, 102: 431–436. 10.1073/pnas.0405193102
Dumic J, Dabelic S, Flogel M: Galectin-3: an open-ended story. Biochim Biophys Acta 2006, 1760: 616–35.
Fabriek B, Dijkstra C, van den Berg T: The macrophage scavenger receptor CD163. Immunobiology 2005, 210: 153–60. 10.1016/j.imbio.2005.05.010
Minami M, Kume N, Shimaoka T, Kita T: Expression of SR-PSOX, a Novel Cell-Surface Scavenger Receptor for Phosphatidylserine and Oxidized LDL in Human Atherosclerotic Lesions. Arteriosclerosis, Thrombosis, and Vascular Biology 2001, 21: 1796–1800. 10.1161/hq1001.096652
Beutler B: Inferences, questions and possibilities in Toll-like receptor signalling. Nature 2004, 430: 257–263. 10.1038/nature02761
Lock K, Zhang J, Lu J, Lee S, Crocker P: Expression of CD33-related siglecs on human mononuclear phagocytes, monocyte-derived dendritic cells and plasmacytoid dendritic cells. Immunobiology 2004, 209: 199–207. 10.1016/j.imbio.2004.04.007
Desjardins M: ER-mediated phagocytosis: a new membrane for new functions. Nat Rev Immunol 2003, 3: 280–291. 10.1038/nri1053
Martin S, Parton R: Lipid droplets: a unified view of a dynamic organelle. Nat Rev Mol Cell Biol 2006, 7: 373–378. 10.1038/nrm1912
Schmitz G, Buechler C: ABCA1: regulation, trafficking and association with heteromeric proteins. Ann Med 2002, 34: 334–47. 10.1080/078538902320772098
Peiser L, Mukhopadhyay S, Gordon S: Scavenger receptors in innate immunity. Curr Opin Immunol 2002, 14: 123–128. 10.1016/S0952-7915(01)00307-7
Branton H, Kopp J: TGF-beta and fibrosis. Microbes Infect 1999, 1: 1349–65. 10.1016/S1286-4579(99)00250-6
Forman HJ, Torres M: Reactive Oxygen Species and Cell Signaling; Respiratory Burst in Macrophage Signaling. American Journal of Respiratory and Critical Care Medicine 2002, 166: S4-S8. 10.1164/rccm.2206007
Nordberg J, Arner ES: Reactive oxygen species, antioxidants, and the mammalian thioredoxin system. Free Radic Biol Med 2001, 31: 1287–1312. 10.1016/S0891-5849(01)00724-9
Wang Y, Qiao M, Mieyal J, Asmis L, Asmis R: Molecular mechanism of glutathione-mediated protection from oxidized low-density lipoprotein-induced cell injury in human macrophages: role of glutathione reductase and glutaredoxin. Free Radic Biol Med 2006, 41: 775–85. 10.1016/j.freeradbiomed.2006.05.029
Wang Y, Zeigler MM, Lam GK, Hunter MG, Eubank TD, Khramtsov VV, Tridandapani S, Sen CK, Marsh CB: The Role of the NADPH Oxidase Complex, p38 MAPK, and Akt in Regulating Human Monocyte/Macrophage Survival. American Journal of Respiratory and Critical Care Medicine 2006, 36: 68–77.
Cathcart MK: Regulation of Superoxide Anion Production by NADPH Oxidase in Monocytes/Macrophages. Arteriosclerosis, Thrombosis, and Vascular Biology 2004, 24: 23–28. 10.1161/01.ATV.0000097769.47306.12
Kustermans G, Benna JE, Piette J, Legrand-Poels S: Perturbation of actin dynamics induces NF-κB activation in myelomonocytic cells through an NADPH oxidase-dependent pathway. Biochem J 2005, 387: 531–540. 10.1042/BJ20041318
Østerud B, Bjørklid E: Role of Monocytes in Atherogenesis. Physiol Rev 2003, 83: 1069–1112.
Park D, Thomsen A, Frevert C, Pham U, Skerrett S, Kiener P, Liles W: Fas (CD95) induces proinflammatory cytokine responses by human monocytes and monocyte-derived macrophages. J Immunol 2003, 170: 6209–6216.
Bishopric NH, Webster KA: Preventing Apoptosis With Thioredoxin. Circ Res 2002, 90: 1237–1239. 10.1161/01.RES.0000025101.04065.83
Masutani HKYYJ, Bai J: Thioredoxin as a neurotrophic cofactor and an important regulator of neuroprotection. Mol Neurobiol 2004, 29: 229–249. 10.1385/MN:29:3:229
Himes S, Sester D, Ravasi T, Cronau S, Sasmono T, Hume D: The JNK are important for development and survival of macrophages. Journal of Immunology 2006, 176: 2219–2228.
Becker S, Warren M, Haskill S: Colony-stimulating factor-induced monocyte survival and differentiation into macrophages in serum-free cultures. The Journal of Immunology 1987, 139: 3703–3709.
Toshiyuki M, Reed JC: Tumor suppressor p53 is a direct transcriptional activator of the human bax gene. Cell 1995, 80: 293–299. 10.1016/0092-8674(95)90412-3
Scorrano L, Oakes SA, Opferman JT, Cheng EH, Sorcinelli MD, Pozzan T, Korsmeyer1dagger SJ: BAX and BAK Regulation of Endoplasmic Reticulum Ca2+: A Control Point for Apoptosis. Science 2003, 300: 135–139. 10.1126/science.1081208
Ostrakhovitch E, Olsson P, Jiang S, Cherian M: Interaction of metallothionein with tumor suppressor p53 protein. FEBS Lett 2006, 580: 1235–1238. 10.1016/j.febslet.2006.01.036
Arnaudeau S, Frieden M, Nakamura K, Castelbou C, Michalak M, Demaurex N: Calreticulin differentially modulates calcium uptake and release in the endoplasmic reticulum and mitochondria. Journal of Biological Chemistry 2002, 277: 46696–46705. 10.1074/jbc.M202395200
Liebermeister W: Linear modes of gene expression determined by independent component analysis. Bioinformatics 2002, 18: 51–60. 10.1093/bioinformatics/18.1.51
Chiappetta P, Roubaud M, Torresani B: Blind source separation and the analysis of microarray data. J Comput Biol 2004, 11: 1090–109. 10.1089/cmb.2004.11.1090
Martinez F, Gordon S, Locati M, Mantovani A: Transcriptional profiling of the human monocyte-to-macrophage differentiation and polarization: new molecules and patterns of gene expression. Journal of Immunology 2006, 177: 7303–7311.
Schmitz G, Grandl M: Role of redox regulation and lipid rafts in macrophages during Ox-LDL-mediated foam cell formation. Antioxid Redox Signal 2007, 9(9):1499–518. 10.1089/ars.2007.1663
Gene Expression Omnibus[http://www.ncbi.nlm.nih.gov/geo/]
J-F Cardoso AS: Blind beamforming for non-Gaussian signals. IEEE Proc 1993, 140: 362–370.
Cardoso JF, Souloumiac A: Jacobi Angles For Simultaneous Diagonalization. SIAM Journal Mat Anal Appl 1996, 17: 161–164. 10.1137/S0895479893259546
Lee SI, Batzoglou S: Application of independent component analysis to microarrays. Genome Biology 2003, 4(11):R76. 10.1186/gb-2003-4-11-r76
Nair NU, Sankaran PG: Characterization of the Pearson family of distributions. IEEE Transactions on Reliability 1991, 40: 75–77. 10.1109/24.75339
Genomatix Software GmbH, Munich, Germany[http://www.genomatix.de]
Intelligent Systems and Bioinformatics Laboratory[http://vortex.cs.wayne.edu/projects.htm]
KEGG: Kyoto Encyclopedia of Genes and Genomes[http://www.genome.jp/kegg/]
The authors thank the anonymous reviewers for their comments during the preparation of this manuscript. The authors acknowledge helpful suggestions of Thomas Langman and Christoph Moehle concerning the microarray data treatment. The authors would also like to thank Kurt Stadlthanner, Florian Bloechl and Peter Gruber for suggestions on the statistical analysis. Microarray experiments were conducted in the KFB (Kompetenzzentrum fuer Fluoreszente Bioanalytik), D-93053 Regensburg. This work has been performed in association with the SFB Transregio 13 and Reform C: Molecular Pathomechanisms of the Metabolic Syndrome.
DL conceived the study, performed bioinformatical analysis and drafted the manuscript. PU helped with the biochemical analysis. MG and EO performed the microarray studies. FT and EWL guided the bioinformatical analysis. GS guided the study and coordinated the project.
Electronic supplementary material
About this article
Cite this article
Lutter, D., Ugocsai, P., Grandl, M. et al. Analyzing M-CSF dependent monocyte/macrophage differentiation: Expression modes and meta-modes derived from an independent component analysis. BMC Bioinformatics 9, 100 (2008). https://doi.org/10.1186/1471-2105-9-100
- Independent Component Analysis
- Gene Expression Signature
- Gene Expression Mode
- Macrophage Differentiation
- Major Histocompatibility Class