Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

Differential RelA- and RelB-dependent gene transcription in LTβR-stimulated mouse embryonic fibroblasts

Agnes Lovas1, Dörte Radke256, Daniela Albrecht3, Z Buket Yilmaz14, Ulrich Möller2, Andreas JR Habenicht5 and Falk Weih1*

Author Affiliations

1 Research Group Immunology, Leibniz Institute for Age Research, Fritz Lipmann Institute, Beutenbergstrasse 11, 07745 Jena, Germany

2 Bioinformatics – Pattern Recognition, Leibniz Institute for Natural Product Research and Infection Biology, Hans-Knöll-Institute, Beutenbergstrasse 11a, 07745 Jena, Germany

3 Molecular and Applied Microbiology, Leibniz Institute for Natural Product Research and Infection Biology, Hans-Knöll-Institute, Beutenbergstrasse 11a, 07745 Jena, Germany

4 Signal Transduction in Tumor Cells, Max Delbrück Center for Molecular Medicine, Robert Rössle Strasse 10, 13092 Berlin-Buch, Germany

5 Institute for Vascular Medicine, Friedrich Schiller University of Jena, Bachstrasse 18, 07743 Jena, Germany

6 Institute for Community Medicine, Ernst Moritz Arndt University Greifswald, Walther Rathenau Strasse 48, 17475 Greifswald, Germany

For all author emails, please log on.

BMC Genomics 2008, 9:606  doi:10.1186/1471-2164-9-606

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/9/606


Received:3 July 2008
Accepted:16 December 2008
Published:16 December 2008

© 2008 Lovas et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Lymphotoxin signaling via the lymphotoxin-β receptor (LTβR) has been implicated in biological processes ranging from development of secondary lymphoid organs, maintenance of spleen architecture, host defense against pathogens, autoimmunity, and lipid homeostasis. The major transcription factor that is activated by LTβR crosslinking is NF-κB. Two signaling pathways have been described, the classical inhibitor of NF-κB α (IκBα)-regulated and the alternative p100-regulated pathway that result in the activation of p50-RelA and p52-RelB NF-κB heterodimers, respectively.

Results

Using microarray analysis, we investigated the transcriptional response downstream of the LTβR in mouse embryonic fibroblasts (MEFs) and its regulation by the RelA and RelB subunits of NF-κB. We describe novel LTβR-responsive genes that were regulated by RelA and/or RelB. The majority of LTβR-regulated genes required the presence of both RelA and RelB, revealing significant crosstalk between the two NF-κB activation pathways. Gene Ontology (GO) analysis confirmed that LTβR-NF-κB target genes are predominantly involved in the regulation of immune responses. However, other biological processes, such as apoptosis/cell death, cell cycle, angiogenesis, and taxis were also regulated by LTβR signaling. Moreover, LTβR activation inhibited expression of a key adipogenic transcription factor, peroxisome proliferator activated receptor-γ (pparg), suggesting that LTβR signaling may interfere with adipogenic differentiation.

Conclusion

Microarray analysis of LTβR-stimulated fibroblasts provided comprehensive insight into the transcriptional response of LTβR signaling and its regulation by the NF-κB family members RelA and RelB.

Background

NF-κB transcription factors are essential for innate and adaptive immunity, cell survival, cellular stress responses, development and maintenance of lymphoid organ structures, and other biological functions [1-3]. The vertebrate NF-κB family includes five structurally related members, the Rel proteins RelA (p65), RelB, cRel, and the NF-κB proteins p50 and p52. Among the Rel/NF-κB family, only RelA, RelB, and cRel contain C-terminal transcriptional activation domains enabling them to directly regulate transcription. The other two members, p50 and p52, are synthesized as p105 and p100 precursors, respectively. The Rel and NF-κB proteins can form different homo- and heterodimers (for example p50-RelA or p52-RelB) that bind to DNA target sites, so-called κB sites. In resting cells, Rel/NF-κB proteins associate with inhibitory κB molecules (IκBs) and are retained in the cytoplasm as inactive forms [4].

Two major NF-κB signaling pathways can be distinguished, the classical or canonical and the alternative or non-canonical pathway. In response to stimulation of transmembrane receptors like tumor necrosis factor receptor (TNFR)-1 or Toll-like receptor (TLR)-4, signaling cascades are initiated that lead to the liberation of Rel/NF-κB complexes from their IκB molecules. As a result, they translocate to the nucleus and regulate transcription of numerous target genes. This classical pathway involves phosphorylation of IκBα by the NEMO (NF-κB essential modulator)/IKKγ- and IKKβ-containing IκB kinase (IKK) complex followed by its ubiquitin-dependent proteasomal degradation. Regulation of gene transcription is predominantly mediated through p50-RelA and p50-cRel heterodimers and target genes are mainly involved in innate immunity, cell survival, and inflammation. A few inducers of NF-κB, like LTβR, are able to trigger an additional, so-called alternative or non-canonical pathway through the activation of NF-κB-inducing kinase (NIK) and IKKα. The alternative pathway governs gene regulation mainly through p52-RelB heterodimers that are generated from the inactive cytoplasmic p100-RelB complex via signal-dependent processing of the p100 inhibitor to p52. This pathway controls genes that are predominantly involved in adaptive immunity and lymphoid organ development [5-11]. Recent findings by Hoffmann and colleagues extend this scenario. They could show that not only RelB- but also RelA-containing complexes can be released from the p100 inhibitor after LTβR stimulation [12-14].

This report focuses on the transcriptional response downstream of the LTβR and its regulation by RelA and RelB. The role of LTβR signaling in development and organization of secondary lymphoid structures is well documented (reviewed in [8,15-17]). We are interested in similarities and differences in RelA and RelB function in lymphoid organ development. However, a major problem is that RelA-deficient (relA-/-) mice are embryonic lethal due to tumor necrosis factor (TNF)-induced hepatocyte apoptosis [18]. Moreover, RelB-deficient (relB-/-) mice display impaired secondary lymphoid organ development and suffer from an autoinflammatory syndrome that also affects organization and function of lymphoid tissues [19,20]. Thus, stromal compartments that display LTβR signaling and thereby have an organizational role in the development of lymphoid organs cannot be used for in vivo gene expression studies from the above animals. Therefore, we applied MEFs established from wild-type (wt), relA-/-, and relB-/- mice as an in vitro model system. Also, there is increasing evidence that LTβR functions beyond lymphoid organs, as it is involved in liver regeneration, hepatitis [21], and hepatic lipid metabolism [22]. We therefore hypothesized that LTβR signaling, via RelA and/or RelB, may participate in physiological processes other than lymphorganogenesis. MEFs with different genotypes (wt, relA-/-, and relB-/-) allowed us to dissect specific RelA and RelB activities in the regulation of gene transcription after LTβR stimulation. In wt MEFs, LTβR signals were predominantly transduced by RelA- and/or RelB-containing dimers. Upon LTβR signaling in relA-/- cells, gene regulatory events were mediated by RelB and vice versa in relB-/- cells, changes in gene expression were mediated by RelA. Using this system, we describe novel LTβR-responsive genes that were regulated solely by RelA or RelB or by both RelA and RelB.

