Long Noncoding RNAs and Messenger RNAs Expression Profiles Potentially Regulated by ZBTB7A in Nasopharyngeal Carcinoma

Liu, Fei; Wei, Jiazhang; Hao, Yanrong; Tang, Fengzhu; Jiao, Wei; Qu, Shenhong; He, Ning; Cai, Yonglin; Lan, Jiao; Yang, Yong; Wang, Yongli; Li, Min; Weng, Jingjin; Li, Bing; Lu, Jinlong; Han, Xing

doi:https://doi.org/10.1155/2019/7246491

BioMed Research International

On this page

Abstract Introduction Methods Results Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 7246491 | https://doi.org/10.1155/2019/7246491

Long Noncoding RNAs and Messenger RNAs Expression Profiles Potentially Regulated by ZBTB7A in Nasopharyngeal Carcinoma

Fei Liu,¹Jiazhang Wei,²Yanrong Hao,³Fengzhu Tang,²Wei Jiao,¹Shenhong Qu,²Ning He,²Yonglin Cai,⁴Jiao Lan,¹Yong Yang,²Yongli Wang,²and Min Li² et al.

Academic Editor: Paul W. Doetsch

Received15 Dec 2018

Revised15 Apr 2019

Accepted20 May 2019

Published11 Jun 2019

Abstract

Our previous studies showed that ZBTB7A played an important role in promoting nasopharyngeal carcinoma (NPC) progression. However, molecular mechanisms of different levels of ZBTB7A are still unclear. It is necessary to search molecular markers which are closely connected with ZBTB7A. We selected NPC sublines CNE2 with stably transfecting empty plasmid (negative control, NC) and short hair RNA (shRNA) plasmid targeting ZBTB7A as research objectives. Microarray was used to screen differentially expressed long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) via shRNA-CNE2 versus NC-CNE2. Quantitative PCR (qPCR) was used to validate lncRNAs and mRNAs from the sublines, chronic rhinitis, and NPC tissues. Bioinformatics was used to analyze regulatory pathways which were connected with ZBTB7A. The 1501 lncRNAs (long noncoding RNAs) and 1275 differentially expressed mRNAs were upregulated or downregulated over 2-fold. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the upregulated or downregulated carbohydrate and lipid metabolisms probably involved in carcinogenicity of shRNA-CNE2 (P-value cut-off was 0.05). In order to find the molecular mechanisms of ZBTB7A, we validated 12 differentially expressed lncRNAs and their nearby mRNAs by qPCR. Most of the differentially expressed mRNAs are closely connected with carbohydrate and lipid metabolisms in multiply cancers. Furthermore, part of them were validated in NPC and rhinitis tissues by qPCR. As a result, NR_047538, ENST00000442852, and fatty acid synthase (FASN) were closely associated with NPC. ZBTB7A had a positive association with NR_047538 and negative associations with ENST00000442852 and FASN. The results probably provide novel candidate biomarkers for NPC progression with different levels of ZBTB7A.

1. Introduction

Nasopharyngeal carcinoma (NPC) is an endemic malignant head and neck tumor in southern China [1, 2]. Epstein-Barr virus (EBV) and tumor metastasis related factors are closely connected with NPC progression [3, 4]. The biomarkers of EBV include Epstein-Bar encoded small nuclear RNA (EBER) in tissues, EBV DNA, and anti-EBV antibodies in circulating plasma or serum. All of them are important predictors in early diagnosis and prognosis of NPC. Furthermore, plasma EBV DNA is a better marker in advanced NPC than others [5–7]. However, a few patients have negative EBV DNA in plasma [8, 9]. The negative result indicates a limitation of EBV DNA. In order to improve detection rate of NPC, it is imperative to search novel candidate biomarkers.

ZBTB7A is also named as Pokemon, FBI-1, OCZF, and LRF. It is a critical transcription factor in the poxvirus and zinc finger/broad complex, tramtrack, and bric-a-brac (POZ/BOZ) and Krüppel (POK) family, which can specifically bind DNA through the Krüppel-like C₂H₂ zinc fingers and repress transcription by interacting transcriptional cofactors with POZ/BTB domain [10]. Takahiro Maeda [11] firstly discovered that ZBTB7A was a proto-oncogene because it peculiarly repressed the transcription of tumor suppressor alternative reading frame (ARF). Overexpression of ZBTB7A promotes lymphomas progression [11]. It also plays an oncogenic role in non-small cell lung cancer [12], ovarian cancer [13], breast cancer [14], hepatocellular carcinoma [15], and osteosarcoma [16]. However, loss of ZBTB7A shows it is a tumor suppressor in prostate cancer [17], melanoma [18]. The dual roles of ZBTB7A controversially appear in a kind of cancer such as colorectal cancer [19, 20] and gastric cancer [21, 22]. Basing on the complex roles in multiply cancers, ZBTB7A is confirmed to be a significant target of prognosis and therapy [11–22]. It is also an important proto-oncogene in NPC [23–25]. ZBTB7A level of NPC tissues was mostly higher than that of chronic rhinitis tissues [23]. Overexpression of ZBTB7A promoted cell vitality, migration, and invasion (Supplementary Materials, Figure S1 and S2). The tumorigenicity of NPC cell lines CNE2 and CNE3 with transiently decreasing ZBTB7A was weaker than that of cells transfecting empty plasmid [24, 25]. However, progression of NPC cell lines CNE2 and 5-8F was compensatorily promoted when ZBTB7A was stably knocked down [24]. The results indicate that ZBTB7A probably connects with complex pathways including oncogene and tumor suppressor genes.

Long noncoding RNAs (lncRNAs) are non-protein-coding transcripts with more than 200 nucleotides in length [26]. Since lncRNA microarray is continuously developed [27], lncRNAs have been demonstrated to involve in biological processes of cancer [28–30]. Some lncRNAs such as antisense noncoding RNA in the INK4 locus (ANRIL), HOXA transcript at the distal tip (HOTTIP), and differentiation antagonizing non-protein-coding RNA (DANCR) deteriorate NPC development [31–33]. Maternally expressed gene 3 (MEG3) and LINC0086 suppress proliferation of NPC cells [34, 35]. Although the connections between ZBTB7A and lncRNAs are reported in osteosarcoma and non-small cell lung cancer [16, 36, 37], their relationships are still mysterious in NPC.

Therefore, it is meaningful to discover the unknown field. We will screen differentially expressed lncRNAs and messenger RNAs (mRNAs) from NPC cell sublines with stably transfecting empty plasmid and short hair RNA (shRNA) plasmid targeting ZBTB7A by lncRNA microarray. Bioinformatics will be used to search potential connections between ZBTB7A and relevant pathways. Part of them will be selected and validated by qPCR in the sublines, chronic rhinitis, and NPC tissues.

2. Methods

2.1. Patient Selection

After being approved by patients, 20 chronic rhinitis and 60 NPC tissues were obtained using biopsy (July 2016-February 2017). All of them were preserved in liquid nitrogen. The tissues were confirmed by pathology from the formalin-fixed paraffin wax-embedded samples. The pathological type of NPC was undifferentiated carcinomas of the nasopharyngeal type (UCNT). This study was approved by the ethics committee of The People’s Hospital of Guangxi Zhuang Autonomous Region.

2.2. Cell Culture

