- Open Access
Overexpression of NELFE contributes to gastric cancer progression via Wnt/β-catenin signaling-mediated activation of CSNK2B expression
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 54 (2021)
Accumulating evidence has highlighted the importance of negative elongation factor complex member E (NELFE) in tumorigenesis. However, the relationship between NELFE and gastric cancer (GC) remains unclear. This study aimed to explore the expression pattern and specific function of NELFE in GC.
NELFE expression was evaluated by immunohistochemistry and qRT-PCR in GC tissues, respectively. Cell proliferation, migration and invasion were measured by CCK-8, colony formation, transwell assays, and nude mice model. Bioinformatics analysis was performed to search potential target genes of NELFE, and a Cignal Finder 10-Pathway Reporter Array was used to explore potential signaling pathways regulated by NELFE. Dual-luciferase reporter assays, qRT-PCR and western blotting were conducted to verify their regulatory relationship. The expression correlations among NELFE, β-catenin and CSNK2B were further explored by immunohistochemistry on consecutive resections.
NELFE was significantly overexpressed in GC tissues both in protein and mRNA level and negatively correlated with the prognosis of GC patients. Gain- and loss-of-function experiments showed that NELFE potentiated GC cell proliferation and metastasis in vitro and in vivo. CSNK2B was identified as a downstream effector of NELFE. Wnt/β-catenin signaling may mediate the regulation of CSNK2B by NELFE. In addition, NELFE, β-catenin and CSNK2B were all remarkably upregulated in tumor tissues compared with adjacent normal tissues, and their expression levels in GC were positively correlated with each other.
Our findings reveal a new NELFE-Wnt/β-catenin-CSNK2B axis to promote GC progression and provide new candidate targets against this disease.
Dynamic regulation of transcription elongation shortly after initiation by RNA polymerase II (RNAPII) plays key roles in the implementation of gene expression widespread in metazoans [1, 2]. This step is modulated by positive and negative transcription elongation factors which are known as P-TEF and N-TEF, respectively [3, 4]. Acting as N-TEFs, the proteins negative elongation factor (NELF) and 5,6-dichloro-1β-D-ribofuranosylbenzimidazole (DRB)-sensitivity-inducing factor (DSIF) cooperatively inhibit elongation by associating with the RNAPII complex, and this effect can be alleviated by P-TEFb-mediated phosphorylation of RNAPII, DSIF and NELF, thus continuing gene transcription [4, 5]. Dysregulation of the elongation step has been found to contribute to human diseases including cancer [6, 7].
NELF has been shown the essential roles in the regulation of promoter-proximal pausing of RNAPII [8, 9]. It includes four subunits, NELFA, B, C/D, and E, and only NELFE (also known as RDBP) has a RNA binding ability. The structure of NELFE comprises of an N-terminal leucine zipper motif, a central domain abundant in Arg-Asp dipeptide repeats (the RD motif) and a C-terminal RNA recognition motif (RNM) containing two highly conserved elements dubbed RNP1 and RNP2, respectively . It has been implicated that NELFE binds to a set of RNA sequences and contributes to the activities of NELF [10, 11]. More importantly, NELFE has been found to be involved in cancer progression in several studies. For instance, oncogenic activation of NELFE promotes hepatocellular carcinoma (HCC) cell proliferation and predicts the strong metastatic potential of HCC cells , whereas knockdown of NELFE results in enhanced cell colony growth of breast cancer cells when treated with estrogen-like reagents . While these findings implicate a critical role of NELFE in cancers, its specific functions and underlying molecular mechanisms are not fully understood.
Gastric cancer (GC) is one of the most lethal human malignancies worldwide . Although advances of modern medical science have improved the survival rates of patients with GC, there is still a dearth of effective biomarkers and effective treatments, which represents a major challenge in the management of GC. In the present study, we focused on NELFE and investigated its role in the development of GC and explore the underlying molecular mechanisms. We showed that NELFE was overexpressed in GC tissues and could predict poor prognosis. The upregulation of NELFE facilitated GC cell proliferation and metastasis in vitro and in vivo. CSNK2B was identified as a downstream effector of NELFE. Further mechanistic study revealed that Wnt/β-catenin signaling was involved in the regulation of CSNK2B by NELFE.
Materials and methods
Tissue samples and immunohistochemistry methods
To analyze the protein expression of NELFE in GC, 2 independent commercial tissue microarrays (TMA) containing 68 primary tumors, 53 adjacent normal tissues and 35 metastatic tumors from patients with GC were purchased from Shanghai Outdo Biotech Co., Ltd., Shanghai, China (#HStm-Ade076Met-01 and #HStm-Ade120lym-01). For further IHC analysis of NELFE, β-catenin and CSNK2B expression and their correlation in GC, consecutive sections of a TMA containing 75 paired GC tumor and adjacent normal tissues were purchased from Shanghai Outdo Biotech Co., Ltd., Shanghai, China (#HStmA150CS02). The clinical characteristics of the patients can be checked at the website: http://www.superchip.com.cn/biology/tissue.html. Standard immunohistochemical staining procedures were performed using a specific antibody against NELFE (#sc-32,912, Santa Cruz Biotechnology, USA) at a dilution of 1:200, and the immunostained sections were scanned using an Aperio ScanScope CS scanner (Aperio, Vista, USA) and visualized via Aperio ImageScope software (Aperio, Vista, USA). The results were analyzed by two pathologists blinded to the clinical information independently based on the percentage of stained cells and staining intensity. Briefly, the percentage of NELFE/β-catenin/CSNK2B-positive cells was classified into 5 groups: < 10% (0), 10–25% (1), 25–50% (2), 50–75 (3), and > 75% (4). The staining intensity was divided into 4 groups: no staining (0), light brown (1), brown (2), and dark brown (3). The overall IHC scores of NELFE/β-catenin/CSNK2B were calculated using the following formula: overall score = percentage score×intensity score. The samples with an overall score ≤6 were defined as weak staining, and > 6 were defined as strong staining. In addition, 28 paired of GC samples and adjacent non-cancer tissues were collected from patients who underwent surgical resection and signed an informed consent in Shanghai East Hospital. After resection, these GC tissue samples were immediately frozen in liquid nitrogen and stored at − 80 °C for RNA extraction. The use of human samples in this study was approved by the ethics committee of Shanghai East Hospital, Tongji University School of Medicine.
Cell culture and reagents
The human GC cell lines AGS, MGC-803, HGC-27, SGC-7901 and BGC-823 were obtained from the Shanghai Cell Bank of the Chinese Academy of Sciences, Shanghai, China. Cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM; Corning, Inc., Corning, NY, USA) supplemented with 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin (M&C Gene Technology Ltd., Beijing, China). The cells were cultured under standard conditions in a humidified incubator with 5%CO2. To activate Wnt signaling, the indicated cells were treated with 5 μM Wnt agonist 1  (#S8178, Selleck Chemicals, Shanghai, China) for 24 h.
Plasmid construction, RNA inference and lentivirus transduction
A full-length cDNA encoding NELFE gene was PCR-amplified cloned into pcDNA 3.1 mammalian expression vector (pNELFE), and empty pcDNA 3.1 vector (pVEC) was used as control group. For transient expression of NELFE, cells were seeded into a 6-well plate at 80% confluence and were transfected with 2.5 μg of pNELFE or pVEC using Lipofectamine 3000 (Invitrogen, Thermo Fisher Scientific, US) according to the manufacturer’s protocol. NELFE or CSNK2B-specific small inference RNAs (siRNAs) and negative control siRNA (siNC) were synthesized by GenePharma (Shanghai, China). The sequences are as follow: siNELFE-1: 5′-GGAAAGGGAAUACUCUCUAdTdT-3′; siNELFE-2: 5′-GAUUCCUUGUGCCUCAUAUdTdT-3′; siCSNK2B-1: 5′-CUCCGUGGCAAUGAAUUCUdTdT-3′; siCSNK2B-2: 5′-GUCAAGACGAUUCGCUGAUdTdT-3′; siNC: 5′-UUCUCCGAACGUGUCACGUdTdT-3′. For siRNA transfection, the indicated cells were seeded into a 6-well plate with a confluence of 20–40%, and the siRNAs were transfected into these cells using Lipofectamine 3000 (Invitrogen, Thermo Fisher Scientific, US) as per the manufacturer’s instructions. NELFE-knockdown lentivirus (shNELFE) was packaged (siNELFE-1 sequence was used) and supplied by Genepharma (Shanghai, China), the NELFE-overexpression lentivirus (LV-NELFE) was constructed and packaged by Tuzhu Biotech (Shanghai, China). For lentivirus transduction, the host cells at a confluence of 50% were transduced with the lentivirus in the presence of polybrene (8 μg/ml) and stably infected cells were selected with puromycin (2 μg/ml).
Microarray data collection and bioinformatic analysis
The microarray data involved in this study were obtained from GEO datasets of NCBI (GEO accession number: GSE19826 and GSE13911), and the expression levels of CSNK2B were analyzed by GEO2R tool online. Kaplan-Meier survival curves were generated using Kaplan–Meier plotter (http://kmplot.com/analysis/) , a comprehensive online platform which is capable to assess the effect of over 54,000 genes (mRNA, miRNA, protein) on survival in 21 cancer types. The co-expressing genes of NELFE in GC were identified using the Coexpedia database (http://www.coexpedia.org/) , and their correlation was further verified by Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn/) . CSNK2B-interacting proteins and their KEGG pathway enrichments were analyzed by STRING database (Version 11.0, https://string-db.org/) . The KEGG pathway map was generated using ggplot2 and Cario packages under R v3.5.2.
Dual luciferase reporter assay
TCF/LEF1-Luc reporter plasmid was obtained from Genomeditech (Shanghai, China). For CSNK2B promoter-reporter plasmid construction, 1200 bp upstream promoter region of CSNK2B gene was cloned into pGL3 luciferase reporter vector. Dual luciferase reporter assays were carried out with a dual luciferase kit (Promega, Wisconsin, USA) as follows: Cells were seeded into a 24-well plate with a confluence of 70–90% prior to transfection. The indicated cells were co-transfected with the TCF/LEF1-Luc reporter plasmids or CSNK2B promoter-reporter plasmids and Renilla luciferase constructs using Lipofectamine 3000 as described above. After incubation for 24 h, cells were lysed using the dual luciferase kit and the luciferase activities were measured on a Glomax-multi Luminometer (Promega, Wisconsin, USA) according to the manufacturer’s protocols. Relative luciferase activity was calculated according to the following formula: luciferase activity = firefly luciferase bioluminescence/ Renilla luciferase bioluminescence.
Pathway reporter array
Cignal™ Finder 10-Pathway Reporter Array (Qiagen, Hilden, Germany) was performed to concurrently assess whether NELFE influenced the activity of 10 different tumor-associated signaling pathways in accordance with the manufacturer’s specifications, and luciferase activities of different pathways were determined using a dual luciferase kit (Promega, Wisconsin, USA) as described above.
Cell proliferation assay
The viability of GC cell lines was assessed using the Cell Counting Kit-8 (CCK-8) assay kit (Dojindo, Kumamoto, Japan). Briefly, cells were seeded into a 96-well plate at a density of 3000 cells per well suspended in 100 μl of complete culture medium. After incubation overnight, cells were treated with CCK-8 solution (10 μl per well) at the scheduled time intervals (0, 24, 48, 72, 96 h) and incubated for 1 h. Subsequently, the optical density (OD) values were measured at a wavelength of 450 nm on an automated microplate reader (SpectraMax M5, Molecular Devices, USA). For colony formation assay, cells were seeded into a 6-well plate with a density of 1000 cells/well in triplicate and cultured under normal conditions for 2–3 weeks. Colonies were fixed with 4% paraformaldehyde and stained with 0.5% crystal violet for 20 min, after which the plates were rinsed, dried, and colonies were photographed and counted using Image J software (v1.52, National Institutes of Health Freeware, USA).
Cell migration assay
Migration assays were performed using a 24-well transwell chamber system (8 μm pore diameter, Costar #3422, Corning, USA) according to the protocol of manufacturers. The indicated cells suspended with 300 μl FBS-free medium were seeded into the upper chamber with a density of 30,000/well, while 600 μl medium containing 10% FBS was added into the lower chamber. After incubating for 24 h at 37 °C, the cells left in the upper surface were removed with a cotton swab, and the migrated cells on the lower surface were fixed with 4% paraformaldehyde and stained with 0.5% crystal violet for 30 min. Cells that migrated through the membranes were photographed and counted under an inverted microscope.
Cell invasion assay
Cell invasion was measured with BD BioCoat Matrigel invasion chambers (BD Biosciences, #354480, USA). Cells suspended in 400 μl of serum-free DMEM were seeded into the upper chamber at a density of 5× 104 cells/well, while 800 μl of DMEM containing 10% FBS was added to the lower chamber as chemoattractant. After incubation for 24 h at 37 °C, cells on the upper side were wiped off with a cotton swab, and invaded cells on the lower side of the membrane were fixed with 4% paraformaldehyde and stained with 0.5% crystal violet for 30 min. The cells were photographed and counted under an inverted microscope at × 100 magnification.
Expression and subcellular location of NELFE and CSNK2B were detected by immunofluorescence assay. The indicated cells were planted onto glass coverslips and cultured under normal conditions. After fixation with 4% paraformaldehyde and permeabilized with 0.1% NP40 in PBS, the coverslips were blocked with 50% horse serum in PBS for 1 h at 4 °C and subjected to incubation of anti-NELFE (1:50) and anti-CSNK2B (1:50) at 4 °C overnight. After incubated with fluorescein-conjugated secondary antibody (Invitrogen, USA) for 1 h, the coverslips were mounted with 4,6-Diamidino-2-Phenylindole (DAPI)-containing mounting medium (Sigma, USA). Immunofluorescence staining was visualized under a Leica SP8 confocal laser scanning microscope (Leica Biosystems, Germany).
Western blot analysis
Total proteins were extracted from GC cell lines using RIPA lysis buffer (#20–188, Merck Millipore, MA, USA) with 1× protease and phosphatase inhibitor cocktails (Sigma-Aldrich, MO, USA). Total protein lysates were obtained by centrifuged at 4 °C and 12,000 rpm for 15 min, and protein concentration was determined by using a BCA Protein Assay Kit (#T9300A, Takara Bio, Kyoto, Japan). 4× SDS loading buffer (#KGP101X, Keygentec Inc., Shanghai, China) was added into the protein lysates and the samples were boiled for 5 min prior to electrophoresis. Aliquots of the proteins were separated by 10% SDS-PAGE and transferred onto a nitrocellulose membrane (0.45 μm, Bio-Rad, CA, USA), followed by block of non-specific binding with 5% non-fat milk in PBS containing 0.05% Tween-20 (PBST) for 1 h at room temperature. Then membranes were subjected to incubation with specific primary antibodies at 4 °C overnight and secondary antibodies at room temperature for 1 h (avoid light). Immunoblots were visualized by the Odyssey Digital Infrared Imaging System (LiCoR Biosciences, NE, USA), and β-actin was used for internal reference. Primary antibodies used in the present study are as follow: anti-NELFE (1:200, # sc-32,912, Santa Cruz Biotechnology, CA, USA), anti-CSNK2B (1:1000, # 20234–1-AP, Proteintech Group, Wuhan, China), anti-β-catenin (1:1500, #8480, Cell Signaling Technology, MA, USA) and anti-β-actin (1:5000, #60008–1-Ig, Proteintech Group, Wuhan, China). Goat anti-rabbit-DyLight 800 (1:1000, #SA5–35571, Thermo Fisher Scientific, MA, USA) and anti-mouse-DyLight 800 (1:1000, #SA5–35521, Thermo Fisher Scientific, MA, USA) served as secondary antibodies.
Quantitative Real Time PCR (qRT-PCR) analysis
Total RNA was isolated using TRI Reagent (Sigma-Aldrich, MO, USA) and then reverse transcribed into cDNA by PrimeScript RT reagent Kit with gDNA Eraser (Takara Bio, Kyoto, Japan) according to the manufacturers’ protocols. Synthesized cDNA was subjected to qPCR analysis using TB Green Premix Ex Taq II (Takara Bio, Kyoto, Japan) on an ABI QuantStudio 6 Flex Real-time PCR system (Applied Biosystems, CA, USA) in accordance with the manufacturer’s instructions. β-actin was used as the internal reference gene and relative gene expression level was calculated using the comparative cycle threshold (ddCt) method. Primer sequences for qPCR analysis were synthesized HuaGene Biotech (Shanghai, China) and listed in Supplementary Table 1.
Four to six weeks male athymic BALB/c nude mice were obtained from SLAC Laboratory Animal Company (Shanghai, China) and were raised under pathogen-free conditions with access to distilled food and water. To establish a subcutaneous xenograft model, 2× 106 of MGC-803/BGC-823 cells suspended in 100 μl of PBS were inoculated subcutaneously into left (LV-VEC/shNC) or right (LV-NELFE/shNELFE) flank axillary region of nude mice (5 mice per group). The size of tumors was measured every 7 days, and tumor volume was calculated using the following formula: Tumor volume (mm3) = length × width2/2. Tumor bearing mice were euthanatized on day 28 after cell injection, histologically intact tumors were removed, weighted and photographed. For the lung metastasis model, 3× 106 the above-mentioned GC cells in 100 μl of PBS were injected intravenously into the nude mice (6 mice per group) through tail vein. On day 42, the nude mice were sacrificed, lung tissues were collected and the number of metastatic lung nodules was counted. The lung tissues were immersed in 4% paraformaldehyde for subsequent histological examination, and potential lung metastatic lesions were finally confirmed by hematoxylin and eosin (H&E) staining method. All animal experiments were approved by the Medical Ethics Committee of Shanghai East Hospital.
All in vitro and in vivo experiments were performed independently at least three times in triplicate, and continuous variables were presented as mean ± standard deviation (SD). Statistical analysis was performed by using GraphPad Prism software 7.0. The student’s t-test or one-way analysis of variance (ANOVA) was used to compare the mean difference between binary or multiple variables, the χ2 test was used to explore the relationships between NELFE/β-catenin/CSNK2B expression and clinicopathological parameters of GC patients. Pearson correlation coefficient analysis was used for comparing the correlation among the expression of three genes. A P value of < 0.05 was considered statistically significant.
NELFE is significantly upregulated in GC tissues
To investigate the expression pattern of NELFE in GC, we performed immunohistochemistry (IHC) analysis using commercial tissue microarrays (TMAs) containing 68 primary gastric tumor tissues, 35 metastatic tumor tissues and 53 adjacent normal tissues derived from 102 GC patients, and representative IHC images of NELFE were shown in Fig. 1a. Compared with adjacent normal tissues, elevated NELFE expression was observed in both primary and metastatic GC tissues (Fig. 1b). Intriguingly, metastatic tumor tissues exhibited significantly higher NELFE expression than primary tumor tissues (P < 0.05, Fig. 1b), suggesting that NELFE might be a potential biomarker to distinguish metastatic gastric malignancies from primary gastric malignancies. Besides, additional tumor tissues and matched adjacent normal tissues from 28 GC patients were collected and subjected to quantitative real-time PCR (qRT-PCR) analysis. The result showed that the mRNA level of NELFE was also remarkably elevated in GC tissues (P < 0.0001, Fig. 1d), among which 18 cases (64.29%) exhibited an increase of more than 1.5-fold change (Fig. 1c). In order to confirm the reliability of above data, we analyzed the expression profiles of NELFE in human gastric cancer samples from public Gene Expression Omnibus (GEO) datasets (GEO accession number: GSE19826 and GSE13911). Consistently, NELFE expression in GC tumor samples was higher than that in normal tissues (Fig. 1d). To evaluate the prognostic significance of elevated NELFE expression in GC, survival analysis was conducted on an online Kaplan-Meier Plotter database using published microarray data from GC patients. As shown in Fig. 1e, high NELFE expression was negatively correlated with overall survival (P < 0.05) and first progression survival (P < 0.05) rates of GC patients. By immunofluorescence assay, we observed that the subcellular location of NELFE in MGC-803 cells was mostly in nucleus and partly in cytoplasm (Fig. 1f). In addition, the expression pattern of NELFE in various GC cell lines was shown in Fig. 1g. Taken together, these data suggested that NELFE is overexpressed in GC and may serve as a potential biomarker for the progression of GC.
NELFE promotes GC cell proliferation in vitro and in vivo
We next assessed the biological functions of NELFE in GC. To achieve NELFE overexpression or knockdown in GC cell lines, NELFE-overexpression plasmid (pNELFE) was transfected into MGC-803 and HGC-27 cells with a relatively low NELFE expression, and specific siRNAs against NELFE (siNELFE-1 and siNELFE-2) were constructed and transfected into BGC-823 and AGS cells with high levels of NELFE expression, after which western blot analyses confirmed the transfection efficiencies (Fig. 2a). Subsequently, cell viability was measured using CCK-8 assays. The growth curves indicated that NELFE overexpression facilitated the proliferation of MGC-803 and HGC-27 cells (Fig. 2b), whereas BGC-823 and AGS cells with NELFE silencing exhibited reduced proliferation rates as compared to their control groups (Fig. 2c). The effects of NELFE on long-term proliferation abilities of GC cell lines were investigated using colony formation assays after the indicated cells were transduced with lentivirus for stable overexpression or silencing of NELFE. As expected, exogenous NELFE augmented colony formation of MGC-803 and HGC-27 cells (Fig. 2d), while reduced NELFE expression dramatically inhibited the colony growth of BGC-823 and AGS cells (Fig. 2e). To verify whether these in vitro findings were relevant to GC tumor growth in vivo, a xenograft model of nude mice was generated by subcutaneous injection of GC cells with stable expression or knockdown of NELFE. Consistent with the in vitro observations, forced expression of NELFE substantially increased tumor size and tumor weight (Fig. 2f), while knockdown of NELFE led to opposite results (Fig. 2g). These findings revealed that NELFE plays important roles in GC cell growth in vitro and in vivo.
Abnormal upregulation of NELFE is essential for GC cell migration, invasion and metastasis
Considering the IHC results that NELFE expression was significantly upregulated in metastatic tumor tissues in comparison with primary tumor tissues, we assumed that NELFE plays a role in tumor metastasis of GC. To verify this hypothesis, transwell migration and invasion assays were carried out in GC cells following transfection with pNELFE or siNELFE-1/2. The results showed that NELFE overexpression drastically increased the number of migratory cells (Fig. 3a), while siRNAs-mediated silencing of NELFE resulted in opposite results (Fig. 3b). As the same time, similar phenotypes were also observed in cell invasion assays (Supplementary Fig. 1a and b). Moreover, a lung metastasis model of athymic nude mice was generated using GC cells with stable overexpression or knockdown of NELFE. As expected, NELFE overexpression increased pulmonary metastatic nodules of MGC-803 cells (Fig. 3c), whereas NELFE knockdown in BGC-823 cells reduced pulmonary metastatic nodules compared with the control group (Fig. 3d). Collectively, the above findings suggested that NELFE acts as a facilitating factor for metastasis of GC.
Oncogenic functions of NELFE in GC is mediated by CSNK2B
To explore potential mechanisms possible for the cancer-promoting functions of NELFE, we performed gene co-expression analysis using data from Coexpedia database, and 12 co-expressed genes of NELFE were identified in GC (Fig. 4a). For further validation, Pearson correlation analysis was conducted to compare the relationships between NEFLE and these genes using an online Cancer Genome Atlas (TCGA) data visualization web tool GEPIA. Among the 12 genes, casein kinase 2 beta (CSNK2B, r = 0.77, P < 0.01) and BAG cochaperone 6 (BAG6, r = 0.70, P < 0.01) showed the highest correlation with NELFE in GC (Fig. 4b). Subsequently, we performed qRT-PCR analysis to determine their mRNA expression levels after knockdown of NELFE in AGS cells. Notably, a concomitant decrease in CSNK2B mRNA was observed, while BAG6 expression showed no significant change in AGS-shNELFE cells as compared to the control cells (Fig. 4c), indicating that CSNK2B was regulated by NELFE. Further western blot results confirmed the positive effects of NELFE on CSNK2B protein expression (Fig. 4d). CSNK2B is a regulatory subunit of casin kinase 2 (CK2) and plays a key role in regulating its kinase activity . Previous studies have shown that CSNK2B is an oncogene in different human malignancies by directly trans-activating NF-κB signaling pathway [21, 22].
On the basis of these findings, we next synthesized two siRNAs against CSNK2B (siCSNK2B-1 and siCSNK2B-2) to test whether NELFE exerts its function dependent on CSNK2B expression. Western blot results showed that both siCSNK2B-1 and siCSNK2B-2 successfully silenced CSNK2B expression in MGC-803 cells, while siCSNK2B-2 exhibited a better knockdown efficiency than siCSNK2B-1 (Fig. 4e), thus siCSNK2B-2 was selected for further cell biology experiments. CCK-8 assays were carried out after transfection of siCSNK2B-2 into MGC-803 cells with NELFE overexpression. In line with previous reports [22,23,24], knockdown of CSNK2B significantly arrested cell proliferation, and more importantly, it also strongly weakened the positive effects of NELFE on proliferation ability of MGC-803 cells (Fig. 4f). Similar results were also observed in colony formation and transwell migration assays (Fig. 4g, h), hinting us that NELFE enhanced GC cell proliferation and migration through activating CSNK2B expression. Accordingly, it can be deduced that CSNK2B is a downstream factor of NELFE and mediates positive effects of NELFE on GC progression.
NELFE enhances CSNK2B gene expression via activation of Wnt/β-catenin pathway
In an attempt to better understand how NELFE regulates CSNK2B expression in GC, we used the Cignal Finder 10-pathway reporter array to screen NELFE-associated signaling pathways in GC cells. Among the 10 signaling pathways directly related with tumor progression, Wnt/β-catenin pathway was strongly activated by NELFE overexpression (Fig. 5a). Next, a TCF/LEF1 luciferase reporter assay was performed to validate the effect of NELFE on the activity of canonical Wnt/β-catenin signaling pathway. Consistent with the data described above, higher levels of Wnt-dependent activities were observed in cells with NELFE overexpression, whereas silenced NELFE led to contrary results (Fig. 5b). Furthermore, upon NELFE overexpression or knockdown in GC cell lines, we examined the protein and mRNA level of β-catenin, a major component of Wnt signaling pathway , respectively. Interestingly, the protein level of β-catenin was enhanced by NELFE (Fig. 5c), while no significant change in mRNA level was observed after NELFE overexpression or knockdown (Fig. 5d), implying the regulation of β-catenin by NELFE is not at the transcriptional level. To explore whether CSNK2B is associated with Wnt/β-catenin signaling pathway, CSNK2B-interacting proteins were identified using the STRING database and further KEGG pathway analysis was performed. Notably, CSNK2B-interacting proteins were highly enriched in Wnt/β-catenin signaling pathway (Fig. 5e), suggesting a relationship between CSNK2B and Wnt signaling. Thereby, we aimed to explore whether Wnt/β-catenin signaling pathway could regulate CSNK2B expression in GC cells. By comparing sequences, 2 potential TCF/LEF response elements with the consensus sequence CANNTG (E-box motifs) were recognized in the promoter region of CSNK2B (Fig. 5f), indicating that CSNK2B might be a target gene of TCF/LEF transcription factors. To verify this possibility, CSNK2B promoter-luciferase reporter plasmid was constructed and a Wnt/β-catenin pathway activator (Wnt agonist 1) was used for dual-luciferase reporter assays. In agreement with our expectations, treatment with Wnt agonist 1 dramatically increased the activity of CSNK2B promoter in MGC-803 and AGS cells (Fig. 5g), indicating that Wnt/β-catenin pathway mediated transcriptional activation of CSNK2B gene. Further immunofluorescence analysis showed a significant enhancement of CSNK2B staining intensity after Wnt agonist 1 treatment (Supplementary Fig. 1c). Combined with the above analyses, we questioned whether NELFE regulated CSNK2B through activation of Wnt/β-catenin pathway. Dual luciferase assays demonstrated that NELFE inhibition significantly reduced CSNK2B promoter activity (P < 0.01), which could be rescued by treatment of Wnt agonist 1 (Fig. 5h). Collectively, these data strongly suggested that Wnt/β-catenin pathway mediate the activation of CSNK2B expression induced by NELFE.
Increased expression of NELFE, β-catenin and CSNK2B in GC tissues are positively correlated with each other
To investigate possible clinical relevance of these 3 genes in GC, we performed an IHC analysis using consecutive sections on a TMA containing 75 paired tumor/adjacent normal tissues from GC patients. Consistent with the IHC data described above, there was a significant increase in NELFE expression in tumor tissues than in adjacent normal tissues (P < 0.01, Fig. 6a), and the same trend was observed in β-catenin expression (P < 0.001, Fig. 6b), which has been well understood in previous studies [26, 27]. Similarly, high CSNK2B expression was also found in tumor tissues when compared with adjacent normal tissues (P < 0.05, Fig. 6c). Further analysis of the same locations on consecutive tumor tissues sections indicated a consistent trend among the expression levels of all three genes in GC (Fig. 6d). Moreover, Pearson correlation analyses were applied based on the IHC scores of NELFE, β-catenin and CSNK2B (Fig. 6e). As expected, both β-catenin (R = 0.5072, P < 0.0001) and CSNK2B (R = 0.5292, P < 0.0001) expression showed a strong correlation with NELFE expression in GC, and both of them were positively associated with each other (R = 0.6143, P < 0.0001). Thus, these findings further supported the regulatory relationship among NELFE, β-catenin and CSNK2B in GC.
It is well established that tumorigenesis is a complex process commonly accompanied with aberrant activation of oncogenic genes . NELFE, a member of RNA-binding proteins (RBP) which is involved in regulating transcription elongation and RNA homeostasis, has received increasing attention in recent years in the development of human malignancies for its key roles in dysregulated transcriptional events [29,30,31]. Nevertheless up to date, the relationship between NELFE and GC progression remains largely unknown. In the present study, we for the first time investigated the expression features of NELFE in GC. We found NELFE was dramatically overexpressed in GC tumor tissues than in adjacent normal tissues by IHC and qPCR analysis, and metastatic tumors exhibited an increased NELFE expression than primary tumors, indicating the participation of NELFE in GC progression. Gene expression data from GEO database further confirmed these results, and a negative correlation of NELFE expression with the overall and first progression survival rates of GC patients were found based on the online survival database Kaplan-Meier plotter. Through gain- or loss-of-function experiments, we showed that NELFE facilitated GC cell proliferation, migration and invasion in vitro, and further animal experiments confirmed this conclusion in vivo. These observations unambiguously supported the cancer-promoting role of NELFE in GC, consistent with previous studies on NELFE in HCC or pancreatic cancer [12, 30].
In mechanistic studies, we performed comprehensive analysis combining the bioinformatics data and our experimental evidence to elucidate the exact mechanisms by which NELFE exerts its effects during GC tumorigenesis, and CSNK2B, a regulatory subunit of casein kinase II (CK2), was identified as a downstream effector of NELFE. Studies have demonstrated that CSNK2B is essential for the kinase activity of CK2, and overexpression of CSNK2B has been reported to promote HCC tumorigenesis via trans-activation of NF-κB [20, 21]. Although there is currently no available evidence regarding the roles of CSNK2B in GC, our cell biology experiments confirmed that CSNK2B facilitated GC cell proliferation and migration, and it also mediated the pro-oncogenic functions of NELFE in GC.
Another important finding of our work is that Wnt/β-catenin signaling pathway was uncovered to mediate the regulation of NELFE on CSNK2B expression. As we known, activation of β-catenin is then able to translocate into the nucleus and form a complex with T-cell factor/lymphoid enhancer factor (TCF/LEF) proteins to induce the transcription of genes involved in cell growth and proliferation, such as c-myc and cyclin D1 . Consistently, dual-luciferase assays showed that CSNK2B gene was transcriptionally activated by TCF/LEF, which means Wnt/β-catenin signaling pathway is a positive regulator of CSNK2B expression. Therefore, these findings identified a NELFE/β-catenin/CSNK2B axis in promoting the initiation and development of GC. However, the exact mechanism for how NELFE enhanced the protein expression of β-catenin, still needs more in-depth studies. In addition, it is worth noting that a recent study reported that NELFE is able to decrease the stabilization of N-Myc downstream-regulated gene 2 (NDRG2), which is involved in the inhibition of the Wnt/β-catenin signaling pathway , thus providing a possible explanation for the effects of NELFE on β-catenin expression. Certainly, more experimental evidence is needed to verify this hypothesis. Finally, our IHC analysis on consecutive sections of a TMA confirmed that NELFE, β-catenin and CSNK2B expression were all highly expressed in tumor tissues as compared to adjacent normal tissues, and all these three genes expression in GC were positively correlated with each other, which provides further support to the above experimental findings.
This study revealed that NELFE functions as an oncogene in GC and a new oncogenic NELFE/β-catenin/CSNK2B signaling axis was identified, which may help to develop new diagnosis and treatment strategies against GC (Fig. 7).
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Negative elongation factor complex member E
Quantitative real-time PCR
RNA polymerase II
Negative elongation factor
RNA recognition motif
Cell Counting Kit-8
PBS containing 0.05% Tween-20
Hematoxylin and eosin
One-way analysis of variance
Gene Expression Omnibus
BAG cochaperone 6
Casein kinase 2 beta
Cancer Genome Atlas
Casein kinase II
T-cell factor/lymphoid enhancer factor
N-Myc downstream-regulated gene 2
Chen FX, Smith ER, Shilatifard A. Born to run: control of transcription elongation by RNA polymerase II. Nat Rev Mol Cell Biol. 2018;19(7):464–78.
Uptain SM, Kane CM, Chamberlin MJ. Basic mechanisms of transcript elongation and its regulation. Annu Rev Biochem. 1997;66:117–72.
Peterlin BM, Price DH. Controlling the elongation phase of transcription with P-TEFb. Mol Cell. 2006;23(3):297–305.
Yamaguchi Y, Takagi T, Wada T, Yano K, Furuya A, Sugimoto S, et al. NELF, a multisubunit complex containing RD, cooperates with DSIF to repress RNA polymerase II elongation. Cell. 1999;97(1):41–51.
Yamaguchi Y, Wada T, Watanabe D, Takagi T, Hasegawa J, Handa H. Structure and function of the human transcription elongation factor DSIF. J Biol Chem. 1999;274(12):8085–92.
Lin C, Smith ER, Takahashi H, Lai KC, Martin-Brown S, Florens L, et al. AFF4, a component of the ELL/P-TEFb elongation complex and a shared subunit of MLL chimeras, can link transcription elongation to leukemia. Mol Cell. 2010;37(3):429–37.
Liang K, Smith ER, Aoi Y, Stoltz KL, Katagi H, Woodfin AR, et al. Targeting Processive transcription elongation via SEC disruption for MYC-induced Cancer therapy. Cell. 2018;175(3):766–79 e17.
Vos SM, Farnung L, Urlaub H, Cramer P. Structure of paused transcription complex pol II-DSIF-NELF. Nature. 2018;560(7720):601–6.
Narita T, Yamaguchi Y, Yano K, Sugimoto S, Chanarat S, Wada T, et al. Human transcription elongation factor NELF: identification of novel subunits and reconstitution of the functionally active complex. Mol Cell Biol. 2003;23(6):1863–73.
Yamaguchi Y, Inukai N, Narita T, Wada T, Handa H. Evidence that negative elongation factor represses transcription elongation through binding to a DRB sensitivity-inducing factor/RNA polymerase II complex and RNA. Mol Cell Biol. 2002;22(9):2918–27.
Rao JN, Neumann L, Wenzel S, Schweimer K, Rosch P, Wohrl BM. Structural studies on the RNA-recognition motif of NELF E, a cellular negative transcription elongation factor involved in the regulation of HIV transcription. Biochem J. 2006;400(3):449–56.
Iida M, Iizuka N, Tsunedomi R, Tsutsui M, Yoshida S, Maeda Y, et al. Overexpression of the RD RNA binding protein in hepatitis C virus-related hepatocellular carcinoma. Oncol Rep. 2012;28(2):728–34.
Aiyar SE, Sun JL, Blair AL, Moskaluk CA, Lu YZ, Ye QN, et al. Attenuation of estrogen receptor alpha-mediated transcription through estrogen-stimulated recruitment of a negative elongation factor. Genes Dev. 2004;18(17):2134–46.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Liu J, Wu X, Mitchell B, Kintner C, Ding S, Schultz PG. A small-molecule agonist of the Wnt signaling pathway. Angew Chem. 2005;44(13):1987–90.
Szasz AM, Lanczky A, Nagy A, Forster S, Hark K, Green JE, et al. Cross-validation of survival associated biomarkers in gastric cancer using transcriptomic data of 1,065 patients. Oncotarget. 2016;7(31):49322–33.
Yang S, Kim CY, Hwang S, Kim E, Kim H, Shim H, et al. COEXPEDIA: exploring biomedical hypotheses via co-expressions associated with medical subject headings (MeSH). Nucleic Acids Res. 2017;45(D1):D389–D96.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–W102.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–D13.
Litchfield DW. Protein kinase CK2: structure, regulation and role in cellular decisions of life and death. Biochem J. 2003;369(Pt 1):1–15.
Yu M, Yeh J, Van Waes C. Protein kinase casein kinase 2 mediates inhibitor-kappaB kinase and aberrant nuclear factor-kappaB activation by serum factor(s) in head and neck squamous carcinoma cells. Cancer Res. 2006;66(13):6722–31.
Xiao Y, Huang S, Qiu F, Ding X, Sun Y, Wei C, et al. Tumor necrosis factor alpha-induced protein 1 as a novel tumor suppressor through selective downregulation of CSNK2B blocks nuclear factor-kappaB activation in hepatocellular carcinoma. EBioMedicine. 2020;51:102603.
Chua MMJ, Lee M, Dominguez I. Cancer-type dependent expression of CK2 transcripts. PLoS One. 2017;12(12):e0188854.
Wang W, Wang X, Li C, Chen T, Zhang N, Liang Y, et al. Huaier suppresses breast Cancer progression via linc00339/miR-4656/CSNK2B signaling pathway. Front Oncol. 2019;9:1195.
Nusse R. Wnt signaling in disease and in development. Cell Res. 2005;15(1):28–32.
Wu C, Zhuang Y, Jiang S, Liu S, Zhou J, Wu J, et al. Interaction between Wnt/beta-catenin pathway and microRNAs regulates epithelial-mesenchymal transition in gastric cancer (review). Int J Oncol. 2016;48(6):2236–46.
Song Y, Li ZX, Liu X, Wang R, Li LW, Zhang Q. The Wnt/beta-catenin and PI3K/Akt signaling pathways promote EMT in gastric cancer by epigenetic regulation via H3 lysine 27 acetylation. Tumour Biol. 2017;39(7):1010428317712617.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Dang H, Takai A, Forgues M, Pomyen Y, Mou H, Xue W, et al. Oncogenic activation of the RNA binding protein NELFE and MYC signaling in hepatocellular carcinoma. Cancer Cell. 2017;32(1):101–14 e8.
Han L, Zan Y, Huang C, Zhang S. NELFE promoted pancreatic cancer metastasis and the epithelialtomesenchymal transition by decreasing the stabilization of NDRG2 mRNA. Int J Oncol. 2019;55(6):1313–23.
Dang H, Pomyen Y, Martin SP, Dominguez DA, Yim SY, Lee JS, et al. NELFE-dependent MYC signature identifies a unique Cancer subtype in hepatocellular carcinoma. Sci Rep. 2019;9(1):3369.
Ashihara E, Takada T, Maekawa T. Targeting the canonical Wnt/beta-catenin pathway in hematological malignancies. Cancer Sci. 2015;106(6):665–71.
This study was supported by grants from the National Key Research and Development Program of China (Code: 2017YFC1308900), the National Natural Science Foundation of China (Code: 81772568 and 81772567), Key Disciplines Group Construction Project of Pudong Health Bureau of Shanghai (Code: PWZxq2017–13), and the Outstanding Clinical Discipline Project of Shanghai Pudong (Code: PWYgy2018–02).
Ethics approval and consent to participate
This study has been approved by the ethics committee of Shanghai East Hospital, Tongji University School of Medicine, China.
Consent for publication
The authors declare there is no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primer sequences for qRT-PCR analysis.
a-b Cell invasion assays were carried out at 24 h after transfection of NELFE-expressing plasmids (a) or siRNAs (b), and representative images (upper) and quantification analysis (bottom) were shown. Magnification: 100×. *P < 0.05, **P < 0.01. c Representative images of CSNK2B immunofluorescence staining (red) in MGC-803 cells after Wnt agonist 1 treatment (5 μM). Magnification: 40×.
About this article
Cite this article
Yu, S., Li, L., Cai, H. et al. Overexpression of NELFE contributes to gastric cancer progression via Wnt/β-catenin signaling-mediated activation of CSNK2B expression. J Exp Clin Cancer Res 40, 54 (2021). https://doi.org/10.1186/s13046-021-01848-3
- Gastric cancer