Results and discussion

LTβR stimulation of MEFs

For LTβR stimulation, MEFs of each genotype were either left untreated or were treated with agonistic anti-LTβR monoclonal antibody (mAb) for 2.5 or 10 h. For each treatment group, cells from four experiments were pooled. Nuclear protein extracts were used in electrophoretic mobility shift assays (EMSAs) to verify proper LTβR signaling (Figure 1). In wt cells, LTβR signaling resulted in modest induction of κB-binding complexes at the early time point (2.5 h) but strong induction after 10 h of stimulation. Dissection of these complexes with supershifting antibodies revealed that the faster migrating complex contained RelB and the slower migrating complex contained RelA. As expected, in wt cells both RelA and RelB complexes were activated in response to LTβR signaling, whereas in relA-/- cells only RelB- and in relB-/- cells only RelA-containing κB-binding complexes were induced (Figure 1). Recently, slow and relatively weak DNA-binding of NF-κB complexes in response to LTβR ligation was reported [12]. The plateau was reached between 10 and 15 h of LTβR stimulation corresponding to a 2- to 3-fold induction of NF-κB DNA binding. Our results are in agreement with these observations: for each genotype the strongest induction of κB-binding complexes was observed at 10 h. For gene expression profiling we therefore used total RNA isolated from untreated (0 h) and 10 h agonistic anti-LTβR mAb treated wt, relA-/-, and relB-/- MEFs, assuming that stronger DNA-binding activity reflects stronger gene expression changes controlled by NF-κB transcription factor complexes.

thumbnailFigure 1. Confirmation of LTβR stimulation: induction of RelA- and RelB-containing DNA-binding complexes. Wild-type, relA-/-, and relB-/- MEFs were treated with agonistic anti-LTβR mAb for the indicated times and subsequently nuclear extracts were prepared and analyzed by EMSA for NF-κB DNA-binding activity using an Igκ oligo. Specific Igκ DNA-binding complexes are indicated by arrow (RelA-containing dimers) and arrowhead (RelB-containing dimers). Non-specific DNA binding complexes (ns, lower lane) serve as loading control. Supershift analysis was performed using pre-immune serum (pre-imm. serum), anti-RelA antibody (α-RelA Ab), and anti-RelB antibody (α-RelB Ab). Supershifted complexes are indicated by asterisk.

Global gene expression in response to LTβR stimulation in MEFs

To identify RelA- and RelB-regulated genes after LTβR stimulation, we carried out microarray analysis using total RNA from the experiment described above hybridized to CodeLink UniSet Mouse 20K I bioarrays. For statistical analysis, different genotypes were analyzed separately and significantly differentially expressed genes between time points 0 h and 10 h were identified (p < 0.05). The fold change (FC) threshold was determined from the minimal detectable fold change (MDFC) calculated by the CodeLink Expression Analysis v4.1 software (wt: 1.48; relA-/-: 1.54; relB-/-: 1.36). In response to LTβR stimulation, a total of 528 genes were regulated in wt cells. In line with the moderate NF-κB activation seen in the EMSAs the observed gene regulation was also modest: gene expression changes were in the range of +5-fold (induction) and -5-fold (repression). We assigned the 528 LTβR-responsive genes to 4 categories: genes that were significantly regulated (i) only in wt cells (category I, n = 366), (ii) in wt and relA-/- cells (category II, n = 30), (iii) in wt and relB-/- cells (category III, n = 102), and (iv) genes that were significantly regulated in all 3 genotypes (category IV, n = 30) (Figure 2A; for the list of LTβR-responsive genes in wt cells see Additional file 1).

thumbnailFigure 2. LTβR-responsive genes can be allocated into distinct categories. (A) Venn-diagram of significantly (p < 0.05) regulated genes. (B) Schematic depiction of gene expression patterns. The four main categories in (A) can be segregated into further subcategories, depending on whether their genes were upregulated or downregulated. The arrows in the plots show the direction of gene expression changes from non-induced (0 h) to the 10 h induced state in response to LTβR stimulation. The first arrow describes gene expression behavior in wild-type, the second in relA-/-, and the third in relB-/- cells. Horizontal arrows indicate lack of change or statistically insignificant change in gene expression. Arrows pointing upwards or downwards indicate significant positive or negative regulation, respectively.

Additional file 1. Total LTβR transcriptome in wt cells. List of the 528 genes that were LTβR responsive in wt cells (10 h), regardless whether they were also regulated in relA-/- or relB-/- cells, or not.

Format: XLS Size: 88KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

The genes in these four categories could be segregated into further subcategories, which helped us to assign regulatory mechanisms underlying the expression patterns of individual genes (see schematic depiction of gene expression behavior in Figure 2B and lists of genes belonging to different subcategories in Additional files 2, 3, 4, 5).

Additional file 2. LTβR-responsive genes in wt cells. List of genes that were significantly regulated in wt cells, but not in relA-/- or relB-/- cells (10 h; upregulation, cat I/1, n = 161; downregulation, cat I/2, n = 205).

Format: XLS Size: 80KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 3. LTβR-responsive genes in wt and relA-/- cells. List of genes that were significantly regulated in wt and in relA-/- cells (10 h; upregulation, cat II/1, n = 13; downregulation, cat II/2, n = 17).

Format: XLS Size: 15KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 4. LTβR-responsive genes in wt and relB-/- cells. List of genes that were significantly regulated in wt and in relB-/- cells (10 h; upregulation, cat III/1, n = 54; downregulation, cat III/2, n = 43; cat III/3, n = 3; cat III/4, n = 2).

Format: XLS Size: 29KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 5. LTβR-responsive genes in wt, relA-/- and relB-/- cells. List of genes that were significantly regulated in each of the genotypes (10 h; upregulation, cat IV/1, n = 20; downregulation, cat IV/2, n = 10).

Format: XLS Size: 16KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Category (cat) I genes were significantly regulated only in wt cells in response to LTβR stimulation. This group of genes required both RelA and RelB for their LTβR-dependent activation (cat I/1, n = 161) or repression (cat I/2, n = 205). Therefore, expression of these genes did not significantly change in either of the mutant cell lines in response to agonistic anti-LTβR mAb treatment (Figure 2B, Additional file 2).

Category II genes were significantly regulated in wt and relA-/- cells upon LTβR ligation. Genes upregulated (cat II/1, n = 13) or downregulated (cat II/2, n = 17) in both wt and relA-/- cells, but not significantly regulated in relB-/- cells, were considered to be RelB target genes in response to LTβR signaling. Other theoretical patterns could also be appointed to category II, but we did not find any example in our analysis for these subcategories (cat II/3, n = 0 and cat II/4, n = 0) (Figure 2B, Additional file 3).

