Journal of Oncology
Volume 2019 (2019), Article ID 1785042, 15 pages
https://doi.org/10.1155/2019/1785042

Research Article

Low GAS5 Levels as a Predictor of Poor Survival in Patients with Lower-Grade Gliomas

Yanfang Wang,¹ Shan Xin,¹ Kai Zhang,² Run Shi,¹ and Xuanwen Bao3,4

¹Ludwig-Maximilians-Universität München (LMU), 80539 Munich, Germany
²Department of Cardiology, Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, 310016 Hangzhou, China
³Institute of Radiation Biology, Helmholtz Center Munich, German Research Center for Environmental Health, 85764 Neuherberg, Germany
⁴Technical University Munich (TUM), 80333 Munich, Germany

Correspondence should be addressed to Run Shi; shirun@outlook.com and Xuanwen Bao; xuanwen.bao@helmholtz-muenchen.de

Received 1 October 2018; Revised 18 December 2018; Accepted 3 January 2019; Published 3 February 2019

Academic Editor: Srikumar P. Chellappan

Copyright © 2019 Yanfang Wang 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.

Abstract.

Introduction. Gliomas are infiltrative neoplasms of a highly invasive nature. Different stages of gliomas feature distinct genomic, genetic, and epigenetic changes. The long noncoding RNA Growth Arrest Specific Transcript 5 (GAS5) is an identified tumour suppressor involved in several cancers. However, the underlying roles of the GAS5 gene in lower-grade glioma (LGG) patients are not clear. Methods. Via bioinformatic analysis based on TCGA-LGG and TCGA-GBM data, we explored the mechanisms of GAS5 expression in LGG (grades II and III) and high-grade glioma (glioblastoma multiforme, grade IV). The log-rank test and multivariate Cox analysis were performed to find the association between GAS5 and overall survival (OS) in LGG patients. Weighted gene coexpression network analysis (WGCNA) and RNA-Seq analysis were applied to find the key gene network associated with GAS5. Results. We found that GAS5 expression was downregulated in both LGG and glioblastoma multiforme (GBM) compared with normal brain tissue. Low methylation in the GAS5 promoter region was detected in both LGG and GBM tissues. The amplification type was the predominant type of GAS5 gene alteration in both LGG and GBM. High GAS5 expression was more associated with long overall survival (OS) in LGG patients than in GBM patients. The multivariate survival analysis of GAS5 and clinical and molecular characteristics in LGG patients further confirmed the association between GAS5 and OS in LGG patients. We then developed a nomogram for clinical use. WGCNA and RNA-Seq analysis indicated that ribosomal biogenesis and translation initiation were the predominant events regulated by GAS5 in LGG patients. Conclusion. Taken together, these results demonstrate that GAS5 expression is associated with OS in LGG patients and that its underlying roles involve the regulation of ribosomal biogenesis and translation initiation, which may aid in identifying a new target for the treatment of LGG.

1. Introduction

Gliomas are commonly classified as low-grade glioma (LGG, grades II and III) or glioblastoma multiforme (GBM, grade IV). LGGs are infiltrative neoplasms that arise most often in the cerebral hemispheres of adults, including oligodendrogliomas, oligoastrocytomas, and astrocytomas [1]. A subset of LGGs progress into GBM within a few months, and others may stay stable for several years, causing large variations in median survival [2, 3]. Hence, individual treatments should be performed based on the identification of histologic class, grade, and reactions to chemotherapy and radiotherapy [4]. GBM is the most common cause of death among children with central nervous system (CNS) neoplasms, and no effective therapies currently exist [5]. In most cases, GBM will recur after surgical resection, resulting in a poor prognosis. Hundreds of molecular alterations exist in LGG and GBM, making them react differently to chemotherapy, radiotherapy, and surgical resection. Identifying the key genes or genome alterations that drive the progression of gliomas will contribute to the understanding of the molecular mechanism behind gliomas and help to improve the effects of therapy.

