Skip to main content
Download PDF

Expression profile analysis of antisense long non-coding RNA identifies WDFY3-AS2 as a prognostic biomarker in diffuse glioma

Cancer Cell International volume 18, Article number: 107 (2018) Cite this article

Abstract

Background

Increasing evidence has shown that long non-coding RNAs (lncRNAs) are important prognostic biomarkers and epigenetic regulators with critical roles in cancer initiation and progression. However, the expression and clinical prognostic value of antisense lncRNAs in diffuse glioma patients remain unknown.

Methods

Here, we profiled differentially expressed antisense lncRNAs in glioma using RNA sequencing data from Chinese Glioma Genome Atlas database. Cox regression was performed to evaluate the prognostic value. Gene oncology (GO) and gene set enrichment analysis (GSEA) were used for functional analysis of antisense LncRNAs.

Results

Expression profiling identified 169 aberrantly expressed antisense lncRNAs between lower grade glioma (LGG) (grade II and III) and glioblastoma multiforme (GBM), 113 antisense lncRNAs between LGG IDH-wt and IDH-mut samples, and 70 antisense lncRNAs between GBM IDH-wt and IDH-mut samples, respectively. Among them, three antisense lncRNAs (WDFY3-AS2, MCM3AP-AS1 and LBX2-AS1) were significantly associated with prognosis and malignant progression of patients. WDFY3-AS2, the top one of downregulated antisense lncRNAs in GBM with fold change of 0.441 (P < 0.001), showed specific decreased expression in classical, mesenchymal, LGG IDH-wt, GBM IDH-wt or MGMT promoter unmethylated stratified patients. Chi square test found that WDFY3-AS2 was significantly associated with the clinical and molecular features of glioma. Univariate and multivariate Cox regression analysis indicated that WDFY3-AS2 was independently correlated with overall survival (OS) of patients. Kaplan–Meier analysis found that patients with high WDFY3-AS2 expression had longer OS than the low expression ones in the stratified cohorts. Additionally, GO and GSEA showed that gene sets correlated with WDFY3-AS2 expression were involved in regulation of synaptic transmission, glutamate receptor and TNF signaling pathway.

Conclusion

Our findings provided convincing evidence that WDFY3-AS2 is a novel valuable prognostic biomarker for diffuse glioma.

Background

Gliomas, the common primary brain tumor, mostly remain incurable, despite aggressive therapies with surgery, chemo and radiotherapy [1]. Glioblastoma multiforme (GBM) (IDH wild type) is the most aggressive glioma in adults characterized by glioma cell invasiveness and robust neovascularization, with a median survival of 15 months and 5-year survival rate less than 3% following standard treatment [2]. Due to the infiltrative growth and inherent resistance to chemo and radiotherapy, it is urgent to identify new targets to increase glioma patients’ survival. Although recent studies have found several genetic mutations and deregulated signaling pathways, these findings are still insufficient to uncover the molecular basis of gliomagenesis [3]. Since the coding genes only account for 1–2% of whole human genome [4], noncoding RNAs are emerging as new therapeutic targets in glioma.

Long non-coding RNAs (lncRNAs) are defined as transcripts longer than 200 nucleotides without protein-coding potential [5]. Antisense lncRNAs are a class of lncRNAs relative to their complementary protein-coding genes, which are orientated in an antisense direction [6]. Recent studies have demonstrated that antisense lncRNAs are involved in multiple biological processes, including cell proliferation, survival, migration, invasion and apoptosis through epigenetic modification and chromatin remodeling [7]. Moreover, accumulating reports found that dysregulated expression of antisense lncRNAs plays a crucial role in development and progression of various cancers. For example, it was reported that antisense lncRNA HAND2-AS1 was significantly down-regulated in endometrioid endometrial carcinoma, and HAND2-AS1 overexpression inhibited invasion and metastasis through inactivating neuromedin U [8]. Zhu et al. found that upregulation of lncRNA ZEB1-AS1 promoted tumor metastasis and predicted poor prognosis in hepatocellular carcinoma [9]. In glioma, Mineo et al. reported that lncRNA HIFA-AS2, upregulated in mesenchymal glioma stem cells, facilitated the maintenance of self-renewal in hypoxia niches [10]. Han et al. found that HOXA11-AS, a novel cell cycle-associated lncRNA, could serve as a biomarker and therapeutic target for glioma patients [11]. Our previous study reported that upregulation of lncRNA HOXA-AS3 promoted tumor progression and predicted poor prognosis in glioma [12]. Hence, understanding the role and prognostic value of cancer associated antisense lncRNAs would reveal new perspectives on the molecular basis of gliomagenesis.