Genes belonging to category III were significantly regulated in wt and relB-/- cells in response to LTβR stimulation. Genes upregulated (cat III/1, n = 54) or downregulated (cat III/2, n = 43) in both wt and relB-/- cells, but not significantly regulated in relA-/- cells, were considered to be RelA target genes in response to LTβR signaling. Negligible numbers of genes in category III could also be allocated to cat III/3 and III/4 (n = 3 and n = 2, respectively) (Figure 2B, Additional file 4). However, these genes were not further analyzed. The significantly larger number of RelA- (cat III) compared to RelB-regulated genes (cat II; Figure 2A) is likely to be a consequence of the stronger LTβR-induced DNA binding of RelA compared to RelB complexes (Figure 1).

Category IV genes were significantly regulated in each of the genotypes in response to LTβR ligation. Although eight theoretically possible gene expression behaviors exist, we only found genes that belonged to two easily explainable scenarios: genes were either upregulated (cat IV/1, n = 20), or downregulated (cat IV/2, n = 10) in each genotype upon LTβR signaling (Figure 2B, Additional file 5). Most likely, both RelA and RelB contributed redundantly to their regulation or alternatively, a third factor/pathway controlled these genes in response to LTβR stimulation. JNK (c-Jun N-terminal kinase) is a possible candidate for such a third pathway, as there are indications that LTβR stimulation leads to activation of JNK. However, the experimental setup in those studies was different from ours as LTβR-overexpressing HEK293 cells [23] or treatment of MEFs with the LTβR agonist LIGHT (lymphotoxin-related inducible ligand that competes for glycoprotein D binding to herpesvirus entry mediator on T cells) [24] were studied.

FC values observed in the three cell lines at 10 h compared to 0 h are displayed in a heatmap that also reflects the four categories and their subcategories (Figure 3, for a zoomable/enlarged version of FC heatmaps supplied with gene symbols and GenBank Accession Numbers see Additional file 6).

thumbnailFigure 3. Fold change heatmaps. Heatmaps displaying the fold change values observed in the three different cell lines at 10 h compared to 0 h. The color code indicates the fold change values between -2.5-fold downregulation (green) and +2.5-fold upregulation (red). Fold change of -2.5 and below are depicted in the brightest green and fold change of +2.5 and above are shown in the brightest red. Black indicates no change in gene expression. Each horizontal line on the heatmap corresponds to one gene. Genes are arranged by their subcategory (see bars on the left) and main categories are divided by a horizontal white line.

Additional file 6. Zoomable/enlarged version of fold change heatmaps. Heatmaps displaying the fold change values observed in the three different cell lines at 10 h compared to 0 h. For figure legend see Figure 3. Gene symbols and GenBank Accession Numbers (in brackets) are also displayed.

Format: PDF Size: 124KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Interestingly, in the two subcategories with the largest number of genes both RelA and RelB together were required for LTβR-induced gene regulation (161 cat I/1 genes for their activation and 205 cat I/2 genes for their repression). In case one of the transcription factors was missing the other one was not able to ensure regulation alone, suggesting significant crosstalk between the two NF-κB activation pathways. In response to LTβR stimulation, sequential engagement of the classical and alternative pathway was suggested, resulting in initial DNA binding by RelA followed by RelB complexes [7,9]. These findings may suggest a scenario where RelA binds first to the DNA in the promoter of category I genes, loosens up chromatin, thereby enabling subsequent DNA binding and gene regulatory action by RelB [25]. Alternatively, since relB is an NF-κB target gene [26] RelA may ensure sufficiently high expression of RelB and in the absence of RelA the reduced RelB levels cannot mediate proper regulation of certain LTβR target genes. This possibility is supported by the observation that in the absence of RelA both RelB protein levels and binding of RelB to κB sites were reduced (Figure 1 and data not shown) [13].

Meta analysis of LTβR-dependent transcriptomes

LTβR signaling is best known in the context of secondary lymphoid organ development and a recent expression profiling study described LTβR-dependent transcriptomes in lymph nodes and follicular dendritic cells (FDCs) [27]. However, increasing evidence suggests that LTβR also plays a role in non-lymphoid organs such as epithelial tissues during embryonic development [28] and adult liver [21,22].

To interpret our results in the light of other studies investigating LTβR signaling, we compared our LTβR-responsive genes with two recently published LTβR-dependent transcriptomes. Huber et al. identified transcripts in murine mesenteric lymph nodes affected in vivo by the administration of a soluble LTβR-Ig decoy receptor which blocks LTβR signaling [27]. A gene cluster of 80 unique transcripts that showed decreased expression after LTβR blockade was further analyzed. Twelve genes in this cluster were also associated with germinal centers (GCs)/FDC. A few common genes were found between our analysis and the LTβR-dependent transcriptomes described by Huber et al. Dclk1 and enpp2 (doublecortin-like kinase 1; GenBank Accession Number: NM_019978 and ectonucleotide pyrophosphatase/phosphodiesterase 2 or autotaxin; GenBank Accession Number: NM_015744) expression was moderately decreased 3 d after LTβR blockade (FC: 0.70× and 0.66×, respectively) [27]. In our hands, both genes were upregulated in response to LTβR stimulation in a RelA-dependent manner (cat III/1, for enpp2 see also Table 12). Enpp2 was also found to be associated with GC/FDC in mesenteric lymph nodes [27]. Moreover, Enpp2 (also called autotaxin) has been recently described as a new molecule in lymphocyte homing through high endothelial venules (HEVs) [29]. Collectively, these findings suggest that LTβR, in addition to its well-described effect on the HEV differentiation program [30], further contributes via RelA-dependent upregulation of enpp2 to lymphocyte homing through HEVs. Unfortunately, we could not detect further genes with a similar regulation pattern in our and Huber and colleagues' studies. This lack of overlap could be the consequence of several reasons: (i) different modes of function and kinetics of antagonistic LTβR-Ig vs agonistic anti-LTβR mAb application, (ii) incubation time (3 d treatment with LTβR-Ig vs 10 h treatment with agonistic anti-LTβR mAb), or (iii) in vivo collection of different cell types influenced by the treatment vs in vitro cell culture system using MEFs.

Lo et al. described a hepatic gene expression profile of wt vs lck-LIGHT transgenic mice (overexpressing the LTβR ligand LIGHT on the surface of T lymphocytes) [22]. A group of significantly regulated genes (n = 19) involved in lipid and cholesterol metabolism was identified. The gene that displayed the highest level of regulation (23-fold repression in transgenic vs wt mice) encodes for hepatic lipase, a key enzyme in lipid metabolism. We did not observe repression of hepatic lipase in our experiments, most probably due to its restricted expression on the surface of hepatocytes. However, we found another gene belonging to the lipid/cholesterol metabolism-related group described by Lo and colleagues. Ralgds (ral guanine nucleotide dissociation stimulator, GenBank Accession Number: NM_009058) expression was increased in the liver of transgenic mice and also upregulated in our LTβR stimulation experiments, belonging to the RelA-responsive genes (cat III/1, Table 12).