Long noncoding RNAs (lncRNAs) are non-protein-coding transcripts that play essential roles in cellular regulation at various levels and in diverse biological functions, including chromatin modification, transcriptional regulation, cell differentiation, immune responses, and epigenetic regulation [6–8]. GAS5 produces a spliced lncRNA, which sensitizes cells to apoptosis by suppressing glucocorticoid-mediated induction of several responsive genes [9]. By binding to the DNA-binding domain of the glucocorticoid receptor, GAS5 acts as a decoy glucocorticoid response element (GRE), thus competing with DNA GREs to bind to the glucocorticoid receptor [9]. By analysing GAS5 level and clinical parameters, it has been found that GAS5 correlates with tumour progression and poor prognosis in many tumour types [10–13]. Some publications have also indicated that GAS5 is a predictor of survival in GBM patients [14, 15]. One study revealed an association between serum GAS5 level and the OS of GBM patients [14]. However, the potential roles of GAS5 involvement in gliomas and the regulatory mechanism of GAS5 expression in gliomas are not clearly identified. Besides, the difference of GAS5 expression and the regulatory mechanism of GAS5 expression in LGG and GBM patients have not been shown yet. In this study, via bioinformatic analysis, we explored the regulatory mechanisms of GAS5 expression in gliomas and compared its prognostic value in LGG and GBM. In addition, based on TCGA-LGG RNA-Seq datasets, we also applied WGCNA and enrichment analysis to identify the pathways highly associated with GAS5 in LGG and the network of genes highly associated with GAS5.

2. Materials and Methods

2.1. Clinical and Omic Data Download