CNE2 and CNE3 were preserved in Research Center of Medical Sciences, The People’s Hospital of Guangxi Zhuang Autonomous Region (Nanning, Guangxi, China). CNE2 was established from a patient suffering from NPC in Guangdong Province. The histological type of NPC was a poorly differentiated squamous carcinoma [38]. CNE3 was established from a patient suffering from liver metastatic carcinoma of primary NPC in Guangxi Province. The histological type of liver metastatic carcinoma and CNE3 was a poorly differentiated adenosquamous carcinoma [39], but that of CNE3 turned poorly differentiated adenocarcinoma after 20 years [40]. As a team member of Professor Yi Zeng, Wei Jiao has cultured and detected the cell lines in our hospital [23–25, 41, 42]. 5-8F and 6-10B were kindly provided by Professor Musheng Zeng (State Key Laboratory of Oncology in South China, Sun Yat-Sen University Cancer Center, Guangzhou, Guangdong, China). They were grown in RPMI Medium 1640 basic with 10% fetal bovine serum (Gibco; Thermo Fisher Scientific, Inc., Waltham, MA, USA). NP69 was kindly provided by Professor Sai-Wah Tsao (School of Biomedical Sciences, University of Hong Kong, Hong Kong SAR). NP69 was in Keratinocyte-SFM with 5% Bovine Pituitary Extract and Recombinant Epidermal Growth Factor (Gibco). All of them were in supplementary material. Based on the prior study [24], we used CNE2 sublines with stably transfecting empty plasmid and shRNA plasmid targeting ZBTB7A to screen differentially expressed lncRNAs and mRNAs. The sublines were named as negative control-CNE2 (NC-CNE2) and shRNA-CNE2, of which tumorigenicity was stronger than that of NC-5-8F and NC-6-10B [24].

2.3. LncRNA Microarray

shRNA-CNE2 and NC-CNE2 were used in lncRNA microarray. Both of them were acquired from three different passages. Total RNA was extracted by TRIZOL (Invitrogen, Carlsbad, CA, USA). RNA quantity and quality were measured by NanoDrop ND-1000 (Thermo Fisher Scientific, Inc.). RNA integrity was assessed by standard denaturing agarose gel electrophoresis. The Arraystar Human lncRNA Microarray V3.0 was designed for the global profiling of human lncRNAs and protein-coding mRNAs. 30,586 lncRNAs and 26,109 mRNAs could be detected by third-generation lncRNA microarray.

2.4. RNA Labeling and Array Hybridization

Sample labeling and array hybridization were performed according to the Agilent One-Color Microarray-Based Gene Expression Analysis protocol (Agilent Technology, Santa Clara, CA 95051, USA) with minor modifications. Briefly, mRNA was purified from total RNA after removal of ribosome RNA (rRNA; mRNA-ONLY™ Eukaryotic mRNA Isolation Kit; Epicentre, Madison, WI 53719, USA). Each sample was amplified and transcribed into fluorescent complementary RNA (cRNA) along the entire length of the mRNAs without 3′ bias utilizing a random priming method (Arraystar Flash RNA Labeling Kit; Arraystar, Inc., Rockville, MD 20850, USA). The labeled cRNAs were purified by RNeasy Mini Kit (Qiagen, Düsseldorf, Nordrhein-Westfalen, Germany). The concentration and specific activity of the labeled cRNAs (pmol Cy3/μg cRNA) were measured by NanoDrop ND-1000. 1 μg of each labeled cRNA was fragmented by adding 5 μl 10 x Blocking Agent and 1 μl of 25 x Fragmentation Buffer, then the mixture was heated at 60°C for 30 min, finally 25 μl 2 x GE hybridization buffer was added to dilute the labeled cRNA. 50 μl of hybridization solution was dispensed into the gasket slide and assembled to the lncRNA expression microarray slide. The slides were incubated for 17 h at 65°C in an Agilent Hybridization Oven. The hybridized arrays were washed, fixed, and scanned with using the Agilent DNA Microarray Scanner (part number G2505C; Agilent Technology).

2.5. Data Analysis

The Agilent Feature Extraction software (version 11.0.1.1) was used to analyze acquired array images. Quantile normalization and subsequent data processing were performed with using the GeneSpring GX v12.1 software package (Agilent Technologies). Differentially expressed lncRNAs and mRNAs between the two groups were identified using fold-change >2 as the cut-off. Hierarchical clustering and combined analysis were performed using homemade scripts. Gene ontology (GO) and pathway analysis were performed in the standard enrichment computation method.

2.6. Validation of the Differentially Expressed LncRNAs and MRNAs by QPCR

The total RNA of NC-CNE2 and shRNA-CNE2 cells, chronic rhinitis, and NPC tissues was extracted by TRIZOL (Invitrogen) and reversely transcribed by SuperScript ™ III Reverse Transcriptase Kit (Invitrogen). 12 lncRNAs and their nearby mRNAs were selected basing on data analysis. All of them were differentially expressed. The quantitative PCR (qPCR) was executed with 2 × PCR master mix (Arraystar) by ViiA 7 qPCR System (Applied Biosystems; Thermo Fisher Scientific, Inc.). Specific primers of the lncRNAs and nearby mRNAs were designed by Primer 5.0 (Table 1). The reaction program consisted of an initial denaturation step at 95°C for 10 minutes; denaturation at 95°C for 10 seconds, and annealing at 60°C for 60 seconds for 40 cycles; the dissociation stage at 95°C for 10 seconds, 60°C for 1 minute, 95°C for 15 seconds, and 60°C for 15 seconds. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as an internal control. The experiments were performed in triplicate. The fold-change of the lncRNAs and mRNAs were calculated by the 2 Delta-delta cycle threshold (Ct) method.

2.7. Statistical Analysis

Statistical analysis was performed by SPSS 13.0 and sigmaPlot 12.5 (SPSS Inc., Chicago, IL, USA), GraphPad Prism 5.0 (GraphPad Software, San Diego, CA, USA). The results of the assays were presented as mean ± Standard Deviation (SD). Statistical differences between two groups were evaluated with independent samples t-test. In GO analysis, Fisher’s exact test was used to search more overlap between the differentially expressed list and the GO annotation list than that of being expected by chance. In scatter plot of tissues by qPCR, Mann Whitney U test was used if the variances were significantly different through F test. P<0.05 was considered to be statistically significant.

3. Results

3.1. Differentially Expressed LncRNAs and MRNAs in NPC Cells with Low Expression of ZBTB7A

After quantile normalization of the raw data, the expression profiles of 23947 lncRNAs and 21058 mRNAs were obtained from shRNA-CNE2 and NC-CNE2 (Figure 1). Cluster analysis showed some lncRNAs and mRNAs were differentially expressed via shRNA-CNE2 versus NC-CNE2 (Figures 1(a) and 1(b)). Box plot showed the distribution of normalized intensity values of test and control groups was generally symmetrical (Figures 2(c) and 2(d)). Scatter plot figuratively showed relative folds of 1501 lncRNAs (738 upregulated and 763 downregulated) and 1276 mRNAs (679 upregulated and 597 downregulated), which were over 2-fold. The number of differentially expressed lncRNAs and mRNAs was listed. We divided them into 3 groups according to different fold-change. There were 2-5-fold-change, 5-10-fold-change, and more than 10-fold-change. 77 lncRNAs and 44 mRNAs were over 10-fold-change (Table 2).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1