Gene Ontology (GO) enrichment analysis

Our goal was not only to define the LTβR-dependent transcriptome in MEFs, but also to assign regulatory mechanisms to LTβR signaling, i.e. to examine which part of the LTβR transcriptome is regulated by RelA, RelB, or both. We started out with GO enrichment analysis of significantly regulated genes to identify biological processes, molecular functions, and cellular components putatively regulated in the categories described above. Compared to molecular functions and cellular components, GO analysis of biological processes yielded the most conclusive results.

First, GO analysis was performed on the total LTβR transcriptome in wt cells to see how LTβR signaling influences biological processes in these fibroblasts, regardless whether these genes were also regulated in relA-/- or relB-/- cells (Category: Total wild-type, Table 1). For interpretation of our data we chose GO terms with p < 0.01. As lower limit, we did not consider GO terms with less than 3 annotated genes in the list of differentially regulated genes since they are too specific. As upper limit we did not use GO terms represented by more than 600 genes on the microarray since they are too general. Among the considered GO terms we found that apoptosis/cell death (A/CD)- and cell cycle (CCY)-related processes were overrepresented. We also found that genes annotated with "response to biotic stimulus", "immune system process" (immune related (IR) features) as well as "blood vessel morphogenesis" and "angiogenesis" (blood vessel development related (BR) features) were enriched. Collectively, these data indicate that LTβR signaling largely influences cell survival/cell proliferation features. Moreover, it has an impact on immune responses and blood vessel development/angiogenesis related processes. Since these GO terms were found in LTβR-stimulated "non-immune" fibroblasts it is likely that LTβR signaling regulates similar biological processes in stromal cells of secondary lymphoid tissues governing lymphorganogenesis and maintaining lymphoid tissue architecture.

Table 1. Gene Ontology analysis of total LTβR transcriptome in wild-type cells

Next, we carried out GO analysis for the four main categories and for all subcategories with at least 20 genes. Interpretation of the data was performed applying the same criteria as above. GO analysis of category I genes revealed those biological processes that were overrepresented only in LTβR-stimulated wt cells, i.e. in the presence of both RelA and RelB (Table 2). Amongst these processes, CCY-related terms dominated. Subsequently, we analyzed cat I/1 (containing genes that were upregulated exclusively in wt cells) and found enrichment of IR- and cell/biological adhesion (important events in immune cell migration)-related terms on the list of biological processes (Table 3). This finding indicates that in the absence of RelA or RelB a considerable portion of LTβR-stimulated immune response-related events cannot be carried out; fibroblasts need both molecules to execute these processes. In cat I/2 (containing genes that are downregulated exclusively in wt cells) we found enrichment of CCY-related terms on the list of overrepresented biological processes (Table 4). This finding indicates that in wt cells an important action of RelA and RelB is to downregulate numerous genes that are implicated in cell cycle regulation in response to LTβR signaling.

Table 2. Gene Ontology analysis of category I

Table 3. Gene Ontology analysis of category I/1

Table 4. Gene Ontology analysis of category I/2

Since cat II/1 and II/2 had only few genes (n = 13 and n = 17, respectively), investigation of GO terms for these groups of genes was not meaningful. GO analysis of the main category II (containing genes that were regulated – either up or down – in wt and relA-/- cells, n = 30) revealed only one enriched GO term, the cell cycle (Table 5). Thus, in response to LTβR signaling a characteristic feature of RelB was to influence cell cycle-related events.

Table 5. Gene Ontology analysis of category II

Category III contains genes that were regulated – either up or down – in wt and relB-/- cells in response to LTβR stimulation. Among enriched biological processes, the new and in previous categories not yet observed theme taxis and response to external/chemical stimulus (T) dominated, but A/CD-related events also appeared (Table 6). As expected, the theme IR was also represented among the enriched biological processes. This shows that RelA is not only a signal transducer for immune responses and apoptosis/cell death, but also has an impact on the transcription of taxis- and stimulus-responsive genes following LTβR ligation. Among the enriched biological processes of cat III/1 we observed again overrepresentation of T and IR processes (Table 7), revealing that in response to LTβR signaling RelA strongly influenced these events via upregulation of several genes. In cat III/2 we found genes that were repressed by RelA. In this subcategory RelA on one hand regulated several BR events. On the other hand, it turned out to be a negative regulator of genes involved in ion homeostasis (ION) downstream of the LTβR (Table 8).

Table 6. Gene Ontology analysis of category III

Table 7. Gene Ontology analysis of category III/1

Table 8. Gene Ontology analysis of category III/2

Category IV contains genes that were regulated – either up or down – in each of the cell types in response to LTβR stimulation (Table 9). IR processes were overrepresented, but the terms related to hematopoietic or lymphoid organ development (LY) and taxis (T) were also present on the list of enriched biological processes. Unfortunately, we could not analyze cat IV/2, as it comprises too few genes (n = 10). Cat IV/1 contains 20 genes that were upregulated, irrespective of the genotype (Table 10). These genes primarily belong to IR and T. Possibly, RelA and RelB redundantly regulate these events or alternatively a RelA- and RelB-independent third factor/pathway (e.g. JNK) controls these biological processes following LTβR ligation. Table 11 shows a summary of our GO analysis.

Table 9. Gene Ontology analysis of category IV

Table 10. Gene Ontology analysis of category IV/1

Table 11. Summary of Gene Ontology analysis results

Verification of microarray results by qRT-PCR

The changes in mRNA levels of several known as well as novel LTβR-responsive genes on the microarray were confirmed by quantitative real-time reverse-transcription-PCR (qRT-PCR), using RNA from three independent LTβR stimulation experiments (Table 12). In agreement with previous reports, we also found induction of nfkb2 [5,6], ccl2/mcp1 [6], and ikba expression [31] in LTβR-stimulated wt fibroblasts. In addition, our data indicate that both RelA and RelB redundantly contributed to the proper regulation of these genes in response to LTβR stimulation. However, we did not observe LTβR-dependent upregulation of lymphorganogenic chemokines as described by others. Ccl21, ccl19, cxcl13, and cxcl12 were shown to be LTβR-induced genes in spleen 8 h after peritoneal injection of an agonistic anti-LTβR mAb [5]. Possibly, cell context-specific signaling accounts for the difference observed between splenocytes and established 3T3 fibroblasts used in our experiments. Basak et al. observed modest upregulation of cxcl13 and ccl21 in established wt 3T3 fibroblasts after 24 h treatment with agonistic anti-LTβR mAb [13]. To reduce indirect gene regulatory effects due to rather long stimulation we activated LTβR signaling only for 10 h, where modulation of these chemokines was not observed.