In this study, we profiled differentially expressed antisense lncRNAs in diffuse glioma. Combined with prognosis and tumor progression, we identified a new antisense lncRNA WDFY3-AS2 with length of 3383 nt which is located in chromosome 4q21.23. WDFY3-AS2 downregulation was closely correlated with tumor grade and poor prognosis in patients. Bioinformatic analysis predicted that WDFY3-AS2 was involved in synaptic transmission, glutamate receptor and TNF signaling pathway. These results indicated that WDFY3-AS2 could be a potential prognostic biomarker for diffuse glioma.

Methods

Patients and tissues

309 diffuse glioma samples from the CGGA database were included in this study. Each sample was diagnosed by two neuropathologists according to the World Health Organization (WHO) classification guidelines of central nervous system tumors [13]. The study was approved by the ethics committee of Tiantan Hospital, and the written informed consent was obtained from all patients.

Datasets

The RNA sequencing data and corresponding clinical information (age, gender, histology, methylguanine methyltransferase (MGMT) promoter status, isocitrate dehydrogenase (IDH) mutation status and 1p/19q status) were downloaded from CGGA database (http://www.cgga.org.cn) [14]. The characteristics of glioma patients are summarized in Table 1. IDH mutation and MGMT promoter status were determined by DNA pyrosequencing as described in previous study [15]. 1p/19q status was determined using dual-color fluorescence in situ hybridization. According to the standard quality control criteria, the raw RNA sequencing data was first processed to define improper reads which would be removed. Then, the gene expression was estimated using the reads per kilobase transcriptome per million reads (RPKM). The calculated gene expression data was further used to compare the expression difference among samples [15].

Table 1 Clinical characteristics of glioma patients
Full size table

Gene oncology (GO), Kyoto encyclopedia of genes and genomes (KEGG) and gene set enrichment analysis (GSEA)

GO was applied to analyze the main function of differentially expression genes. KEGG was performed to analyze the pathway enrichment (http://david.ncifcrf.gov/). GSEA was carried out to identify gene sets of statistical difference between two groups by using GSEA v3 software (http://www.brodinstitute.org/gsea/index.jsp) [16].

Statistical analysis

Patients were divided into high-expression and low-expression groups based on WDFY3-AS2 level using the median value. Kaplan–Meier with 2-sided log-rank test was used to evaluate the OS differences between these two groups. Chi square test was performed to detect difference of the pathologic features between two groups of patients. Univariate and multivariate Cox regression analysis was conducted to identify independent prognostic factors. Significance analysis of microarray (SAM) was calculated the differential expression genes with R package “samr”. ROC curve analysis was used to predict OS with R package “pROC”. All statistical analyses were conducted using SPSS or R software. P < 0.05 was considered significant.

Results

Antisense lncRNAs profile in diffuse glioma using CGGA RNA sequencing data

The RNA sequencing data of 309 glioma samples was downloaded from CGGA database, and the clinicopathological features of the selected patients are listed in Table 1. Out of the whole transcriptome mRNAs, 420 antisense lncRNAs were identified. SAM with R package “samr” was performed to detect the differentially expressed antisense lncRNAs between LGG and GBM. Totally, 70 antisense lncRNAs were downregulated in GBM samples compared with LGG samples, while 99 antisense lncRNAs were upregulated in GBM samples. Volcano plot was used for showing the antisense lncRNAs of differential expression (Fig. 1a). Hierarchical clustering analysis conducted with R package “pheatmap” showed systematic variations in expression of these 169 antisense lncRNAs among samples (Fig. 1d). Considering the IDH mutation status, differentially expressed antisense lncRNAs within LGG and GBM samples stratified by IDH status were profiled. 113 antisense lncRNAs were abnormally expressed between LGG IDH-wt and IDH-mut samples (Fig. 1b, e), and 70 antisense lncRNAs were found between GBM IDH-wt and IDH-mut samples (Fig. 1c, f). These results suggested that the expression profiles of antisense lncRNAs were significantly different in lower grade gliomas and GBM.

Fig. 1
figure 1

Expression analysis of antisense lncRNAs in diffuse glioma. a–c Volcano plots show the expression of 420 antisense lncRNAs between LGG and GBM, LGG IDH-mut and IDH-wt, GBM IDH-mut and IDH-wt samples. d–f Heat maps show the differential expression of antisense lncRNAs between LGG and GBM, LGG IDH-mut and IDH-wt, GBM IDH-mut and IDH-wt samples, respectively

Full size image

Screening of prognosis and tumor progression-related antisense lncRNAs

In order to find out antisense lncRNAs involved in prognosis and malignant progression in glioma, Cox regression analysis was performed using R package “survival”. Based on the Wald P values, 248 antisense lncRNAs were correlated with patients’ overall survival. Meanwhile, SAM was used to assess the differentially expressed antisense lncRNAs between tumor grades. Among samples of different grade, 47 (II vs III), 168 (II vs IV) and 55 (III vs IV) abnormally expressed lncRNAs were obtained respectively (P < 0.05), of which 6 overlapping antisense lncRNAs associated with grade progression were identified (Additional file 1: Figure S1). Combining prognosis, tumor progression and IDH status, three antisense lncRNAs were ultimately identified by computing the intersection. Among them, LBX2-AS1 was significantly upregulated, while the other 2 lncRNAs (WDFY3-AS2 and MCM3AP-AS1) were significantly downregulated as tumor grade increased. The Venn diagram presented the workflow used to identify antisense lncRNAs associated with prognosis and tumor progression (Fig. 2a). Heat map showed the expression profile of these three antisense lncRNAs (Fig. 2b).

Fig. 2
figure 2

Screening of prognosis and tumor progression-related antisense lncRNAs. a Venn diagram presents the workflow used to identify antisense lncRNAs associated with prognosis and tumor progression. b Heat map shows the expression profile of three antisense lncRNAs (WDFY3-AS2, MCM3AP-AS1 and LBX2-AS1)

Full size image

WDFY3-AS2 preferred downregulation in malignant glioma

Among the obtained three antisense lncRNAs, WDFY3-AS2 was one of the most significantly downregulated lncRNAs in GBM with fold change of 0.441 (P < 0.001). we further detected its expression in patients stratified by grade, subtype, IDH, 1p/19q and MGMT promoter status. As shown in Fig. 3, WDFY3-AS2 expression was downregulated along with histological grades, and the decreased expression was also observed in classical, mesenchymal, MGMT promoter unmethylated, LGG IDH-wt or GBM IDH-wt stratified patients. Additionally, the statistical difference of these pathologic features between high-expression and low-expression group of WDFY3-AS2 (based on the medial value) was assessed with Chi square test. Except gender (P = 0.102), most of features were found different between these two groups (Additional file 1: Table S1). These findings indicated a significant correlation between WDFY3-AS2 expression and clinical factors.

Fig. 3
figure 3

WDFY3-AS2 expression in stratified patients. a WDFY3-AS2 expression among different pathologic grades. b WDFY3-AS2 expression among different subtypes of patients. c WDFY3-AS2 expression between MGMT promoter methylated and unmethylated patients. d WDFY3-AS2 expression between LGG IDH-mut and IDH-wt patients. e WDFY3-AS2 expression between GBM IDH-mut and IDH-wt patients

Full size image

WDFY3-AS2 showed strong prognostic power

We further evaluated the prognostic value of WDFY3-AS2 for diffuse glioma patients. By using the expression of WDFY3-AS2 as a continuous covariate, univariate and multivariate Cox regression analysis were performed to determine the prognostic significance. The results showed that WDFY3-AS2 expression could serve as an independent prognostic factor in glioma patients (95% CI 0.423–0.771, P < 0.001) (Table 2). Furthermore, we assessed the predictive accuracy through computing the AUC (area under the curve) of WDFY3-AS2 and pathologic features with ROC curve. As shown in Fig. 4, the AUC of WDFY3-AS2 (0.796) was much higher than that of age (0.611), grade (0.764), subtype (0.676), IDH status (0.706) and MGMT promoter status (0.629). These data demonstrated the powerful ability of WDFY3-AS2 expression for predicting overall survival of diffuse glioma patients.

Table 2 Univariate and multivariate Cox regression analysis of clinical pathologic features for OS in CGGA cohort
Full size table
Fig. 4
figure 4

ROC analysis of age, grade, subtype, IDH, MGMT and WDFY3-AS2 for overall survival. AUC area under curve, OS overall survival

Full size image

WDFY3-AS2 downregulation was associated with poor prognosis in diffuse glioma

To further explore the relationship between WDFY3-AS2 expression and patients’ outcome, we also assessed the correlation between WDFY3-AS2 expression and overall survival. Based on the median value of WDFY3-AS2 expression, patients were divided into high expression and low expression groups. Kaplan–Meier analysis found that patients in low expression group had shorter OS than ones in high expression group (P < 0.001, Fig. 5a). Moreover, the similar trend was observed in patients stratified by grade, IDH status and 1p/19q status although no significant difference was found in IDH-wt, IDH-mut and 1p/19q codeleted patients (probably due to the small sample size) (Fig. 5b–h). Our results indicated that low level of WDFY3-AS2 was significantly associated with unfavorable prognosis in patients with diffuse glioma.

Fig. 5
figure 5

Survival analysis in patients stratified by grade, IDH and 1p/19q status based on WDFY3-AS2 expression. a–c Kaplan-Meier analysis of  WDFY3-AS2 in all grade, LGG and GBM patients. d–f Kaplan-Meier analysis of WDFY3-AS2 in LGG stratified by IDH and 1p/19q status. g, h Kaplan-Meier analysis of WDFY3-AS2 in GBM IDH-wt and IDH-mut patients 

Full size image

Functional and pathway analysis of antisense lncRNA WDFY3-AS2

To explore the function of WDFY3-AS2 in glioma, the differentially expressed genes between high and low expression groups of patients were identified by R package “samr”. GO and KEGG analysis were performed based on the top 2000 genes positively and negatively associated with WDFY3-AS2. GO enrichment found that the top involved biological processes were chemical synaptic transmission, positive regulation of GTPase activity, neurotransmitter secretion, cell division, extracellular matrix organization and mitotic nuclear division (Fig. 6a). Pathway analysis showed the top significant pathways included insulin secretion, cAMP signaling pathway, phosphatidylinositol signaling system, ECM-receptor interaction, proteasome and cell cycle (Fig. 6b). In addition, GSEA revealed that low expression of WDFY3-AS2 was associated with synaptic transmission, glutamate receptor signaling, cAMP mediated signaling and l-amino acid transport (Fig. 6c), while high expression was enriched in tumor necrosis factor mediated signaling and macromolecule metabolic processing (Fig. 6d).

Fig. 6
figure 6

Functional and pathway analysis of antisense lncRNA WDFY3-AS2 in glioma. a GO annotations based on the top 2000 genes positively and negatively associated with WDFY3-AS2. b KEGG pathway analysis based on the top 2000 genes positively and negatively associated with WDFY3-AS2. c, d GSEA analysis based on the median value of WDFY3-AS2 expression level

Full size image

Discussion

Recently, large numbers of lncRNAs were found to be aberrantly expressed and involved in tumor initiation and progression in various cancers. Such lncRNAs could function as oncogenes and tumor suppressor genes and their expression could correlate with good or poor prognosis, making them valuable prognostic biomarkers. lncRNA HOTAIR has been used as a prognostic marker in different cancer types. High HOTAIR expression was significantly associated with poor overall survival and could serve as an independent prognostic factor [17]. It had been shown that MALAT1 was upregulated in many cancers and might act as a biomarker to predict survival in lung cancer [18]. lncRNA MEG3 was downregulated in a variety of primary cancers and found to be a prognostic factor for patients with glioma [19]. It was reported that lncRNA GAS5 expression was significantly reduced in liver cancer and could be an independent prognostic factor for patients [20]. In our study, we profiled differentially expressed antisense lncRNAs in glioma and identified a new antisense lncRNA WDFY3-AS2. Univariate and multivariate COX regression analysis found that WDFY3-AS2 downregulation was independently correlated with overall survival in patients, which indicated that WDFY3-AS2 could be a valuable prognostic biomarker for glioma.

As we known, increasing evidence has demonstrated that antisense lncRNAs are involved in regulation of a large range of biological processes through various mechanisms, such as transcriptional, post-transcriptional regulation and epigenetic modification. For example, antisense lncRNAs may regulate gene transcription in cis or trans through recruiting chromatin-modifying enzymes. Zhang et al. reported that lncRNA EZR-AS1 could enhance EZR transcription and expression by recruiting SMYD3 (a histone H3K4-methyltransferase) to its promoter region in esophageal squamous cell carcinoma [21]. lncRNA AGAP2-AS1 could bind with EZH2 and LSD1 and recruit them to KLF2 and LATS2 promoter regions to repress their transcription in non-small-cell lung cancer [22]. Antisense lncRNAs could function as competing endogenous RNAs (ceRNAs) to regulate targeted genes by binding with miRNAs. Wang et al. found lncRNA HOXD-AS1 promoted tumor progression by regulating the expression of frizzled family receptor 4 (FZD4) through competitively binding to miR-608 in ovarian cancer [23]. lncRNA TP73-As1 modulated cell proliferation through miR-200a dependent HMGB1/RAGE regulation in hepatocellular carcinoma [24]. In addition, antisense lncRNAs could also act as a molecular partner. Ding et al. reported that HNF1A-AS1 functioned as phosphatase activator through direct interaction with SHP-1 in hepatocellular carcinoma [25]. WDFY3, known as ALFY, encodes a phosphatidylinositol 3-phosphate-binding protein that function as a master conductor for aggregate clearance by autophagy [26]. Tschan reported that WDFY3 is critical for the granulocytic differentiation of AML cells [27]. Since antisense lncRNAs generally function as regulators of their complementary protein-coding genes [21], we first detected the expression correlation between WDFY3-AS2 and WDFY3 to evaluate the possible function and mechanism of WDFY3-AS2 in glioma. As shown in Additional file 1: Figure S1, WDFY3-AS2 expression was positively correlated with WDFY3 level in glioma (Pearson correlation coefficient r2 = 0.6418, P < 0.001), indicating a possible regulatory relationship between WDFY3-AS2 and WDFY3. Moreover, we also predicted the miRNAs that may interact with WDFY3-AS2 using LncBase v2 (https://omictools.com/lncrna2target-tool). mir-221, mir-26a, mir-135a, mir-9 and mir-139 showed the potential interaction with WDFY3-AS2. Among these miRNAs, mir-221 and mir-26a were overexpressed and could promote the malignant progression in glioma [28,29,30]. Therefore, we speculated that WDFY3-AS2 might modulate glioma cell malignant phenotype through interacting with mir-221 or mir-26a. GO and GSEA revealed that gene sets correlated with WDFY-AS2 expression were involved in regulation of synaptic transmission, glutamate receptor and TNF signaling pathways. Glutamate receptor signaling pathway plays important role in regulating glioma cell proliferation and invasion [31, 32]. Ciceroni found that type-3 glutamate receptors regulate chemoresistance in glioma stem cells, and their levels are inversely related to survival in patients with gliomas [33]. TNF signaling pathway is also an important regulator of the progression of glioma. Glioma cells tend to develop a resistance mechanism opposing to the TRAIL-induced apoptosis by overexpressing a wide variety of antiapoptotic proteins [34]. Deshayes reported that abnormal production of the TNF-homologue APRIL increases the proliferation of human malignant glioblastoma cell lines via a specific receptor [35]. Therefore, we proposed to determine whether WDFY-AS2 affected glioma cell behaviors through regulating glutamate receptor or TNF signaling pathway. Overall, further studies are required to clarify the function and molecular mechanism of WDFY3-AS2 in glioma.

Conclusion

In summary, we profiled abnormally expressed antisense lncRNAs in diffuse glioma and identified a new antisense lncRNA WDFY3-AS2 whose expression was closely correlated with tumor grade and poor prognosis in patients. Bioinformatic analysis predicted that WDFY3-AS2 was involved in synaptic transmission, glutamate receptor and TNF signaling pathway. Our data suggested that WDFY3-AS2 might be a potential prognostic biomarker for diffuse glioma.

References

  1. Ricard D, Idbaih A, Ducray F, Lahutte M, Hoang-Xuan K, Delattre JY. Primary brain tumours in adults. Lancet. 2012;379:1984–96.

    Article  PubMed  Google Scholar 

  2. Sadetzki S, Zach L, Chetrit A, Nass D, Hoffmann C, Ram Z, Zaaroor M, Umansky F, Rappaport ZH, Cohen A, Wald U, Rothman S, Hadani M. Epidemiology of gliomas in Israel: a nationwide study. Neuroepidemiology. 2008;31:264–9.

    Article  PubMed  Google Scholar 

  3. Ceccarelli M, Barthel FP, Malta TM, Sabedot TS, Salama SR, Murray BA, et al. Molecular profiling reveals biologically discrete subsets and pathways of progression in diffuse glioma. Cell. 2016;164:550–63.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Dunham I, Kundaje A, Aldred SF, Collins PJ, Davis CA, Doyle F, Epstein CB, Frietze S, Harrow J, Kaul R, Khatun J, Lajoie BR, Landt SG, et al. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.

    Article  CAS  Google Scholar 

  5. Ulitsky I, Bartel DP. lincRNAs: genomics, evolution, and mechanisms. Cell. 2013;154:26–46.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Isin M, Dalay N. LncRNAs and neoplasia. Clin Chim Acta. 2015;15:280–8.

    Article  CAS  Google Scholar 

  7. Batista PJ, Chang HY. Long noncoding RNAs: cellular address codes in development and disease. Cell. 2013;152:1298–307.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Yang X, Wang CC, Lee WYW, Trovik J, Chung TKH, Kwong J. Long non-coding RNA HAND2-AS1 inhibits invasion and metastasis in endometrioid endometrial carcinoma through inactivating neuromedin U. Cancer Lett. 2018;413:23–34.

    Article  PubMed  CAS  Google Scholar 

  9. Li T, Xie J, Shen C, Cheng D, Shi Y, Wu Z, Deng X, Chen H, Shen B, Peng C, Li H, Zhan Q, Zhu Z. Upregulation of long noncoding RNA ZEB1-AS1 promotes tumor metastasis and predicts poor prognosis in hepatocellular carcinoma. Oncogene. 2016;35:1575–84.

    Article  PubMed  CAS  Google Scholar 

  10. Mineo M, Ricklefs F, Rooj AK, Lyons SM, Ivanov P, Ansari KI, Nakano I, Chiocca EA, Godlewski J, Bronisz A. The long non-coding RNA-HIF1A-AS2 facilitates the maintenance of mesenchymal glioblastoma stem-like cells in hypoxic niches. Cell Rep. 2016;15:2500–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Wang Q, Zhang J, Liu Y, Zhang W, Zhou J, Duan R, Pu P, Kang C, Han L. A novel cell cycle-associated lncRNA, HOXA11-AS, is transcribed from the 5-prime end of the HOXA transcript and is a biomarker of progression in glioma. Cancer Lett. 2016;373:251–9.

    Article  PubMed  CAS  Google Scholar 

  12. Wu F, Zhang C, Cai J, Yang F, Liang T, Yan X, Wang H, Wang W, Chen J, Jiang T. Upregulation of long noncoding RNA HOXA-AS3 promotes tumor progression and predicts poor prognosis in glioma. Oncotarget. 2017;8:53110–23.

    PubMed  PubMed Central  Google Scholar 

  13. Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, Ohgaki H, Wiestler OD, et al. The 2016 World Health Organization classification of tumors of the central nervous system: a summary. Acta Neuropathol. 2016;131:803–20.

    Article  PubMed  Google Scholar 

  14. Zhang C, Cheng W, Ren XF, et al. Tumor purity as an underlying key factor in glioma. Clin Cancer Res. 2017;23:6279–91.

    Article  PubMed  CAS  Google Scholar 

  15. Yan W, Zhang W, You G, et al. Molecular classification of gliomas based on whole genome gene expression: a systematic report of 225 samples from the Chinese Glioma Cooperative Group. Neuro Oncol. 2012;14:1432–40.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Subramanian A, Tamayo P, Mootha VK, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102:15545–50.

    Article  PubMed  CAS  Google Scholar 

  17. Zhang S, Chen S, Yang G, Gu F, Li M, Zhong B, et al. Long noncoding RNA HOTAIR as an independent prognostic marker in cancer: a meta-analysis. PLoS ONE. 2014;9:e105538.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Müller-Tidow C, Diederichs S, Thomas M, Serve H. Genome-wide screening for prognosis-predicting genes in early-stage non-small-cell lung cancer. Lung Cancer. 2004;45:S145–50.

    Article  PubMed  Google Scholar 

  19. Zhou Y, Zhang X, Klibanski A. MEG3 noncoding RNA: a tumor suppressor. J Mol Endocrinol. 2012;48:R45–53.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Tu ZQ, Li RJ, Mei JZ, Li XH. Down-regulation of long non-coding RNA GAS5 is associated with the prognosis of hepatocellular carcinoma. Int J Clin Exp Pathol. 2014;7:4303–9.

    PubMed  PubMed Central  Google Scholar 

  21. Zhang XD, Huang GW, Xie YH, He JZ, Guo JC, Xu XE, Liao LD, Xie YM, Song YM, Li EM, Xu LY. The interaction of lncRNA EZR-AS1 with SMYD3 maintains overexpression of EZR in ESCC cells. Nucleic Acids Res. 2018;46:1793–809.

    Article  PubMed  Google Scholar 

  22. Li W, Sun M, Zang C, Ma P, He J, Zhang M, Huang Z, Ding Y, Shu Y. Upregulated long non-coding RNA AGAP2-AS1 represses LATS2 and KLF2 expression through interacting with EZH2 and LSD1 in non-small-cell lung cancer cells. Cell Death Dis. 2016;7:e2225.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Wang Y, Zhang W, Wang Y, Wang S. HOXD-AS1 promotes cell proliferation, migration and invasion through miR-608/FZD4 axis in ovarian cancer. Am J Cancer Res. 2018;8:170–82.

    PubMed  PubMed Central  Google Scholar 

  24. Li S, Huang Y, Huang Y, Fu Y, Tang D, Kang R, Zhou R, Fan XG. The long non-coding RNA TP73-AS1 modulates HCC cell proliferation through miR-200a-dependent HMGB1/RAGE regulation. J Exp Clin Cancer Res. 2017;36:51.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Ding CH, Yin C, Chen SJ, Wen LZ, Ding K, Lei SJ, Liu JP, Wang J, et al. The HNF1α-regulated lncRNA HNF1A-AS1 reverses the malignancy of hepatocellular carcinoma by enhancing the phosphatase activity of SHP-1. Mol Cancer. 2018;17:63.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Isakson P, Holland P, Simonsen A. The role of ALFY in selective autophagy. Cell Death Differ. 2013;20:12–20.

    Article  PubMed  CAS  Google Scholar 

  27. Schläfli AM, Isakson P, Garattini E, Simonsen A, Tschan MP. The autophagy scaffold protein ALFY is critical for the granulocytic differentiation of AML cells. Sci Rep. 2017;7:12980.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Quintavalle C, Garofalo M, Zanca C, Romano G, Iaboni M, De Caro MD, Martinez-Montero JC, et al. miR-221/222 overexpression in human glioblastoma increases invasiveness by targeting the protein phosphate PTPμ. Oncogene. 2012;31:858–68.

    Article  PubMed  CAS  Google Scholar 

  29. Zhang J, Han L, Ge Y, Zhou X, Zhang A, Zhang C, Zhong Y, You Y, Pu P, Kang C. miR-221/222 promote malignant progression of glioma through activation of the Akt pathway. Int J Oncol. 2010;36:913–20.

    PubMed  CAS  Google Scholar 

  30. Qian X, Zhao P, Li W, Shi ZM, Wang L, Xu Q, Wang M, et al. MicroRNA-26a promotes tumor growth and angiogenesis in glioma by directly targeting prohibitin. CNS Neurosci Ther. 2013;19:804–12.

    Article  PubMed  CAS  Google Scholar 

  31. Yoshida Y, Tsuzuki K, Ishiuchi S, Ozawa S. Serum-dependence of AMPA receptor-mediated proliferation in glioma cells. Pathol Int. 2006;56:262–71.

    Article  PubMed  CAS  Google Scholar 

  32. Ishiuchi S, Yoshida Y, Sugawara K, Aihara M, Ohtani T, Watanabe T, et al. Ca2+ -permeable AMPA receptors regulate growth of human glioblastoma via Akt activation. J Neurosci. 2007;27:7987–8001.

    Article  PubMed  CAS  Google Scholar 

  33. Ciceroni C, Bonelli M, Mastrantoni E, Niccolini C, Laurenza M, Larocca LM, et al. Type-3 metabotropic glutamate receptors regulate chemoresistance in glioma stem cells, and their levels are inversely related to survival in patients with malignant gliomas. Cell Death Differ. 2013;20:396–407.

    Article  PubMed  CAS  Google Scholar 

  34. Woo JS, Kim SM, Jeong CH, Ryu CH, Jeun SS. Lipoxygenase inhibitor MK886 potentiates TRAIL-induced apoptosis through CHOP- and p38 MAPK-mediated up-regulation of death receptor 5 in malignant glioma. Biochem Biophys Res Commun. 2013;431:354–9.

    Article  PubMed  CAS  Google Scholar 

  35. Deshayes F, Laprée G, Portier A, Richard Y, Pencalet P, et al. Abnormal production of the TNF-homologue APRIL increases the proliferation of human malignant glioblastoma cell lines via a specific receptor. Oncogene. 2004;23:3005–12.

    Article  PubMed  CAS  Google Scholar 

Download references

Authors’ contributions

FW designed the study and wrote the manuscript. RC and ZZ performed the gene analysis. KW, ZW, GL, HJ, RH and KZ collected the clinical data. All authors read and approved the final manuscript.

Acknowledgements

The authors accomplishing this work represent the Chinese Glioma Cooperative Group (CGCG).

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data used in this study were available from CGGA database (http://www.cgga.org.cn).

Consent for publication

Not applicable.

Ethics approval and consent to participate

The study was approved by the ethics committee of Tiantan Hospital, and the written informed consent was obtained from all patients.

Funding

This work was supported by National Natural Science Foundation of China (NSFC)/Research Grants Council (RGC) Joint Research Scheme (81761168038), the National Key Research and Development Plan (No. 2016YFC0902500) and National Natural Science Foundation of China (81672479, 81773208).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Author notes

  1. Fan Wu and Zheng Zhao contributed equally to this work

Authors and Affiliations

  1. Department of Molecular Neuropathology, Beijing Neurosurgical Institute, Capital Medical University, Beijing, China

    Fan Wu, Zheng Zhao, Ruichao Chai, Yuqing Liu, Kuanyu Wang, Zhiliang Wang, Guanzhang Li, Ruoyu Huang, Haoyu Jiang & Kenan Zhang

  2. Department of Neurosurgery, Beijing Tiantan Hospital, Capital Medical University, Beijing, China

    Fan Wu, Zheng Zhao, Ruichao Chai, Yuqing Liu, Kuanyu Wang, Zhiliang Wang, Guanzhang Li, Ruoyu Huang, Haoyu Jiang & Kenan Zhang

  3. Chinese Glioma Genome Atlas Network (CGGA) and Asian Glioma Genome Atlas Network (AGGA), Beijing, China

    Fan Wu, Zheng Zhao, Ruichao Chai, Yuqing Liu, Kuanyu Wang, Zhiliang Wang, Guanzhang Li, Ruoyu Huang, Haoyu Jiang & Kenan Zhang

  4. No. 6, Tiantan Xili, Dongcheng District, Beijing, 100050, China

    Fan Wu

Authors
  1. Fan Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Zheng Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ruichao Chai
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yuqing Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Kuanyu Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Zhiliang Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Guanzhang Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Ruoyu Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Haoyu Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Kenan Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Fan Wu.

Additional file

Additional file 1.

Additional figures and table.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, F., Zhao, Z., Chai, R. et al. Expression profile analysis of antisense long non-coding RNA identifies WDFY3-AS2 as a prognostic biomarker in diffuse glioma. Cancer Cell Int 18, 107 (2018). https://doi.org/10.1186/s12935-018-0603-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-018-0603-2

Keywords

Cancer Cell International

ISSN: 1475-2867

Contact us