Induction of NIK Transcription Directly Correlates with Invasion
Consistent with our previous findings11, we observed that the invasiveness of the human-derived glioblastoma (GBM) cell lines BT25, BT114, and BT116 was significantly induced by treatment with TWEAK, whereas treatment with TNFα did not stimulate invasion (Figure 1A). Transcriptome analysis of BT25 cells treated with TWEAK or TNFα, which preferentially activate the canonical or noncanonical NF-κB pathways, respectively 13,14, revealed that the expression of NIK (MAP3K14) directly correlated with glioma cell invasion and was highly induced in response to TWEAK treatment but not TNFα treatment (Figure 1B, C). Moreover, we observed TWEAK-specific upregulation of tumor necrosis factor receptor (TNFR)-associated factor 1 (TRAF1), a signaling adapter that interacts with and stabilizes NIK for activation of the noncanonical NF-κB pathway 15. Elevated levels of integrin β3 (ITGB3), integrin subunit alpha 11 (ITGA11), matrix metalloproteinase 9 (MMP9), and Fms-related tyrosine kinase (FLT1) were also observed, consistent with increased invasive and migratory potential 16. In addition to the induction of ITGB3, we observed that TWEAK treatment also elevated the expression of ITGA11, which has previously been described in glioma invasion 17. Furthermore, only TWEAK-treated glioma cells exhibited increased expression of integrins, as well as MMP9, all of which are associated cancer markers 18–21. Ingenuity Pathway Analysis (IPA) of groups of genes belonging to specific diseases and functions revealed that TWEAK treatment elevated the overall expression of cancer pathways, including tumor formation, invasion, and metastasis (Figure 1D). Additionally, RNA-seq analysis demonstrated that relative to other MAP3 kinases and upstream regulatory kinases in the NF-κB signaling pathways, TWEAK treatment singularly induced NIK/MAP3K14 gene expression (Figure 1E). For example, while qPCR analysis showed high expression of MAP3K14 as well as MAP3K8 in TWEAK-treated cells, MAP3K8 was also induced by both TNF and TWEAK. Moreover, analysis of cells actively undergoing invasion in collagen matrices revealed that NIK/MAP3K14 was the most highly upregulated gene compared to other kinases (Figure 1F). We also observed TWEAK-induction of NIK mRNA and invasion in other GBM cell lines, as well as mouse embryonic fibroblasts (MEFs) (Supp. Figure 1A, B, C).
After observing elevated NIK transcription and invasion upon TWEAK treatment, we investigated potential paracrine effects among glioma cells during invasion. We found that BT116 cells underwent increased invasion when treated with conditioned media (Control CM) from highly invasive BT25 cells compared with unconditioned media. This increased invasion was dependent on NIK, as it was not observed when BT116 cells were cultured with conditioned media from NIK knockout BT25 cells (NIK KO CM), while conditioned media from NIK KO cells rescued with ectopic expression of murine NIK (NIK KO-mNIK CM) restored the increased invasion (Supp. Figure 1D-E). Consistent with the results from media treatments, direct coculture of BT25 and BT114 cells increased the cell invasion/migration of the BT114 cells compared to cells cultured alone (Supp. Figure 1F). These results demonstrate that transcriptional induction of NIK is associated with elevated glioma invasiveness that is propagated by NIK-dependent paracrine signaling.
NIK Expression is Upregulated in Highly Invasive Glioma Cells and Promotes Collective Invasion
To monitor the induction of NIK transcription during invasion in vivo, we generated BT116 glioma cells stably expressing red fluorescent protein (RFP) under the promoter of NIK (pNIK-RFP). Analysis of invading BT116 pNIK-RFP cells revealed a general induction of NIK expression in invading cells (red signal from monolayer pseudocolored white) (Figure 2A), with the farthest invading cells exhibiting the highest RFP intensity or highest expression of NIK (Figure 2B). Treatment of BT116 pNIK-RFP cells with TWEAK further increased NIK expression (RFP signal-red/yellow), with RFP-positive cells invading farther under either untreated (NT) or TWEAK-treated conditions (Figure 2C).
To evaluate collective invasion and leader-follower phenotypes in a 3-dimensional view, we performed invasion assays of glioma tumor spheres, which better mimic cell-cell and cell-matrix interactions. Imaging of collagen-embedded BT116 pNIK-RFP spheres revealed an increase in pNIK-RFP expression upon TWEAK treatment (Figure 2D), with the highest pNIK-RFP expression observed among the most invasive cells at the leading edge of the sphere (Supp. Fig 2A, Figure 2D, E). Consistent with single-cell invasion assays, the induction of NIK expression directly correlated with the migration and dispersion of cells from glioma spheres, which was enhanced upon TWEAK treatment (Figure 2F) compared to the invasion of similar-sized spheres from the untreated group (Figure 2G). Cells with a high pNIK-RFP signal were observed among cells undergoing collective invasion at the sphere peripheries (Figure 2H). These results demonstrate that cells with the highest TWEAK-induced NIK expression directly correlated with the most invasive glioma cells, facilitating collective invasion, consistent with increased invasion and metastasis gene signatures.
Inhibition of NIK Activity Reduces Glioma Cell Invasion
Next, we evaluated whether NIK catalytic activity was required to promote GBM invasion. Treatment of cells with mangiferin, a natural inhibitor of NIK 22–24, significantly attenuated TWEAK-stimulated invasion to levels comparable with untreated or TNF-treated conditions (Figure 3A, B). Mangiferin was verified to inhibit activation of the noncanonical NF-κB pathway by TWEAK treatment in glioma cells, as seen with a reduction in p100-p52 processing and RelB (Figure 3C). Mangiferin also inhibited TWEAK-induced NIK transcription (Supp. Fig. 2B). Although mangiferin inhibited NIK and the noncanonical NF-κB pathway, it did not affect cell proliferation (Supp. Figure 2C). These data suggest therapeutic potential of inhibiting NIK to attenuate glioma invasion.
E2F4 and E2F5 Regulate NIK Gene Expression
The mechanisms underlying the regulation of NIK gene expression are not fully understood. A recent study identified MAP3K14 (NIK) as a significantly upregulated gene in human osteosarcoma cells with E2F activation25, and analysis of the NIK promoter revealed the presence of E2F binding sites (Supp. Fig. 3D). Thus, we next investigated whether E2F transcription factors played a role in early NIK gene expression in response to TWEAK treatment in glioma cells. We observed that TWEAK treatment increased nuclear E2F4 and E2F5 protein levels in glioma cells, with E2F4 having greater nuclear translocation in BT25 cells and E2F5 in BT114 cells (Figure 4A). We also found that overexpression of E2F5, but not E2F1, increased NIK transcript levels (Supp. Figure 3A). E2F4 and E2F5 single and double knockout cell lines generated by CRISPR-Cas9 gene editing exhibited a significant reduction in NIK protein levels when treated with TWEAK and MG132 to inhibit proteosome-dependent degradation (Supp. Figure 3B, Figure 4B). Moreover, E2F4/5 double knockout cells (E2F DKO) exhibited reduced p52/RelB nuclear translocation, demonstrating impaired activation of NIK-driven noncanonical NF-κB signaling (Supp. Figure 3C).
Next, we investigated whether E2F proteins directly regulate NIK transcription. Chromatin immunoprecipitation (ChIP) analyses demonstrated that antibodies specific to E2F4 and E2F5, but not IgG, were bound to certain regions of the NIK promoter (Figure 4C). Additionally, qPCR analysis of E2F-DKO cells showed significantly attenuated induction of NIK mRNA expression even with TWEAK treatment (Figure 4D). Reduced NIK expression in E2F DKO cells proved to have functional consequences, as these cells were poorly invasive, even after TWEAK treatment (Figure 4E), which was restored with ectopic expression of NIK (NIK OE) in E2F DKO cells (Figure 4F). Overall, these data establish a novel role for E2F regulation of NIK transcription in a stimulated state and thus affect the cell’s overall ability to invade (Supp. Fig. 3F).