RNA-Seq and DNA methylation datasets from patients with primary LGG or GBM were downloaded from the UCSC Xena browser (https://xenabrowser.net/) [16]. GAS5 mRNA expressions in different types of solid tumours and in corresponding normal tissues were obtained from TCGA datasets and GTEx datasets [17]. A box plot was created using ggplot2 [18].

2.2. WGCNA of RNA-Seq Data from LGG Patients and Network Visualization

The WGCNA R package was used to evaluate the correlation of GAS5 expression and module membership by the ‘p.weighted’ function [19]. An adjacency matrix was generated to evaluate the weighted coexpression values between all the subjects in the probe set. The coexpression similarity was defined as the absolute value of the correlation coefficient between the profiles of nodes and :

where and were expression value of for genes and and represented Pearson’s correlation coefficient of gene and gene .

A weighed network adjacency was defined by raising the coexpression similarity to a power :

with [20]. We selected the power of = 4 and scale free = 0.95 as the soft-thresholding parameter to ensure a signed scale-free coexpression gene network. The salmon and blue modules, which had most significant adjusted p-values, were selected. Genes involved in the salmon and blue modules and their weights as calculated by WGCNA were visualized by Cytoscape 3.4.0 [21]. GO analysis was performed using the clusterProfiler package [22] and Metascape (http://metascape.org/gp/index.html) [23]. Gene set enrichment analysis (GSEA) was applied by the GSEA software from Broad Institute [24].

2.3. Statistical Analysis

Statistical analysis was performed using R (R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL http://www.R-project.org/). The association between GAS5 expression and the clinicopathological features was evaluated by the different mathematical methods listed in the table [25]. The log-rank test was performed to assess the difference between the Kaplan-Meier curves. The fold change and Q value for samples with high and low GAS5 expression were calculated using the limma package [26]. A Q value < 0.05 was considered statistically significant.

3. Results

3.1. GAS5 Was Upregulated in Both LGG and GBM Compared with Normal Brain Tissues

GAS5 expression levels were characterized in several types of solid tumours, including LGG and GBM (Figure 1). The results indicated that GAS5 was downregulated in both LGG and GBM tissue compared with the normal brain tissue. Then, we compared the expression level of GAS5 between LGG and GBM (Figures 2(a), 2(b), and S1a). No significant difference was identified between the LGG and GBM datasets. In addition, grade II, grade III, and grade IV gliomas had similar expression levels of GAS5, which were all higher than that of normal brain tissues. The DNA methylation level of GAS5 was analysed by DNA methylation 450K chips, and the six probes were chosen from the CpG island in the promoter region of GAS5. All probes showed low β values (<0.1), and there were no differences in β values from any of the six probes between LGG and GBM tissues, indicating low methylation of the GAS5 promoter in both LGG and GBM patients (Figure 2(c)). The association between GAS5 methylation and GAS5 expression was tested, revealing no significant relationship in either GBM or LGG tissues (Figures S1B and S1C). Then, to investigate the driver of high GAS5 expression in LGG and GBM tissues, we also examined copy number alterations (CNAs) in LGG and GBM (Figures 2(d) and 2(e)). The results showed that the amplification-type alteration contributed to the high expression of GAS5 in both LGG and GBM tissues (P = 0.044 and P = 0.047).

Figure 1: GAS5 mRNA expression in different types of solid tumours and corresponding normal tissues. The RNA-Seq data of normal brain tissue were obtained from GTEx dataset.

(a)

(b)

(c)

(d)

(e)

Figure 2: GAS5 expression, DNA methylation, and copy number alteration (CNA) in LGG and in GBM. (a) Heatmap of GAS5 mRNA expression, copy number alteration, and DNA methylation in patients with primary LGG or GBM. Data were obtained from TCGA-LGG and TCGA-GBM. (b) Violin plots of GAS5 expression in LGG and GBM tissues. (c) Bar plots of GAS5 methylation level in LGG and in GBM tissues; six probes were chosen from the CpG island before the transcription start site of GAS5. (d-e) Box plots of GAS5 expression in LGG (d) and GBM (e) tissues with indications of genetic status. Data were obtained from TCGA-LGG and TCGA-GBM.

3.2. Low GAS5 Expression Was a Predictor of Poor Survival in LGG

The associations between GAS5 expression and the pathological parameters of patients with primary LGG and GBM were summarized in Table 1. In patients with LGG, higher GAS5 expression was significantly correlated with higher tumour purity (Table 1 and Figure S2). Nevertheless, we could not find a similar tendency in patients with GBM. Furthermore, in LGG patients, higher expression of GAS5 was associated with longer OS according to the log-rank test (Figure 3(a)). However, the level of GAS5 expression did not show a strong association with the OS of GBM patients (Figure 3(b)).

Table 1: The associations between GAS5 expression and pathological parameters in patients.

Attributes 
GAS expression and …
GBM
LGG

Statistic
P-value
FDR (BH)
Statistic
P-value
FDR (BH)

Ethnicity (Wilcox test)
0.096
0.620
0.843
0.002
0.487
0.487

Histological type (Kruskal-Wallis test)
2.467
0.291
0.843
3.918
0.141
0.197

Race (Kruskal-Wallis test)
1.110
0.574
0.843
4.180
0.243
0.283

Radiation therapy (Wilcox test)
-0.003
0.843
0.843
-0.025
0.019
0.045

Tumour purity (Spearman correlation)
-0.023
0.785
0.843
0.300
1.775e-11
1.243e-10

Age (Spearman correlation)
-0.064
0.446
0.843
-0.093
0.041
0.073

(a)

(b)

Figure 3: The associations between GAS5 expression and OS in LGG (a) and GBM (b) patients.

3.3. GAS5 Expression in LGG Patients with Different Clinical and Molecular Characteristics

To further identify the relationship between GAS5 expression and overall survival of LGG patients, we first examined GAS5 expression in LGG patients with different IDH1 status (Figure 4(a)), histologic grades (Figure 4(b)), tumour sizes (Figure 4(c)), and histologic types (Figure 4(d)). The results showed that the expression of GAS5 in IDH1 mutant patients was significantly higher than that in IDH1 wild-type patients. In contrast, histologic grade, tumour size, and histological type had no association with GAS5 expression. A Kaplan-Meier survival analysis stratified by clinicopathological risk factors was performed to verify the effect of GAS5 on survival. In patients with large or small tumours of histologic grade II or III, GAS5 showed a good predictive effect for survival (Figures 5(a)–5(d)).

(a)

(b)

(c)

(d)

Figure 4: GAS5 expression in LGG patients with different IDH1 statuses (a); histologic grade (b); tumour sizes (c); and histological type (d).

(a)

(b)

(c)

(d)

Figure 5: Kaplan-Meier survival analysis for LGG patients according to GAS5 expression levels, stratified by clinicopathological risk factors. (a, b) Tumour size. (c, d) Tumour grade. P values were calculated using the log-rank test.

3.4. Multivariate Survival Analysis of GAS5 with Clinical and Molecular Characteristics in LGG Patients

We applied the Cox regression model to conduct a multivariable survival analysis and used Cox regression coefficients to generate a nomogram (Figure 6(a)). In the multivariable survival analysis, we included histologic grade, gender, histological type, IDH mutation status, age, and GAS5 expression level as risk factors. The nomogram predicts LGG patients’ overall 3-year and 5-year survival probability. The results indicated that the expression level of GAS5 was highly associated with overall survival in LGG patients. Calibration plots showed that the nomograms performed well compared with an ideal model (Figure 6(b)). Then, the risk scores from Cox regression were assessed with a time-dependent receiver operating characteristic (ROC) curve (Figure 6(c)). The area under the curve (AUC) was 0.897 and 0.825 for 3 years and 5 years, respectively. The patients were separated into 2 groups, namely, those with high and low risk scores according to the Cox regression model. The Kaplan–Meier curve, which was generated by the log-rank test, showed good performance in predicting the OS of LGG patients (Figure 6(d)).

(a)

(b)

(c)

(d)

Figure 6: Multivariate survival analysis of GAS5 and clinical and molecular characteristics in LGG patients. Data are AUC (95% CI) or p values. ROC: receiver operating characteristic. AUC: area under the curve. (a) The nomogram to predict overall survival in LGG patients. (b) Plots depict the calibration of each model in terms of agreement between predicted and observed 3-year and 5-year outcomes. Model performance is shown by the plot, relative to the 45-degree line, which represents perfect prediction. (c) Time-dependent ROC curves indicate the prognostic accuracy of the nomogram over 3 and 5 years. We used AUCs at 3 and 5 years to assess prognostic accuracy. (d) Kaplan-Meier survival analysis was conducted for LGG patients in the high-risk-score group and the low-risk-score group from a Cox regression model, and p values were calculated using the log-rank test.

3.5. WGCNA of RNA-Seq in LGG

To identify the gene network that was highly associated with GAS5 expression, we further analysed LGG RNA-Seq data by WGCNA. WGCNA has been widely applied to study coexpressed gene networks by scale-free weighted network construction and module detection. In our analysis, 43 modules were detected; the clustering of genes is shown as a dendrogram (Figures 7(a) and S3). The module similarity was quantified and is shown as a heatmap below (Figure 7(b)). The weighted network of identified genes from RNA-Seq is shown as a heatmap below, which depicts the topological overlap matrix among all genes (Figure 7(c)). Then, 26 modules were identified to be related with GAS5 expression by the heatmap of associations between modules and GAS5 expression (Figure 7(d)).

(a)

(b)

(c)

(d)

Figure 7: WGCNA of RNA-Seq data from LGG patients and GBM patients. (a) Dendrogram of modules identified by WGCNA. (b) Heatmap of module-GAS5 expression associations. (c) Topological overlap matrix among detected probes from RNA-Seq. (d) Heatmap for module similarity. (e-f) Scatterplot of gene significance for GAS5 expression versus module membership in the salmon (e) and blue modules (f).

3.6. Enrichment Analysis of Genes Highly Associated with GAS5 Expression in LGG

Based on the results above, the salmon and blue modules, which had the most significant relationship with GAS5 expression, were selected. The scatterplot below illustrates genetic significance between GAS5 expression (the correlation of genes with GAS5 expression) and module membership (the correlation of genes with clusters) for the salmon module (Figure 8(a)) and the blue module (Figure 8(b)). The results indicated that genes that had high correlations with the salmon and blue modules were also strongly associated with GAS5 expression. The genes from the salmon and blue modules with coefficients higher than 0.5 were selected and subjected to GO enrichment analysis (Figures 8(c) and 9). rRNA processing, ribonucleoprotein complex biogenesis, translation factor activity, ribosomal subunit biogenesis, and other terms were enriched in the analysis, indicating the key roles of GAS5 in the regulation of ribosomes in LGG tissues (Figures 9(a), 9(c), and 9(e)). The genes involved in GO analysis and their fold changes between the high GAS5 expression group (the 50 samples with the highest GAS5 expression) and the low GAS5 expression group (the 50 samples with the lowest GAS5 expression) are shown on the circular plots (Figures 9(b), 9(d), and 9(f)). RPL18A, RPL23, RPL24, PRL32, PRS12, and other ribosome proteins were highly expressed along with high GAS5 expression, resulting in the upregulation of rRNA processing and ribosome biogenesis. Finally, the network connections in the salmon and blue modules with weight parameters higher than 0.19 from WGCNA were loaded into Cytoscape and analysed. A variety of ribosome proteins, such as RPL10, RPL12, RPL17, and RPS24, were included in the network (Figures 10(a) and 10(b)). High interconnectivity within the gene cluster is depicted by thick edges. The genes most strongly related to GAS5 were analysed in the TCGA dataset and validated with a CGGA (Chinese Glioma Genome Atlas) dataset (Figures 10(c) and 10(d)). The most related genes included RPL6, RPL37A, RPL37, RPS24, and C20orf199, which were included in the salmon module. The results further confirmed the relationship between salmon module and GAS5.

(a)

(b)

(c)

Figure 8: The relationship between gene significance and module membership for the salmon and blue modules. (a-b) Scatterplot of gene significance for GAS5 expression versus module membership for the salmon (a) and blue modules (b). Red indicates a negative correlation with GAS5, and blue indicates a positive correlation with GAS5. (c) The network of gene ontology analysis for biological processes, cellular components, and molecular functions based on genes in the salmon and blue modules.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 9: Enrichment analysis of the genes in the salmon and blue modules. (a-b) Gene ontology analysis for biological processes and the gene regulation involved in biological processes. (c-d) Gene ontology analysis for cellular components and the gene regulation involved in cellular components. (e-f) Gene ontology analysis for molecular functions and the gene regulation involved in molecular functions.

(a)

(b)

(c)

(d)

Figure 10: Visualization of genes in the salmon and blue modules with weights higher than the threshold (weight > 0.19) and correlations between GAS5 and the most closely related molecules. Node sizes indicate betweenness centrality, which reflects the amount of control that this node exerts over the interactions of other nodes in the network. Edge sizes depict the weight of the connection between genes. (a) Genes in the salmon module with weights higher than 0.19. (b) Genes in the blue module with weights higher than 0.19. (c) The correlations between GAS5 and the most closely related genes in the TCGA dataset. (d) The validation of the association in the CGGA dataset.

3.7. RNA-Seq Analysis of the High and Low GAS5 Expression Groups

To validate the results of WGCNA, we performed GO analysis and GSEA in the high and low GAS5 expression groups. The workflow is shown in Figure 11(a). The volcano plot indicates the fold changes and Q values of differentially expressed genes (Figure 11(b)). A heatmap of the 200 differentially expressed genes with the highest significance levels between the high and low GAS5 expression groups is shown in Figure 11(c). Then, differentially expressed genes were analysed by GO analysis (Figure 11(d)). The results indicated a high degree of association between GAS5 expression and translation initiation. GSEA was performed using all genes from GAS5 high-expressed group and GAS5 low-expressed group, indicating translation initiation, ribosomal biogenesis, and ribonucleoprotein complex biogenesis were significantly upregulated (Figures 11(e), 11(f), and 11(g)). An MA plot was generated to show the differentially expressed genes with high significance levels between the high and low GAS5 expression groups (Figure S4).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 11: RNA-Seq analysis of LGG tissues with the highest and lowest GAS5 expression. (a) Workflow for RNA-Seq and downstream analysis based on differentially expressed genes and all genes in RNA-Seq. (b) Volcano plot for differentially expressed genes in the groups with high and low expression of GAS5. (c) Heatmap of the 200 differentially expressed genes with the highest significance levels between the groups with high and low expression of GAS5. (d) Heatmap of gene ontology analysis for biological processes. (e-g) Gene set enrichment analysis for the groups with high and low expression of GAS5.

4. Discussion

LncRNAs are non-protein-coding transcripts longer than 200 nucleotides and have been demonstrated to play critical roles in diverse cellular processes and the regulation of biological functions [7, 27, 28]. LncRNAs regulate gene expression at the transcriptional, posttranscriptional, and epigenetic levels. The expression levels of lncRNAs are frequently deregulated in cancer [29–32]. For instance, GAS5 has been reported as a crucial tumour suppressor in a variety of human cancers. GAS5 expression is significantly decreased in prostate cancer cells compared with prostate epithelial cells [33]. Ectopic expression of GAS5 induces cell-cycle arrest in the G0–G1 phase by increasing the activity of the P promoter. In addition, GAS5 interacts with E2F1 and enhances the binding of E2F1 to the P promoter in prostate cancer [33]. GAS5 expression is also markedly downregulated in gastric cancer tissues and influences gastric cancer cell proliferation via regulating the expression of E2F1 and P21 [11]. One study revealed that GAS5 was significantly downregulated in lung adenocarcinoma tissues compared with paired adjacent nontumorous tissue samples. The overexpression of GAS5 varied inversely with the expression of EGFR pathway and IGF-1R proteins in lung adenocarcinoma tissues [12].

Recent studies have implied that lncRNAs participate in the development and malignancy of GBM by regulating glioma stem cell self-renewal, proliferation, differentiation, therapeutic response, and resistance [34–37]. In this study, we focused on GAS5, one crucial noncoding gene in glioma progression [37–39]. GAS5 plays a critical role in the control of mammalian apoptosis and acts as a competitive glucocorticoid response element to mediate cell population growth by suppressing glucocorticoid receptors [9]. GAS5 inhibits the association of glucocorticoid receptors with their DNA recognition sequences by binding to their DNA-binding domain through its double-stranded RNA glucocorticoid receptor element mimic [9]. Thus, GAS5 is upregulated in growth-arrested cells and sensitizes mammalian cells to apoptosis [40]. In this study, we focused on the mRNAs coexpressed with GAS5 in gliomas and how GAS5 controlled the survival of glioma patients. By characterizing the levels of GAS5 expression based on data from TCGA-LGG and TCGA-GBM patients, we demonstrated a significant downregulation of GAS5 expression in LGG and GBM compared with normal brain tissue. One previous study revealed that GAS5 expression decreased as the histologic grade of glioma increased [38]. However, by analysing TCGA-LGG and GBM RNA-Seq datasets, we arrived at the conclusion that histologic grade did not influence GAS5 expression. The discrepancy may be due to the sample numbers. In this study, we analysed the expression levels of GAS5 in 702 samples, far more than Zhao’s study used. One previous report revealed that hypermethylation of CpG islands in promoter regions contributes to the downregulation of GAS5 expression [10]. The methylation level of the GAS5 gene was detected by the six probes targeting the CpG island before the transcription start site of GAS5. Both LGG and GBM tissues had low GAS5 methylation levels. However, there was no significant association between the methylation level and the expression level of GAS5. Amplification-type alteration was associated with increased GAS5 expression in both LGG and GBM tissues. GAS5 expression played different roles in affecting the OS of LGG and GBM patients. GAS5 expression was strongly associated with OS in LGG patients. We also examined GAS5 expression in LGG patients with different molecular characteristics and histologic grades. GAS5 expression was higher in IDH1 mutant patients than in patients with wild-type IDH1. To confirm the association between GAS5 and OS in LGG patients, we performed Kaplan-Meier survival analysis on LGG patients according to their GAS5 expression levels, stratified by clinicopathological characteristics. Kaplan-Meier survival analysis showed high association between GAS5 expression level and OS in both grade II and grade III gliomas. We also performed multivariate survival analysis of GAS5 with clinical and molecular characteristics in LGG patients. The Cox HR model indicated that GAS5 was strongly associated with the overall survival of patients. Based on the results discussed above, a nomogram was developed to predict the 3-year and 5-year survival probabilities of LGG patients. Time-dependent ROCs and log-rank tests further confirmed the results.

To understand the molecular mechanism behind the effect of GAS5 in LGG tissues, we performed WGCNA and downstream enrichment analysis. The results indicated that the genes in the salmon and blue modules played important roles in ribosome function, which was consistent with previous reports on the involvement of GAS5 in ribosomal RNA biosynthesis [41]. In addition, the genes highly correlated with GAS5 expression showed intensive connection when visualized by Cytoscape. A variety of the hub genes were ribosomal proteins, e.g., RPL13, RPL37, PRL38, and PRL10A, indicating the strong association between GAS5 and ribosomal function. The results were consistent with a previous study, which revealed that GAS5 may act as a “ribo-repressor” of glucocorticoid receptors to influence cell survival and metabolic activity during starvation by modulating the transcriptional activity of the glucocorticoid receptor [40]. To further confirm the results from WGCNA, we also performed RNA-Seq analysis and downstream GO analysis and GSEA. Translation initiation was the most predominant result from GO, implying that GAS5 may be an important regulator of transcription initiation. One study demonstrated that GAS5 was enriched with eIF4E, which had two RNA binding motifs for GAS5 [42]. The deletion of either motif inhibited the binding of GAS5 to eIF4E. In addition, GSEA results also indicated upregulation of translation initiation and ribosomal biogenesis in group with high GAS5 expression, which was consistent with GO analysis and previous research [43].

Altogether, the data presented here suggest that GAS5 plays essential roles in the physiological process of LGG. The mechanism may include ribosome biogenesis and regulation of translation initiation. In conclusion, based on molecular and clinical analysis, our study demonstrated that GAS5 is linked tightly to ribosomal biogenesis and translation initiation in LGG tissues and acts as a predictor of survival in patients with early-stage gliomas.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

This article does not contain any studies with human participants performed by any of the authors.

Conflicts of Interest

We declare that we have no financial or personal relationships with other people or organizations that could inappropriately influence our work.

Authors’ Contributions

Xuanwen Bao and Run Shi conceived and designed the experiments. Xuanwen Bao, Shan Xin, and Kai Zhang analysed the data. Xuanwen Bao and Yanfang Wang wrote the paper. All authors read and approved the final manuscript. Yanfang Wang and Shan Xin contributed equally to this work.

Acknowledgments

We would like to thank Dr. Omid Azimzadeh, Dr. Natasa Anastasov, and Dr. Michael Rosemann for helpful discussions and suggestions.

Supplementary Materials

Figure S1: the GAS5 expression in different stages of glioma (a) and the correlation of GAS5 methylation level and GAS5 expression level in GBM (b) and LGG (c) tissues. Figure S2: the associations between GAS5 expression level and tumour purity in LGG (a) and GBM (b) patients. Figure S3: the dendrogram of modules identified by WGCNA. Figure S4: MA plot of differentially expressed genes with high significance levels in the groups with high and low GAS5 expression groups in LGG tissues. (Supplementary Materials)

References

Q. T. Ostrom, H. Gittleman, P. Farah et al., “CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2006–2010,” Neuro-Oncology, vol. 15, supplement 2, pp. ii1–ii56, 2013.

Network CGAR, “Comprehensive, Integrative Genomic Analysis of Diffuse Lower-Grade Gliomas,” The New England Journal of Medicine, vol. 372, no. 26, pp. 2481–2498, 2015.

M. J. Van Den Bent, “Practice changing mature results of RTOG study 9802: Another positive PCV trial makes adjuvant chemotherapy part of standard of care in low-grade glioma,” Neuro-Oncology, vol. 16, no. 12, pp. 1570–1574, 2014.

A. S. Jakola, K. S. Myrmel, R. Kloster et al., “Comparison of a strategy favoring early surgical resection vs a strategy favoring watchful waiting in low-grade gliomas,” Journal of the American Medical Association, vol. 308, no. 18, pp. 1881–1888, 2012.

E. Bouffet, V. Larouche, B. B. Campbell et al., “Immune checkpoint inhibition for hypermutant glioblastoma multiforme resulting from germline biallelic mismatch repair deficiency,” Journal of Clinical Oncology, vol. 34, no. 19, pp. 2206–2211, 2016.

P. Han and C.-P. Chang, “Long non-coding RNA and chromatin remodeling,” RNA Biology, vol. 12, no. 10, pp. 1094–1098, 2015.

J. T. Lee, “Epigenetic regulation by long noncoding RNAs,” Science, vol. 338, no. 6113, pp. 1435–1439, 2012.

S. Carpenter, D. Aiello, M. K. Atianand et al., “A long noncoding RNA mediates both activation and repression of immune response genes,” Science, vol. 341, no. 6147, pp. 789–792, 2013.

T. Kino, D. E. Hurt, T. Ichijo, N. Nader, and G. P. Chrousos, “Noncoding RNA Gas5 is a growth arrest- and starvation-associated repressor of the glucocorticoid receptor,” Science Signaling, vol. 3, no. 107, p. ra8, 2010.

X. Shi, M. Sun, H. Liu et al., “A critical role for the long non-coding RNA GAS5 in proliferation and apoptosis in non-small-cell lung cancer,” Molecular Carcinogenesis, vol. 54, pp. E1–E12, 2015.

M. Sun, F.-Y. Jin, R. Xia et al., “Decreased expression of long noncoding RNA GAS5 indicates a poor prognosis and promotes cell proliferation in gastric cancer,” BMC Cancer, vol. 14, no. 1, article 319, 2014.

S. Dong, X. Qu, W. Li et al., “The long non-coding RNA, GAS5, enhances gefitinib-induced cell death in innate EGFR tyrosine kinase inhibitor-resistant lung adenocarcinoma cells with wide-type EGFR via downregulation of the IGF-1R expression,” Journal of Hematology & Oncology, vol. 8, no. 1, Article ID 43, 2015.

Z.-Q. Tu, R.-J. Li, J.-Z. Mei, and X.-H. Li, “Down-regulation of long non-coding RNA GAS5 is associated with the prognosis of hepatocellular carcinoma,” International Journal of Clinical and Experimental Pathology, vol. 7, no. 7, pp. 4303–4309, 2014.

J. Shen, T. R. Hodges, R. Song et al., “Serum HOTAIR and GAS5 levels as predictors of survival in patients with glioblastoma,” Molecular Carcinogenesis, vol. 57, no. 1, pp. 137–141, 2018.

X.-Q. Zhang, S. Sun, K.-F. Lam et al., “A long non-coding RNA signature in glioblastoma multiforme predicts survival,” Neurobiology of Disease, vol. 58, pp. 123–131, 2013.

M. S. Cline, B. Craft, T. Swatloski et al., “Exploring TCGA pan-cancer data at the UCSC cancer genomics browser,” Scientific Reports, vol. 3, Article ID 2652, 2013.

G. Consortium, “The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans,” Science, vol. 348, pp. 648–660, 2015.

L. Wilkinson, ggplot2: Elegant Graphics for Data Analysis, Springer, 2016.

P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, article 559, 2008.

S. Horvath, B. Zhang, M. Carlson et al., “Analysis of oncogenic signaling networks in glioblastoma identifies ASPM as a molecular target,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 103, no. 46, pp. 17402–17407, 2006.

P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software Environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.

G. Yu, L.-G. Wang, Y. Han, and Q.-Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.

S. Tripathi, M. O. Pohl, Y. Zhou et al., “Meta-and orthogonal integration of influenza, OMICs , data defines a role for UBR4 in virus budding,” Cell Host & Microbe, vol. 18, pp. 723–735, 2015.

A. Subramanian, P. Tamayo, V. K. Mootha et al., “Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 102, no. 43, pp. 15545–15550, 2005.

S. V. Vasaikar, P. Straub, J. Wang, and B. Zhang, “LinkedOmics: Analyzing multi-omics data within and across 32 cancer types,” Nucleic Acids Research, vol. 46, no. 1, pp. D956–D963, 2018.

G. K. Smyth, Limma: Linear Models for Microarray Data. Bioinformatics and Computational Biology Solutions Using R and Bioconductor, Springer, 2005.

A. Fatica and I. Bozzoni, “Long non-coding RNAs: new players in cell differentiation and development,” Nature Reviews Genetics, vol. 15, no. 1, pp. 7–21, 2014.

P. J. Batista and H. Y. Chang, “Long noncoding RNAs: cellular address codes in development and disease,” Cell, vol. 152, no. 6, pp. 1298–1307, 2013.

A. M. Schmitt and H. Y. Chang, “Long noncoding RNAs in cancer pathways,” Cancer Cell, vol. 29, no. 4, pp. 452–463, 2016.

A. Lin, C. Li, Z. Xing et al., “The LINK-A lncRNA activates normoxic HIF1α signalling in triple-negative breast cancer,” Nature Cell Biology, vol. 18, no. 2, pp. 213–224, 2016.

W. Nie, H.-J. Ge, X.-Q. Yang et al., “LncRNA-UCA1 exerts oncogenic functions in non-small cell lung cancer by targeting miR-193a-3p,” Cancer Letters, vol. 371, no. 1, pp. 99–106, 2016.

A. Bhan and S. S. Mandal, “LncRNA HOTAIR: A master regulator of chromatin dynamics and cancer,” Biochimica et Biophysica Acta (BBA) - Reviews on Cancer, vol. 1856, no. 1, pp. 151–164, 2015.

G. Luo, D. Liu, C. Huang et al., “LncRNA GAS5 inhibits cellular proliferation by targeting P27Kip1,” Molecular Cancer Research, vol. 15, no. 7, pp. 789–799, 2017.

E.-B. Bian, J. Li, Y.-S. Xie, G. Zong, J. Li, and B. Zhao, “LncRNAs: New players in gliomas, with special emphasis on the interaction of lncRNAs with EZH2,” Journal of Cellular Physiology, vol. 230, no. 3, pp. 496–503, 2015.

P. Wang, Z. Ren, and P. Sun, “Overexpression of the long non-coding RNA MEG3 impairs in vitro glioma cell proliferation,” Journal of Cellular Biochemistry, vol. 113, no. 6, pp. 1868–1874, 2012.

L. Han, K. Zhang, Z. Shi et al., “LncRNA profile of glioblastoma reveals the potential role of lncRNAs in contributing to glioblastoma pathogenesis,” International Journal of Oncology, vol. 40, no. 6, pp. 2004–2012, 2012.

X. Zhao, Y. Liu, J. Zheng et al., “GAS5 suppresses malignancy of human glioma stem cells via a miR-196a-5p/FOXO1 feedback loop,” Biochimica et Biophysica Acta (BBA) - Molecular Cell Research, vol. 1864, no. 10, pp. 1605–1617, 2017.

X. Zhao, P. Wang, J. Liu et al., “Gas5 Exerts Tumor-suppressive Functions in Human Glioma Cells by Targeting miR-222,” Molecular Therapy, vol. 23, no. 12, pp. 1899–1911, 2015.

Q. Liu, S. Sun, W. Yu et al., “Altered expression of long non-coding RNAs during genotoxic stress-induced cell death in human glioma cells,” Journal of Neuro-Oncology, vol. 122, no. 2, pp. 283–292, 2015.

M. Mourtada-Maarabouni, M. R. Pickard, V. L. Hedge, F. Farzaneh, and G. T. Williams, “GAS5, a non-protein-coding RNA, controls apoptosis and is downregulated in breast cancer,” Oncogene, vol. 28, no. 2, pp. 195–208, 2009.

C. M. Smith and J. A. Steitz, “Classification of gas5 as a multi-small-nucleolar-RNA (snoRNA) host gene and a member of the 5'-terminal oligopyrimidine gene family reveals common features of snoRNA host genes,” Molecular and Cellular Biology, vol. 18, no. 12, pp. 6897–6909, 1998.

G. Hu, Z. Lou, and M. Gupta, “The long non-coding RNA GAS5 cooperates with the eukaryotic translation initiation factor 4E to regulate c-Myc translation,” PLoS ONE, vol. 9, no. 9, Article ID e107016, 2014.

W. Sun, Y. Yang, C. Xu, and J. Guo, “Regulatory mechanisms of long noncoding RNAs on gene expression in cancers,” Cancer Genetics, vol. 216-217, pp. 105–110, 2017.