Research Paper Advance Articles
RP11-552D4.1: a novel m6a-related LncRNA associated with immune status in glioblastoma
- 1 Department of Neurosurgery, Shanghai Pudong New Area People’s Hospital, Shanghai, China
- 2 Key Molecular Lab, Shanghai Pudong New Area People’s Hospital, Shanghai, China
Received: January 22, 2022 Accepted: July 5, 2022 Published: July 18, 2022https://doi.org/10.18632/aging.204177
How to Cite
Copyright: © 2022 Zheng et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Glioblastoma (GBM) is the most malignant form of brain cancer in the world. Nevertheless, the survival rate of patients with GBM is extremely low. N6-methyladenosine (m6A) and long noncoding RNAs (lncRNAs) conduct important biological functions in patients’ survival status and the immunotherapeutic response. Here, m6A-related lncRNAs were identified by a co-expression method. Univariate and multivariate Cox regression together with LASSO were applied to establish the risk model. Kaplan-Meier and ROC analysis were applied to evaluate the prediction power of this risk model. Finally, the related immune profiling and chemical sensitivity targets were also investigated. The risk model holding three m6A-related lncRNAs was confirmed as an independent predictor for the prognosis. Furthermore, we found the risk model based on m6A-related lncRNAs is associated with the immune status, immunosuppressive biomarkers, and chemo-sensitivity in GBM patients. The RP11-552D4.1 is found to facilitate neuronal proliferation. This risk model consisted of m6A-related lncRNAs may be available for the clinical interventions in GBM patients.
Glioma is one of the common malignant brain tumors, which accounts for 30% of adult brain cancer, and its incidence rate is rising every year [1, 2]. The most malignant pathological type is glioblastoma (GBM), which accounts for about 70–80% of adult glioma . Compared with low-grade glioma (LGG) and astrocytoma, GBM demonstrates a very poor outcome and is more likely to metastasize . When GBM appears to metastasize, its outcome is less favorable . The curative treatment of early GBM is total resection. Nevertheless, one-third of patients are under a recurrence after brain surgery , even treated together with gene therapy and immune therapy, and the overall survival of GBM patients is still very shorter .
Several lncRNAs in the brain are specifically related to neuronal and synaptic function . LncRNAs constitute a large class of post-transcriptional regulators, some of which can act as ceRNAs to inhibit microRNAs (miRNAs) in the brain [8, 9]. For instance, lnc -NKILA acts as a miR-195 sponge, leading to increased expression of miR-195 target gene NLRX1, particularly in neocortical and hippocampal neurons .
Epigenetic modification in GBM plays a key role in tumor initiation and progression [11, 12]. M6A is a dominant RNA modification on both coding and non-coding RNAs. A series of studies have demonstrated that m6A methylation has a critical effect in the tumor initiation and progression [13–15]. M6A-related genes can be divided into writers, erasers and readers [16, 17]. We previously reported that lncRNA-ATB may participate in tumorigenesis through targeting NF-κB and P38/MAPK ; On the other hand, lncRNA-ATB could also activate the astrocyte via miR-204-3p  and miR-200a  to promote the progression of glioma. However, lncRNA MEG3 can be regulated by DNA methylation and further inhibit glioma growth . The double-faced role of lncRNAs in glioma might be due to its m6A modification style. Few studies have investigated the combination of m6A related lncRNAs in cancers, especially in GBM. Therefore, we constructed a GBM risk model to study the prognostic utility of m6A related genes in GBM.
Materials and Methods
GBM data was available from the Cancer Genome Atlas (TCGA) on Oct 01, 2020, which consists of clinical data (survival status, age, sex, race), mutation files of the SNP, and gene expression profiles. Clinical information was available on 169 TCGA GBM samples and five control samples. Chinese Glioma Genome Atlas (CGGA) data are used for verification of the interested genes. Clinical information of GBM patients was downloaded from the TCGA website as well. Only the clinical information of patients with confirmed GBM pathology was included for final bioinformatic analysis. The candidate m6A gene sets were searched from “c2.cp.kegg. v7.0. symbols” of Gene Set Enrichment Analysis (GSEA). The correlation between lncRNAs and m6A genes was analyzed with co-expression methods. The inclusion criteria of m6A-related lncRNAs were based on Pearson coefficient more than 0.4 and p-value less than 0.001. The “Limma” R package was used to identify the m6A-related lncRNAs, and logFC >1 and false discovery rate (FDR) <0.05 was used to identify differently- expressed m6A-related lncRNA (DEmlncRNA). The flow chart for the current study is briefly summarized in Figure 1.
Figure 1. Flow chart of the study. The RNA-seq data for GBM were obtained from TCGA and m6A-related immune lncRNAs were intersected among immune genes, 19 m6A genes, and 494 differently-expressed lncRNAs. Adding the clinical information, univariate and multivariate COX regression and LASSO were applied to identify one m6A-regulated ir-lncRNA which is correlated with the survival status of GBM patients. Then, the ir-lncRNAs were analyzed for immune status, immune checkpoint inhibitors, and chemical sensitivity for GBM.
Survival status correlation
We only included patients with GBM holding effective survival data from 2 days to 3881 days. Both univariate and multiple Cox model was used to assess the relationship between m6A-related lncRNAs and the survival status of GBM patient (P < 0.05). Risk score = Coefficient of mRNA × lncRNAs. The least absolute shrinkage and selection operator (LASSO) method is further applied to identify the exact coefficient of each DEmlncRNA based on the survival status. R packages “survival” and “survminer” were used to compare the survival status between the high-and low-risk group.
Prognostic signature construction
The survivalROC package was used to draw the Receiver Operating Characteristic (ROC) curve, in order to evaluate the prediction ability of the risk model. By integrating clinical features into Cox models, a prognostic signature was identified as dependent or independent factors. In addition, we drew the one-, two- and three-year ROC curves and obtained the highest area under curve (AUC) value.
Tumor-infiltrating cells (TILs) in GBM
We next investigated the association between risk group and TILs, with different immune infiltration analysis algorithm. Wilcoxon signed-rank test was used to compare the number of TILs between high- and low-risk groups. The correlation between riskScores and TILs was investigated by the Spearman method. A lollipop plot shows the correlation coefficient.
Expression of immune-checkpoint inhibitor (ICI) markers
The correlation between riskScores and the ICI markers was investigated by Wilcoxon signed-rank test and further visualized with ggstatsplot package.
The association between the risk model and the clinical therapy
In order to evaluate the clinical value of the risk model in GBM, the half maximal inhibitory concentration (IC50) of commonly used chemotherapy drugs for GBM patients was investigated. Most medicines were evaluated with pRRophetic package. Wilcoxon signed-rank test was carried out to compare the different IC50 between the high- and low-risk scores.
Reagents and antibodies
Mouse monoclonal to MAP2 (ab11267, Abcam) and rabbit monoclonal to β-tubulin (ab201831, Abcam) were utilized to evaluate the cell proliferation for IFC.
Primary cortical neuronal cultures and plasmids treated IFC
Primary culture of cortical neurons was previously reported . Plasmid-mediated lncRNA overexpression and knockdown vector were obtained from Yazai biotech (Shanghai, China), siRNAs targeting RP11-552D4.1 were bought from Thermo Fisher Scientific (n505669, Shanghai, China) and control plasmid was ordered from GeneCopoeia (Rockville, MD, USA). The immunofluorescence (IFC) method was previously reported as well. In this study, we probed the cultured neurons with MAP2 and β-tubulin overnight and the secondary antibody was HRP.
Each group of cells was adjusted at 1,000 cells per well. 10 mL of CCK-8 solution (Beyotime Biotechnology, Shanghai, China) was added to the cell dish after 24 hours, and a blank control has only CCK-8 solution. Absorbance (OD) value of each well was read for each well at 450 nm and tested every 24 hours for three days.
Differential expression analysis of m6A-related lncRNAs
The gene transfer format (GTF) files from Ensemble were used to differentiate the noncoding and coding RNAs, and m6A-related lncRNAs were identified based on a co-expression form. In total, 460 m6A-related lncRNAs were identified, and among them, 152 were DEmlncRNAs. Among which, 127 were up-regulated and 25 were down-regulated (Figure 2A, 2B).
Figure 2. Establishment of a risk assessment model based on DEmlncRNA differentially expressed m6A-related lncRNAs (DEmlncRNAs). The heatmap (A) and a volcano plot (B) were displayed. Establishment of the LASSO regression (C and D). 3 DEmlncRNA were shown by a forest map with univariate and multivariate Cox analysis (E and F).
DEmlncRNA screening and risk model construction
152 DEmlncRNAs were further identified and LASSO analysis was applied to find differentially-expressed lncRNAs. Figure 2C–2F shows the DElncRNAs. At the same time, the maximum AUC value was at 0.679, and the optimal DEmlncRNAs were identified (Figure 3A). We also calculated the Akaike information criterion (AIC) to look for the maximum inflection point as the cut-off point in the 3-year ROC curve (Figure 3B). ROC curves for one, two, and three years were drawn respectively, and all AUC values were more than 0.65 (Figure 3C). In addition, the one-year AUC value was much better against common clinical features. (Figure 3D).
Figure 3. Establishment of a risk assessment model based on DEmlncRNA. The curve of each AUC value generated by ROCs of DEmlncRNAs was drawn, and the highest point of AUC was determined; The maximum inflection point is the cut-off point acquired from AIC (A and B). The 1-year, 2-year, and 3-year ROC of the optimal model showed that all AUC values exceeded 0.65 (C). Compared with other common clinical features, the 1-year ROC curves showed the superiority of the risk score against age and gender (D).
Risk model prediction in the clinical features
Seventy cases and 89 cases were classified into high- and low-risk groups based on the cut-off value. Figure 4A and 4B show the riskScores and survival rate of each patient. We found the outcome in patients from the high-risk group was poorer than that from the low-risk group, shown by KM analysis (p = 0.003) (Figure 4C). Then, a Chi-square test was applied to investigate the relationship between the risk model in GBM and clinical features. The scatter plots examined through Wilcoxon signed-rank test displayed age and gender (Figure 4D, 4E) were not significantly correlated with the risk model. Consequently, univariate and multiple Cox model showed that riskScores is an independent predicting factor (p < 0.001, HR = 1.927, 95% CI [1.395–2.660]) (Figure 5A), and riskScores (p < 0.001, HR = 1.952, 95% CI [1.395–2.733]) (Figure 5B). The m6A-lncRNA coexpression network was visualized using the Sankey diagram in Figure 5C, and 152 m6A-related lncRNAs were discerned as m6A-related lncRNAs.
Figure 4. Prognostic prediction of the risk assessment model. Risk scores (A) and survival outcome (B) of each case. Kaplan-Meier survival curve of high-risk group and low-risk group (C). The scatter plot showed that age (D) and sex (E) were not significantly correlated with a risk score.
The risk model is correlated with immune status in GBM
As lncRNAs play a role in the immune infiltration of GBM, we next explored whether this model correlates with the immune infiltration status. Spearman correlation was performed, and the lollipop diagram showed their relationships (Figure 6A). We found that a positive correlation between the high-risk group and tumor-infiltrating immune cells, such as neutrophils and monocytes (Figure 6D, 6E), was negatively correlated with CD8+T cells and CD4+ cells (Figure 6B, 6C). As ICIs have an important role in treating solid tumors, however, their roles in GBM remain elusive, we further investigated the association between the risk model and ICI markers. Our results demonstrated that high-risk group was correlated with the expression level of HAVCR2 (p < 0.001) and CTLA4 (p < 0.001) (Figure 7A, 7B). There was no significant difference between the risk scores and ICI related-genes, such as LAG3, CD276, PDCD1, and TNFRSF18 (All p > 0.05, Figure 7C–7F).
Figure 6. Estimate tumor-infiltrating cells through risk assessment model. (A) Association between the risk scores and TILs. High-risk group has higher neutrophils and monocytes and a lower proportion of CD4+ and CD8+ T cells (B–E).
Figure 7. Estimate immunosuppressive molecules through risk assessment model. High-risk scores were associated with the expression of HAVCR2 (A) and CTLA4 (B). There was no statistical difference between the risk scores and immune-related genes, such as PDCD1 (C), LAG3 (D), CD274 (E) and TNFRSF18 (F).
Correlational assessment of risk model and chemotherapy drugs
Besides the relationship with ICIs, we further explored the relationship between the risk and efficacy of commonly used chemotherapy drugs for GBM. The results showed that the high-risk group had a higher half inhibitory concentration (IC50) of chemotherapy drugs for metformin (p < 0.05), while a lower IC50 for AG.014699, PARP inhibitor (p < 0.001), indicating that this risk model could be used as a potential predictor for the chemotherapy sensitivity (Figure 8A–8D).
Validation of the independent risk gene RP11-552D4.1 in CGGA
Because the expression level of RP11-552D4.1 is an independent factor in the survival prediction for GBM (Figure 2F, with a p-value = 0.038), we next investigated its expression in Chinese Genome Atlas (CGGA), we first explored the expression of RP11-552D4.1 in different subtypes of CGGA and found it was lower in the Proneural group. The AUC of ROC in predicting RP11-552D4.1 for CGGA subtype (Proneural) is 0.637. In addition, the expression of RP11-552D4.1 is lower in IDH-1 mutant type compared to IDH-1 wildtype groups. As for the survival status, we found the expression of RP11-552D4.1 is positively associated with the survival status in WHO IV, while negatively correlated with the OS in WHO II, although it did not reach the statistical difference (p = 0.4 and 0.072, respectively). This might be due to the smaller sample size in CGGA compared to the TCGA dataset Supplementary Figure 1. Consistently, in the TCGA dataset, GBM patients with higher expression of RP11-552D4.1 has a longer survival period, while LGG patients with higher expression of this gene shows a shorter survival period (Supplementary Figures 2 and 3).
RP11-552D4.1 facilitates neuronal survival
From both TCGA and CGGA dataset, a similar trend between the expression of RP11-552D4.1 and survival outcome was identified, we next investigated the effect of RP11-552D4.1 on primary cultured neurons. Si- RP11-552D4.1 was used to knock down the expression of RP11-552D4.1 and si-NC was used as a control. CCK-8 assays indicated that primary cultured neurons transfected with si-RP11-552D4.1 showed decreased cell proliferation (Figure 9A and 9B). Similarly, transfection with si- RP11-552D4.1 also significantly reduced the expression of MAP2 and β-tubulin in cell culture (*P < 0.05, Figure 9C–9E). In contrast, overexpressed RP11-552D4.1 could upregulate the CCK level of primary neurons and increase the expression of MAP2 and β-tubulin as well (Figure 10).
Figure 9. RP11-552D4.1 facilitates neuronal proliferation. (A) The expression of RP11-552D4.1 in si-RP11 was verified. (B) After primary cortical neurons were transfected with si-NC or si-RP11-552D4.1, a CCK-8 assay was applied to assess the cell proliferation. (C) Effect of si-RP11-552D4.1 on neuronal proliferation (Green: MAP2, Red: β-tubulin, Merged with DAPI). (D, E) Quantification of the mean density of MAP2 and β-tubulin. The scale bar is 100 μm. *P < 0.05 (compared to the si-NC group).
Figure 10. Si-RP11-552D4.1 inhibits neuronal proliferation. (A) Expression of RP11-552D4.1 in overexpressed RP11-552D4.1 was verified. (B) Primary cortical neurons were transfected with NC or RP11-552D4.1 plasmid. (C) Effect of RP11-552D4.1 on neuronal proliferation (Green: MAP2, Red: β-tubulin, Merged with DAPI). (D, E) Quantification of the mean density of MAP2 and β-tubulin. The scale bar is 100 μm. *P < 0.05 (compared to the NC group).
M6A-related lncRNA is a recent hot topic in the neuroscience field. Current studies have investigated the correlation between the m6A-related lncRNA and the risk of cancer. However, its role in glioma, even GBM is missing. So, whether the expression level of m6A-related lncRNAs has a relationship with the outcome of GBM patients’ needs to be addressed. The lncRNAs included in this model are related to m6A-genes based on co-expressed methods, so these lncRNAs may be regulated by m6A genes.
A recent review paper summarized m6A RNA methylation in brain tumor and reported that several m6A molecules differentially expressed in glioma from the CGGA dataset . However, most m6A molecules directly regulated by miRNAs or indirectly adjusted by lncRNAs or circRNAs with a competitive endogenous RNA mechanism. There were very few reports regarding the lncRNAs directly modified by m6A methylation. Meanwhile, another bioinformatical study focusing on the m6A regulators in GBM identified that the risk score consists of HNRNPC, ALKBH5 and FTO could predict the survival outcome in GBM patients independently . Again, this study did not further investigate the target genes of these m6A regulators. Therefore, we identified the prognosis-related DEmlncRNAs to construct a reliable model to predict the OS in GBM patients. Univariate and multivariate Cox models demonstrated that this risk model could predict the prognostic outcome of GBM patients, and it is independent of other clinical indices, such as patients’ age and gender. The risk score obtained from the model was more accurate against other clinical features in predicting survival status based on ROC analysis.
In addition, the chemotherapy sensitivity for GBM was investigated. The difference between risk groups and immune status and ICI genes were also investigated, which suggested this risk model was able to classify GBM patients who might be suitable for immunotherapy.
The response of ICI blockers is associated with TILs [25–27]. Our study applied customized immunol gene sets to study the correlation between risk factors and TILs. Our results demonstrated that DEmlncRNAs show a positive relationship with tumor-infiltrating immune cells such as neutrophils and monocytes, while negatively correlated with CD4+ and CD8+T cells. In addition, our study indicates that the high-risk group in GBM was associated with chemotherapy sensitivity such as PARP inhibitor, instead of metformin, which provides potential targets for immunotherapy in GBM.
RP11-552D4.1 is previously found to decrease in heart failure , which indicates that it might exert its function on cellular proliferation. In our study, we investigated the effect of RP11-552D4.1 on the primary cultured neurons and found the knockdown of RP11-552D4.1 could inhibit the proliferation of neurons while overexpressed RP11-552D4.1 have the opposite effect, which indicates that RP11-552D4.1 is able to facilitate neuronal proliferation. And this is consistent with our findings that RP11-552D4.1 could become an independent risk factor for glioblastoma. From the Ensemble website, we could find the Ensembl ID of RP11-552D4.1 is ENSG00000228058 and the official symbol is LINC10736.
After literature review, LINC01736 is previously found to act as a protective factor in THCA, although with a high expression . We, therefore, investigated the expression of related prognosis in TCGA and found the expression of LINC01736 is all increased in GBM, LGG, and THCA, however, it is positively correlated with the survival status in LGG, while negatively associated with the OS in THCA and GBM (the latter did not reach the statistical difference) from the KM plot (Supplementary Figure 2). Then, we investigated the LINC01736-related immune cell infiltration with the ssGSEA method . We found LINC01736 is positively associated with the NK CD56bright enrichment in THCA and GBM (the latter did not reach the statistical difference), and negatively correlated with the NK CD56bright enrichment in LGG (Supplementary Figure 3), which suggests that LINC01736 exerts a different role in NK cell recruitment in different tumors. This explained why the expression of LINC01736 is increased in GBM and THCA, and seems to be a protectant factor in these tumors.
We have to commit that there are several limitations in this study. Firstly, the data was from public datasets and the sensitivity of drugs commonly used in the treatment of GBM, such as temozolomide has not been investigated; Secondly, the established risk model needs external verification, such as in-vivo study because the expression level of GBM shows heterogeneity, which results in the unreliability of the algorithm model. Although the lack of in-vivo validation, we applied COX and LASSO methods to establish the risk model with an independent prediction ability and identified one novel m6A-related lncRNA and further in-vivo studies are needed to validate its role in anti-GBM.
In summary, we used m6A-related lncRNAs to construct a risk model that exactly predicts the outcomes in GBM patients. Furthermore, this model reflects the dysregulated immune infiltration status in GBM and provides potential targets for immunotherapy.
Ping Zheng conceived and designed the study. Ping Zheng, Xiaoxue Zhang and Dabin Ren performed the data analysis. Dabin Ren did the cell study. Ping Zheng and Yisong Zhang wrote the manuscript. All authors read and approved the final version.
Conflicts of Interest
The authors declare no conflicts of interest related to this study.
This work was supported by Natural Science Foundation of China (Grant Number:81701231), and also funded by The Featured Clinical Discipline Project of Shanghai Pudong (PWYst2018-01) and Key Discipline Group Construction Project of Shanghai Pudong (PWZxq2017-02).
- 1. Dong Z, Cui H. The Emerging Roles of RNA Modifications in Glioblastoma. Cancers (Basel). 2020; 12:736. https://doi.org/10.3390/cancers12030736 [PubMed]
- 2. Rajesh Y, Pal I, Banik P, Chakraborty S, Borkar SA, Dey G, Mukherjee A, Mandal M. Insights into molecular therapy of glioma: current challenges and next generation blueprint. Acta Pharmacol Sin. 2017; 38:591–613. https://doi.org/10.1038/aps.2016.167 [PubMed]
- 3. Sun J, Ma Q, Li B, Wang C, Mo L, Zhang X, Tang F, Wang Q, Yan X, Yao X, Wu Q, Shu C, Xiong J, et al. RPN2 is targeted by miR-181c and mediates glioma progression and temozolomide sensitivity via the wnt/β-catenin signaling pathway. Cell Death Dis. 2020; 11:890. https://doi.org/10.1038/s41419-020-03113-5 [PubMed]
- 4. Mohme M, Maire CL, Schliffke S, Joosse SA, Alawi M, Matschke J, Schüller U, Dierlamm J, Martens T, Pantel K, Riethdorf S, Lamszus K, Westphal M. Molecular profiling of an osseous metastasis in glioblastoma during checkpoint inhibition: potential mechanisms of immune escape. Acta Neuropathol Commun. 2020; 8:28. https://doi.org/10.1186/s40478-020-00906-9 [PubMed]
- 5. Schipmann S, Müther M, Stögbauer L, Zimmer S, Brokinkel B, Holling M, Grauer O, Suero Molina E, Warneke N, Stummer W. Combination of ALA-induced fluorescence-guided resection and intraoperative open photodynamic therapy for recurrent glioblastoma: case series on a promising dual strategy for local tumor control. J Neurosurg. 2021; 134:426–36. https://doi.org/10.3171/2019.11.JNS192443 [PubMed]
- 6. Pierscianek D, Ahmadipour Y, Kaier K, Darkwah Oppong M, Michel A, Kebir S, Stuschke M, Glas M, Sure U, Jabbarli R. The SHORT Score for Preoperative Assessment of the Risk for Short-Term Survival in Glioblastoma. World Neurosurg. 2020; 138:e370–80. https://doi.org/10.1016/j.wneu.2020.02.131 [PubMed]
- 7. Zhang L, Wang H. Long Non-coding RNA in CNS Injuries: A New Target for Therapeutic Intervention. Mol Ther Nucleic Acids. 2019; 17:754–66. https://doi.org/10.1016/j.omtn.2019.07.013 [PubMed]
- 8. Shan W, Chen W, Zhao X, Pei A, Chen M, Yu Y, Zheng Y, Zhu S. Long noncoding RNA TUG1 contributes to cerebral ischaemia/reperfusion injury by sponging mir-145 to up-regulate AQP4 expression. J Cell Mol Med. 2020; 24:250–9. https://doi.org/10.1111/jcmm.14712 [PubMed]
- 9. Zhang X, Liu Z, Shu Q, Yuan S, Xing Z, Song J. LncRNA SNHG6 functions as a ceRNA to regulate neuronal cell apoptosis by modulating miR-181c-5p/BIM signalling in ischaemic stroke. J Cell Mol Med. 2019; 23:6120–30. https://doi.org/10.1111/jcmm.14480 [PubMed]
- 10. He B, Chen W, Zeng J, Tong W, Zheng P. Long noncoding RNA NKILA transferred by astrocyte-derived extracellular vesicles protects against neuronal injury by upregulating NLRX1 through binding to mir-195 in traumatic brain injury. Aging (Albany NY). 2021; 13:8127–45. https://doi.org/10.18632/aging.202618 [PubMed]
- 11. Crake RLI, Burgess ER, Royds JA, Phillips E, Vissers MCM, Dachs GU. The Role of 2-Oxoglutarate Dependent Dioxygenases in Gliomas and Glioblastomas: A Review of Epigenetic Reprogramming and Hypoxic Response. Front Oncol. 2021; 11:619300. https://doi.org/10.3389/fonc.2021.619300 [PubMed]
- 12. Kunadis E, Lakiotaki E, Korkolopoulou P, Piperi C. Targeting post-translational histone modifying enzymes in glioblastoma. Pharmacol Ther. 2021; 220:107721. https://doi.org/10.1016/j.pharmthera.2020.107721 [PubMed]
- 13. Zhang Y, Xu S, Xu G, Gao Y, Li S, Zhang K, Tian Z, Guo J, Li X, Xu J, Li Y. Dynamic Expression of m6A Regulators During Multiple Human Tissue Development and Cancers. Front Cell Dev Biol. 2021; 8:629030. https://doi.org/10.3389/fcell.2020.629030 [PubMed]
- 14. Liu F, Su X. Effects of m6A modifications on signaling pathways in human cancer (Review). Oncol Rep. 2021; 45:36. https://doi.org/10.3892/or.2021.7987 [PubMed]
- 15. Zhu W, Wang JZ, Wei JF, Lu C. Role of m6A methyltransferase component VIRMA in multiple human cancers (Review). Cancer Cell Int. 2021; 21:172. https://doi.org/10.1186/s12935-021-01868-1 [PubMed]
- 16. Tian S, Lai J, Yu T, Li Q, Chen Q. Regulation of Gene Expression Associated With the N6-Methyladenosine (m6A) Enzyme System and Its Significance in Cancer. Front Oncol. 2021; 10:623634. https://doi.org/10.3389/fonc.2020.623634 [PubMed]
- 17. Yang Y, Hsu PJ, Chen YS, Yang YG. Dynamic transcriptomic m6A decoration: writers, erasers, readers and functions in RNA metabolism. Cell Res. 2018; 28:616–24. https://doi.org/10.1038/s41422-018-0040-8 [PubMed]
- 18. Tang F, Wang H, Chen E, Bian E, Xu Y, Ji X, Yang Z, Hua X, Zhang Y, Zhao B. LncRNA-ATB promotes TGF-β-induced glioma cells invasion through NF-κB and P38/MAPK pathway. J Cell Physiol. 2019; 234:23302–14. https://doi.org/10.1002/jcp.28898 [PubMed]
- 19. Bian EB, Chen EF, Xu YD, Yang ZH, Tang F, Ma CC, Wang HL, Zhao B. Exosomal lncRNA-ATB activates astrocytes that promote glioma cell invasion. Int J Oncol. 2019; 54:713–21. https://doi.org/10.3892/ijo.2018.4644 [PubMed]
- 20. Ma CC, Xiong Z, Zhu GN, Wang C, Zong G, Wang HL, Bian EB, Zhao B. Long non-coding RNA ATB promotes glioma malignancy by negatively regulating miR-200a. J Exp Clin Cancer Res. 2016; 35:90. https://doi.org/10.1186/s13046-016-0367-2 [PubMed]
- 21. Li J, Bian EB, He XJ, Ma CC, Zong G, Wang HL, Zhao B. Epigenetic repression of long non-coding RNA MEG3 mediated by DNMT1 represses the p53 pathway in gliomas. Int J Oncol. 2016; 48:723–33. https://doi.org/10.3892/ijo.2015.3285 [PubMed]
- 22. Chen W, He B, Tong W, Zeng J, Zheng P. Astrocytic Insulin-Like Growth Factor-1 Protects Neurons Against Excitotoxicity. Front Cell Neurosci. 2019; 13:298. https://doi.org/10.3389/fncel.2019.00298 [PubMed]
- 23. Jiang X, Liu B, Nie Z, Duan L, Xiong Q, Jin Z, Yang C, Chen Y. The role of m6A modification in the biological functions and diseases. Signal Transduct Target Ther. 2021; 6:74. https://doi.org/10.1038/s41392-020-00450-x [PubMed]
- 24. Du J, Hou K, Mi S, Ji H, Ma S, Ba Y, Hu S, Xie R, Chen L. Malignant Evaluation and Clinical Prognostic Values of m6A RNA Methylation Regulators in Glioblastoma. Front Oncol. 2020; 10:208. https://doi.org/10.3389/fonc.2020.00208 [PubMed]
- 25. Xu F, Huang X, Li Y, Chen Y, Lin L. m6A-related lncRNAs are potential biomarkers for predicting prognoses and immune responses in patients with LUAD. Mol Ther Nucleic Acids. 2021; 24:780–91. https://doi.org/10.1016/j.omtn.2021.04.003 [PubMed]
- 26. Wang H, Meng Q, Ma B. Characterization of the Prognostic m6A-Related lncRNA Signature in Gastric Cancer. Front Oncol. 2021; 11:630260. https://doi.org/10.3389/fonc.2021.630260 [PubMed]
- 27. Zhang H, Zhao L, Li S, Wang J, Feng C, Li T, Du X. N6-Methylandenosine-Related lncRNAs in Tumor Microenvironment Are Potential Prognostic Biomarkers in Colon Cancer. Front Oncol. 2021; 11:697949. https://doi.org/10.3389/fonc.2021.697949 [PubMed]
- 28. Di Salvo TG, Guo Y, Su YR, Clark T, Brittain E, Absi T, Maltais S, Hemnes A. Right ventricular long noncoding RNA expression in human heart failure. Pulm Circ. 2015; 5:135–61. https://doi.org/10.1086/679721 [PubMed]
- 29. Zhang Y, Jin T, Shen H, Yan J, Guan M, Jin X. Identification of Long Non-Coding RNA Expression Profiles and Co-Expression Genes in Thyroid Carcinoma Based on The Cancer Genome Atlas (TCGA) Database. Med Sci Monit. 2019; 25:9752–69. https://doi.org/10.12659/MSM.917845 [PubMed]
- 30. Wang M, Zhou Z, Zheng J, Xiao W, Zhu J, Zhang C, Jiang X. Identification and Validation of a Prognostic Immune-Related Alternative Splicing Events Signature for Glioma. Front Oncol. 2021; 11:650153. https://doi.org/10.3389/fonc.2021.650153 [PubMed]