RNA expression profiles in shRNA-CNE2 and NC-CNE2 cells. Hierarchical clustering indicated (a) lncRNA and (b) mRNA profiles. Box plot of (c) lncRNA or (d) mRNA profiles was a traditional method for visualizing separately the distribution of dataset. Scatter plot of (e) lncRNA or (f) mRNA profiles was a convenient visualization for assessing the variation of chips. T1-T3 meant test groups which were stably transfected with shRNA plasmid targeting ZBTB7A. C1-C3 meant negative control groups which were stably transfected with blank plasmid. ‘Red’ showed relatively high expression. ‘Green’ showed relatively low expression. The values of the x- and y-axes were the average normalized signal values of NC-CNE2 and shRNA-CNE2 in the scatter plot ( scaled).

(a)

(b)

Figure 2

GO (gene ontology) analysis. (a) GO upregulated. (b) GO downregulated. The P value denoted the significance of GO terms enrichment in the DE mRNAs. The lower the p-value, the more significant the GO Term (P value ≤0.05 was recommended). DE, differentially expressed; BP, biological process; CC, cellular component. MF, molecular function. It meant oxidoreductase activity, acting on the CH-CH group of donors, NAD or NADP as acceptor. It meant oxidoreductase activity, acting on the CH-CH group of donors. It meant oxidoreductase activity, acting on the aldehyde or oxo group of donors.

3.2. GO Analysis

GO is a system of standard classification of gene function, including biological process (BP), cellular component (CC), and molecular function (MF). The results showed the significant GO terms of differentially expressed genes. We found that the highest enriched GO terms of upregulated mRNAs were response to stress (GO: BP), cytoplasmic part (GO: CC), and oxidoreductase activity (GO: MF) (Figure 2(a)). The highest enriched GO terms of downregulated mRNAs were response to stress (GO: BP), cytoplasm (GO: CC), and catalytic activity (GO: BF) (Figure 2(b)). The results indicated some differentially expressed lncRNAs and relevant mRNAs probably involved in the activities.

3.3. KEGG Pathway Analysis

Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database resource for knowing high-level biological functions and utilities. KEGG pathways were constructed to better know the biological pathways and search the molecular mechanisms of NPC progression. The results showed significant pathways of differential genes. KEGG pathways connected with upregulated (Figures 3(a) and 3(b), left) mRNAs and downregulated (Figures 3(a) and 3(b), right) differentially distributed mRNAs. The top upregulated pathway was steroid hormone biosynthesis. Steroid is a kind of lipid. The top downregulated pathway was protein digestion and absorption. The third downregulated pathway was peroxisome proliferators-activated receptor (PPAR) signaling pathway, which included important lipid metabolism. 463 differentially expressed mRNAs were closely associated with carbohydrate and lipid metabolisms. The pathways of lipid metabolism probably were more meaningful than those of carbohydrate metabolism because of more enrichment scores in lipid pathways. In upregulated pathways, the scores of steroid hormone biosynthesis and steroid biosynthesis were 5.04301 and 3.129842. In downregulated pathways, the score of PPAR signaling pathway was 2.035416 (all p<0.05) (Table 3).

(a)

(b)

Figure 3

KEGG pathways in shRNA-CNE2 compared to NC-CNE2. (a) The most significantly enriched KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were illustrated by bubbles. The y-axis indicated the enrichment factors, which refers to the upregulated (left) and downregulated (right) differentially expressed mRNA number to the total mRNAs in a certain pathway. The size of bubble indicated mean number of mRNA enriched in a given pathway. The colour of bubble indicated P value (adjusted p value). (b) Bar plots represented enriched KEGG pathways associated with upregulated (left) and downregulated (right) DE mRNAs. The y-axis indicated average change fold of the implicated mRNAs. Sig, significant; DE, differentially expressed.

3.4. Signal Pathway Analysis

Steroid hormone biosynthesis was markedly activated in NPC cells with low expression of ZBTB7A. Aldo-keto reductase family 1 member C3 (AKR1C3) was markedly upregulated in the pathway. PPAR signaling pathway was markedly suppressed in the cells. Fatty acid binding protein 4 (FABP4) was markedly downregulated. The differentially expressed mRNAs of the pathways were associated with lipid metabolism, such as AKR1C3 and FABP4 (Figure 4).

(a)

(b)

Figure 4

The differentially expressed mRNAs in important ranked pathways. (a) Upregulated mRNAs in steroid hormone biosynthesis. 2.8.2.2, SULT2B1 (sulfotransferase family 2B member 1); 1.3.1.22, SRD5A3 (steroid 5 alpha-reductase 3); 1.1.1.213, AKR1C2 (aldo-keto reductase family 1 member C2); 1.1.1.50, AKR1C4; 1.1.1.149, AKR1C1; 1.1.1.51, AKR1C3; 1.14.14.1, CYP3A5 (cytochrome P450 family 3 subfamily A member 5); 1.14.13.-, CYP3A4 (cytochrome P450 family 3 subfamily A member 4); 2.4.1.17, UGT1A6, (UDP glucuronosyltransferase family 1 member A6). (b) Downregulated mRNAs in PPAR signaling pathway. PPAR, peroxisome proliferators-activated receptor; FABP, fatty acid binding protein; aP2, FABP4; HMGCS2, 3-hydroxy-3-methylglutaryl-CoA synthase 2; Apo-CIII, apolipoprotein C3; MCAD, ACADM (acyl-CoA dehydrogenase medium chain); GyK, GK (glycerol kinase). The regulated mRNAs in shRNA-CNE2 compared to NC-CNE2 were marked in coloured boxes. Red box indicated an upregulated mRNA; yellow box indicated a downregulated mRNA; green box indicated that no significant difference was observed.

3.5. QPCR Validation

lncRNAs mainly include antisense lncRNAs, enhancer lncRNAs, and long intergenic noncoding RNAs (lincRNAs). They have different biological functions through regulating nearby mRNAs [43–45]. In order to validate authenticity and reliability of lncRNA microarray, 12 differentially expressed lncRNAs and their nearby mRNAs were detected by qPCR. The differentially expressed mRNAs are closely connected with carbohydrate and/or lipid metabolisms of cancers, which means the lncRNAs probably interact with cancers via the mRNAs. The results of qPCR were consistent with those of lncRNA microarray through shRNA-CNE2 versus NC-CNE2. All of them have shown the same upregulation or downregulation trend (Figure 5). 6 lncRNAs (NR_033967, ENST00000398216, uc001enh.1, TCONS_00019671, TCONS_00025256, TCONS_00029013) and 1 mRNA (integrin subunit alpha 5, ITGA5) would not be used in next assay because the results of qPCR were not considered to be statistically significant (P>0.05). The other differentially expressed lncRNAs (NR_047538, ENST00000585189, ENST00000442852, TCONS_00020439, TCONS_00029159, TCONS_00013537) and mRNAs (cyclin-dependent kinase inhibitor 2B (CDKN2B), solute carrier family 2 member 1 (SLC2A1), aldolase C (ALDOC), interferon-stimulated gene 15 (ISG15), immediate early response 3 (IER3), thioredoxin-interacting protein (TXNIP), lysine acetyltransferase 5 (KAT5), FASN (fatty acid synthase), dual specificity tyrosine phosphorylation regulated kinase 1A (DYRK1A), salt inducible kinase 1 (SIK1), and paraoxonase 3 (PON3)) were continuingly validated by qPCR in chronic rhinitis and NPC tissues (Table 4). The results showed that NR_047538, ENST00000442852, fatty acid synthase (FASN), and ZBTB7A were closely connected with NPC (Figure 6). ZBTB7A had a positive association with NR_047538 and negative associations with ENST00000442852 and FASN (Figure 7).