Table 12. Verification of microarray results by qRT-PCR

Table 13. LTβR responsive qRT-PCR verified genes in literature

Importantly, we verified novel LTβR-responsive genes and appointed regulatory molecules to them. For a complete list of verified genes see Table 12. Here, some of those verified genes are discussed in more detail.

GO analysis revealed that LTβR stimulation resulted in the regulation of IR processes (Table 11). Except category "Total wild-type", where we could not assign regulatory molecules, in all categories where IR processes were enriched, RelA alone or together with RelB acted as a positive factor. Cx3cl1 (chemokine C-X3-C motif ligand 1/fractalkine) is one of the IR genes in cat I/1. Several studies document that NF-κB upregulates cx3cl1, e.g. in rat aortic endothelial cells upon interleukin-1β (IL-1β), TNF, and lipopolysaccharide treatment [32] or in human coronary artery smooth muscle cells [33]. The latter work shows that atherogenic lipids induce adhesion of artery smooth muscle cells to macrophages via the upregulation of cx3cl1 in a TNF/NF-κB-dependent manner. In our experiments this gene was upregulated in response to LTβR stimulation dependent on RelA and RelB. This data suggests that LTβR, via employing RelA and RelB together, may act as a proatherogenic factor.

IR- and T-related processes were also enriched in cat III and cat III/1 according to the GO analysis. Cd74/ii (invariant polypeptide of major histocompatibility complex, class II antigen-associated) and cxcl10/ip10 (chemokine C-X-C motif ligand 10/interferon-inducible protein-10) are two genes in cat III/1 and assigned to IR and T. CD74/Ii is involved in antigen processing and presentation and CXCL10 is chemotactic for monocytes and T cells. Moreover, expression of CXCL10, along with two other CXCR3-binding chemokines CXCL9 and CXCL11, can be induced in carcinoma cells by LTβR agonists. These chemokines function as potent chemoattractants for activated T, NK, and dendritic cells, which may contribute to antitumor immune responses [34]. In our experiments, expression of cd74/ii and cxcl10/ip10 was upregulated by LTβR signaling in wt and relB-/- cells. Thus, LTβR signaling via RelA may (i) attract T lymphocytes and promote antigen presentation by dendritic cells in the context of MHC class II and (ii) facilitate antitumor responses against cancer cells.

As indicated by GO analysis, IR- and T-related biological processes were significantly regulated in cat IV and cat IV/1. Amongst others, genes encoding proteins that participate in innate immune responses, like ccl7/mcp3, are also represented in these groups. Ccl7/mcp3 encodes the proinflammatory chemokine C-C motif ligand 7/monocyte chemotactic protein-3. Expression of ccl7/mcp3 was upregulated by LTβR signaling in each of the genotypes, indicating redundant positive regulation by RelA and RelB or upregulation via another RelA- and RelB-independent pathway.

Collectively, positive regulation of the expression of proinflammatory chemokines like cx3cl1, cxcl10, ccl7 (but also others, see Table 12) by LTβR suggests that LTβR signaling, besides regulating development and organization of secondary lymphoid structures, also participates in innate/inflammatory immune responses and for that primarily RelA action seems to be necessary.

Moreover, we found that LTβR signaling functions beyond the regulation of immune responses and organization of lymphoid structures. PPARγ (peroxisome proliferator activated receptor γ) is a key-regulatory transcription factor in the process of adipocyte differentiation and activation of PPARγ promotes the storage of fat [35]. The work of Fu and colleagues suggests that LTβR affects lipid homeostasis by downregulating hepatic lipase expression [22]. Hepatic lipase is expressed on the surface of hepatocytes in the liver. It promotes receptor-mediated uptake of plasma lipoproteins that harbor triglycerides and cholesterol and specifically catalyzes hydrolysis of triglycerides, actions that are suppressed when LTβR signaling is switched on. Expression of pparg was negatively affected by LTβR signaling in wt and relA-/- but not in relB-/- cells (belonging to cat II/2 genes), indicating that this gene was downregulated by RelB in response to LTβR stimulation. Our finding is a further indication that LTβR signaling represses lipogenesis and it may do so via RelB. It has been shown that ligand-induced transactivation by PPARγ is suppressed by IL-1 and TNF and that this suppression is mediated through NF-κB (p50-RelA) [36]. However, unlike suppression of PPARγ by p50-RelA, where this heterodimer blocks PPARγ binding to DNA by forming a complex with PPARγ and its co-activator PGC-2, LTβR-mediated suppression of pparg occurred via transcriptional repression executed by RelB. Further experiments are required to find out whether RelB directly or indirectly mediates repression of pparg transcription in response to LTβR signaling. The repressive effect of LTβR signaling on adipogenesis has been confirmed in MEFs that were induced for adipogenic differentiation. LTβR stimulation resulted in attenuated lipid droplet accumulation as well as in reduced pparg and adipogenic marker gene (fabp4/ap2) expression under conditions that promote differentiation into adipocytes (unpublished results).

Conclusion

This study is the first systematic dissection of the RelA- and RelB-driven transcriptome response downstream of the LTβR. We confirmed previously described LTβR-regulated genes. More importantly, we identified novel LTβR-responsive genes and assigned underlying regulatory mechanisms executed by RelA and/or RelB to them (Table 13). We found that the majority of LTβR-regulated genes required the presence of both RelA and RelB, suggesting significant crosstalk between the two NF-κB activation pathways. Gene Ontology analysis confirmed that LTβR-NF-κB target genes were predominantly involved in the regulation of immune responses. However, other biological processes such as apoptosis/cell death, cell cycle, angiogenesis, and taxis were also regulated by LTβR signaling. Furthermore, we show that LTβR stimulation downregulated expression of the gene encoding PPARγ, suggesting that LTβR signaling may repress adipogenic differentiation by attenuating the levels of this key adipogenic transcription factor. Our findings are significant since they indicate a role for LTβR signaling beyond immune responses and lymphoid organ development and assign underlying gene expression regulatory mechanisms to the LTβR transcriptome.

Methods

Cell culture

Mouse embryonic 3T3 fibroblasts (wild-type, relA-/-, and relB-/-; kind gift from A. Hoffmann) were cultured at 37°C in Dulbecco's modified Eagle's medium (GIBCO/Invitrogen, Karlsruhe, Germany) supplemented with 10% heat-inactivated bovine calf serum (Perbio Science, Bonn, Germany), penicillin (100 U/ml), streptomycin (100 μg/ml), and Glutamax I (2 mM) (GIBCO/Invitrogen) and treated with agonistic anti-LTβR mAb (1 μg/ml, clone AC.H6; kind gift from J. Browning and P. Rennert).

EMSA

Preparation of nuclear extracts and EMSAs were essentially performed as previously described [37]. Nuclear and cytoplasmic fractions were prepared according to standard procedures [38].

RNA isolation

