Research article | Open | Published:
Identification of novel transcriptional regulators involved in macrophage differentiation and activation in U937 cells
BMC Immunologyvolume 10, Article number: 18 (2009)
Monocytes and macrophages play essential role in innate immunity. Understanding the underlying mechanism of macrophage differentiation and the identification of regulatory mechanisms will help to find new strategies to prevent their harmful effects in chronic inflammatory diseases and sepsis.
Maturation of blood monocytes into tissue macrophages and subsequent inflammatory response was mimicked in U937 cells of human histocytic lymphoma origin. Whole genome array analysis was employed to evaluate gene expression profile to identify underlying transcriptional networks implicated during the processes of differentiation and inflammation. In addition to already known transcription factors (i.e. MAFB, EGR, IRF, BCL6, NFkB, AP1, Nur77), gene expression analysis further revealed novel genes (i.e. MEF2, BRI, HLX, HDAC5, H2AV, TCF7L2, NFIL3) previously uncharacterized to be involved in the differentiation process. A total of 58 selected genes representing cytokines, chemokines, surface antigens, signaling molecules and transcription factors were validated by real time PCR and compared to primary monocyte-derived macrophages. Beside the verification of several new genes, the comparison reveals individual heterogeneity of blood donors.
Up regulation of MEF2 family, HDACs, and H2AV during cell differentiation and inflammation sheds new lights onto regulation events on transcriptional and epigenetic level controlling these processes. Data generated will serve as a source for further investigation of macrophages differentiation pathways and related biological responses.
The mammalian innate immune system, comprised of macrophages, provides a front line of defense against pathogens. Infection or dysfunction of macrophages results in clinical situations like septic shock or chronic inflammatory disease such as atherosclerosis and rheumatoid arthritis. Under normal steady-state conditions, monocytes migrate randomly from blood to various organs or body cavities and differentiate into macrophages through coordinate expression of numerous genes. During a local infection, migration of blood monocytes into the inflamed tissues is accelerated, and subsequent differentiation into macrophage occurs rapidly . Upon the entry of microorganisms, conserved structures found in a broad range of pathogens are recognized by toll-like receptors (TLRs) of macrophages. Phagocytic macrophages are activated to produce inflammatory cytokines (i.e. TNFalpha) and chemokines (i.e. IP10). Signaling pathways triggered by individual TLRs involve recruitment of adaptor molecules (i.e. MyD88) and autocrine production of IFNbeta. Two distinct cascades, MyD88 dependent and MyD88 independent (IFNbeta mediated), lead to activation of immuno-regulatory transcription factors (i.e. IRFs, STATs, NFkB, AP1) to modulate inflammatory gene expression .
Transcriptional regulators play important roles during developmental processes, and were described in a variety of well-studied differentiation systems. During the process of hematopoietic differentiation of monocytes, functions of transcription factors encompass up regulation of genes specific to the lineage and repression of lineage-inappropriate genes by forming distinct protein complexes . Similarly, cell-type-specific gene expression during macrophage differentiation and function is a result of complex transcriptional regulation by the interplay between endogenous transcription factors and external signals, and also by the networks of individual transcription factors. Whatever other mechanisms can modulate their activity, expression of specific transcription factors during differentiation process is likely to be regulated transcriptionally. Several transcription factors including MAFB, EGRs, NFkB, AP1, IRFs and BCL6 that control various stages of macrophage development or function have already been identified [4, 5]. However, only a genome-wide expression study will be able to identify all transcriptional regulators, gene regulatory networks and help to understand molecular mechanisms underlying relevant immune disease.
Furthermore, the role of chromatin structure in gene regulation of immune system has been increasingly studied. In activated macrophages stimulated by LPS, the transcriptional response has been attributed to the expression and function of chromatin remodeling genes, mainly related to histone modifications [6–8]. Timely coordinated relationship regarding transition of chromatin structure, nucleosome remodeling and transcription factors are known to be important . Under such a circumstance, it is of utmost importance to study genome-wide gene expression changes in model systems that make it feasible to continue investigating and dissecting underlying molecular mechanisms.
The use of a model system circumvents many problems including available cell numbers, high cost and tedious procedure, generally encountered while having to work with primary blood cells. However, the validity of any single myeloid model system has never been fully verified. Among several myeloid cell lines (HL60, THP-1, Mononomac-1, U937), U937 cells and THP1 cells are the most widely used models for investigating monocytic differentiation and subsequent biological functions of differentiated cells [10, 11]. In contrast to THP-1 cells at less mature stage due to their blood leukemic origin, U937 cells of histocytic lymphoma origin are arrested in a more advanced stage of differentiation (promonocyte/monocyte). Various stimuli (i.e. DMSO, Retinoic acid, VitD3, PMA, IFNgamma), either alone or in combination, have been attributed to induce their terminal differentiation into monocytes or mature macrophages. Upon differentiation, U937 cells acquire a large repertoire of macrophage function through the concerted expression of numerous genes. Differentiated U937 cells can be further stimulated with LPS (E. coli Lipopolysaccharide) to mimic inflammatory response of activated macrophages.
In this study, we employed the U937 model system combined with a cRNA hybridization-based whole genome array to analyze genome-wide expression patterns during monocyte-macrophage differentiation and LPS-responsive activation of differentiated cells. The identification of novel transcriptional regulators intensifies directions for further investigation of the molecular mechanism involving common regulatory networks and epigenetic mechanisms.
U937 cells were differentiated into monocyte and macrophage-like cells
U937 cells have been widely used as a model to investigate a variety of biological processes related to monocyte and macrophage function. Here, we used either the hormonally active form of Vitamin D3 (1, 25-dihydroxyvitamin D3) to differentiate U937 cells toward monocytes, or PMA (phorbol 12-myristate 13-acetate) to induce differentiation into a macrophage-like phenotype. Both PMA and VitD3 were known to induce cell cycle arrest prior to proceeding into differentiation stage. PMA-treated cells exhibit apparent growth arrest already 6 hrs after the addition of PMA and almost 100% adherence of cells after 24 hrs. Cells then tend to loose adhesiveness gradually afterwards. Growth arrest phenotype was much less prominent with VitD3-treated cells that continued moderate proliferation even after addition of the reagent (data not shown). Cells differentiated for 24 hrs by each inducer were used for further analysis and LPS stimulation.
Fluorescence activated cell-sorting (FACS) was carried out to monitor two surface proteins characteristic of differentiated U937 cells (data not shown): (1) CD14, serving as a marker for both monocyte and macrophage, is involved in LPS recognition and facilitating LPS contact with TLR4 receptor; (2) CD11b, serving as a macrophage marker, is involved in phagocytosis of bacteria. CD14 expression was increased after VitD3 differentiation, yielding the characteristics of monocyte-like cell type. This result is also in accordance with mRNA expression, and is certainly indicative of enhanced sensitivity of these cells to LPS stimulation . CD11b expression was increased only after PMA inducement, indicating macrophage-like differentiation, meanwhile no CD14 induction was observed with these cells.
According to well-known molecular mechanisms underlying macrophage activation, TNFalpha is a gene product of MyD88-dependent pathway (common to both monocytes and macrophages) while IFNbeta is an autocrine product destined to launch MyD88-independent pathway specific for macrophages . In the case of sepsis, persisting and systemic release of TNFalpha by activated macrophages can cause organ failure and death as a result of septic shock. Inevitably, it would require tight regulation. In order to demonstrate successful differentiation leading to activated macrophages, mRNA accumulation of TNFalpha and IFNbeta was analyzed by real-time PCR. The result from stimulation of PMA and VitD3 differentiated cells with 10 μg/ml LPS is shown in Figure 1. Immediate expression (2 hrs) of TNFalpha was common to all samples (control, VitD3, PMA). With control U937 cells and VitD3-differentiated cells, decreasing production of TNFalpha mRNA was observed across the later period of stimulation (12 ~22 hrs). PMA differentiated cells shows a dramatically decreased TNFalpha expression after 12 hrs but increased again at 22 hrs. This distinct biphasic behavior is a manifestation that macrophage-like cells possess an intrinsic 'automatic-bridle' system to modulate production of this cytokine in a time dependent manner. Stimulation at the lower concentration of LPS (2 μg/ml) showed the same pattern except that the biphasic behavior was extended up to 48 hours. Time dependent mRNA accumulation of IFNbeta was also characteristic for PMA differentiated cells representative of macrophages. Initially, a moderate increase of mRNA was observed (at 2 hrs). Afterwards, further increase was prominently induced (10 hrs). By the time of 22 hrs, the mRNA level was maintained at the level comparable to initial stage. On the contrary, control U937 cells and VitD3 differentiated cells manifest much less prominent accumulation of IFNbeta after both 2 and 10 hrs, and are completely devoid of it by the time of 22 hrs. It indicates that these cells are inefficient (during earlier period) and finally lacking (during later period) in IFNbeta induced MyD88-independent pathway.
Differential gene expression of U937 cells manifests functional macrophage phenotype
Whole transcriptome expression profiling was performed with two biological replicates and two technical replicates from each of seven different samples (control, PMA 6 hrs, PMA 12 hrs, PMA 24 hrs, PMA 32 hrs, LPS, VitD3) while satisfying each set of replicates on separate slides. To evaluate the consistency of data across biological or technical replicates, correlation coefficients were calculated. In all cases, the correlation between biological replicates and also technical replicates were higher than 0.98 (see Additional file 1). During PMA differentiation, there was a shift of correlation coefficients from 0.961 (after 6 hrs) to 0.919 (after 32 hrs), indicating 3.9% to 8.1% of gene features contained on the array were differentially regulated. There appeared only 0.6% (0.925 ~0.919) of change between 24 hrs and 32 hrs of PMA differentiation. This shows that the cells were already reached a state of full differentiation after 24 hrs. It is noteworthy that additional 3.8% of features were promptly induced within 2 hrs of subsequent LPS stimulation.
Gene expression of differentiated or stimulated cells was evaluated in comparison to undifferentiated U937 cells. Differential gene expression was considered 'significant' when the normalized intensity ratio of sample vs control was >1.9 (up) or <-2 (down), respectively (Figure 2A–D). At 6 hrs of PMA differentiation, 324 genes were up regulated and 169 genes down regulated. After 12 hrs, 337 genes were additionally up regulated while 126 genes were down regulated. By the time of 24 hrs, total 595 genes were up regulated and 324 genes down regulated. After 2 hrs of subsequent LPS stimulation, 278 additional genes were up regulated while 107 additional genes were down regulated. With VitD3 differentiation, only 35 genes were significantly up regulated while 15 genes down regulated. The complete list of gene expression profiling is available in Additional file 2.
Top 50 most highly up regulated genes, sorted by 24 hrs of PMA differentiation, are shown in Table 1. Among them are several genes well known to be related with differentiation process as well as monocyte and macrophage functions. This involves proteins known to be involved in the differentiation process (CDK inhibitor p21), cytokine transporter like A2M, chemokines like CCL3 and IP10 as well as matrix metalloproteinase protein such as MMP7 and MMP9 involved in the proteolysis of structural and adhesive matrixes.
Temporal pattern of array gene expression can be categorized
Large scale comparisons of differential gene expression profiles between different experimental conditions (for example, different time points of differentiation and activation) have potential to identify co-regulated and/or inter-related genes. Such categories reflect their functional relationships and can help to identify common transcriptional control mechanisms. Transcription factors in such groups of genes may be responsible for regulation of the genes in the same group. Beside genes not affected, three distinct groups could be categorized. Nine examples of time dependent gene expression were illustrated in Figure 3. IL6 (cytokine), CXCL2 (chemokine) and Nur77 (transcription factor) represent a group of genes that are highly up regulated only after LPS-responsive activation. NCF1, A2M, MMP9, all of which related to macrophage function, and transcription factor MEF2A represent a second group of genes whose expression was constantly up regulated during differentiation period and also activation. CTSG and TNFSF13B, represent a third group of genes specifically down regulated during macrophage differentiation. Both proteins are related to functions in neutrophils and B cells respectively.
Validity of expression analysis was verified by real timePCR
Among 58 selected genes validated by real time PCR (Table 2), were cytokines and chemokines (TNFalpha, IL1A, IL1B, IL6, IL8, CCL3, CCL4, CXCL2, IP10), inflammatory effector molecules (MMP7, MMP9, PLAU, C3, A2M, NCF1, LAMB3, LL37, CD38), surface antigens (CD11b, CD14, CD69, CD83), signaling molecules (MyD88, TLR2, SGK, JAK3) and transcription factors.
The fold change of greater than 20 from the array correspond to a factor of greater than 100 RQ fold (relative quantification, in which RQ = 2-ddCt) based on real time PCR (>ddCt -7). Likewise, array fold change between 10 and 20 usually fell within the range of 10 to 50 RQ fold while low fold change around 2 or 3 were equal to various range from less than 5 up to 30 RQ fold. Taking an example of LPS stimulated sample, confirmation rate of up or down regulation by real time PCR result was 100%. With PMA samples, a few cases of inconsistency, as with IL8, IL1A and CCL4, were observed
Identification of transcriptional regulators
Transcriptional regulators already known to be involved in monocytic-lineage differentiation (i.e. EGR1, EGR2 and MAFB) are constantly up regulated during the differentiation process and subsequent activation (see Additional file 3). Previously unknown factors (i.e. BRI, HLX1, TCF7L2, MEF2 family) were also shown to be significantly and constantly up regulated. Of the most interesting new identified transcription factors are proteins belonging to the MEF2 (myocyte enhancer factor 2) family. Moreover, a target gene (i.e. Nur77) of MEF2 during T cell apoptosis and postsynaptic differentiation were also up regulated. This is in line with previous reports [13–15]. Among the population of transacting proteins identified, several factors acting as transcriptional repressors (ID3, BCL6, NFIL3) and histone modifying proteins (HDAC5, HDAC7/9) were verified. Histone deacetylases (HDACs) were previously implicated as key molecules to repress production of pro-inflammatory gene expression, thus prevent excessive inflammatory responses . KLF2 and KLF4 have been recently described as critical regulators controlling monocyte differentiation and proinflammatory activation [17, 18]. Both factors are validated to be up regulated by real time PCR (Table 2). This together with up regulation of KLF6 suggests importance of kruppel-like family of transcription factors (KLFs) in the differentiation and function of macrophages. Upon LPS stimulation, NUR77, FOS, JUN, STAT3, IRF family, NFkB family and SMAD family show significant up regulation, thus confirm their classical role in macrophage function of inflammatory response. One of the main physiological roles of NFkB family is in the development and functioning of the immune system by regulating transcription of cytokines and antimicrobial effectors as well as genes that regulate cellular differentiation, survival and proliferation . Consistently, RELB, cREL, NFkB1 and NFkB2 are significantly induced upon LPS stimulation, thus confirming the validity of U937 model system.
Cross validation between U937 cells and primary blood monocytes
Various macrophage populations from different tissues of the body exhibit different morphological and functional phenotypes. Two cytokines, macrophage colony-stimulating factor (M-CSF) and granulocyte-macrophage colony-stimulating factor (GM-CSF), are known to promote differentiation of monocytes into macrophages with distinct phenotypes existing in vivo. It is known that M-CSF differentiated macrophages are very similar to peritoneal macrophages [20, 21]. As a counterpart of PMA-differentiated U937 cells, M-CSF differentiated macrophages were generated from primary blood monocytes, and further activated with LPS (see materials and methods). Real time PCR shows that both CD14 and CD11b were up regulated with M-CSF differentiated macrophages (Figure 4). Additionally, monocyte-derived macrophages manifest down regulation of TNFalpha in contrast to PMA-differentiated U937 cells where up regulation was shown.
Expression data of 24 top regulated transcriptional regulators from U937 cells including activators, repressors, histone modifying enzymes and histone variants were subjected to cross validation in comparison to primary human monocytes derived macrophages (M-CSF/macrophages). Two distinct groups of genes could be identified (Table 3). Within the first group, 54.2% (n = 13) of the genes (MEF2 family, HLX1, BRI, HDAC5, H2AV, BCL6, EGR1, ID3, IRF9, FOS, Nur77) were differentially expressed during differentiation and/or activation in both U937 cells and primary M-CSF/macrophages. Interestingly, both cell types manifest moderate up regulation (2 ~4 RQ fold) of H2AV. The expression consistency of BCL6, EGR1 and FOS depended upon the stage of cells (either differentiated or activated). Despite the generally known concept that Nur77 is induced only after LPS stimulation, one individual donor manifest up regulation already upon macrophage differentiation. Nevertheless, these factors are most likely to be true regulators in the differentiation and function of monocytic cells. A second group of genes (TCF7L2, HDAC7, EGR2, KLF2, KLF4, FOXA2, OCT2, HES1, PCAF, IRF7, MEF2B) showed a very opposite or no differential gene expression (45.8%; n = 11). Among these factors, TCF7L2 is noteworthy. It is a ubiquitous transcription factor in Wnt signaling pathway, but has not been connected to any immune cell gene transcription up to date.
Cellular localization of transcription factors
As a step toward verifying actual involvement of selected transcription factors in activating gene program, protein localization was traced and quantified in cytosolic and nucleic fraction by western blot. In order to probe 15 proteins among top regulated transcription factors, twenty different antibodies were tested by ELISA and western blot. Western Blot could detect proteins in cytosolic and nuclear fraction with a clear tendency toward nucleic fraction. Three antibodies were determined to meet the criteria (showing one specific band) that allow their use for protein quantification. Figure 5 shows the results of western blot with these three antibodies. MEF2 antibody recognizes three members (MEF2A, 2C and 2D) simultaneously. One specific band of about 55 kDa in size (white bar in Figure 5, corresponding to the size of MEF2A) was detected in cytosolic fraction while two bands of similar size and intensity (~50 kDa and ~60 kDa, corresponding to MEF2C and 2D) were detected in nucleic fraction (black bar in Fig. 5 represent upper size of ~60 kDa). While expression of cytosolic proteins appears to be decreased, nucleic proteins show 3 ~4 fold induction during macrophage differentiation and subsequent activation by LPS. Up regulation of proteins were comparable to mRNA induction of 3 ~9 RQ fold (ddCt value -1.6 ~-3.2, Table 2). Both NFkB1/p50 and TCF7L2 antibodies produced specific bands with nucleic fractions only. Corresponding to mRNA expression of 5.3 RQ fold (ddCt value -2.4 for LPS sample) and 8.6 RQ fold (ddCt value -3.1 for PMA sample), protein expression of TCF7L2 was clearly induced (up to 5 fold for PMA sample) while LPS sample shows moderate induction of less than 2 fold. Prominent induction of NFkB1/p50 protein after LPS stimulation is also in line with its classical role during inflammatory responses. With all antibodies tested by ELISA, fluorescence detection of antibody binding to protein substrate was observable only with nucleic fraction (data not shown).
U937 cells as valid model system
The relatively small number of genes up and down regulated after VitD3 differentiation is a strong indication that U937 cells used in this study are nearly proximal to functional monocytes (Figure 2). Despite both CD14 and CD11b are induced in monocyte-derived macrophages, U937 cells-derived macrophages manifest up regulation of CD11b alone (Table 2 and Figure 4). However, lack of CD14 induction should not impair the essential function representing macrophages. CD14 binds LPS to facilitate recognition by TLR4, but does not participate directly in signaling. TLR4 alone can function as a singular signaling receptor for LPS. Moreover, LPS signaling in U937 cells can bypass the necessity of CD14 component of receptor complex due to high concentration of LPS (2 μg/ml or 10 μg/ml) used in this study. The abundance of TNFalpha mRNA prior to differentiation was higher with primary monocytes (data not shown) when compared to U937 cells as manifested by initial Ct value. Taken together, the lack of CD14 expression and the induction of TNFalpha in PMA-differentiated U937 cells can be attributed to the effect of the artificial differentiation system. The MyD88-independent pathway is responsible for sustenance of the proinflammatory program with delayed kinetics of NFkB activation, in the meantime the MyD88-dependent pathway triggers response to initial challenge [22, 23]. It is likely that the biphasic behavior of TNFalpha production (up regulation at early and late time points) observed with macrophage-like cells (Figure 1) reflects cross talks and temporal segregation between these two pathways. Up regulation of RB together with p21 upon PMA differentiation is in line with proliferation arrest prior to differentiation process of cells. Both JAK3 and STAT3 (Additional file 3) are known to play a vital role in monocytic differentiation [24, 25]. JAK3 (tyrosine kinase) mediates STAT signaling through direct tyrosine phosphorylation of STAT proteins followed by STAT dimerization. Down regulation of MPO, an effector gene of neutrophil function is certainly indicative of macrophage specific gene repression. The presence of Nur77 in the category of LPS inducible genes (Figure 3) is directly corresponding to previous reports in which the role of NR4A orphan receptors (Nur77, Nurr1, Nor1) was described as key regulators of chronic inflammatory diseases thus in the innate immune response [26, 27]. Taken together, U937 cells serve a good model system with propensity of differentiation and activation pathways intrinsic to primary monocytes and macrophages.
Transcriptional regulators involved in macrophages differentiation and activation
Gene expression analysis of transcriptional regulators in both U937 cells and primary macrophages revealed several potentially novel factors. MEF2 family comprises 4 members (MEF2A, B, C, D) that are encoded by different genes. They are described as key regulators in muscle cells and related cell types as in the case of cardiac development . Beside from T cell selection and function [29–31], they are fairly unknown in the context of the immune system. Previously, MEF2C was proved to be activated by p38- mediated phosphorylation to induce transcription of JUN which in turn regulates cytokine gene expression during inflammatory response . More recently, it has been suggested that the MEF2 family is a versatile regulator in diverse cellular systems [33–35]. However, no significant enrichment of MEF2 binding site was seen by promoter analysis of differentially expressed genes (data not shown). Possibly, not all MEF2 binding sites are yet known to be available for complete analysis, thus prevent the final conclusion. On the other hand, MEF2 might act also through differential protein modification and indirect binding to modulate gene expression. Specificity of MEF2 and co-regulator complexes is likely to confer differentiation and immune regulatory function in macrophages. Interestingly, HDAC5 and HDAC9, histone deacetylases, known to form repression complexes with MEF2 family members in many cellular systems are concomitantly up regulated. Both histone deacetylases belong to class II histone deacetylase known to be involved in phosphorylation dependent nucleo – cytoplasmic trafficking in response to activation signals and the formation of repression complexes [36–39]. During muscle development, expression of MEF2 and HDAC9 was shown to be increased during initial phase but diminished expression of HDAC9 at late stage . A likewise operation during macrophage activation mechanism by which MEF2 up regulates its own repressor HDAC9 can be suggested. Among other factors up regulated in both U937 cells and primary macrophages are HLX1 and BRI. HLX1 was described as a marker of immature hematopoietic cells and also known to be involved in the activation of T lymphocyte or NK cells [41, 42]. BRI was identified as a LPS-inducible gene in murine macrophages in the presence of CSF-1 . However, there has been no report correlating expression and function of HLX1 and BRI in macrophages. Microarray demonstrated their constant up regulation across different time points of PMA differentiation (Table 2). Additionally, RT-PCR result shows constant up regulation of BRI across LPS stimulation period in U937 cells (data not shown). Taken together, it can be strongly suggested they play a role as effecter molecules of transcriptional reprogramming in macrophages.
Inconsistent expression pattern and heterogeneity of transcription factors
There are several factors manifesting inconsistent expression patterns between U937 and primary cell systems. TCF7L2 is one of them. On the contrary to PMA/U937 cells with which TCF7L2 induction of both mRNA and protein were verified (Figure 5), M-CSF/macrophages show no induction of TCF7L2 mRNA. Keeping in mind that it is known as the hottest susceptibility gene for type 2 diabetes while being involved in glucose homeostasis, individual heterogeneity of TCF7L2 expression observed with M-CSF/macrophages in this study might be informative. During the past decade, many evidences pointed out type 2 diabetes as a disease of innate immune system, particularly of macrophages, thus closely link the metabolic syndrome to inflammatory pathway [44–46]. Such an inter-relationship promises high probability of TCF7L2 acting in macrophage biology. Another particularly interesting example of inconsistency is up regulation of KLF2 observed with PMA/U937 cells on the contrary to M-CSF/macrophages that shows down regulation. KLF2 was previously reported to be down regulated upon monocyte differentiation or activation in THP-1 cells, and has been proved to be a negative regulator of proinflammatory activation of monocytes . The basis for this difference between U937 and THP-1 cells may be due to the different origin and maturation stage of cells. U937 cells are of tissue origin (histocytic lymphoma) thus at more mature stage. THP1 cells are of blood leukemic origin at less mature stage. Importantly, up regulation rather than down regulation of KLF2 did not produce the consequences of abrogating inflammatory gene activation in U937 cells. Based on the observation of Kumar and colleagues in which p65 and HDACs repress MEF2 transcriptional activation of KLF2 promoter in endothelial cells, it might be reasonable to assume that altered activity balance of antagonistic factors contribute to different cellular activation . The roles of KLF4 and IRF7 in the differentiation or function of macrophages have been also described previously [18, 48, 49], thus down regulation rather than up regulation with primary macrophages doesn't seem to be immediately accountable. One possible explanation is that an initial amount of protein available in different cell type determine up or down regulation to adjust correct cellular dose. Another explanation is that up regulation can be the effect of U937 cells that are not equal to primary monocytes in the first place. Taken together, the result emphasizes a necessary caution to be exercised when interpreting data obtained from a transformed cell line, alluding differential transcriptional regulatory and signaling pathways between two systems. Nevertheless, the heterogeneity of human individuals observed with several factors (33%; n = 8) including ID3, IRF9, IRF7, TCF7L2, HDAC7/9, OCT2, PCAF and MEF2B raises a possibility that functional capacities of individual macrophages might determine differential susceptibility to and severity of various immune diseases. This enforce the strong argument that investigations using primary cells from different blood donors can greatly impede reliable conclusion of result, thus further emphasize the necessity of a unified model system.
Potential role of chromatin architecture in macrophage activation
Being surmised that increased expression of genes upon differentiation may be already indicative of their serving important roles in macrophage function afterward, prompt induction of additional 278 genes within 2 hrs of LPS stimulation is intriguing (Figure 2). It might have an implication that promoters of these genes acquired poised status ready to be quickly remodeled upon LPS signal. Accordingly, increased expression of variant histone H2AV, which is also known as H2AZ-2, suggests an exciting possibility. As a dynamic regulator relying on its deposit and loss, incorporation of H2AV in exchange with H2A (canonical histone) could serve as a stable epigenetic mark to keep quiescent promoters in a repressed state while creating more permissive chromatin architecture until appropriate activation signal is received [50–53]. Both forms of protein, H2AZ-1 (previously H2AZ) and H2AZ-2 (previously H2AV), may play new or complementary functions . There are probably several types of chromatin remodeling in general. There may be a large scale remodeling early in development and then more regional remodeling as sections of chromatin become active or are shut down. A more localized type of remodeling takes previously prepared regions and opens them up completely upon receipt of more immediate signals for active transcription. Each type may shares some protein factors and also has unique factor involvement. It is tempting to speculate that identified transcription factors are operating to repress and induce target genes by directing decondensed chromatin structure (H2AV incorporation) during macrophage development which proceed rapid nucleosome remodeling (H2AV eviction) upon LPS stimulation . One possible speculation is whether it functions to activate or silence depend upon the nature of the factors recruited to H2AZ containing nucleosomes while turning off a set of genes and turning on a new subset during the differentiation and inflammatory process.
Identification of novel transcriptional regulators
The focus of this study was the identification of transcriptional regulators during macrophage differentiation and activation using U937 cells as a model system. Several novel transcription factors (i.e. MEF2 family, BRI, HLX1, TCF7L2, ID3, NFIL3) were found to be up regulated in addition to those already known (i.e. MAFB, EGR, BCL6, IRF7, NFkB, JUN, FOS, STAT, Nur77) to be involved during macrophage differentiation and/or inflammatory response. We report for the first time both mRNA and protein induction of MEF2 proteins, a MADS-box transcription factor family, during macrophage differentiation and activation in U937 model system. In both U937 cells and primary macrophage cells, three different isoforms of MEF2 (MEF2A, 2C, 2D) appear to be concomitantly induced, thus suggesting functional redundancy of different members. Additionally, MEF2 is directly linking gene activation with epigenetic gene regulation such as histone acetylation and chromatin remodeling . This hypothesis is supported by the fact that interaction partners known to modulate histone modification (HDAC 5 and 9) and chromatin structure (H2AV) are concomitantly expressed.
Validity of U937 cells as a model system
U937 cells used in this study are proved to be very close to functional monocytes as indicated by the low number of differently expressed genes upon VitD3 differentiation, thus can be an ideal model system to investigate gene regulatory mechanisms in differentiating macrophages. Array analysis of PMA differentiated and LPS activated cells revealed a variety of genes intrinsic to macrophages, inflammation and innate immunity. Nevertheless, how faithfully do these cells resemble primary macrophages is a question to be challenged with carefulness because differentiation pathways in U937 cells-derived macrophages and primary macrophages may diverges to some degree considering the malignant origin of U937 cells. This was asserted by cross validation of 24 transcriptional regulators, which proved a 54.2% (n = 13) of tested factors with consistent expression pattern between PMA differentiated U937 cells and M-CSF differentiated monocytes. The observed differences are based on different macrophage maturation thus highlight the necessity of cross validation. Individual heterogeneity of blood donors further emphasizes the advantage of using a unified model system. Taken together, target genes under investigation should be selected carefully keeping in mind that cells are characterized by how they respond to various stimuli as defined by gene expression. This allows PMA differentiated U937 cells to act as representative of macrophages.
As a next step, the role of new identified genes should be verified through functional studies such as RNAi knock down. A further task would be to elucidate how newly identified transcription regulators communicate with factors that modify and remodel chromatin structure of target gene promoters.
U937 cell culture, differentiation, stimulation
U937 cells were maintained in RPMI 1640 medium (Gibco BRL) supplemented with 10% FBS (Biochrome) under humidified air with 5% CO2 at 37°C. Cell viability was estimated regularly by trypan blue dye exclusion. To induce differentiation into monocytes or adherent macrophages, the cells were seeded at an initial density of 2 × 105 cells/ml of VLE-RPMI (Biochrome) and cultured in the presence of 100 nM VitD3 (Calbiochem) or 10 nM PMA (Sigma) for up to 48 hrs at 37°C. Differentiation was monitored by FACS analysis for positive staining of cell surface markers, CD14 and CD11b. Briefly, control and differentiated cells were collected, washed twice in PBS at 4°C, and incubated with 1 μg/ml anti-human FITC-conjugated anti-CD11b and PE-conjugated anti-CD14 (Miltenyi Biotech. Inc.) for 10 min at 4°C and washed once before being resuspended in FACS buffer. For stimulation, differentiated cells were washed, resuspended in VLE-RPMI and stimulated by adding 2 or 10 μg/ml of LPS (Sigma).
Generation of macrophages from primary blood monocytes
Human white blood cells were separated from blood buffy coat of healthy donors by Ficoll-Paque density centrifugation. To isolate CD14+ monocytes, selection was performed using immunomagnetic particles on AutoMACS (Miltenyi Biotech, Germany) according to manufacturer's instruction. Purity of monocyte suspension was confirmed by cell counting (Sysmex-F 820, Germany), and only preparation proven to have a level of contaminating platelets below the detection limit < 1 × 103/μl was selected. To generate monocyte derived macrophages, 1 × 105 CD14+ cells were expanded in 96-well flat bottom plates containing 0.2 ml of RPMI 1640 medium with L-glutamine, 10% FCS (PAA laboratories, Austria) and 0.5% penicillin/streptomycin (10000 U/ml, Bayer, Germany). As a differentiation inducer, 50 ng/ml macrophage-colony stimulating factor (M-CSF, Promokine, Germany) was added and cells were cultured for 4 days at 37°C, 5% CO2 atmosphere. LPS stimulation was done at the final concentration of 10 μg/ml.
RNA isolation, RT-PCR and real-time PCR
Total RNA was extracted from frozen cell pellet by using RNeasy Plus Mini kit (Qiagen) according to the manufacturer's protocol. RNA concentration and purity was measured using NanoDrop 1000A Spectrophotometer. 1 μg of total RNA from each sample was reverse-transcribed to produce cDNA using ImProm-II RT-system (Promega). SyBr-Green (ABgene) real-time PCR was performed with ABI PRISM 7700 to generate 100 ~200 bps of product. All samples were run in duplicate or triplicate and dissociation curve was generated after each run to control formation of primer dimers. GAPDH or betaActin was used as RNA loading control. Normalization of each data point was carried out, in which dCt = Cttest gene - CtGAPDH; ddCt = dCt(test sample) - dCt(U937 control). Primer sequences are listed in Additional file 4.
cRNA Hybridization Array
To ensure quantitative accuracy and reproducibility of expression analysis, two independent rounds of differentiation and stimulation experiment were performed. For each experiment, two biological replicates derived from separate culture plates were run in parallel, and represented onto two separate BeadChip slides. Upon RNA isolation and synthesis of labeled cRNA, an additional sample was represented onto a separate slide thus satisfying a technical replicate. Prior to cRNA hybridization, RNA integrity was evaluated by microfluidics analysis using the Agilent 2100 bioanalyzer. cRNA synthesis and labeling was performed using Illumina TotalPrep RNA Amplification Kit (Ambion) then subsequently hybridized, according to the instruction of manufacturer, onto Illumina Sentrix BeadChip Array, HumanRef-8 containing 24,000 features of oligonucleotide probe representing majority of known human genes. Arrays were scanned using Illumina BeadArray reader (V18.104.22.168). Our experimental design adhered to MIAME (Minimum Information About Microarray Experiments) guidelines.
Differential gene expression analysis
Raw data was acquired as signal intensity of each spot on the array, and normalization was done by variance stabilization method . Signal values of technical or biological replicates were averaged beforehand. Differentially expressed genes were determined based on a comparison of normalized signal intensity ratio between experimental (differentiated and/or stimulated U937 cells) and control (undifferentiated U937 cells) samples, which was expressed as fold change on a spread sheet tables.
Protein localization and quantification
Protein lysates were prepared using Cytosolic/Nucleic fraction kit (BioCat), and concentrations were determined by Bradford (BioRad) assay. 5 to 10 μg of proteins were run on denaturing 10% polyacrylamide gel after 5 min heating in Laemmli buffer, transferred onto polyvinylidene fluoride (PDVF) membranes (Millipore). After blocking in TBST with 1% BSA for one hour at RT, membrane was incubated overnight at 4°C with primary antibodies (1: 500 or 1: 1,000 dilution). Incubation with horseradish peroxidase conjugated anti-rabbit secondary AB (1: 2,000) was done for 1 hr (RT). Equal amount of protein loading was confirmed by betaActin. After ECL detection, signal intensity of specific band was analyzed by AIDA image analysis program. The NFkB p50/105 polyclonal (H-119) and MEF2 polyclonal antibodies (C-21) were purchased from Santa Cruz Biotech. Monoclonal antibody to TCF7L2 (C9B9, #2565) was from Cell Signalling.
Janeway CA, Travers P, Walport M, Shlomchik MJ: Immunobiology, the immune system in health and disease. 2005, Garland Science Publishing/New York
Fujihara M, Muroi M, Tanamoto K, Suzuki T, Azuma H, Ikeda H: Molecular mechanisms of macrophage activation and deactivation by lipopolysaccharide: roles of the receptor complex. Pharmacol Ther. 2003, 100 (2): 171-194. 10.1016/j.pharmthera.2003.08.003.
Dion MF, Kaplan T, Kim M, Buratowski S, Friedman N, Rando OJ: Dynamics of replication-independent histone turnover in budding yeast. Science. 2007, 315 (5817): 1405-1408. 10.1126/science.1134053.
Valledor AF, Borras FE, Cullell-Young M, Celada A: Transcription factors that regulate monocyte/macrophage differentiation. J Leukoc Biol. 1998, 63 (4): 405-417.
Toney LM, Cattoretti G, Graf JA, Merghoub T, Pandolfi PP, Dalla-Favera R, Ye BH, Dent AL: BCL-6 regulates chemokine gene transcription in macrophages. Nat Immunol. 2000, 1 (3): 214-220. 10.1038/79749.
Park GY, Joo M, Pedchenko T, Blackwell TS, Christman JW: Regulation of macrophage cyclooxygenase-2 gene expression by modifications of histone H3. Am J Physiol Lung Cell Mol Physiol. 2004, 286 (5): L956-962. 10.1152/ajplung.00338.2003.
Aung HT, Schroder K, Himes SR, Brion K, van Zuylen W, Trieu A, Suzuki H, Hayashizaki Y, Hume DA, Sweet MJ, et al.: LPS regulates proinflammatory gene expression in macrophages by altering histone deacetylase expression. FASEB J. 2006, 20 (9): 1315-1327. 10.1096/fj.05-5360com.
Foster SL, Hargreaves DC, Medzhitov R: Gene-specific control of inflammation by TLR-induced chromatin modifications. Nature. 2007, 447 (7147): 972-978.
Smale ST, Fisher AG: Chromatin structure and gene regulation in the immune system. Annu Rev Immunol. 2002, 20: 427-462. 10.1146/annurev.immunol.20.100301.064739.
Tsuchiya S, Yamabe M, Yamaguchi Y, Kobayashi Y, Konno T, Tada K: Establishment and characterization of a human acute monocytic leukemia cell line (THP-1). Int J Cancer. 1980, 26 (2): 171-176. 10.1002/ijc.2910260208.
Sundstrom C, Nilsson K: Establishment and characterization of a human histiocytic lymphoma cell line (U-937). Int J Cancer. 1976, 17 (5): 565-577. 10.1002/ijc.2910170504.
Liu HZ, Gong JP, Wu CX, Peng Y, Li XH, You HB: The U937 cell line induced to express CD14 protein by 1,25-dihydroxyvitamin D3 and be sensitive to endotoxin stimulation. Hepatobiliary Pancreat Dis Int. 2005, 4 (1): 84-89.
Woronicz JD, Lina A, Calnan BJ, Szychowski S, Cheng L, Winoto A: Regulation of the Nur77 orphan steroid receptor in activation-induced apoptosis. Mol Cell Biol. 1995, 15 (11): 6364-6376.
Youn HD, Liu JO: Cabin1 represses MEF2-dependent Nur77 expression and T cell apoptosis by controlling association of histone deacetylases and acetylases with MEF2. Immunity. 2000, 13 (1): 85-94. 10.1016/S1074-7613(00)00010-8.
Shalizi A, Gaudilliere B, Yuan Z, Stegmuller J, Shirogane T, Ge Q, Tan Y, Schulman B, Harper JW, Bonni A: A calcium-regulated MEF2 sumoylation switch controls postsynaptic differentiation. Science. 2006, 311 (5763): 1012-1017. 10.1126/science.1122513.
Ito K, Ito M, Elliott WM, Cosio B, Caramori G, Kon OM, Barczyk A, Hayashi S, Adcock IM, Hogg JC, et al.: Decreased histone deacetylase activity in chronic obstructive pulmonary disease. N Engl J Med. 2005, 352 (19): 1967-1976. 10.1056/NEJMoa041892.
Das H, Kumar A, Lin Z, Patino WD, Hwang PM, Feinberg MW, Majumder PK, Jain MK: Kruppel-like factor 2 (KLF2) regulates proinflammatory activation of monocytes. Proc Natl Acad Sci USA. 2006, 103 (17): 6653-6658. 10.1073/pnas.0508235103.
Feinberg MW, Wara AK, Cao Z, Lebedeva MA, Rosenbauer F, Iwasaki H, Hirai H, Katz JP, Haspel RL, Gray S, et al.: The Kruppel-like factor KLF4 is a critical regulator of monocyte differentiation. EMBO J. 2007, 26 (18): 4138-4148. 10.1038/sj.emboj.7601824.
Hayden MS, West AP, Ghosh S: NF-kappaB and the immune response. Oncogene. 2006, 25 (51): 6758-6780. 10.1038/sj.onc.1209943.
Andreesen R, Brugger W, Scheibenbogen C, Kreutz M, Leser HG, Rehm A, Lohr GW: Surface phenotype analysis of human monocyte to macrophage maturation. J Leukoc Biol. 1990, 47 (6): 490-497.
Akagawa KS: Functional heterogeneity of colony-stimulating factor-induced human monocyte-derived macrophages. Int J Hematol. 2002, 76 (1): 27-34. 10.1007/BF02982715.
Kawai T, Adachi O, Ogawa T, Takeda K, Akira S: Unresponsiveness of MyD88-deficient mice to endotoxin. Immunity. 1999, 11 (1): 115-122. 10.1016/S1074-7613(00)80086-2.
Covert MW, Leung TH, Gaston JE, Baltimore D: Achieving stability of lipopolysaccharide-induced NF-kappaB activation. Science. 2005, 309 (5742): 1854-1857. 10.1126/science.1112304.
Mangan JK, Rane SG, Kang AD, Amanullah A, Wong BC, Reddy EP: Mechanisms associated with IL-6-induced up-regulation of Jak3 and its role in monocytic differentiation. Blood. 2004, 103 (11): 4093-4101. 10.1182/blood-2003-06-2165.
Nakajima K, Yamanaka Y, Nakae K, Kojima H, Ichiba M, Kiuchi N, Kitaoka T, Fukada T, Hibi M, Hirano T: A central role for Stat3 in IL-6-induced regulation of growth and differentiation in M1 leukemia cells. EMBO J. 1996, 15 (14): 3651-3658.
Pei L, Castrillo A, Tontonoz P: Regulation of macrophage inflammatory gene expression by the orphan nuclear receptor Nur77. Mol Endocrinol. 2006, 20 (4): 786-794. 10.1210/me.2005-0331.
Pei L, Castrillo A, Chen M, Hoffmann A, Tontonoz P: Induction of NR4A orphan nuclear receptor expression in macrophages in response to inflammatory stimuli. J Biol Chem. 2005, 280 (32): 29256-29262. 10.1074/jbc.M502606200.
Black BL, Olson EN: Transcriptional control of muscle development by myocyte enhancer factor-2 (MEF2) proteins. Annu Rev Cell Dev Biol. 1998, 14: 167-196. 10.1146/annurev.cellbio.14.1.167.
Pan F, Ye Z, Cheng L, Liu JO: Myocyte enhancer factor 2 mediates calcium-dependent transcription of the interleukin-2 gene in T lymphocytes: a calcium signaling module that is distinct from but collaborates with the nuclear factor of activated T cells (NFAT). J Biol Chem. 2004, 279 (15): 14477-14480. 10.1074/jbc.C300487200.
Esau C, Boes M, Youn HD, Tatterson L, Liu JO, Chen J: Deletion of calcineurin and myocyte enhancer factor 2 (MEF2) binding domain of Cabin1 results in enhanced cytokine gene expression in T cells. J Exp Med. 2001, 194 (10): 1449-1459. 10.1084/jem.194.10.1449.
Youn HD, Sun L, Prywes R, Liu JO: Apoptosis of T cells mediated by Ca2+-induced release of the transcription factor MEF2. Science. 1999, 286 (5440): 790-793. 10.1126/science.286.5440.790.
Han J, Jiang Y, Li Z, Kravchenko VV, Ulevitch RJ: Activation of the transcription factor MEF2C by the MAP kinase p38 in inflammation. Nature. 1997, 386 (6622): 296-299. 10.1038/386296a0.
Arnold MA, Kim Y, Czubryt MP, Phan D, McAnally J, Qi X, Shelton JM, Richardson JA, Bassel-Duby R, Olson EN: MEF2C transcription factor controls chondrocyte hypertrophy and bone development. Dev Cell. 2007, 12 (3): 377-389. 10.1016/j.devcel.2007.02.004.
Verzi MP, Agarwal P, Brown C, McCulley DJ, Schwarz JJ, Black BL: The transcription factor MEF2C is required for craniofacial development. Dev Cell. 2007, 12 (4): 645-652. 10.1016/j.devcel.2007.03.007.
Lam BY, Chawla S: MEF2D expression increases during neuronal differentiation of neural progenitor cells and correlates with neurite length. Neurosci Lett. 2007, 427 (3): 153-158. 10.1016/j.neulet.2007.09.030.
McKinsey TA, Zhang CL, Lu J, Olson EN: Signal-dependent nuclear export of a histone deacetylase regulates muscle differentiation. Nature. 2000, 408 (6808): 106-111. 10.1038/35040593.
Zhang CL, McKinsey TA, Chang S, Antos CL, Hill JA, Olson EN: Class II histone deacetylases act as signal-responsive repressors of cardiac hypertrophy. Cell. 2002, 110 (4): 479-488. 10.1016/S0092-8674(02)00861-9.
Vega RB, Harrison BC, Meadows E, Roberts CR, Papst PJ, Olson EN, McKinsey TA: Protein kinases C and D mediate agonist-dependent cardiac hypertrophy through nuclear export of histone deacetylase 5. Mol Cell Biol. 2004, 24 (19): 8374-8385. 10.1128/MCB.24.19.8374-8385.2004.
Dequiedt F, Van Lint J, Lecomte E, Van Duppen V, Seufferlein T, Vandenheede JR, Wattiez R, Kettmann R: Phosphorylation of histone deacetylase 7 by protein kinase D mediates T cell receptor-induced Nur77 expression and apoptosis. J Exp Med. 2005, 201 (5): 793-804. 10.1084/jem.20042034.
Haberland M, Arnold MA, McAnally J, Phan D, Kim Y, Olson EN: Regulation of HDAC9 gene expression by MEF2 establishes a negative-feedback loop in the transcriptional circuitry of muscle differentiation. Mol Cell Biol. 2007, 27 (2): 518-525. 10.1128/MCB.01415-06.
Deguchi Y, Kirschenbaum A, Kehrl JH: A diverged homeobox gene is involved in the proliferation and lineage commitment of human hematopoietic progenitors and highly expressed in acute myelogenous leukemia. Blood. 1992, 79 (11): 2841-2848.
Becknell B, Hughes TL, Freud AG, Blaser BW, Yu J, Trotta R, Mao HC, Caligiuri de Jesus ML, Alghothani M, Benson DM, et al.: Hlx homeobox transcription factor negatively regulates interferon-gamma production in monokine-activated natural killer cells. Blood. 2007, 109 (6): 2481-2487. 10.1182/blood-2006-10-050096.
Sester DP, Trieu A, Brion K, Schroder K, Ravasi T, Robinson JA, McDonald RC, Ripoll V, Wells CA, Suzuki H, et al.: LPS regulates a set of genes in primary murine macrophages by antagonising CSF-1 action. Immunobiology. 2005, 210 (2–4): 97-107. 10.1016/j.imbio.2005.05.004.
Yi F, Brubaker PL, Jin T: TCF-4 mediates cell type-specific regulation of proglucagon gene expression by beta-catenin and glycogen synthase kinase-3beta. J Biol Chem. 2005, 280 (2): 1457-1464. 10.1074/jbc.M411487200.
Hattersley AT: Prime suspect: the TCF7L2 gene and type 2 diabetes risk. J Clin Invest. 2007, 117 (8): 2077-2079. 10.1172/JCI33077.
Fernandez-Real JM, Pickup JC: Innate immunity, insulin resistance and type 2 diabetes. Trends Endocrinol Metab. 2008, 19 (1): 10-16. 10.1016/j.tem.2007.10.004.
Kumar A, Lin Z, SenBanerjee S, Jain MK: Tumor necrosis factor alpha-mediated reduction of KLF2 is due to inhibition of MEF2 by NF-kappaB and histone deacetylases. Mol Cell Biol. 2005, 25 (14): 5893-5903. 10.1128/MCB.25.14.5893-5903.2005.
Kawai T, Sato S, Ishii KJ, Coban C, Hemmi H, Yamamoto M, Terai K, Matsuda M, Inoue J, Uematsu S, et al.: Interferon-alpha induction through Toll-like receptors involves a direct interaction of IRF7 with MyD88 and TRAF6. Nat Immunol. 2004, 5 (10): 1061-1068. 10.1038/ni1118.
Feinberg MW, Cao Z, Wara AK, Lebedeva MA, Senbanerjee S, Jain MK: Kruppel-like factor 4 is a mediator of proinflammatory signaling in macrophages. J Biol Chem. 2005, 280 (46): 38247-38258. 10.1074/jbc.M509378200.
Zhang H, Roberts DN, Cairns BR: Genome-wide dynamics of Htz1, a histone H2A variant that poises repressed/basal promoters for activation through histone loss. Cell. 2005, 123 (2): 219-231. 10.1016/j.cell.2005.08.036.
Farris SD, Rubio ED, Moon JJ, Gombert WM, Nelson BH, Krumm A: Transcription-induced chromatin remodeling at the c-myc gene involves the local exchange of histone H2A.Z. J Biol Chem. 2005, 280 (26): 25298-25303. 10.1074/jbc.M501784200.
Li B, Pattenden SG, Lee D, Gutierrez J, Chen J, Seidel C, Gerton J, Workman JL: Preferential occupancy of histone variant H2AZ at inactive promoters influences local histone modifications and chromatin remodeling. Proc Natl Acad Sci USA. 2005, 102 (51): 18385-18390. 10.1073/pnas.0507975102.
Raisner RM, Hartley PD, Meneghini MD, Bao MZ, Liu CL, Schreiber SL, Rando OJ, Madhani HD: Histone variant H2A.Z marks the 5' ends of both active and inactive genes in euchromatin. Cell. 2005, 123 (2): 233-248. 10.1016/j.cell.2005.10.002.
Eirin-Lopez JM, Gonzalez-Romero R, Dryhurst D, Ishibashi T, Ausio J: The evolutionary differentiation of two histone H2A.Z variants in chordates (H2A.Z-1 and H2A.Z-2) is mediated by a stepwise mutation process that affects three amino acid residues. BMC evolutionary biology. 2009, 9: 31-10.1186/1471-2148-9-31.
Gevry N, Chan HM, Laflamme L, Livingston DM, Gaudreau L: p21 transcription is regulated by differential localization of histone H2A.Z. Genes Dev. 2007, 21 (15): 1869-1881. 10.1101/gad.1545707.
Ohkawa Y, Yoshimura S, Higashi C, Marfella CG, Dacwag CS, Tachibana T, Imbalzano AN: Myogenin and the SWI/SNF ATPase Brg1 maintain myogenic gene expression at different stages of skeletal myogenesis. The Journal of biological chemistry. 2007, 282 (9): 6564-6570. 10.1074/jbc.M608898200.
Huber W, von Heydebreck A, Sultmann H, Poustka A, Vingron M: Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics. 2002, 18 (Suppl 1): S96-104.
This work was supported by the Max Planck Society. We would like to thank Nils Rademacher for help with FACS analysis, Pamela Kepper for help with Illumina hybridization and documentation, Christina Lang for M-CSF/macrophage differentiation. We are also grateful to Dr. Ulrike Dorman, Dr. Wasco Wruck, Dr. Andreas Dahl, Wibke Schwarzer, Mario Drungowski and Daniela Koester for valuable discussions and help.
YB and HS were primarily responsible for the design, coordination, conduct of the study and drafted the manuscript. SH were responsible for array normalization and bioinformatic analysis. HH and GB provided U937 cells and primary macrophages. MHS contributed to experimental design of U937 cell system. SS and HL contributed editorial revisions, provided comments and discussion. All authors read and approved the final manuscript.