(a)

(b)

(c)

Figure 5

Validation of lncRNA microarray by qPCR through shRNA-CNE2 versus NC-CNE2. (a) 12 differentially expressed lncRNAs and (b), (c) nearby mRNAs were validated by qPCR through shRNA-CNE2 versus NC-CNE2. The heights of the columns in the charts represent the mean expression values of fold-change. CDKN2B, cyclin-dependent kinase inhibitor 2B; SLC2A1, solute carrier family 2 member 1, also named glucose transporter 1 (GLUT1); ALDOC (aldolase C); ISG15 (interferon-stimulated gene 15); IER3 (immediate early response 3); TXNIP (thioredoxin-interacting protein); KAT5 (lysine acetyltransferase 5); ITGA5 (integrin subunit alpha 5); FASN (fatty acid synthase); DYRK1A (dual specificity tyrosine phosphorylation regulated kinase 1A); SIK1 (salt inducible kinase 1); PON3 (paraoxonase 3).

(a)

(b)

(c)

(a)

(b)

(c)

4. Discussion

lncRNA microarrays have been used in studies of multiple cancers [27, 46, 47]. Differentially expressed lncRNAs and mRNAs can be found by the method, which reveals some lncRNAs are important regulatory factors in melanoma, renal tumor, hepatocellular carcinoma, and colorectal cancer [48–51]. lncRNA ANRIL, also called NR_047538, promotes NPC progression via reprogramming cell glucose metabolism [31]. lncRNA HOTAIR promotes tumorigenesis of NPC via upregulating FASN, which is a pivotal enzyme of lipid metabolism [52]. Therefore, the technique is effectively applied in NPC study of exploring possible connections between ZBTB7A and differentially expressed lncRNAs.

Comparing to protein levels of ZBTB7A of NPC cell lines CNE2, CNE3, 5-8F, and 6-10B (Supplementary Materials, Figure S1(a)), we selected CNE2 and 5-8F with high expression of ZBTB7A as the objectives of stably knocking down ZBTB7A [24]. Then we selected CNE3 and 6-10B with low expression of ZBTB7A as the objectives of transfecting with ZBTB7A plasmid. The cells were transfected with plasmids and selected by G418. The best concentrations were 200 μg/ml (6-10B) and 230 μg/ml (CNE3). The concentration was, respectively, reduced to 100 μg/ml and 115 μg/ml after 10 days. The clone cells were stably cultured during 15 passages. The cells stably transfected with plasmids were named NC-6-10B, ZBTB7A-6-10B, NC-CNE3, and ZBTB7A-CNE3. The protein levels of ZBTB7A-6-10B and ZBTB7A-CNE3 were higher than those of control groups (Supplementary Materials, Figure S1 (b)). ZBTB7A-6-10B and ZBTB7A-CNE3 separately had stronger tumorigenicity than NC-6-10B and NC-CNE3 by 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide (MTT) and transwell assays (Supplementary Materials, Figure S2). The results of transiently and stably knocking down ZBTB7A in NPC cell line CNE2 were opposite [24]. We hypothesize that some important pathways vicariously maintain carcinogenicity of NPC cells. Therefore, it is a good method to search unknown pathways of connecting with ZBTB7A by lncRNA microarray.

CNE2 cells have specific characteristic of NPC because they are susceptible to EBV infection in vitro, while Hela cells are resistant [53]. They also have specific biomarkers and pathways of NPC [27, 54, 55]. Considering the differentially stable carcinogenicity, NC-CNE2 and shRNA-CNE2 were suitable models for screening differentially expressed lncRNAs and mRNAs NPC by lncRNA microarray.

As a result, some differentially expressed mRNAs indicated that ZBTB7A may connect with oncogenic and oncosuppressive pathways (Supplementary Materials, Figure S3 and Table S1). Loss of ZBTB7A caused upregulation of some oncogenes such as matrix metalloproteinase-10 (MMP10) [56] and lysyl oxidase-like 2 (LOXL2) [57], and downregulation of some tumor suppressor genes such as suppressor of cytokine signalling (SOCS3) [58] and fibulin 1 (FBLN1) [59]. The results indicate that ZBTB7A possibly plays an oncosuppressive role. However, loss of ZBTB7A also caused downregulation of some oncogenes such as myosin light chain kinase (MYLK) [60] and chloride intracellular channel 5 (CLIC5) [61] and upregulation of some tumor suppressor genes such as interleukin-24 (IL24) [62] and arachidonate 15-lipoxygenase, type B (ALOX15B) [63]. The results indicate that ZBTB7A possibly plays an oncogenic role. Interestingly, we find that some genes have opposite effects in carcinogenesis, such as cellular retinoic acid binding protein 2 (CRABP2) [64, 65], epithelial membrane protein 3 (EMP3) [66, 67], transglutaminase 2 (TGM2) [68, 69], Yes associated protein 1 (YAP1) [70, 71], CDKN2A interacting protein (CDKN2AIP) [72], and BRCA2 and CDKN1A interacting protein (BCCIP) [73]. The genes can promote or inhibit cells growth in different conditions or stages of cancer (Supplementary Materials, Figure S3 and Table S1).

We went on disclosing that carbohydrate and lipid metabolisms probably involved in NPC progression through KEGG analysis. In order to search the molecular mechanisms of ZBTB7A, we selected 12 differentially expressed lncRNAs and their nearby mRNAs to validate. 12 differentially expressed mRNAs included FASN, CDKN2B, SLC2A1, ALDOC, ISG15, IER3, TXNIP, KAT5, ITGA5, DYRK1A, SIK1, and PON3. Basing on the bioinformatics analysis, the lncRNAs and mRNAs are possibly associated with cancer cells progression via carbohydrate or lipid pathway. In order to exclude the contaminated possibility by Hela cell, parts of them were validated in NPC and rhinitis tissues by qPCR. lncRNA NR_047538, lncRNA ENST00000442852, and fatty acid synthase (FASN) were closely associated with NPC. ZBTB7A had a positive association with NR_047538, while having negative associations with ENST00000442852 and FASN. The results provide novel candidate biomarkers for NPC progression with different levels of ZBTB7A.

LncRNA NR_047538 and FSAN were important markers and regulatory factor in NPC [31, 74]. LncRNA NR_047538 can promote NPC progression by upregulating GLUT1 expression [31]. It is also associated with lipid pathways [75]. Sex determining region Y-box 2 gene (SOX2) induces proliferation of NPC cells through activating lncRNA NR_047538 [76]. Downregulation of lncRNA NR_047538 inhibits NPC tumorigenicity and enhances the efforts of radiotherapy and chemotherapy via regulating microRNA 125a and let-7a [77, 78]. FASN is an important enzyme of lipogenesis. It promotes NPC progression through remarkably providing endogenous fatty acid [52]. Interestingly, we found that sterol regulatory element binding transcription factor 1 (SREBF1) also was upregulated by lncRNA microarray and qPCR (Supplementary Materials, Figure S4 and Table S1). SREBF1 is the transcript of sterol regulatory element binding protein 1 (SREBP1). Activated SREBP1 enters nucleus and transcribes the genes of lipid metabolism such as FASN [79]. EBV-encoded latent membrane protein 1 (LMP1) induces cell proliferation and NPC metastasis via activating SREBP1 and its downstream FASN [80].