Total cellular RNA was isolated using the RNeasy Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. Possible contamination by genomic DNA was removed by DNaseI treatment using the RNase-Free DNase Set (Qiagen). Quality of RNA samples was checked by spectrophotometry and agarose gel electrophoresis. RNAs (2 μg total RNA per sample) were used for cRNA preparation for microarrays only when the ratio A260:A280 was 1.8–2.1 and the RNA was intact.

Microarrays

Microarray analysis was performed using CodeLink UniSet Mouse 20K I bioarrays (GE Healthcare, Munich, Germany), a one-color system where for each of the investigated 19,801 transcripts there is one 30–mer oligo probe spotted per slide. For gene expression profiling, untreated (0 h) and 10 h agonistic anti-LTβR mAb treated wt, relA-/-, and relB-/- MEFs were used. For every treatment group, cells from 4 experiments were pooled, total RNA isolated, cRNA prepared and hybridized onto the bioarrays in technical triplicates. cRNA target preparation, bioarray hybridization and detection were carried out according to the manufacturer's protocol provided with the CodeLink Expression Assay Reagent Kit. For scanning microarrays, a GenePix 4000B Array Scanner and GenePix Pro 4.0 software (Axon Instruments Inc./Molecular Devices, Munich, Germany) were employed according to settings suggested by the protocol provided with the CodeLink Expression Assay Reagent Kit. Microarray data have been deposited in NCBIs GEO http://www.ncbi.nlm.nih.gov/geo/ webcite and are accessible through GEO series accession number GSE11963.

Microarray data preprocessing

Microarray raw data of stimulated and unstimulated MEFs were analyzed using the Codelink™ Expression Analysis v4.1 software (GE Healthcare) and MDFC values were extracted. All subsequent analyses were performed using R and Bioconductor. For the analysis only genes with probe type 'DISCOVERY' were considered (19,801 genes) and all genes flagged MSR (Manufactory Slide Report) in any sample were excluded (leaving 19,580 genes). To remove negative expression values (local background > spot intensity) raw intensities with values < 0.01 were set to 0.01. The raw intensities of each array were scaled to the array median. After logarithmizing the expression values quantile normalization was applied across all arrays.

Differentially expressed genes

Array data for the different genotypes were analyzed separately. A gene was included in the analysis if it was flagged 'G' (good) or 'S' (contains saturated pixels) on at least two arrays in any of the two groups (stimulated or unstimulated). Furthermore, genes selected were required to have a FC higher than or equal to the FC threshold determined from the maximum MDFC in these groups. To identify genes significantly differentially expressed after stimulation, a Student's t-test was performed for the previously filtered genes. The resulting p values were corrected for multiple testing using the method of Benjamini and Hochberg [39]. Allowing a false discovery rate of 5%, a total of 528 genes were identified that were significantly regulated in wt cells (regardless whether they were regulated somewhere else). From these, 366 genes were regulated exclusively in wt, 30 genes in wt and relA-/-, 102 in wt and relB-/- cells and 30 genes in all 3 genotypes.

Functional analysis with GO

Analysis of functional enrichment was performed employing Fisher's exact test. The resulting p values (p < 0.01) were used to rank GO terms according to their significance. Terms with more than 600 genes on the array or less than 3 genes on the list of investigated genes were regarded as too general or too specific, respectively, and excluded from the analysis. Expert knowledge was used to assign broader themes to specific GO categories.

qRT-PCR

For qRT-PCR, first strand cDNA was obtained from 2 μg of total RNA for each treatment group using oligo-dT primers and M-MLV Reverse Transcriptase kit (Promega, Mannheim, Germany) according to manufacturer's protocols. qRT-PCRs were performed in an iCycler Thermal Cycler real-time PCR machine (Bio-Rad Laboratories, Hercules, CA) using SYBR Green I as detector dye and reagents from the Quantace SensiMix DNA Kit (Quantace Ltd., Watford, UK). Primers for qRT-PCRs with Tm of 60°C were designed using Primer3 software (v. 0.4.0; http://frodo.wi.mit.edu webcite) [40]. For individual samples, each gene was tested in triplicates and the mean of the 3 cycle threshold values was used to calculate relative expression levels. For normalization, β-actin was used as an endogenous reference gene to correct for variation in RNA content and variation in the efficiency of the reverse transcription reaction. Statistical analysis of qRT-PCR results from 3 independent LTβR stimulation experiments was performed employing a Welch test. Forward (F) and reverse primers (R) in 5' to 3' orientation were: Nfkb2_F: GCTAATGTGAATGCCCGGAC, Nfkb2_R: CTTTGGGTATCCCTCTCAGGC, Ccl2_F: CCCACTCACCTGCTGCTACT, Ccl2_R: TCTGGACCCATTCCTTCTTG, IκBα_F: TGCACTTGGCAATCATCCAC, IκBα_R: TTCCTCGAAAGTCTCGGAGCT, Ralgds_F: CATCCACCGCCTAAAGAAGA, Ralgds_R: GGGCTCTCCTAGGGTTCATC, Cx3cl1_F: GGCTAAGCCTCAGAGCATTG, Cx3cl1_R: CATTTTCCTCTGGGGTTGA, Pparg_F: TCATGACCAGGGAGTTCCTC, Pparg_R: GGCGGTCTCCACTGAGAATA, Enpp2_F: TGGCTTACGTGACATTGAGG, Enpp2_R: GTCGGTGAGGAAGGATGAAA, Birc3_F: TGACGTGTGTGACACCAATG, Birc3_R: TGAGGTTGCTGCAGTGTTTC, Cxcl10_F: AAGTGCTGCCGTCATTTTCT, Cxcl10_R: GTGGCAATGATCTCAACACG, Irf1_F: ACCCTGGCTAGAGATGCAGA, Irf1_R: TTTGTATCGGCCTGTGTGAA, Cd74_F: ATGACCCAGGACCATGTGAT, Cd74_R: CCAGGGAGTTCTTGCTCATC, Fosl1_F: CAAAATCCCAGAAGGAGACAAG, Fosl1_R: AAAAGGAGTCAGAGAGGGTGTG, Ccl7_F: AATGCATCCACATGCTGCTA, Ccl7_R: ATAGCCTCCTCGACCCACTT, Cxcl1_F: GCTGGGATTCACCTCAAGAA, Cxcl1_R: TGGGGACACCTTTTAGCATC, Id2_F: CCCCAGAACAAGAAGGTGAC, Id2_R: ATTCAGATGCCTGCAAGGAC, β-actin_F: TGGCGCTTTTGACTCAGGA, β-actin_R: GGGAGGGTGAGGGACTTCC

Abbreviations

