Research | Open | Published:
A roadmap of constitutive NF-κB activity in Hodgkin lymphoma: Dominant roles of p50 and p52 revealed by genome-wide analyses
Genome Medicinevolume 8, Article number: 28 (2016)
NF-κB is widely involved in lymphoid malignancies; however, the functional roles and specific transcriptomes of NF-κB dimers with distinct subunit compositions have been unclear.
Using combined ChIP-sequencing and microarray analyses, we determined the cistromes and target gene signatures of canonical and non-canonical NF-κB species in Hodgkin lymphoma (HL) cells.
We found that the various NF-κB subunits are recruited to regions with redundant κB motifs in a large number of genes. Yet canonical and non-canonical NF-κB dimers up- and downregulate gene sets that are both distinct and overlapping, and are associated with diverse biological functions. p50 and p52 are formed through NIK-dependent p105 and p100 precursor processing in HL cells and are the predominant DNA binding subunits. Logistic regression analyses of combinations of the p50, p52, RelA, and RelB subunits in binding regions that have been assigned to genes they regulate reveal a cross-contribution of p52 and p50 to canonical and non-canonical transcriptomes. These analyses also indicate that the subunit occupancy pattern of NF-κB binding regions and their distance from the genes they regulate are determinants of gene activation versus repression. The pathway-specific signatures of activated and repressed genes distinguish HL from other NF-κB-associated lymphoid malignancies and inversely correlate with gene expression patterns in normal germinal center B cells, which are presumed to be the precursors of HL cells.
We provide insights that are relevant for lymphomas with constitutive NF-κB activation and generally for the decoding of the mechanisms of differential gene regulation through canonical and non-canonical NF-κB signaling.
The NF-κB family of transcription factors (TFs) includes RelA (p65), RelB, c-Rel, and the precursor proteins p105 and p100, which are processed to p50 and p52, respectively. These subunits can form various homo- and heterodimers which regulate gene expression by binding to κB sites in the promoters and enhancers of genes in a wide range of mammalian cell types. NF-κB activation can be achieved through the canonical and non-canonical signaling pathways. The canonical pathway involves the degradation of small IκB proteins and the nuclear translocation of p50-RelA. The non-canonical pathway depends on the activation of NF-κB-inducing kinase (NIK) and IκB kinase α (IKKα), resulting in the induction of NF-κB dimers containing p52 and/or RelB [1–3]. Our recent work, however, has shown that non-canonical signaling is more complex. Stimulation through LTβR results in the parallel production of p50 and p52, and both are important to assure the full range of signal-induced gene expression .
Activated NF-κB plays a major role in cell cycle progression and survival, immune responses, and development . To control such a broad range of physiological processes, the NF-κB response relies on diverse and complex mechanisms that guarantee selectivity by limiting the activation of only subsets of potential NF-κB target genes. The regulation of selected sets of genes by distinct NF-κB subunits may depend on chromatin accessibility, differential expression, dimer formation, differential DNA recognition, and transactivation mechanisms . The diverse phenotypes of mice that lack single NF-κB subunits imply that different dimers regulate distinct sets of genes . However, while previous studies have mainly focused on RelA and c-Rel and revealed considerable redundancy in gene expression regulation , little is known about the dependence of the expression of specific genes on canonical versus non-canonical NF-κB dimers. Previous chromatin immunoprecipitation (ChIP)-based genome-wide NF-κB analyses have mostly investigated canonical NF-κB and been focused on RelA [9–11], or analyzed NF-κB subunits in B cells expressing LMP1 . Nonetheless, to our knowledge, the relation between NF-κB cistromes and subunit-dependent gene regulation has not yet been explored in systematic global analyses.
Understanding how particular subunit combinations control target gene expression is especially important for cancer biology, where the development and maintenance of various solid tumors and lymphoid malignancies are driven by NF-κB [13, 14]. Within B-lymphoid malignancies, constitutive NF-κB activation has been identified in classical Hodgkin lymphoma (HL) and - to various degrees - in non-Hodgkin lymphomas including diffuse large B-cell lymphoma, mucosa-associated lymphoid tissue lymphoma, primary effusion lymphoma, adult T-cell lymphoma/leukemia, and multiple myeloma [13, 15]. The deregulation of NF-κB in these malignancies is caused by diverse mechanisms. Its activation in HL is thought to be associated with multiple genetic alterations and deregulated autocrine and/or paracrine signaling, leading to constitutive IKK activity . The presence of activated RelA and of NIK in primary HL biopsies as well as in HL cell lines [17, 18] indicates a parallel constitutive activation of canonical and non-canonical IKK signaling. This makes HL an ideal model system to study the mechanisms of differential gene regulation by these two pathways.
Here we dissected the composition of activated NF-κB subunits, their cistromes, and p50/RelA- and p52/RelB-specific transcriptomes and functions in HL cells. Our results revealed that NIK promotes the formation of both p50 and p52, which account for most NF-κB activity in HL cells. Major gene expression patterns and cell survival were dependent on the non-canonical pathway. However, each of the pathways regulated both unique and overlapping transcriptomes, and these signatures were consistent with the differential gene expression patterns seen in HL biopsies. A regression analysis of binding patterns and regulation indicates a cross-contribution of p50 and p52 to canonical and non-canonical NF-κB pathways. Our study provides the basis for understanding how canonical and non-canonical NF-κB pathways regulate gene expression in tumor cells.
HL cell lines L428 (ACC 197), L1236 (ACC 530), KM-H2 (ACC 8), L591 (ACC 602), L540 (ACC 72), and HDLM-2 (ACC 17) were purchased from DSMZ (www.dsmz.de). All cell lines were cultivated in RPMI medium 1640 (Gibco) supplemented with inactivated FBS, according to the instructions of the provider.
Immunohistochemical detections of p105/p50 and p100/p52 were performed as described . A total of 20 biopsies were used, from which 20 samples were stained for p105/p50 (#3035, Cell Signaling), and 18 were stained for p100/p52 (HPA008422, Sigma).
Analysis of ChIP-seq data
Short reads from the Illumina sequencer were aligned using Bowtie, retaining only uniquely mapping reads. Quality control was performed according to the ENCODE guidelines . All our samples passed the stringent QC thresholds (see ‘Quality control of ChIP-seq libraries’ in Additional file 1: Supplemental Experimental Procedures). Duplicated reads that are likely to be PCR artefacts were removed using samtools. In order to be able to compare sequencing libraries with different total read counts we used random sampling to reduce all libraries to a size of 10 million reads. The genome wide distribution of binding regions was determined using MACS . All ChIP-seq experiments were conducted in two biological replicates, so we investigated the reproducibility using the IDR method. After IDR analysis, reads from biological replicates were merged and jointly analyzed with MACS. Final peak calls were defined as regions with IDR <0.05 and MACS false discovery rate (FDR) <0.05. See also Additional file 1: Supplemental Experimental Procedures for details about peak calling and combinatorial analysis of NF-κB subunits.
Analysis of gene expression data
Experiments for the knockdown of the two subunit combinations (RelA/p50 and RelB/p52) were performed and analyzed independently. Each experiment was performed in three biological replicates. For Experiment #1, each biological replicate was composed of four samples (sample 1 = siControl #1; sample 2 = siControl#2; sample 3 = siRELA #1 + siNFKB1 #1; sample 4 = siRELA #2 + siNFKB1 #2). For Experiment #2, each biological replicate was composed of an additional four samples (sample 5 = siControl #1; sample 6 = siControl#2; sample 7 = siRELB #1 + siNFKB2 #1; sample 8 = siRELB #2 + siNFKB2 #2). For the data analysis, all 24 samples were separated into four groups (Group1 = all replicates from samples 1 and 2, Group 2 = all replicates from samples 3 and 4, Group 3 = all replicates from samples 5 and 6, and Group 4 = all replicates from samples 7 and 8). Differential expression was determined using a two-sample t-test and P values adjusted for multiple testing using the Benjamini–Hochberg method. A gene was called differentially expressed if at least one transcript cluster assigned to the gene was called differential (FDR <0.05) and showed at least 10 % expression difference between Groups 1 and 2 or 3 and 4. See Additional file 1: Supplemental Experimental Procedures for details.
Integration of target genes with gene expression in human lymphomas
Human lymphoma gene expression data were obtained from GEO (GSE12453). Genes with low average expression (Affymetrix hybridization signal <6) were removed. Differential expression was determined between HL and all remaining samples using the same method as described above. For details see Additional file 1.
Predominance of p50 and p52 and non-canonical IKK signaling in the constitutive NF-κB activity in HL
To test the relative contribution of canonical and non-canonical NF-κB dimers in HL cells, we depleted NIK, a central component of non-canonical NF-κB signaling. Notably, NIK depletion not only affected processing of p100 but also of p105 (Fig. 1a). Ablation of NIK caused a severe reduction of C-terminal p100 and p105 phosphorylation, accumulation of the precursors, and decreased generation of their products p52 and p50 (Fig. 1a). This is similar to the NIK-dependent coupled generation of p50 and p52 that we recently demonstrated for non-canonical LTβ signaling in fibroblasts . The p50 and p52 subunits are pivotal to the total constitutive NF-κB activity in HL cells, since RNA-mediated knockdown (KD) of both NFKB1 (encoding p105/p50) and NFKB2 (p100/p52) strongly affected global NF-κB DNA-binding, while depletion of the single NF-κB subunits, including RelA and RelB, did not (Fig. 1a). Consistent with its role in p50 and p52 production, NIK depletion had a similar effect on total NF-κB DNA-binding (Additional file 2: Figure S1A). Of note, in addition to slower migrating heteromeric species in the DNA binding assays, abundant activity migrated at the position of p50 and p52 dimers (Fig. 1a).
We next investigated the cytoplasmic-nuclear distribution of the NF-κB subunits in various HL cell lines (Fig. 1b). A significant fraction of p50 and p52 but only a low amount of RelA was present in the nucleus. Nuclear RelB was ubiquitously detected with variable intensities, whereas differing concentrations of nuclear c-Rel were observed throughout the cell lines. To evaluate the canonical and non-canonical NF-κB dimer composition, p50 and p52 were immunoprecipitated from nuclear extracts (Additional file 2: Figure S1C). We could readily detect p52-RelB and to a lesser extent p50-RelB dimers in all cell lines. However, p50-RelA dimers were present at variable lower amounts, and p50-c-Rel dimers were only abundant in L428 and KM-H2 cells, correlating with their generally restricted nuclear abundance. Furthermore, p50-p52 heterodimers were detected.
Nuclear enrichment of p50 and p52 was also evident in primary HL samples. A prominent nuclear staining was seen in 60 % of patients for p50 and in 83 % for p52 (Fig. 1c and Additional file 3: Table S1). Since the N-terminal epitope-directed antibodies recognize cytoplasmic precursors and the products, an excess of nuclear p50 and p52 is indicated in these cells. Altogether, these results demonstrate that non-canonical, NIK-driven IKK signaling predominantly accounts for constitutive NF-κB activity in HL cells, which is strongly dominated by production of dimers that contain p50 or p52. Compared to this, canonical p50-RelA heterodimers are under-represented.
Constitutive activation of canonical and non-canonical NF-κB species generates largely overlapping cistromes
To investigate selective functional roles of the RelA, RelB, p50, and p52 subunits and of the prototypic p50-RelA and p52-RelB heterodimers in the gene expression repertoire of HL cells, we first performed ChIP-seq analyses of the four subunits to determine the genome wide distribution of their binding sites. This should also reveal if distinct regions and sequences are preferentially bound by canonical and non-canonical NF-κB species. The number of enriched regions greatly varied, up to 20-fold comparing RelA and p50 (Additional file 2: Figure S2A and Additional file 4: Table S2), reflecting the relative nuclear amounts of the respective NF-κB subunit (Fig. 1b). For each of the NF-κB datasets individual regions were selected and validated by qPCR (Additional file 2: Figure S2B). In each dataset the maximum frequency of ChIP-seq regions was observed around the TSS (Additional file 2: Figure S2C). The genomic locations of ChIP-seq regions in relation to annotated genes are representatively shown for p52 and p50. About half of the regions were located within transcribed gene regions, with a high percentage within introns (Additional file 2: Figure S2E).
Comparison of the cistromes obtained for the single NF-κB subunits revealed multiple intersections (Fig. 2a). RelA shares more than 94 % and RelB more than 85 % of their ChIP-seq regions with at least one other NF-κB subunit. Binding regions for p50 and p52 overlapped by about 70 %. Remarkably, a large fraction (41 %) of binding regions recruits p50 or p52 alone. Occupancy profiling reveals that combinatorial NF-κB subunit recruitment can be grouped into eight clusters (Fig. 2b) with a most prominent contribution of p50 and p52, either alone or in combination. Interestingly, RelB is most frequently seen in combination with p50 and p52, rather than with p52 alone, while RelA occupancy is least frequent. In many cases the tag pattern for all NF-κB subunits is very similar and binding regions are shared, even though the enrichment for RelA is weaker than that for p50, p52 and RelB (for example, NFKBIA gene). In other cases, the tag patterns differ and only subsets of subunits are present in selected ChIP-seq regions (CDH1 gene, Fig. 2c).
Redundancy in DNA binding was also indicated by de novo motif searches using MEME  with the ChIP-seq datasets for the different NF-κB subunits. In all datasets, an NF-κB motif that resembled the known JASPAR database entries (MA0101; MA0107; MA0105) was found with very high frequency (Fig. 3a and b). To assess whether p50 and p52 dimers recognize distinct DNA motifs, MEME de novo search was performed with ChIP-seq regions that selectively recruited p50 or p52. Again, the NF-κB motifs found for each factor (Additional file 2: Figure S3) did not grossly differ from motifs found within regions that recruit more than one subunit (Fig. 3b). However, p50-selective motifs contained a higher conservation of a G in the first position compared to p52-selective motifs, indicating a putative contribution of highly symmetric binding sites for p50 homodimer recruitment.
We next searched for co-occurring binding sites of potential cooperating TFs within the subunit-specific ChIP-seq regions. Apart from the expected, statistically significant occurrence of the NF-κB motif in all datasets, we found an enrichment of AP-1, IRF, ETS, FOX, and CTCF motifs (Additional file 5: Table S3). This is in agreement with our finding that protein-occupied AP-1 and IRF motifs occur closely spaced to NF-κB motifs in global DNAse I hypersensitive site (DHS) mapping in HL cells . Furthermore, we showed before that constitutively activated AP-1 and IRF5 in HL cells functionally synergize with NF-κB [19, 23].
The cistromes of NF-κB subunits in HL cells differ significantly from the binding landscapes previously reported for EBV-transformed lymphoblastoid B cells (LCL), where the membrane protein LMP1 constitutively activates the canonical and non-canonical NF-κB pathways . In particular, the relative recruitment of RelA and RelB is notably less in HL cells compared to LCLs, while p50 and p52 are prevalent in HL cells. The ChIP regions of NF-κB subunits in HL cells overlap with those in LCL in average by 54 % (Additional file 6: Table S4). The binding regions that overlap in both cell types are summarized in Additional file 7: Table S5.
We also compared the combinatorial binding pattern of NF-κB subunits in ChIP-seq regions between HL and LCL (Additional file 8: Table S6). Both cell types strongly differed in combinatorial subunit binding, roughly in accordance with the differences in frequencies of subunit recruitment noted above. Interestingly, patterns exclusively containing the prototypic heterodimers p50-RelA or p52-RelB are rare in both cell types, about 2 %. The most frequent subunit combinations in HL cells contain p50 and/or p52, in contrast to RelA and RelB in LCL. Furthermore, as expected, the cistromes of the NF-κB species in HL cells and the DHS sites that we mapped in these cells  largely overlapped. Similarly, binding regions and DHS sites reveal an equivalent extent of overlap in LCL (Additional file 6: Table S4).
Thus, globally, the different NF-κB subunit recruitment in HL and LCL cells is in part dictated by the only partially overlapping open chromatin landscape of the two cell types.
Distinct biological processes are differentially associated with genes that are up- or downregulated by p50/RelA and p52/RelB
To determine the transcriptomes regulated by canonical and/or non-canonical NF-κB dimers, we depleted the subunit combinations p50/RelA and p52/RelB (Additional file 2: Figure S4A). Concentrations of siRNAs and treatment time schedules were optimized for each NF-κB subunit for efficient depletion (Additional file 2: Figure S4B). Samples were collected 1 day after the end of the respective siRNA treatments and one day before the initial onset of apoptosis. We performed global microarray expression analyses and determined direct targets by matching gene regulation with ChIP-seq binding data (Additional file 9: Table S7). Almost fivefold more genes were regulated by p52/RelB than by p50/RelA. While the majority of genes were transcriptionally activated, surprisingly, 22 % of p50/RelA and 36 % of p52/RelB regulated genes were down-regulated (Fig. 4a and b, top panels, Additional file 9: Table S7).
Alterations of gene regulation could have been caused indirectly by depletion of p100 and p105, which can sequester other NF-κB subunits in the cytoplasm and consequent liberation of bound Rel factors. However, no change in the nuclear patterns of RelA, RelB, or c-Rel was observed (Additional file 2: Figure S4C), suggesting that the altered expression of genes is a direct consequence of the KD of NF-κB subunits and not of a loss of the IκB function of the precursors.
For up- and for downregulated genes, the p50/RelA- and p52/RelB-dependent sets in each case partially overlapped, but also revealed genes that were only controlled by either p50/RelA or p52/RelB (Fig. 4a and b, top panels). These datasets contained a considerable number of novel direct NF-κB targets, but also previously known NF-κB target genes, as expected (Additional file 9: Table S7). A number of genes that are differentially up- or downregulated by p50/RelA- and/or p52/RelB have a high biological significance for HL (for a selection, see Fig. 4c) and encode cytokines and chemokines responsible for clinical symptoms and TF networks implicated in this disease .
We next categorized genes directly regulated by p50/RelA and/or p52/RelB according to their activation or repression status and grouped them by related biological processes, using Gene Ontology (GO) terms. The set of genes upregulated by p50/RelA was predominantly enriched in immune system processes. These genes are also important for cell activation, cell communication, and regulation of programmed cell death. The set of genes upregulated by p52/RelB is also enriched in these categories, however, covering a broader spectrum of biological functions. The cell death categories are shared by both subunit combinations (Fig. 4a) and include known NF-κB induced anti-apoptotic genes such as CFLAR and BCL2L1 (Additional file 9: Table S7).
The gene sets that were downregulated by p50/RelA or p52/RelB (Additional file 9: Table S7) differed greatly in the enriched GO terms (Fig. 4b, Additional file 10: Table S8). Genes that were selectively downregulated by p50/RelA, like BCL6, are mainly related to general cellular developmental and differentiation processes. The genes exclusively downregulated by p52/RelB on the other hand, including FOXP1 and LEF1, are highly enriched in terms such as hematopoietic and lymphoid organ development, negative regulation of cell adhesion, and positive regulation of cell differentiation (Fig. 4b and c, Additional file 9: Table S7).
Collectively, these results show that p52/RelB controls the majority of the NF-κB regulated genes and especially represses genes involved in differentiation processes. However, both p50/RelA and p52/RelB control genes that are involved in modulation of immune system responses as well as cell survival. Both pathways are differentially engaged in the regulation of genes with previously documented high or low expression in HL (Fig. 4c).
Correlation of combinatorial binding patterns of NF-κB subunits and gene expression regulation
By fitting logistic regression models (see Additional file 1: Supplemental Experimental Procedures for details) we analyzed if specific combinations of NF-κB subunits that were determined within the individual ChIP-seq regions could be used to predict modes of gene expression regulation by canonical and non-canonical dimers. The best model parameters were estimated, also considering the available L1236-DHS data and the genomic 3D architecture in LCL [19, 24], as outlined in Additional file 2: Figure S5A. We found that binding of both, p50 and p52, either alone or in combination with RelB or RelB plus RelA was significantly associated with regulation (that is, activation or repression) of gene expression by canonical and non-canonical dimers. Specifically, the genes regulated by p52/RelB revealed a high preference for binding of p50 along with p52, with or without RelB and were independent of proximal or distal localization of the region relative to the TSS (Fig. 5a). Genes sensitive to p50/RelA knockdown also revealed the same preference for p50/p52/RelB but the most significant subunit combination was p50/p52/RelA/RelB (Fig. 5a). In this case all regions were restricted to proximal positions. Likewise, proximal binding of all four subunits, thus including RelA, was highly predictive for genes regulated by canonical compared to non-canonical dimers (Additional file 2: Figure S5B). Furthermore, we identified the preferential NF-κB subunit combinations associated with gene activation or repression. p50/p52/RelB and p50/p52/RelA/RelB were significantly associated with activated genes (Additional file 2: Figure S5C). Interestingly, although with lesser confidence, binding of p50/p52 scored highest with repressed genes when in distal positions, but with activation, when at proximal positions (Additional file 2: Figure S5C and D). The compositional and positional correlation of NF-κB subunit binding pattern with gene activation and repression is schematically summarized in Fig. 5b. With exception of p50/p52 (binary code ‘1100’ in Fig. 5a), the identified high confidence binding patterns did not simply correlate with their frequency of occurrence in the ChIP-seq datasets (Additional file 8: Table S6). Importantly, exclusive recruitment of prototypic dimers (p50-RelA or p52-RelB) was not significantly associated with gene regulation (‘1010’ and ‘0101’ in Fig. 5a and Additional file 2: Figures S5B and C).
Taken together, the logistic regression analyses reveal that the specificity of NF-κB subunit-selective gene regulation is not determined by the exclusive recruitment of the targeted prototypic dimers, but by more complex combinations of NF-κB subunits, involving a parallel contribution of p50 and p52. The most active combinations (1111, 1101, and 1100) are part of clusters C4 and C8 (Fig. 2b), revealing strongest and equal p50 and p52 recruitment.
Cell survival of HL cells is predominantly controlled by the non-canonical NF-κB pathway
We next investigated the relative impact of p50-RelA and p52-RelB, both of which activated genes associated with cell death/apoptosis (Fig. 4a), on cell viability. Close monitoring of siRNA-treated HL cells revealed a 30 % to 40 % decrease of cell viability upon p50/RelA KD and up to 95 % reduction in viability upon p52/RelB KD (Fig. 6a and Additional file 2: Figures S6A and B). The cleavage of initiator caspases 8 and 9, indicating apoptosis, were predominantly increased upon KD of non-canonical and canonical subunits, respectively, while effector caspase 3 was equally activated in both groups (Fig. 6b). Interestingly, we observed a correlation between the activation of initiator caspases and the protein amounts of negative regulators of apoptosis such as c-FLIP and Bcl- XL (Fig. 6b). c-FLIP downregulation sensitizes HL cells to death receptor-induced apoptosis . CFLAR mRNA expression was downregulated by the KD of both combinations of NF-κB subunits, but a decrease in protein levels was only observed upon the KD of p52/RelB, which was concomitant to strong caspase 8 activation (see Fig. 6b). Overexpression of Bcl-XL blocks the intrinsic apoptosis pathway in many cell types [26, 27] and its protein level was downregulated by the KD of p52/RelB as well as of p50/RelA (Fig. 6b). To investigate the possible connection between the downregulation of Bcl-XL caused by the KD of NF-κB dimers and apoptosis, we tested the effect of BCL2L1 knockdown. Following this, we observed an up to 30 % decrease in cell viability in HL cell lines, along with activation of caspases 9 and 3 (Additional file 2: Figure S6C and D). This suggests that Bcl-XL blocks the activation of the intrinsic apoptosis pathway in HL cells and is partially responsible for the NF-κB-dependent control of HL cell survival.
Taken together, our results indicate that the activation of p52-RelB is essential for the survival of HL cells through the upregulation of proteins that block the intrinsic and especially the extrinsic apoptosis pathway (for example, c-Flip), whereas p50-RelA also contributes to survival by the upregulation of Bcl-XL and subsequent inhibition of the intrinsic apoptosis pathway.
p50-RelA and p52-RelB gene signatures distinguish classical Hodgkin lymphoma from other lymphomas and normal B cells
In order to compare the gene sets regulated by the p50/RelA and p52/RelB subunit combinations, reflecting canonical and non-canonical NF-κB signaling in HL cells in vitro, with primary human lymphoma, we performed a retrospective analysis of microarray datasets generated from micro-dissected Hodgkin/Reed Sternberg (HRS) cells of classical HL, other malignant B cell lymphomas, and normal B cells at various stages of differentiation [28, 29]. We identified 123 genes that were direct NF-κB targets in L1236 cells and differentially regulated in primary HRS cells when compared to other malignant and normal B cells (Fig. 7; Additional file 11: Table S9). These genes were classified into three groups, genes regulated by canonical, non-canonical, or both NF-κB pathways. Selected examples of regulated genes were validated by RT-qPCR (Additional file 2: Figure S7).
Some genes that were activated by NF-κB in L1236 cells and identified as classical HL signatures were also upregulated in other B-cell lymphomas, where NF-κB activation has been observed previously, such as nodular lymphocyte-predominant HL (NLPHL) , T cell-rich B cell lymphoma (TCRBL)  and DLBCL  (Fig. 7). These genes were downregulated in normal B cells and also in Burkitt’s lymphoma (BL), where NF-κB is not activated [31, 32].
Interestingly, genes that were repressed in L1236 by NF-κB and downregulated in classical HL, were upregulated in normal germinal center (GC) B cells (centroblasts and centrocytes) (Fig. 7). Thus, NF-κB might contribute to the characteristic de-differentiation of HL cells, which are derived from germinal center B cells .
Our results show that gene sets that are differentially up- or downregulated by canonical and non-canonical NF-κB signaling in L1236 cells, which were used as a model system, are differentially expressed in a similar manner in primary HL tumor cells compared to other B cell malignancies.
In various lymphoid malignancies, oncogenic constitutive NF-κB activity is determined to varying extents by canonical and non-canonical signaling, as indicated by mutations in upstream pathway components . However, the specific contributions of constitutive canonical versus non-canonical NF-κB species to cancer biology remain poorly understood. Here, we analyzed the genomic landscapes of p50, p52, RelA, and RelB in HL cells by ChIP-sequencing and combined these data with the genes regulated by these subunits.
We could demonstrate that NF-κB dependent gene expression in HL cells is largely driven by the non-canonical pathway. The central kinase of this pathway, NIK, is known to be overexpressed in HL cell lines and primary cells [18, 33]. We show that NIK triggers not only p100, but also p105 precursor processing to p52 and p50. This finding is consistent with NIK-dependent processing of both precursors in LTβ-stimulated cells, which we have recently characterized . We also noted a surprisingly low translocation of RelA. Interestingly, low levels of nuclear translocation/activation of RelA and the parallel increase of p50 and p52 in HL cells can be observed at late time points upon stimulation with TWEAK or CD40 [34, 35], suggesting that this might be a common feature of cells with stimulated NIK pathways. A similar pattern is also seen in peripheral T cell lymphoma  and multiple myeloma .
Constitutive non-canonical signaling causes abundant nuclear activities of p50 and p52 and each of them is recruited to more than 10,000 binding regions. This number of binding regions for NF-κB subunits in HL cells is comparable to the number determined in prior studies using lymphoblastoid B cells [10, 12]. However, we show that only a small fraction of these genes are transcriptionally regulated. In addition, we observed a lower contribution of RelA recruitment and of p50-RelA to the control of NF-κB dependent gene expression.
We asked for non-redundant functions, that is, a functional ‘addictiveness’ of regulated genes to particular subunit combinations by knocking down p50/RelA or p52/RelB. Several established NF-κB targets, which were previously identified , such as CCR7, CCND2, FAS, JUNB, or VCAM1, were bound by canonical and non-canonical subunits, but not regulated upon knockdown. This fact can be explained by redundant functions of NF-κB subunits or the applied stringency of the parameters used in our gene expression analysis. From a total of 501 activated genes, 63 were activated by both, p50/RelA and p52/RelB, revealing additive or synergistic regulation. This group of target genes comprised BIRC2 (encoding cIAP1), BCL2L1 (Bcl-XL), CFLAR (c-FLIP), CD40, IRF5, STAT5A, CSF2, CCL22, and HEY1, all of which are known to be expressed at high levels in HL cells [19, 38, 39].
Surprisingly, one-third of all genes regulated by NF-κB were suppressed. The observed suppression of genes by both dimer combinations may be caused by subunit exchange mechanisms. p50 or p52 dimers that lack transactivation domains and do not interact with transcriptional co-activators, may prevent access of activating RelA or RelB containing heterodimers. Furthermore, p50 and p52 may recruit repressor complexes, including histone deacetylases [1, 40, 41], in a context-dependent manner to suppress transcription. Indeed, our logistic regression models indicate that a combined recruitment of p50 and p52 significantly correlated with gene expression.
In addition to the separate biological processes identified by GO analyses for canonical and non-canonical subunit combinations, overlapping functions were implicated in cell death regulation. Indeed, we were able to demonstrate that apoptosis was most prominently inhibited by p52/RelB that primarily affected the extrinsic apoptosis pathway, while p50/RelA contributed to cell death protection by interfering with the intrinsic pathway. Thus, HL cells engage in both apoptosis pathways to block cell death, which is similar to mature murine B cells, where both are essential for maintenance and in part are interchangeable .
Another interesting observation was the fact that genes downregulated by NF-κB in HL cells were upregulated in normal germinal center (GC) B cells (Fig. 7). Among these genes were BCL6 and PAX5 as well as other important regulators of B cell development, such as FOXP1 , JAK1 , and PTPN6 , suggesting that NF-κB signaling may contribute to the known de-differentiation process of HL cells, which are derived from GC cells . Activation of the canonical NF-κB subunits is important for the maintenance and differentiation of GC B cells , but the role of the non-canonical subunits in this stage of the B-cell development remains unclear. Our signatures of suppressed genes will thus be an interesting source for future studies of HL lymphomagenesis.
We also show that the differential gene regulation programs of canonical and non-canonical NF-κBs were neither explainable by differences in their DNA sequence elements, nor by exclusive recruitment of dimer species to the genes they regulate. The MEME analysis of the ChIP-seq datasets for the different NF-κB subunits revealed high redundancy in the conserved motifs. Not taking into account a very weak first position in motifs related to p50, p52, and RelB, we have found in this study very similar 11-bp motifs for each factor (p50, p52, RelA, and RelB). Yet, we could notice a more specific C at the final position for p50 and p52. A similar redundancy of motifs was also found in the recruitment patterns in lymphoblastoid cells .
The motifs found, for example, in regions selective for p50 or p52 did not considerably differ from those associated with additional binding of RelA or RelB, which indicates that other constraints such as protein-protein interactions are important in determining recruitment. Accordingly, the differential interaction of NF-κB with other TFs may establish functional gene selectivity. Motifs for bZIP, ETS, FOX, IRF and CTCF families were significantly over-represented in binding regions of all four NF-κB subunits. Interestingly, these motifs belong to a dense cluster of TF binding sites surrounding cohesin anchor sites, which may contribute to the regulation of nearby genes . Indeed, we have shown before that AP1 and IRF5, which are recruited to bZIP and IRF motifs, functionally synergize with NF-κB [19, 23]. This synergism may control chromatin accessibility, as was shown for AP1 in other systems [48, 49].
It remains a difficult task to correlate binding of a given subunit to an individual gene with the regulation of this gene, since regulation may be conferred by other remote regions. We therefore developed logistic regression models to dissect the overall correlation between transcriptional regulation and recruitment of combinations of NF-κB subunits. This analysis provided insight into the functional subunit organization, such as a cross-contribution of p50 to gene regulation by non-canonical p52-containing dimers and vice versa. It ruled out that exclusive recruitment of the prototypic dimers p50-RelA and p52-RelB alone is significantly associated with gene regulation. This analysis also clearly indicated that RelA-containing binding patterns correlated with function only from proximal, but not distal localizations, while p50, p52, and RelB were distance independent. Consequently, the distinct NF-κB species do have different preferences for action through enhancers or other remote positions versus promoters. Thus, the models developed here using binary codes for NF-κB subunits allowed us to correlate gene expression regulation with specific combinatorial subunit occupancy.
In summary, in this study we have dissected the NF-κB subunit dependent gene expression pattern in HL cells, providing datasets that will be the basis for future comparative analyses once equivalent functional investigations for other NF-κB-associated malignancies become available. The understanding of the differential contributions of canonical and non-canonical NF-κB species to transcriptional networks and of their potential synergisms of action on targets will be crucial for establishing refined functional gene signatures. The latter may provide important information for future stratification of NF-κB associated malignancies and for the evaluation of novel NF-κB pathway directed drugs for individualized therapies.
Ethics approval and consent to participate
The use of human material was approved by the Ethik-Kommission of the Charité, Universitätsmedizin Berlin (EA4/104/11) and was performed in accordance with the Declaration of Helsinki. Informed consent was obtained from all subjects, where required.
Availability of data and materials
Geo accession numbers for ChIP-seq and gene expression data are GSE63736 and GSE64232, respectively.
chromatin immunoprecipitation followed by sequencing
DNAse I hypersensitive site
diffuse large B-cell lymphoma
gene expression omnibus
IκB kinase α
inhibitor of κB
lymphoblastoid B cells
multiple EM for motif elicitation
nuclear factor κB
nodular lymphocyte-predominant HL
- p50/RelA or p52/RelB:
slash refers to combined or single function of subunits
- p50-RelA or p52-RelB:
hyphen refers to the indicated heterodimer
quantitative real-time PCR
reverse transcription quantitative PCR
T cell-rich B cell lymphoma
transcription start site
TNF-related weak inducer of apoptosis
Hayden MS, Ghosh S. Shared principles in NF-kappaB signaling. Cell. 2008;132(3):344–62.
Scheidereit C. IkappaB kinase complexes: gateways to NF-kappaB activation and transcription. Oncogene. 2006;25(51):6685–705.
Vallabhapurapu S, Karin M. Regulation and function of NF-kappaB transcription factors in the immune system. Annu Rev Immunol. 2009;27:693–733.
Yilmaz ZB, Kofahl B, Beaudette P, Baum K, Ipenberg I, Weih F, et al. Quantitative dissection and modeling of the NF-kappaB p100-p105 module reveals interdependent precursor proteolysis. Cell Rep. 2014;9(5):1756–69.
Hayden MS, Ghosh S. NF-kappaB, the first quarter-century: remarkable progress and outstanding questions. Genes Dev. 2012;26(3):203–34.
Smale ST. Hierarchies of NF-kappaB target-gene regulation. Nat Immunol. 2011;12(8):689–94.
Gerondakis S, Banerjee A, Grigoriadis G, Vasanthakumar A, Gugasyan R, Sidwell T, et al. NF-kappaB subunit specificity in hemopoiesis. Immunol Rev. 2012;246(1):272–85.
Smale ST. Dimer-specific regulatory mechanisms within the NF-kappaB family of transcription factors. Immunol Rev. 2012;246(1):193–204.
Jin F, Li Y, Dixon JR, Selvaraj S, Ye Z, Lee AY, et al. A high-resolution map of the three-dimensional chromatin interactome in human cells. Nature. 2013;503(7475):290–4.
Kasowski M, Grubert F, Heffelfinger C, Hariharan M, Asabere A, Waszak SM, et al. Variation in transcription factor binding among humans. Science. 2010;328(5975):232–5.
Lim CA, Yao F, Wong JJ, George J, Xu H, Chiu KP, et al. Genome-wide mapping of RELA(p65) binding identifies E2F1 as a transcriptional activator recruited by NF-kappaB upon TLR4 activation. Mol Cell. 2007;27(4):622–35.
Zhao B, Barrera LA, Ersing I, Willox B, Schmidt SC, Greenfeld H, et al. The NF-kappaB genomic landscape in lymphoblastoid B cells. Cell Rep. 2014;8(5):1595–606.
Lim KH, Yang Y, Staudt LM. Pathogenetic importance and therapeutic implications of NF-kappaB in lymphoid malignancies. Immunol Rev. 2012;246(1):359–78.
DiDonato JA, Mercurio F, Karin M. NF-kappaB and the link between inflammation and cancer. Immunol Rev. 2012;246(1):379–400.
Jost PJ, Ruland J. Aberrant NF-kappaB signaling in lymphoma: mechanisms, consequences, and therapeutic implications. Blood. 2007;109(7):2700–7.
Kuppers R. The biology of Hodgkin’s lymphoma. Nat Rev Cancer. 2009;9(1):15–27.
Bargou RC, Emmerich F, Krappmann D, Bommert K, Mapara MY, Arnold W, et al. Constitutive nuclear factor-kappaB-RelA activation is required for proliferation and survival of Hodgkin’s disease tumor cells. J Clin Invest. 1997;100(12):2961–9.
Ranuncolo SM, Pittaluga S, Evbuomwan MO, Jaffe ES, Lewis BA. Hodgkin lymphoma requires stabilized NIK and constitutive RelB expression for survival. Blood. 2012;120(18):3756–63.
Kreher S, Bouhlel MA, Cauchy P, Lamprecht B, Li S, Grau M, et al. Mapping of transcription factor motifs in active chromatin identifies IRF5 as key regulator in classical Hodgkin lymphoma. Proc Natl Acad Sci U S A. 2014;111(42):E4513–22.
Landt SG, Marinov GK, Kundaje A, Kheradpour P, Pauli F, Batzoglou S, et al. ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. Genome Res. 2012;22(9):1813–31.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, et al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(Web Server issue):W202–8.
Mathas S, Hinz M, Anagnostopoulos I, Krappmann D, Lietz A, Jundt F, et al. Aberrantly expressed c-Jun and JunB are a hallmark of Hodgkin lymphoma cells, stimulate proliferation and synergize with NF-kappa B. EMBO J. 2002;21(15):4104–13.
Rao SS, Huntley MH, Durand NC, Stamenova EK, Bochkov ID, Robinson JT, et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014;159(7):1665–80.
Mathas S, Lietz A, Anagnostopoulos I, Hummel F, Wiesner B, Janz M, et al. c-FLIP mediates resistance of Hodgkin/Reed-Sternberg cells to death receptor-induced apoptosis. J Exp Med. 2004;199(8):1041–52.
Boise LH, Gonzalez-Garcia M, Postema CE, Ding L, Lindsten T, Turka LA, et al. bcl-x, a bcl-2-related gene that functions as a dominant regulator of apoptotic cell death. Cell. 1993;74(4):597–608.
Roth KA, Kuan C, Haydar TF, D'Sa-Eipper C, Shindler KS, Zheng TS, et al. Epistatic and independent functions of caspase-3 and Bcl-X(L) in developmental programmed cell death. Proc Natl Acad Sci U S A. 2000;97(1):466–71.
Brune V, Tiacci E, Pfeil I, Doring C, Eckerle S, van Noesel CJ, et al. Origin and pathogenesis of nodular lymphocyte-predominant Hodgkin lymphoma as revealed by global gene expression analysis. J Exp Med. 2008;205(10):2251–68.
Tiacci E, Doring C, Brune V, van Noesel CJ, Klapper W, Mechtersheimer G, et al. Analyzing primary Hodgkin and Reed-Sternberg cells to capture the molecular and cellular pathogenesis of classical Hodgkin lymphoma. Blood. 2012;120(23):4609–20.
Davis RE, Brown KD, Siebenlist U, Staudt LM. Constitutive nuclear factor kappaB activity is required for survival of activated B cell-like diffuse large B cell lymphoma cells. J Exp Med. 2001;194(12):1861–74.
Dave SS, Fu K, Wright GW, Lam LT, Kluin P, Boerma EJ, et al. Molecular diagnosis of Burkitt’s lymphoma. N Engl J Med. 2006;354(23):2431–42.
Hummel M, Bentink S, Berger H, Klapper W, Wessendorf S, Barth TF, et al. A biologic definition of Burkitt’s lymphoma from transcriptional and genomic profiling. N Engl J Med. 2006;354(23):2419–30.
Otto C, Giefing M, Massow A, Vater I, Gesk S, Schlesner M, et al. Genetic lesions of the TRAF3 and MAP3K14 genes in classical Hodgkin lymphoma. Br J Haematol. 2012;157(6):702–8.
Saitoh T, Nakayama M, Nakano H, Yagita H, Yamamoto N, Yamaoka S. TWEAK induces NF-kappaB2 p100 processing and long lasting NF-kappaB activation. J Biol Chem. 2003;278(38):36005–12.
Coope HJ, Atkinson PG, Huhse B, Belich M, Janzen J, Holman MJ, et al. CD40 regulates the processing of NF-kappaB2 p100 to p52. EMBO J. 2002;21(20):5375–85.
Odqvist L, Sanchez-Beato M, Montes-Moreno S, Martin-Sanchez E, Pajares R, Sanchez-Verde L, et al. NIK controls classical and alternative NF-kappaB activation and is necessary for the survival of human T-cell lymphoma cells. Clin Cancer Res. 2013;19(9):2319–30.
Thu YM, Richmond A. NF-kappaB inducing kinase: a key regulator in the immune system and in cancer. Cytokine Growth Factor Rev. 2010;21(4):213–26.
Hinz M, Lemke P, Anagnostopoulos I, Hacker C, Krappmann D, Mathas S, et al. Nuclear factor kappaB-dependent gene expression profiling of Hodgkin’s disease tumor cells, pathogenetic significance, and link to constitutive signal transducer and activator of transcription 5a activity. J Exp Med. 2002;196(5):605–17.
Kochert K, Ullrich K, Kreher S, Aster JC, Kitagawa M, Johrens K, et al. High-level expression of Mastermind-like 2 contributes to aberrant activation of the NOTCH signaling pathway in human lymphomas. Oncogene. 2011;30(15):1831–40.
Zhong H, May MJ, Jimi E, Ghosh S. The phosphorylation status of nuclear NF-kappa B determines its association with CBP/p300 or HDAC-1. Mol Cell. 2002;9(3):625–36.
Kim T, Yoon J, Cho H, Lee WB, Kim J, Song YH, et al. Downregulation of lipopolysaccharide response in Drosophila by negative crosstalk between the AP1 and NF-kappaB signaling modules. Nat Immunol. 2005;6(2):211–8.
Kaileh M, Sen R. NF-kappaB function in B lymphocytes. Immunol Rev. 2012;246(1):254–71.
Hu H, Wang B, Borde M, Nardone J, Maika S, Allred L, et al. Foxp1 is an essential transcriptional regulator of B cell development. Nat Immunol. 2006;7(8):819–26.
Rodig SJ, Meraz MA, White JM, Lampe PA, Riley JK, Arthur CD, et al. Disruption of the Jak1 gene demonstrates obligatory and nonredundant roles of the Jaks in cytokine-induced biologic responses. Cell. 1998;93(3):373–83.
Pao LI, Lam KP, Henderson JM, Kutok JL, Alimzhanov M, Nitschke L, et al. B cell-specific deletion of protein-tyrosine phosphatase Shp1 promotes B-1a cell development and causes systemic autoimmunity. Immunity. 2007;27(1):35–48.
Heise N, De Silva NS, Silva K, Carette A, Simonetti G, Pasparakis M, et al. Germinal center B cell maintenance and differentiation are controlled by distinct NF-kappaB transcription factor subunits. J Exp Med. 2014;211(10):2103–18.
Yan J, Enge M, Whitington T, Dave K, Liu J, Sur I, et al. Transcription factor binding in human cells occurs in dense clusters formed around cohesin anchor sites. Cell. 2013;154(4):801–13.
Biddie SC, John S, Sabo PJ, Thurman RE, Johnson TA, Schiltz RL, et al. Transcription factor AP1 potentiates chromatin accessibility and glucocorticoid receptor binding. Mol Cell. 2011;43(1):145–55.
Sheffield NC, Thurman RE, Song L, Safi A, Stamatoyannopoulos JA, Lenhard B, et al. Patterns of regulatory activity across diverse human cell types predict tissue identity, transcription factor binding, and long-range interactions. Genome Res. 2013;23(5):777–88.
Vogel MJ, Xie L, Guan H, Tooze RM, Maier T, Kostezka U, et al. FOXO1 repression contributes to block of plasma cell differentiation in classical Hodgkin lymphoma. Blood. 2014;124(20):3118–29.
We thank Benjamin Gewurz, Harvard Medical School, Boston, for discussion and for providing LCL peak calling data. We thank Mirjam Feldkamp, Wei Chen lab, and Gabriele Born for ChIP-sequencing and Affymetrix chip hybridization, respectively. We further acknowledge Sabine Jungmann and Inge Krahn for expert technical assistance and Michael Hinz, Z. Buket Yilmaz, and Ruth Schmidt-Ullrich for critical comments on the manuscript. This work was supported in part by BMBF (MSBN) to NH, MA, and CS, by DFG grant TRR54 and BMBF CancerSys project ProSiTu to CS and by Berliner Krebshilfe to KAPO and CS.
The authors declare that they have no competing interests.
RNAi system establishment, biochemical analysis of NF-κB subunits and kinases, and IP experiments were performed by KAPO. Immunohistochemistry was performed by MHu. EK performed ChIP experiments and prepared ChIP-seq libraries. ChIP-seq analysis was performed by GP and MHe together with EK. Genomic location analysis of ChIP seq regions was done by EMM and analyzed together with MAN. Clover and MEME analysis was performed by JFF and analyzed together with MAN. Validation of ChIP/microarray data by RT-qPCR and primer design was done by EK. Microarray experiments were performed by KAPO and analyzed by MHe and KAPO. Integration of ChIP-seq and microarray data and development of logistic regression models were performed by MHe. Cell death analyses were done by KAPO. Integration of datasets with gene expression data of primary samples and GO analysis were performed by MHe, together with KAPO. NH, SM, and CS analyzed data. KAPO, EK, and CS wrote the manuscript with contributions and feedback of all authors. All authors read and approved the final manuscript.
Supplemental Experimental Procedures shows a detailed version of materials and methods. (DOCX 60 kb)
shows co-immunoprecipitation of p50 and p52 with NF-κB subunits. Figure S2 shows the number of ChIP-seq regions bound by the individual NF-κB subunits, the number of genes assigned to ChIP-seq regions, the validation of ChIP-seq regions by qPCR, the frequency of ChIP-seq regions for NF-κB subunits relative to the transcription start site (TSS), p52 recruitment to selected ChIP-seq regions upon treatment of L1236 cells with NFKB2-siRNAs and control-siRNAs, and the genomic locations of p50 and p52-ChIP-seq regions with respect to annotated genes. Figure S3 shows p50- and p52-selective motifs detected by using MEME. Figure S4 shows the optimization and specificity controls for gene knockdown experiments. Figure S5 shows the prediction of gene expression patterns from combinatorial binding of NF-κB subunits. Figure S6 shows cell death induction by depletion of p50/RelA, p52/RelB or Bcl-XL. Figure S7 shows RT-qPCR analyses of genes that were differentially expressed in classical HL compared to other lymphoid malignancies and normal B cells. (PDF 3236 kb)
shows immunohistochemical analyses of HL biopsy samples with p100/p52 and p105/p50 antibodies. (DOCX 12kb)
shows ChIP-seq regions for RelA, RelB, p50, and p52 and assigned genes. (XLSX 5263 kb)
shows motifs over-represented in NF-κB ChIP-seq regions found by CLOVER analysis. (XLSX 11 kb)
shows the overlap of NF-κB subunit binding and DNAse I hypersensitive sites (DHS) in L1236 and GM12878 cells. (DOCX 13 kb)
shows the overlap of DHS and ChIP-seq regions for RelA, RelB, p50, and p52 in L1236 and GM12878 cells. (XLSX 1227 kb)
shows the comparison of combinatorial subunit binding between L1236 and GM12878 cells. (DOCX 13 kb)
shows the gene level summary of gene expression results. (XLSX 3346 kb)
shows the functional enrichment analysis of p50/RelA and p52/RelB targets using ‘Gene Ontology (GO) for biological processes’ analyses. (XLSX 6646 kb)
shows p50/RelA and p52/RelB targets that are differentially expressed in primary samples of lymphoid malignancies and normal B lymphocytes. (XLSX 40 kb)
About this article
- Transcription factor
- ChIP sequencing
- gene expression
- consensus sequence
- B lymphocytes
- cell death