The potential biomarker lncRNA ENST00000442852 has not been reported. We speculate it may be connected with the nearby mRNA IER3, which is also called immediate early response gene X-1 (IEX-1). It is immediately regulated by transcriptional factors, inflammatory cytokines, growth factors, and so on [81]. It is connected with poor or good prognosis in different types of cancers [82, 83]. Therefore, the functions and characteristics of lncRNA ENST00000442852 will be explored in NPC.

Although SLC2A1, CDKN2B, and ISG15 are closely connected with NPC progression in other studies [31, 84, 85], our results do not show the positive results. The quantity of 20 chronic rhinitis and 60 NPC tissues is limited. The genotyping of Guangxi population probably is different from that of other region. The protein expression of the genes should be validated by Western blot or immunohistochemistry in future because of the difference between mRNA and protein expression of the same gene.

5. Conclusions

Above all, the results showed the changes of the differentially expressed lncRNAs and mRNAs with stable loss of ZBTB7A expression in shRNA-CNE2 cells. The results may provide potential biomarkers for NPC progression with different levels of ZBTB7A, such as lncRNA NR_047538, lncRNA ENST00000442852, and FASN. In the future, we will deeply dig the connections between the lncRNAs/mRNAs and ZBTB7A in NPC.

Data Availability

The data used to support the findings of this study are included in the link of https://drive.google.com/drive/folders/12uRcwskYElRlXr4j_rdVRpO4n61yNFLR?usp=sharing.

Conflicts of Interest

The authors have no conflicts of interests to declare.

Acknowledgments

The authors are very grateful to Professor Musheng Zeng from Sun Yat-Sen University for generously providing the cell lines 5-8F and 6-10B and Professor Sai-Wah Tsao from University of Hong Kong for kindly providing the cell line NP69. The study was supported by the Natural Science Foundation of Guangxi Province (2015GXNSFAA139166, 2016GXNSFBA380144, 2016GXNSFCB380003).

Supplementary Materials

Supplementary 1. Figure S1. The protein levels of ZBTB7A in NPC cell lines.

Supplementary 2. Figure S2. The vitality, migration, and invasion of 6-10B and CNE3 stably transfected cells with ZBTB7A overexpression.

Supplementary 3. Figure S3. Screening and validation of 14 differentially expressed mRNAs by microarray and qPCR.

Supplementary 4. Figure S4. Screening and validation of differentially expressed SREBF1 by microarray and qPCR.

Supplementary 5. Table S1. List of specific primers of 15 differentially expressed mRNA designed utilizing primer 5.0.

References