LTβR: lymphotoxin-β receptor; IκBα: inhibitor of NF-κB α; MEF: mouse embryonic fibroblasts; GO: Gene Ontology; PPARγ/pparg: peroxisome proliferator activated receptor-γ; TNFR1: tumor necrosis factor receptor 1; TLR4: Toll-like receptor 4; NEMO: NF-κB essential modulator; IKK: IκB kinase; NIK: NF-κB-inducing kinase; relA-/-: RelA-deficient; TNF: tumor necrosis factor; relB-/-: RelB-deficient; wt: wild-type; mAb: monoclonal antibody; EMSA: electrophoretic mobility shift assay; FC: fold change; MDFC: minimal detectable fold change; cat: category; JNK: c-Jun N-terminal kinase; LIGHT: lymphotoxin-related inducible ligand that competes for glycoprotein D binding to herpesvirus entry mediator on T cells; FDC: follicular dendritic cell; GC: germinal center; HEV: high endothelial venule; A/CD: apoptosis/cell death; CCY: cell cycle; IR: immune related; BR: blood vessel development related; T: taxis, response to external/chemical stimulus; ION: ion homeostasis; LY: hematopoietic or lymphoid organ developmental processes; qRT-PCR: quantitative real-time reverse-transcription PCR; IL-1: interleukin-1; MSR: Manufactory Slide Report; SD: standard deviation.

Authors' contributions

AL: carried out the molecular genetic studies, analyzed and interpreted data, drafted manuscript. DR: carried out the bioinformatic and statistic analysis, participated in study design, analyzed and interpreted data. DA: participated in the bioinformatic and statistic analysis, analyzed and interpreted data. ZBY: initiated and participated in the molecular genetic studies. UM: participated in the bioinformatic and statistic analysis. AJRH: supported bioinformatic analysis. FW: conceived the study, participated in its design and coordination, interpreted data, helped to write the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Jeffrey Browning and Paul Rennert for agonistic anti-LTβR mAb and Alexander Hoffmann for mouse embryonic 3T3 fibroblasts (wild-type, relA-/-, and relB-/-). We gratefully acknowledge Heike Mondrzak, Ulrike Schure, Melissa Lehmann, Kerstin Andreas, and Sandra Westhaus for technical help with the qRT-PCR. We are indebted to Hans Peter Saluz for providing us with the possibility to use the GenePix 4000B Array Scanner. We also thank Marc Riemann for valuable discussions on this manuscript.

