Pan-cancer characterization of expression and clinical relevance of m6A-related tissue-elevated long non-coding RNAs
Molecular Cancer volume 20, Article number: 31 (2021)
N6-methyladenosine (m6A) has become a critical internal RNA modification, and it plays important roles in the development and progression of cancer . m6A has also been found in diverse non-coding RNAs, such as microRNAs and long noncoding RNAs (lncRNAs) . LncRNAs comprise a large class of RNA transcripts and are critical regulators of gene expression. The regulatory effectiveness of lncRNAs is closely associated with spatial expression, whose dysregulation often influences cancer development and progression . For these reasons, global characterization of lncRNA spatial expression across tissues or cancers could improve our understanding of lncRNA functions. Recently, LncRNA Spatial Atlas (LncSpA) and landscape of m6A have been proposed as valuable resources to understand lncRNA and m6A regulatory functions across different tissues [4, 5]. However, we still lack understanding of the distribution and functions of m6A modification in lncRNAs, particularly the tissue-elevated (TE) lncRNAs.
In this study, we aimed to systematically characterize the distribution and clinical relevance of m6A-related TE lncRNAs across tissues and cancer types. We found that TE lncRNAs were found to be regulated by m6A modification across tissues, particular brain tissues. We also investigated the correlation between expression of m6A regulators and TE lncRNAs, and found that numbers of m6A-related TE lncRNAs were associated with expression of m6A regulators. We assessed the clinical prognostic values of m6A-regulated TE lncRNAs. We identified several m6A-related TE lncRNAs as potentially useful markers for prognostic stratification. Our analysis highlights the importance of m6A modification in the regulation of lncRNA expression and helps bridge the knowledge gap between lncRNA expression and phenotypes.
TE lncRNAs are associated with m6A modification across tissues
We first retrieved the TE lncRNAs from 38 normal tissues from LncSpA in 4 data resources (Fig. 1a), including Human Body Map (HBM2.0), Human Protein Atlas (HPA), the Genotype-Tissue Expression (GTEx), and the Function Annotation Of The Mammalian Genome (FANTOM) project. In total, 9837, 13,337, 10,718, and 74,767 TE lncRNAs were obtained from GTEx, HPA, HBM2.0, and FANTOM5, respectively. Higher numbers of TE lncRNAs were found in tissues of the brain and testis tissues than in other tissues (Fig. 1b and Additional file 1: Table S1). Next, we mapped all the m6A modification peaks to lncRNAs and identified approximately 511–1600 lncRNAs regulated by m6A across tissues (Fig. 1c, Additional file 2: Figure S1 and Additional file 3: Table S2). We next assessed the proportion of m6A-modified TE lncRNAs among human tissues. We found that brain tissues had the highest proportion of TE lncRNAs with m6A modifications (Fig. 1d and Additional file 2: Figure S2). Approximately 14.89–19.20% TE lncRNAs were m6A-modified in brain tissues than in other tissues in four data resources. Although there were higher numbers of TE lncRNAs in testis tissues, the proportion of m6A-modified TE lncRNAs was small (Fig. 1d).
Next, we compared the overlap of m6A modified TE lncRNAs among tissues from different data resources. The Simpson index was calculated for two tissues from different sources. High correlations were observed for the same tissues across different sources (Fig. 1e), suggesting that m6A modified TE lncRNAs were conserved across different resources. To investigate potential tissue specificity of the m6A-modified TE lncRNAs, we calculated the percentage of m6A-modified TE lncRNAs and non-TE lncRNAs in each tissue. There were no significant differences observed for the two lncRNA categories for the most tissues, which is consistent with the observations in protein coding genes . However, the proportion of m6A-modified TE lncRNAs is significantly higher than that of non-TE lncRNAs in brain tissues (Fig. 1f and Additional file 2: Figure S3). We explored the number of m6A peaks for lncRNAs across tissues. We found that the majority of m6A peaks in lncRNAs were in brain tissues (Additional file 2: Figure S4). Collectively, these results indicated that TE lncRNAs are associated with m6A modification across tissues and are more prone to be regulated by m6A in brain tissues than in other tissues.
Co-expression network of TE lncRNAs and m6A regulators
The regulatory effects of m6A modification are primarily determined by regulators, including readers, writers, and erasers . The extent to which variation in m6A modification of TE lncRNAs may be attributed to the expression of m6A regulators remains unknown. Thus, we next sought to analyze the correlation between the expressions of m6A-modified TE lncRNAs and regulators. In total, we identified 4862 correlations among 860 TE lncRNAs and 20 m6A regulators in 4 resources (Additional file 2: Figure S5A and Additional file 4: Table S3). Numbers of TE lncRNAs were associated with expression of m6A regulators in all four sources, including AC091878.1, LINC00854 and AC007879.5 (Additional file 2: Figure S5B). In contrast, we calculated the number of TE lncRNAs correlated with each m6A regulators. Higher numbers of TE lncRNAs were found to be correlated with the expression of IGF2BP1, METTL3 and VIRMA (Additional file 2: Figure S6).
Notably, we identified several TE lncRNA-regulator pairs that had been verified in literature. We took PVT1 as an example and found its expression to be significantly correlated with YTHDF2 (Additional file 2: Figure S5C, R = 0.64, P = 0.0003). Evidence has shown that YTHDF2 and PVT1 interact and that YTHDF2 plays critical roles in the stability of PVT1 . Another example is SOX2-OT, which has been reported to play an oncogenic role in cancer. It was identified as a TE lncRNA in brain tissues from all four sources. We found its expression to be significantly closely correlated with HNRNPA2B1 (Additional file 2: Figure S5D, R = 0.56, P = 0.0005). It has been shown that SOX2-OT can regulate cancer proliferation and metastasis through the miR-146b-5p/HNRNPA2B1 pathway. We also found a significant correlation between KCNK15 − AS1 and ALKBH5 (Additional file 2: Figure S5E, R = 0.48, P = 0.0108). ALKBH5 had been demonstrated to inhibit cancer motility by demethylating lncRNA KCNK15-AS1 . Together, all these results suggest that m6A modification of TE lncRNAs is partially regulated by the expression of m6A regulators.
Association of m6A-modified TE lncRNAs with tumor prognosis
LncRNA has been identified as a biomarker suitable for the classification of cancer patients. We next investigated the relationship between expression of m6A modified TE lncRNAs and patient survival. We first manually mapped the m6A modification in human tissues to cancer types and identified 104–621 TE lncRNAs in 16 cancers (Fig. 2a). Cancers with similar tissue of origin were clustered together based on the overlap of TE lncRNAs, such as LGG and GBM, COAD, and READ. In addition, numbers of m6A-modified TE lncRNAs were identified across cancer types, ranging from 3 to 105 (Fig. 2a and Additional file 2: Figure S7).
We next explored the differences in survival between patients with high- and low-levels of lncRNA expression and identified 83 protective and 18 risky m6A-modified TE lncRNAs across cancer types (Fig. 2a and Additional file 5: Table S4). Moreover, we identified 28 m6A-modified TE lncRNAs that had significantly higher expression in cancer patients than in healthy controls and 8 m6A-modified TE lncRNAs that had significantly lower expression (Fig. 2a and Additional file 5: Table S4). There were two m6A-modified TE lncRNAs (F11-AS1 and LINC01018) showing significantly lower expression in hepatocellular carcinoma patients than in controls, and these lower expressions were associated with worse survival rates (Fig. 2b-e). F11-AS1 can inhibit HBV-related hepatocellular carcinoma progression by regulating NR1I3 via binding to microRNA-211-5p. LINC01018 has a novel tumor suppressor role in hepatocellular carcinoma by sponging miR-182-5p [9, 10]. We also found lower expression of m6A-modified MIR325HG to be correlated with worse patient survival in both LGG and GBM (Fig. 2f-g). These results suggest that these TE lncRNAs could be potentially tumor suppressors in cancer.
We next tried to determine the functions of F11-AS1, LINC01018 and MIR325HG. We performed Gene Set Enrichment Analysis (GSEA) on cancer patients. We found that these m6A-modified lncRNAs were involved in a number of cancer hallmark-related functions (Additional file 2: Figure S8 and Additional file 6: Table S5), such as DNA repair and epithelial mesenchymal transition pathways (Additional file 2: Figure S9). Taken together, all these results suggest a connection between m6A modified TE lncRNAs and the risk of diseases.
We have shown the prevalence of m6A modification in TE lncRNAs across tissues and cancer types. The expression levels of m6A-modified TE lncRNAs were significantly closely associated with the activity of m6A regulators. Several studies have also shown that m6Am can regulate the expression of noncoding RNAs. Thus, it would also be interesting to integrate such m6A and m6Am data to identify potential lncRNA biomarkers in cancer. In summary, our work reveals the landscape of m6A-modified TE lncRNAs and provides a valuable resource for functional studies of m6A and lncRNA functions in the future.
Availability of data and materials
The gene expression profiles and clinical data can be found at the GDC portal (https://portal.gdc.cancer.gov/). The TE lncRNAs across tissues were obtained from LncSpA (http://bio-bigdata.hrbmu.edu.cn/LncSpA/). Software and resources used for the analyses are described in each method section. All results generated in this study can be found in supplementary tables.
methylation of N6 adenosine
Brain low grade glioma
Long non-coding RNA
LncRNA Spatial Atlas of expression
Gene Set Enrichment Analysis
The Genotype-Tissue Expression
Human Protein Atlas
Human body epigenome maps
Functional Annotation of the Mamaliam Genome
Wang S, Chai P, Jia R, Jia R. Novel insights on m (6) a RNA methylation in tumorigenesis: a double-edged sword. Mol Cancer. 2018;17:101.
Chen Y, Lin Y, Shu Y, He J, Gao W. Interaction between N (6)-methyladenosine (m (6) a) modification and noncoding RNAs in cancer. Mol Cancer. 2020;19:94.
Hon CC, Ramilowski JA, Harshbarger J, Bertin N, Rackham OJ, Gough J, Denisenko E, Schmeier S, Poulsen TM, Severin J, Lizio M, Kawaji H, Kasukawa T, Itoh M, Burroughs AM, Noma S, Djebali S, Alam T, Medvedeva YA, Testa AC, Lipovich L, Yip CW, Abugessaisa I, Mendez M, Hasegawa A, Tang D, Lassmann T, Heutink P, Babina M, Wells CA, Kojima S, Nakamura Y, Suzuki H, Daub CO, de Hoon MJ, Arner E, Hayashizaki Y, Carninci P, Forrest AR. An atlas of human long non-coding RNAs with accurate 5′ ends. Nature. 2017;543:199–204.
Lv D, Xu K, Jin X, Li J, Shi Y, Zhang M, Jin X, Li Y, Xu J, Li X. LncSpA: LncRNA spatial atlas of expression across Normal and Cancer tissues. Cancer Res. 2020;80:2067–71.
Liu J, Li K, Cai J, Zhang M, Zhang X, Xiong X, Meng H, Xu X, Huang Z, Peng J, Fan J, Yi C. Landscape and regulation of m (6) a and m (6) am Methylome across human and mouse tissues. Mol Cell. 2020;77:426–440 e426.
Li Y, Xiao J, Bai J, Tian Y, Qu Y, Chen X, Wang Q, Li X, Zhang Y, Xu J. Molecular characterization and clinical relevance of m (6) a regulators across 33 cancer types. Mol Cancer. 2019;18:137.
Chen S, Zhou L, Wang Y. ALKBH5-mediated m (6) a demethylation of lncRNA PVT1 plays an oncogenic role in osteosarcoma. Cancer Cell Int. 2020;20:34.
He Y, Hu H, Wang Y, Yuan H, Lu Z, Wu P, Liu D, Tian L, Yin J, Jiang K, Miao Y. ALKBH5 inhibits pancreatic Cancer motility by decreasing long non-coding RNA KCNK15-AS1 methylation. Cell Physiol Biochem. 2018;48:838–46.
Wang S, Xu M, Sun Z, Yu X, Deng Y, Chang H. LINC01018 confers a novel tumor suppressor role in hepatocellular carcinoma through sponging microRNA-182-5p. Am J Physiol Gastrointest Liver Physiol. 2019;317:G116–26.
Deng Y, Wei Z, Huang M, Xu G, Wei W, Peng B, Nong S, Qin H. Long non-coding RNA F11-AS1 inhibits HBV-related hepatocellular carcinoma progression by regulating NR1I3 via binding to microRNA-211-5p. J Cell Mol Med. 2020;24:1848–65.
The results here are in whole based upon data generated by the TCGA Research Network, GTEx, Human Body Map, Human Protein Atlas and FANTOM projects.
This work was supported by the National Natural Science Foundation of China (31871338, 32060152, 32070673, 31970646, 61873075); Major science and technology projects of Hainan Province (ZDKJ202003); Hainan Provincial Natural Science Foundation of China (820MS053); the National Key R&D Program of China (2018YFC2000100); Natural Science Foundation for Distinguished Young Scholars of Heilongjiang Province (JQ2019C004) and Heilongjiang Touyan Innovation Team Program. The funders played no roles in the study design, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
Patient data we used were acquired by publicly available datasets that were collected with patients’ informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Number of TE lncRNAs, m6A-modified lncRNAs and the proportion of m6A-modified TE lncRNAs across tissues in four resources.
Supplemental materials and methods, and supplemental figure S1-S9.Figure S1. Numbers of m6A-regulated lncRNAs across tissues. Figure S2. Numbers of m6A-regulated TE lncRNAs across tissues in four data resources. Figure S3. Distribution of odds ratios for comparison between TE lncRNAs and non-TE lncRNAs across tissues in four data resources. Figure S4. Number of m6A peaks correlated with TE lncRNAs in four resources. Figure S5.Co-expression between m6A regulators and m6A-modified TE lncRNAs. A, River plot showing the expression correlation between m6A modified TE lncRNAs and m6A regulators. B, Bar plots showing the number of m6A regulators correlated with each m6A modified TE lncRNA. Color indicated the different data resources. C-E, Scatter plots showing the correlation between the expression of lncRNAs and m6A regulators. C for PVT1 and YTHDF2; D for SOX2-OT and HNRNPA2B1; E for KCNK15-AS1 and ALKBH5. Figure S6. Numbers of TE lncRNAs correlated with m6A regulators. Figure S7. Number of m6A modified TE lncRNAs across cancer types. Figure S8. GSEA for m6A modified lncRNAs in HCC, LGG and GBM. A, F11-AS1 in HCC; B, LINC01018 in HCC; C, MIR325HG in LGG, D, MIR325HG in GBM. Figure S9. GSEA figures for m6A-modified lncRNAs in HCC, LGG and GBM. A, F11-AS1 enriched in DNA repair pathway in HCC; B, LINC01018 enriched in DNA repair pathway in HCC; C, MIR325HG enriched in EMT pathway in LGG, D, MIR325HG enriched in EMT pathway in GBM.
List of m6A-regulated TE lncRNAs across tissues and cancers.
Expression correlation between m6A regulators and lncRNAs.
Clinical association of m6A-modified TE lncRNAs.
GSEA results for four m6A-modified TE lncRNAs.
About this article
Cite this article
Xu, K., Cai, Y., Zhang, M. et al. Pan-cancer characterization of expression and clinical relevance of m6A-related tissue-elevated long non-coding RNAs. Mol Cancer 20, 31 (2021). https://doi.org/10.1186/s12943-021-01324-8