W. Zhao, H. Lei, X. Zhu, L. Li, S. Qu, and X. Liang, “Investigation of long-term survival outcomes and failure patterns of patients with nasopharyngeal carcinoma receiving intensity-modulated radiotherapy: A retrospective analysis,” Oncotarget , vol. 7, no. 52, pp. 86914–86925, 2016.
View at: Google Scholar
W. Ye, E. T. Chang, Z. Liu et al., “Development of a population-based cancer case-control study in southern china,” Oncotarget, vol. 8, no. 50, pp. 87073–87085, 2017.
View at: Publisher Site | Google Scholar
W.-W. Sung, P.-R. Chen, M.-H. Liao, and J.-W. Lee, “Enhanced aerobic glycolysis of nasopharyngeal carcinoma cells by Epstein-Barr virus latent membrane protein 1,” Experimental Cell Research, vol. 359, no. 1, pp. 94–100, 2017.
View at: Publisher Site | Google Scholar
K. Wang, Y. Ge, C. Ni et al., “Epstein-Barr virus-induced up-regulation of TCAB1 is involved in the DNA damage response in nasopharyngeal carcinoma,” Scientific Reports, vol. 7, no. 1, p. 3218, 2017.
View at: Publisher Site | Google Scholar
J. C. Lin, W. Y. Wang, and K. Y. Chen, “Quantification of plasma Epstein-Barr virus DNA in patients with advanced nasopharyngeal carcinoma,” The New England Journal of Medicine, vol. 350, pp. 2461–2470, 2004.
View at: Publisher Site | Google Scholar
X. Hou, L. Zhang, C. Zhao et al., “Prognostic impact of plasma Epstein-Barr virus DNA concentration on distant metastasis in nasopharyngeal carcinoma,” Ai Zheng, vol. 25, no. 7, pp. 785–792, 2006.
View at: Google Scholar
V. E. Gurtsevitch, N. B. Senyuta, A. V. Ignatova et al., “Epstein-barr virus biomarkers for nasopharyngeal carcinoma in non-endemic regions,” Journal of General Virology, vol. 98, no. 8, Article ID 000889, pp. 2118–2127, 2017.
View at: Publisher Site | Google Scholar
M. Chen, L. Yin, J. Wu et al., “Impact of plasma epstein-barr virus-DNA and tumor volume on prognosis of locally advanced nasopharyngeal carcinoma,” BioMed Research International, vol. 2015, Article ID 617949, 5 pages, 2015.
View at: Publisher Site | Google Scholar
L. Lu, J. Li, C. Zhao et al., “Prognostic efficacy of combining tumor volume with Epstein-Barr virus DNA in patients treated with intensity-modulated radiotherapy for nasopharyngeal carcinoma,” Oral Oncology, vol. 60, pp. 18–24, 2016.
View at: Publisher Site | Google Scholar
J. A. Costoya, “Functional analysis of the role of POK transcriptional repressors,” Briefings in Functional Genomics & Proteomics, vol. 6, no. 1, pp. 8–18, 2007.
View at: Publisher Site | Google Scholar
T. Maeda, R. M. Hobbs, T. Morghoub et al., “Role of the proto-oncogene Pokemon in cellular transformation and ARF repression,” Nature, vol. 433, no. 7023, pp. 278–285, 2005.
View at: Publisher Site | Google Scholar
Z.-H. Zhao, S.-F. Wang, L. Yu et al., “Overexpression of Pokemon in non-small cell lung cancer and foreshowing tumor biological behavior as well as clinical results,” Lung Cancer, vol. 62, no. 1, pp. 113–119, 2008.
View at: Publisher Site | Google Scholar
L. Jiang, M. K. Siu, O. G. Wong et al., “Overexpression of proto-oncogene FBI-1 activates membrane type 1-matrix metalloproteinase in association with adverse outcome in ovarian cancers,” Molecular Cancer, vol. 9, no. 1, p. 318, 2010.
View at: Publisher Site | Google Scholar
H. Qu, D. Qu, F. Chen, Z. Zhang, B. Liu, and H. Liu, “ZBTB7 overexpression contributes to malignancy in breast cancer,” Cancer Investigation, vol. 28, no. 6, pp. 672–678, 2010.
View at: Publisher Site | Google Scholar
F. Fang, L. Yang, Y. Tao, and W. Qin, “FBI-1 promotes cell proliferation and enhances resistance to chemotherapy of hepatocellular carcinoma in vitro and in vivo,” Cancer, vol. 118, no. 1, pp. 134–146, 2012.
View at: Publisher Site | Google Scholar
L. Zhang, Y. Wang, L. Zhang et al., “ZBTB7A, a miR-663a target gene, protects osteosarcoma from endoplasmic reticulum stress-induced apoptosis by suppressing LncRNA GAS5 expression,” Cancer Letters, vol. 448, pp. 105–116, 2019.
View at: Publisher Site | Google Scholar
G. Wang, A. Lunardi, J. Zhang et al., “Zbtb7a suppresses prostate cancer through repression of a Sox9-dependent pathway for cellular senescence bypass and tumor invasion,” Nature Genetics, vol. 45, no. 7, pp. 739–746, 2013.
View at: Publisher Site | Google Scholar
X.-S. Liu, M. D. Genet, J. E. Haines et al., “Zbtb7a suppresses melanoma metastasis by transcriptionally repressing mcam,” Molecular Cancer Research, vol. 13, no. 8, pp. 1206–1217, 2015.
View at: Publisher Site | Google Scholar
Y. Zhao, Y. Yao, L. Li et al., “Pokemon enhances proliferation, cell cycle progression and anti-apoptosis activity of colorectal cancer independently of p14ARF–MDM2–p53 pathway,” Medical Oncology, vol. 31, no. 12, p. 288, 2014.
View at: Publisher Site | Google Scholar
X. Liu, J. E. Haines, E. K. Mehanna et al., “ZBTB7A acts as a tumor suppressor through the transcriptional repression of glycolysis,” Genes & Development, vol. 28, no. 17, pp. 1917–1928, 2014.
View at: Publisher Site | Google Scholar
D.-B. Shi, Y.-W. Wang, A.-Y. Xing et al., “C/EBPα-induced miR-100 expression suppresses tumor metastasis and growth by targeting ZBTB7A in gastric cancer,” Cancer Letters, vol. 369, no. 2, pp. 376–385, 2015.
View at: Publisher Site | Google Scholar
G. Sun, B. Peng, Q. Xie, J. Ruan, and X. Liang, “Upregulation of ZBTB7A exhibits a tumor suppressive role in gastric cancer cells,” Molecular Medicine Reports, vol. 17, no. 2, pp. 2635–2641, 2018.
View at: Google Scholar
W. Jiao, F. Liu, F.-Z. Tang et al., “Expression of the pokemon proto-oncogene in nasopharyngeal carcinoma cell lines and tissues,” Asian Pacific Journal of Cancer Prevention, vol. 14, no. 11, pp. 6315–6319, 2013.
View at: Publisher Site | Google Scholar
F. Liu, F. Tang, J. Lan et al., “Stable knockdown of ZBTB7A promotes cell proliferation and progression in nasopharyngeal carcinoma,” Tumori, vol. 104, no. 1, pp. 37–42, 2018.
View at: Publisher Site | Google Scholar
F. Liu, J. Lan, W. Jiao et al., “Differences in Zbtb7a expression cause heterogeneous changes in human nasopharyngeal carcinoma CNE3 sublines,” Oncology Letters, vol. 14, no. 3, pp. 2669–2676, 2017.
View at: Publisher Site | Google Scholar
K. C. Wang and H. Y. Chang, “Molecular mechanisms of long noncoding RNAs,” Molecular Cell, vol. 43, no. 6, pp. 904–914, 2011.
View at: Publisher Site | Google Scholar
X. Wen, X. Tang, Y. Li et al., “Microarray expression profiling of long non-coding RNAs involved in nasopharyngeal carcinoma metastasis,” International Journal of Molecular Sciences, vol. 17, no. 11, p. 1956, 2016.
View at: Publisher Site | Google Scholar
Y. Nie, X. Liu, S. Qu, E. Song, H. Zou, and C. Gong, “Long non-coding RNA HOTAIR is an independent prognostic marker for nasopharyngeal carcinoma progression and survival,” Cancer Science, vol. 104, no. 4, pp. 458–464, 2013.
View at: Publisher Site | Google Scholar
H. Bo, Z. Gong, W. Zhang et al., “Upregulated long non-coding RNA AFAP1-AS1 expression is associated with progression and poor prognosis of nasopharyngeal carcinoma,” Oncotarget , vol. 6, no. 24, pp. 20404–20418, 2015.
View at: Google Scholar
L. Yang, Y. Tang, Y. He et al., “High Expression of LINC01420 indicates an unfavorable prognosis and modulates cell migration and invasion in nasopharyngeal carcinoma,” Journal of Cancer, vol. 8, no. 1, pp. 97–103, 2017.
View at: Publisher Site | Google Scholar
Z. W. Zou, C. Ma, L. Medoro et al., “LncRNA ANRIL is up-regulated in nasopharyngeal carcinoma and promotes the cancer progression via increasing proliferation, reprograming cell glucose metabolism and inducing sidepopulation stem-like cancer cells,” Oncotarget , vol. 7, no. 38, pp. 61741–61754, 2016.
View at: Google Scholar
M. Shen, M. Li, and J. Liu, “Long noncoding RNA HOTTIP promotes nasopharyngeal cancer cell proliferation, migration, and invasion by inhibiting miR-4301,” Medical Science Monitor, vol. 25, pp. 778–785, 2019.
View at: Publisher Site | Google Scholar
Y. Hao, H. Zhao, X. Jin et al., “Long noncoding RNA DANCR promotes nasopharyngeal carcinoma cell proliferation and migration,” Molecular Medicine Reports, 2019.
View at: Publisher Site | Google Scholar
W.-P. Chak, R. W.-M. Lung, J. H.-M. Tong et al., “Downregulation of long non-coding RNA MEG3 in nasopharyngeal carcinoma,” Molecular Carcinogenesis, vol. 56, no. 3, pp. 1041–1054, 2017.
View at: Publisher Site | Google Scholar
J. Guo, J. Ma, G. Zhao et al., “Long noncoding RNA LINC0086 functions as a tumor suppressor in nasopharyngeal carcinoma by targeting miR-214,” Oncology Research : Featuring Preclinical and Clinical Cancer Therapeutics, vol. 25, no. 7, pp. 1189–1197, 2017.
View at: Publisher Site | Google Scholar
L. Zhang, Y. Wang, X. Li et al., “ZBTB7A enhances osteosarcoma chemoresistance by transcriptionally repressing lncRNALINC00473-IL24 activity,” Neoplasia (United States), vol. 19, no. 11, pp. 908–918, 2017.
View at: Google Scholar
Z. Zhao, J. Wang, S. Wang, H. Chang, T. Zhang, and J. Qu, “LncRNA CCAT2 promotes tumorigenesis by over-expressed Pokemon in non-small cell lung cancer,” Biomedicine & Pharmacotherapy, vol. 87, pp. 692–697, 2017.
View at: Publisher Site | Google Scholar
S. Y. Gu, W. P. Tang, Y. Zeng et al., “An epithelial cell line established from poorly differentiated nasopharyngeal carcinoma,” Ai Zheng, vol. 2, no. 2, pp. 70–72, 1983.
View at: Google Scholar
W. Jiao, “Establishment of a human epithelial cell line of nasopharyngeal carcinoma – CNE3 and its biological characteristics,” Journal of Guangxi Medical University, vol. 12, no. 2, pp. 187–190, 1995.
View at: Google Scholar
F. Liu, W. Jiao, X.-L. Mo et al., “Molecular pathological study of the human nasopharyngeal carcinoma CNE3 cell line,” Oncology Letters, vol. 6, no. 4, pp. 980–984, 2013.
View at: Publisher Site | Google Scholar
W. Chen, Y. Lee, H. Wang et al., “Suppression of human nasopharyngeal carcinoma cell growth in nude mice by the wild-type p53 gene,” Journal of Cancer Research and Clinical Oncology, vol. 119, no. 1, pp. 46–48, 1992.
View at: Publisher Site | Google Scholar
Z.-P. Teng, T. Ooka, D. P. Huang, and Y. Zeng, “Detection of Epstein-Barr Virus DNA in well and poorly differentiated nasopharyngeal carcinoma cell lines,” Virus Genes, vol. 13, no. 1, pp. 53–60, 1996.
View at: Publisher Site | Google Scholar
M. A. Faghihi and C. Wahlestedt, “Regulatory roles of natural antisense transcripts,” Nature Reviews Molecular Cell Biology, vol. 10, no. 9, pp. 637–643, 2009.
View at: Publisher Site | Google Scholar
U. A. Ørom, T. Derrien, M. Beringer et al., “Long noncoding RNAs with enhancer-like function in human cells,” Cell, vol. 143, no. 1, pp. 46–58, 2010.
View at: Publisher Site | Google Scholar
M. N. Cabili, C. Trapnell, L. Goff et al., “Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses,” Genes & Development, vol. 25, no. 18, pp. 1915–1927, 2011.
View at: Publisher Site | Google Scholar
Q. Wang, C. Li, P. Tang, R. Ji, S. Chen, and J. Wen, “A minimal lncRNA-mRNA signature predicts sensitivity to neoadjuvant chemotherapy in triple-negative breast cancer,” Cellular Physiology and Biochemistry, vol. 48, no. 6, pp. 2539–2548, 2018.
View at: Publisher Site | Google Scholar
K. Su, Q. Zhao, A. Bian, C. Wang, Y. Cai, and Y. Zhang, “A novel positive feedback regulation between long noncoding RNA UICC and IL-6/STAT3 signaling promotes cervical cancer progression,” American Journal of Cancer Research, vol. 8, no. 7, pp. 1176–1189, 2018.
View at: Google Scholar
W. Luan, Z. Zhou, X. Ni et al., “Long non-coding RNA H19 promotes glucose metabolism and cell growth in malignant melanoma via miR-106a-5p/E2F3 axis,” Journal of Cancer Research and Clinical Oncology, vol. 144, no. 3, pp. 531–542, 2018.
View at: Publisher Site | Google Scholar
Z. Xiao, L. Han, H. Lee et al., “Energy stress-induced lncRNA FILNC1 represses c-Myc-mediated energy metabolism and inhibits renal tumor development,” Nature Communications, vol. 8, no. 1, 2017.
View at: Publisher Site | Google Scholar
X. Liu, Y. Liang, R. Song et al., “Long non-coding RNA NEAT1-modulated abnormal lipolysis via ATGL drives hepatocellular carcinoma proliferation,” Molecular Cancer, vol. 17, no. 1, p. 90, 2018.
View at: Publisher Site | Google Scholar
L. L. Christensen, K. True, M. P. Hamilton et al., “SNHG16 is regulated by the Wnt pathway in colorectal cancer and affects genes involved in lipid metabolism,” Molecular Oncology, vol. 10, no. 8, pp. 1266–1282, 2016.
View at: Publisher Site | Google Scholar
D. D. Ma, L. L. Yuan, and L. Q. Lin, “LncRNA HOTAIR contributes to the tumorigenesis of nasopharyngeal carcinoma via up-regulating FASN,” European Review for Medical and Pharmacological Sciences, vol. 21, no. 22, pp. 5143–5152, 2017.
View at: Google Scholar
H. Zhang, Y. Li, H. Wang et al., “Author Correction: Ephrin receptor A2 is an epithelial cell receptor for Epstein–Barr virus entry,” Nature Microbiology, vol. 3, no. 2, pp. 1–8, 2018.
View at: Publisher Site | Google Scholar
D. Zong, L. Yin, Q. Zhong et al., “ZNF488 enhances the invasion and tumorigenesis in nasopharyngeal carcinoma via the Wnt signaling pathway involving epithelial mesenchymal transition,” Cancer Research and Treatment, vol. 48, no. 1, pp. 334–344, 2016.
View at: Publisher Site | Google Scholar
X.-K. Qi, H.-Q. Han, H.-J. Zhang et al., “OVOL2 links stemness and metastasis via fine-tuning epithelial-mesenchymal transition in nasopharyngeal carcinoma,” Theranostics, vol. 8, no. 8, pp. 2202–2216, 2018.
View at: Publisher Site | Google Scholar
T. Mariya, Y. Hirohashi, T. Torigoe et al., “Matrix metalloproteinase-10 regulates stemness of ovarian cancer stem-like cells by activation of canonical Wnt signaling and can be a target of chemotherapy-resistant ovarian cancer,” Oncotarget , vol. 7, no. 18, pp. 26806–26822, 2016.
View at: Google Scholar
A. Martin, F. Salvador, G. Moreno-Bueno et al., “Lysyl oxidase-like 2 represses Notch1 expression in the skin to promote squamous cell carcinoma progression,” EMBO Journal, vol. 34, no. 8, pp. 1090–1109, 2015.
View at: Publisher Site | Google Scholar
K. Yuan, Y. Lei, H. Chen et al., “HBV-induced ROS accumulation promotes hepatocarcinogenesis through Snail-mediated epigenetic silencing of SOCS3,” Cell Death & Differentiation, vol. 23, no. 4, pp. 616–627, 2016.
View at: Publisher Site | Google Scholar
Z. Xu, H. Chen, D. Liu, and J. Huo, “Fibulin-1 Is downregulated through promoter hypermethylation in colorectal cancer,” Medicine, vol. 94, no. 13, p. e663, 2015.
View at: Publisher Site | Google Scholar
V. Sundararajan, N. Gengenbacher, M. P. Stemmler, J. A. Kleemann, T. Brabletz, and S. Brabletz, “The ZEB1/miR-200c feedback loop regulates invasion via actin interacting proteins MYLK and TKS5,” Oncotarget , vol. 6, no. 29, pp. 27083–27096, 2015.
View at: Google Scholar
T. N. J. Flores-Téllez, T. V. Lopez, V. R. Vásquez Garzón, and S. Villa-Treviño, “Co-expression of Ezrin-CLIC5-Podocalyxin is associated with migration and invasiveness in hepatocellular carcinoma,” PLoS ONE, vol. 10, no. 7, Article ID e0131605, 2015.
View at: Google Scholar
J. Panneerselvam, J. Jin, M. Shanker et al., “IL-24 inhibits lung cancer cell migration and invasion by disrupting The SDF-1/CXCR4 signaling axis,” PLoS ONE, vol. 10, no. 3, p. e0122439, 2015.
View at: Publisher Site | Google Scholar
Y. Liu, C. Chen, Z. Xu et al., “Deletions linked to TP53 loss drive cancer through p53-independent mechanisms,” Nature, vol. 531, no. 7595, pp. 471–475, 2016.
View at: Publisher Site | Google Scholar
R.-Z. Liu, S. Li, E. Garcia et al., “Association between cytoplasmic CRABP2, altered retinoic acid signaling, and poor prognosis in glioblastoma,” Glia, vol. 64, no. 6, pp. 963–976, 2016.
View at: Google Scholar
Q. Yang, R. Wang, W. Xiao, F. Sun, H. Yuan, and Q. Pan, “Cellular retinoic acid binding protein 2 is strikingly downregulated in human esophageal squamous cell carcinoma and functions as a tumor suppressor,” PLoS ONE, vol. 11, no. 2, Article ID e0148381, 2016.
View at: Google Scholar
Y.-H. Hsieh, S.-C. Hsieh, C.-H. Lee et al., “Targeting EMP3 suppresses proliferation and invasion of hepatocellular carcinoma cells through inactivation of PI3K/Akt pathway,” Oncotarget , vol. 6, no. 33, pp. 34859–34874, 2015.
View at: Google Scholar
Q. Xue, Y. Zhou, C. Wan et al., “Epithelial membrane protein 3 is frequently shown as promoter methylation and functions as a tumor suppressor gene in non-small cell lung cancer,” Experimental and Molecular Pathology, vol. 95, no. 3, pp. 313–318, 2013.
View at: Publisher Site | Google Scholar
S. Y. Yeo, Y. Itahana, A. K. Guo et al., “Transglutaminase 2 contributes to a TP53-induced autophagy program to prevent oncogenic transformation,” eLife, vol. 5, 2016.
View at: Publisher Site | Google Scholar
D. T. Leicht, T. Kausar, Z. Wang et al., “TGM2: A cell surface marker in esophageal adenocarcinomas,” Journal of Thoracic Oncology, vol. 9, no. 6, pp. 872–881, 2014.
View at: Publisher Site | Google Scholar
W. Lehmann, D. Mossmann, J. Kleemann et al., “ZEB1 turns into a transcriptional activator by interacting with YAP1 in aggressive cancer types,” Nature Communications, vol. 7, Article ID 10498, 2016.
View at: Publisher Site | Google Scholar
M. Yuan, V. Tomlinson, R. Lara et al., “Yes-associated protein (YAP) functions as a tumor suppressor in breast,” Cell Death & Differentiation, vol. 15, no. 11, pp. 1752–1759, 2008.
View at: Publisher Site | Google Scholar
C. T. Cheung, R. Singh, R. S. Kalra, S. C. Kaul, and R. Wadhwa, “Collaborator of ARF (CARF) regulates proliferative fate of human cells by dose-dependent regulation of DNA damage signaling,” The Journal of Biological Chemistry, vol. 289, no. 26, pp. 18258–18269, 2014.
View at: Publisher Site | Google Scholar
Y.-Y. Huang, L. Dai, D. Gaines et al., “BCCIP suppresses tumor initiation but is required for tumor progression,” Cancer Research, vol. 73, no. 23, pp. 7122–7133, 2013.
View at: Publisher Site | Google Scholar
Y.-C. Kao, S.-W. Lee, L.-C. Lin et al., “Fatty acid synthase overexpression confers an independent prognosticator and associates with radiation resistance in nasopharyngeal carcinoma,” Tumor Biology, vol. 34, no. 2, pp. 759–768, 2013.
View at: Publisher Site | Google Scholar
G. Bochenek, R. Häsler, N.-E. E. Mokhtari et al., “The large non-coding RNA ANRIL, which is associated with atherosclerosis, periodontitis and several forms of cancer, regulates ADIPOR1, VAMP3 and C11ORF10,” Human Molecular Genetics, vol. 22, no. 22, Article ID ddt299, pp. 4516–4527, 2013.
View at: Publisher Site | Google Scholar
J. Wu, J. Tang, J. Li, and X. Li, “Upregulation of SOX2-activated lncRNA ANRIL promotes nasopharyngeal carcinoma cell growth,” Scientific Reports, vol. 8, no. 1, p. 3333, 2018.
View at: Publisher Site | Google Scholar
X. Hu, H. Jiang, and X. Jiang, “Downregulation of lncRNA ANRIL inhibits proliferation, induces apoptosis, and enhances radiosensitivity in nasopharyngeal carcinoma cells through regulating miR-125a,” Cancer Biology & Therapy, vol. 18, no. 5, pp. 331–338, 2017.
View at: Publisher Site | Google Scholar
Y. Wang, N. Cheng, and J. Luo, “Downregulation of lncRNA ANRIL represses tumorigenicity and enhances cisplatin-induced cytotoxicity via regulating microRNA let-7a in nasopharyngeal carcinoma,” Journal of Biochemical and Molecular Toxicology, vol. 31, no. 7, p. e21904, 2017.
View at: Publisher Site | Google Scholar
D. Guo, E. H. Bell, P. Mischel, and A. Chakravarti, “Targeting SREBP-1-driven lipid metabolism to treat cancer,” Current Pharmaceutical Design, vol. 20, no. 15, pp. 2619–2626, 2014.
View at: Publisher Site | Google Scholar
A. K.-F. Lo, R. W.-M. Lung, C. W. Dawson et al., “Activation of sterol regulatory element-binding protein 1 (SREBP1)-mediated lipogenesis by the Epstein–Barr virus-encoded latent membrane protein 1 (LMP1) promotes cell proliferation and progression of nasopharyngeal carcinoma,” The Journal of Pathology, vol. 246, no. 2, pp. 180–190, 2018.
View at: Publisher Site | Google Scholar
M. X. Wu, I. V. Ustyugova, L. Han, and O. E. Akilov, “Immediate early response gene X-1, a potential prognostic biomarker in cancers,” Expert Opinion on Therapeutic Targets, vol. 17, no. 5, pp. 593–606, 2013.
View at: Publisher Site | Google Scholar
J. Ye, Y. Zhang, Z. Cai et al., “Increased expression of immediate early response gene 3 protein promotes aggressive progression and predicts poor prognosis in human bladder cancer,” BMC Urology, vol. 18, no. 1, 2018.
View at: Publisher Site | Google Scholar
L. Han, L. Geng, X. Liu, H. Shi, W. He, and M. X. Wu, “Clinical significance of IEX-1 expression in ovarian carcinoma,” Ultrastructural Pathology, vol. 35, no. 6, pp. 260–266, 2011.
View at: Publisher Site | Google Scholar
H. Y. Wang, F. Li, and N. Liu, “Prognostic implications of a molecular classifier derived from whole-exome sequencing in nasopharyngeal carcinoma,” Cancer Medicine, 2019.
View at: Publisher Site | Google Scholar
R.-H. Chen, Y. Du, P. Han et al., “ISG15 predicts poor prognosis and promotes cancer stem cell phenotype in nasopharyngeal carcinoma,” Oncotarget , vol. 7, no. 13, pp. 16910–16922, 2016.
View at: Google Scholar

Copyright

Copyright © 2019 Fei Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

916

Downloads

1381

Citations