References

  1. Ghosh S, Karin M: Missing pieces in the NF-kappaB puzzle.

    Cell 2002, 109(Suppl):S81-96. PubMed Abstract | Publisher Full Text OpenURL

  2. Hayden MS, Ghosh S: Signaling to NF-kappaB.

    Genes Dev 2004, 18:2195-2224. PubMed Abstract | Publisher Full Text OpenURL

  3. Hayden MS, Ghosh S: Shared principles in NF-kappaB signaling.

    Cell 2008, 132:344-362. PubMed Abstract | Publisher Full Text OpenURL

  4. Hoffmann A, Natoli G, Ghosh G: Transcriptional regulation via the NF-kappaB signaling module.

    Oncogene 2006, 25:6706-6716. PubMed Abstract | Publisher Full Text OpenURL

  5. Dejardin E, Droin NM, Delhase M, Haas E, Cao Y, Makris C, Li ZW, Karin M, Ware CF, Green DR: The lymphotoxin-beta receptor induces different patterns of gene expression via two NF-kappaB pathways.

    Immunity 2002, 17:525-535. PubMed Abstract | Publisher Full Text OpenURL

  6. Derudder E, Dejardin E, Pritchard LL, Green DR, Korner M, Baud V: RelB/p50 dimers are differentially regulated by tumor necrosis factor-alpha and lymphotoxin-beta receptor activation: critical roles for p100.

    J Biol Chem 2003, 278:23278-23284. PubMed Abstract | Publisher Full Text OpenURL

  7. Muller JR, Siebenlist U: Lymphotoxin beta receptor induces sequential activation of distinct NF-kappa B factors via separate signaling pathways.

    J Biol Chem 2003, 278:12006-12012. PubMed Abstract | Publisher Full Text OpenURL

  8. Weih F, Caamano J: Regulation of secondary lymphoid organ development by the nuclear factor-kappaB signal transduction pathway.

    Immunol Rev 2003, 195:91-105. PubMed Abstract | Publisher Full Text OpenURL

  9. Yilmaz ZB, Weih DS, Sivakumar V, Weih F: RelB is required for Peyer's patch development: differential regulation of p52-RelB by lymphotoxin and TNF.

    EMBO J 2003, 22:121-130. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Bonizzi G, Karin M: The two NF-kappaB activation pathways and their role in innate and adaptive immunity.

    Trends Immunol 2004, 25:280-288. PubMed Abstract | Publisher Full Text OpenURL

  11. Dejardin E: The alternative NF-kappaB pathway from biochemistry to biology: pitfalls and promises for future drug development.

    Biochem Pharmacol 2006, 72:1161-1179. PubMed Abstract | Publisher Full Text OpenURL

  12. Basak S, Kim H, Kearns JD, Tergaonkar V, O'Dea E, Werner SL, Benedict CA, Ware CF, Ghosh G, Verma IM, Hoffmann A: A fourth IkappaB protein within the NF-kappaB signaling module.

    Cell 2007, 128:369-381. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Basak S, Shih VF, Hoffmann A: Generation and activation of multiple dimeric transcription factors within the NF-kappaB signaling system.

    Mol Cell Biol 2008, 28:3139-3150. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Basak S, Hoffmann A: Crosstalk via the NF-kappaB signaling system.

    Cytokine Growth Factor Rev 2008, 19:187-197. PubMed Abstract | Publisher Full Text OpenURL

  15. Mebius RE: Organogenesis of lymphoid tissues.

    Nat Rev Immunol 2003, 3:292-303. PubMed Abstract | Publisher Full Text OpenURL

  16. Schneider K, Potter KG, Ware CF: Lymphotoxin and LIGHT signaling pathways and target genes.

    Immunol Rev 2004, 202:49-66. PubMed Abstract | Publisher Full Text OpenURL

  17. Ware CF: Network communications: lymphotoxins, LIGHT, and TNF.

    Annu Rev Immunol 2005, 23:787-819. PubMed Abstract | Publisher Full Text OpenURL

  18. Beg AA, Sha WC, Bronson RT, Ghosh S, Baltimore D: Embryonic lethality and liver degeneration in mice lacking the RelA component of NF-kappa B.

    Nature 1995, 376:167-170. PubMed Abstract | Publisher Full Text OpenURL

  19. Weih F, Carrasco D, Durham SK, Barton DS, Rizzo CA, Ryseck RP, Lira SA, Bravo R: Multiorgan inflammation and hematopoietic abnormalities in mice with a targeted disruption of RelB, a member of the NF-kappa B/Rel family.

    Cell 1995, 80:331-340. PubMed Abstract | Publisher Full Text OpenURL

  20. Weih DS, Yilmaz ZB, Weih F: Essential role of RelB in germinal center and marginal zone formation and proper expression of homing chemokines.

    J Immunol 2001, 167:1909-1919. PubMed Abstract | Publisher Full Text OpenURL

  21. Anders RA, Subudhi SK, Wang J, Pfeffer K, Fu YX: Contribution of the lymphotoxin beta receptor to liver regeneration.

    J Immunol 2005, 175:1295-1300. PubMed Abstract | Publisher Full Text OpenURL

  22. Lo JC, Wang Y, Tumanov AV, Bamji M, Yao Z, Reardon CA, Getz GS, Fu YX: Lymphotoxin beta receptor-dependent control of lipid homeostasis.

    Science 2007, 316:285-288. PubMed Abstract | Publisher Full Text OpenURL

  23. Chang YH, Hsieh SL, Chen MC, Lin WW: Lymphotoxin beta receptor induces interleukin 8 gene expression via NF-kappaB and AP-1 activation.

    Exp Cell Res 2002, 278:166-174. PubMed Abstract | Publisher Full Text OpenURL

  24. Kim YS, Nedospasov SA, Liu ZG: TRAF2 plays a key, nonredundant role in LIGHT-lymphotoxin beta receptor signaling.

    Mol Cell Biol 2005, 25:2130-2137. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Saccani S, Pantano S, Natoli G: Modulation of NF-kappaB activity by exchange of dimers.

    Mol Cell 2003, 11:1563-1574. PubMed Abstract | Publisher Full Text OpenURL

  26. Bren GD, Solan NJ, Miyoshi H, Pennington KN, Pobst LJ, Paya CV: Transcription of the RelB gene is regulated by NF-kappaB.

    Oncogene 2001, 20:7722-7733. PubMed Abstract | Publisher Full Text OpenURL

  27. Huber C, Thielen C, Seeger H, Schwarz P, Montrasio F, Wilson MR, Heinen E, Fu YX, Miele G, Aguzzi A: Lymphotoxin-beta receptor-dependent genes in lymph node and follicular dendritic cell transcriptomes.

    J Immunol 2005, 174:5526-5536. PubMed Abstract | Publisher Full Text OpenURL

  28. Browning JL, French LE: Visualization of lymphotoxin-beta and lymphotoxin-beta receptor expression in mouse embryos.

    J Immunol 2002, 168:5079-5087. PubMed Abstract | Publisher Full Text OpenURL

  29. Kanda H, Newton R, Klein R, Morita Y, Gunn MD, Rosen SD: Autotaxin, an ectoenzyme that produces lysophosphatidic acid, promotes the entry of lymphocytes into secondary lymphoid organs.

    Nat Immunol 2008, 9:415-423. PubMed Abstract | Publisher Full Text OpenURL

  30. Browning JL, Allaire N, Ngam-Ek A, Notidis E, Hunt J, Perrin S, Fava RA: Lymphotoxin-beta receptor signaling is required for the homeostatic control of HEV differentiation and function.

    Immunity 2005, 23:539-550. PubMed Abstract | Publisher Full Text OpenURL

  31. Bonizzi G, Bebien M, Otero DC, Johnson-Vroom KE, Cao Y, Vu D, Jegga AG, Aronow BJ, Ghosh G, Rickert RC, Karin M: Activation of IKKalpha target genes depends on recognition of specific kappaB binding sites by RelB:p52 dimers.

    EMBO J 2004, 23:4202-4210. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Garcia GE, Xia Y, Chen S, Wang Y, Ye RD, Harrison JK, Bacon KB, Zerwes HG, Feng L: NF-kappaB-dependent fractalkine induction in rat aortic endothelial cells stimulated by IL-1beta, TNF-alpha, and LPS.

    J Leukoc Biol 2000, 67:577-584. PubMed Abstract | Publisher Full Text OpenURL

  33. Barlic J, Zhang Y, Murphy PM: Atherogenic lipids induce adhesion of human coronary artery smooth muscle cells to macrophages by up-regulating chemokine CX3CL1 on smooth muscle cells in a TNFalpha-NFkappaB-dependent manner.

    J Biol Chem 2007, 282:19167-19176. PubMed Abstract | Publisher Full Text OpenURL

  34. Lukashev M, LePage D, Wilson C, Bailly V, Garber E, Lukashin A, Ngam-ek A, Zeng W, Allaire N, Perrin S, Xu X, Szeliga K, Wortham K, Kelly R, Bottiglio C, Ding J, Griffith L, Heaney G, Silverio E, Yang W, Jarpe M, Fawell S, Reff M, Carmillo A, Miatkowski K, Amatucci J, Crowell T, Prentice H, Meier W, Violette SM, Mackay F, Yang D, Hoffman R, Browning JL: Targeting the lymphotoxin-beta receptor with agonist antibodies as a potential cancer therapy.

    Cancer Res 2006, 66:9617-9624. PubMed Abstract | Publisher Full Text OpenURL

  35. Kubota N, Terauchi Y, Miki H, Tamemoto H, Yamauchi T, Komeda K, Satoh S, Nakano R, Ishii C, Sugiyama T, Eto K, Tsubamoto Y, Okuno A, Murakami K, Sekihara H, Hasegawa G, Naito M, Toyoshima Y, Tanaka S, Shiota K, Kitamura T, Fujita T, Ezaki O, Aizawa S, Nagai R, Tobe K, Kimura S, Kadowaki T: PPAR gamma mediates high-fat diet-induced adipocyte hypertrophy and insulin resistance.

    Mol Cell 1999, 4:597-609. PubMed Abstract | Publisher Full Text OpenURL

  36. Suzawa M, Takada I, Yanagisawa J, Ohtake F, Ogawa S, Yamauchi T, Kadowaki T, Takeuchi Y, Shibuya H, Gotoh Y, Matsumoto K, Kato S: Cytokines suppress adipogenesis and PPAR-gamma function through the TAK1/TAB1/NIK cascade.

    Nat Cell Biol 2003, 5:224-230. PubMed Abstract | Publisher Full Text OpenURL

  37. Weih F, Carrasco D, Bravo R: Constitutive and inducible Rel/NF-kappa B activities in mouse thymus and spleen.

    Oncogene 1994, 9:3289-3297. PubMed Abstract OpenURL

  38. Schreiber E, Matthias P, Muller MM, Schaffner W: Rapid detection of octamer binding proteins with 'mini-extracts', prepared from a small number of cells.

    Nucleic Acids Res 1989, 17:6419. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing.

    J R Stat Soc Ser B 1995, 57:289-300. OpenURL

  40. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers.

    Methods Mol Biol 2000, 132:365-386. PubMed Abstract OpenURL

  41. Kutsch S, Degrandi D, Pfeffer K: Immediate lymphotoxin beta receptor-mediated transcriptional response in host defense against L. monocytogenes.

    Immunobiology 2008, 213:353-366. PubMed Abstract | Publisher Full Text OpenURL