Identification of N6-Methyladenosine-Associated Long Non-coding RNAs for Immunotherapeutic Response and Prognosis in Patients With Pancreatic Cancer

Non-coding RNA-Mediated N6-Methyladenosine (m6A) deposition: A pivotal regulator of cancer, impacting key signaling pathways in carcinogenesis and therapy response

The emergence of RNA modifications has recently been considered as critical post-transcriptional regulations which governed gene expression. N6-methyladenosine (m6A) modification is the most abundant type of RNA modification which is mediated by three distinct classes of proteins called m6A writers, readers, and erasers. Accumulating evidence has been made in understanding the role of m6A modification of non-coding RNAs (ncRNAs) in cancer. Importantly, aberrant expression of ncRNAs and m6A regulators has been elucidated in various cancers. As the key role of ncRNAs in regulation of cancer hallmarks is well accepted now, it could be accepted that m6A modification of ncRNAs could affect cancer progression. The present review intended to discuss the latest knowledge and importance of m6A epigenetic regulation of ncRNAs including mircoRNAs, long non-coding RNAs, and circular RNAs, and their interaction in the context of cancer. Moreover, the current insight into the underlying mechanisms of therapy resistance and also immune response and escape mediated by m6A regulators and ncRNAs are discussed.

Non-coding RNA in tumor-infiltrating regulatory T cells formation and associated immunotherapy

Cancer immunotherapy has exhibited promising antitumor effects in various tumors. Infiltrated regulatory T cells (Tregs) in the tumor microenvironment (TME) restrict protective immune surveillance, impede effective antitumor immune responses, and contribute to the formation of an immunosuppressive microenvironment. Selective depletion or functional attenuation of tumor-infiltrating Tregs, while eliciting effective T-cell responses, represents a potential approach for anti-tumor immunity. Furthermore, it does not disrupt the Treg-dependent immune homeostasis in healthy organs and does not induce autoimmunity. Yet, the shared cell surface molecules and signaling pathways between Tregs and multiple immune cell types pose challenges in this process. Noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), regulate both cancer and immune cells and thus can potentially improve antitumor responses. Here, we review recent advances in research of tumor-infiltrating Tregs, with a focus on the functional roles of immune checkpoint and inhibitory Tregs receptors and the regulatory mechanisms of ncRNAs in Treg plasticity and functionality.

A prognostic model based on tumor microenvironment-related lncRNAs predicts therapy response in pancreatic cancer

Pancreatic cancer is an aggressive malignant tumor with high mortality and a low survival rate. The immune and stromal cells that infiltrate in the tumor microenvironment (TME) significantly impact immunotherapy and drug responses. Therefore, we identify the TME-related lncRNAs to develop a prognostic model for predicting the therapy efficacy in pancreatic cancer patients. Firstly, we identified differentially expressed genes (DEGs) for weighted gene co-expression network analysis (WGCNA) to identify the TME-related module eigengenes. According to the module eigengenes, the TME-related prognostic lncRNAs were screened through the univariate Cox, least absolute shrinkage and selection operator (LASSO), and multivariate Cox analyses to construct a prognostic risk score (RS) model. Next, the predictive power of this model was evaluated by the time-dependent receiver operating characteristic (ROC) curve and Kaplan-Meier analyses. In addition, functional enrichment, immune cell infiltration, and somatic mutation analyses were performed. Finally, tumor immune dysfunction and exclusion (TIDE) score and drug sensitivity analyses were applied to predict therapy response. In this study, 11 TME-related prognostic lncRNAs were identified to develop the prognostic RS model. According to the RS, the low-risk patients had a better prognosis, lower rates of somatic mutation, lower TIDE scores, and higher sensitivity to gemcitabine and paclitaxel compared to high-risk patients. The findings above suggested that low-risk patients may benefit more from immunotherapy, and high-risk patients may benefit more from chemotherapy. Within this study, we established a prognostic RS model based on 11 TME-related lncRNAs, which may help improve clinical decision-making.

Identification and validation of an m6A-related gene signature to predict prognosis and evaluate immune features of breast cancer

Breast cancer is the most prevalent cancer, and it is accompanied by high heterogeneity. N6-methyladenosine (m6A) modification significantly contributes to breast cancer tumorigenesis and progression. However, how m6A-related genes affect the clinical outcomes and tumor immune microenvironment (TIME) of breast cancer is largely unknown. Our study developed an m6A-related gene signature on the basis of The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The m6A-related gene signature was constructed using univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses. Breast cancer patients were classified into low- and high-risk groups depending on the median risk score. The reliability and efficiency of the signature were validated using Kaplan-Meier analysis, receiver operating characteristic (ROC) curves, and principal component analysis (PCA). The risk score was validated as an independent indicator associated with overall survival, and a nomogram model was created to estimate the overall survival of patients with breast cancer. Functional annotation suggested that the risk score had a strong relationship with immune-related pathways. Different proportions of immune cell infiltration between the two groups were evaluated using various algorithms. The high-risk group had higher immune checkpoint expression levels. We discovered that one of the 6 prognostic genes, TMEM71, was downregulated in breast cancer tissues. In vitro experiments indicated that overexpression of TMEM71 suppressed breast cancer cell proliferation and migration. In conclusion, the m6A-related gene signature may be a sensitive biomarker for overall survival prediction and guide the individualized treatment for breast cancer patients.

Identification of m6A-related lncRNAs for thyroid cancer recurrence

Background

Although the prognosis of thyroid cancer (THCA) is generally good, many patients have a high risk of recurrence after treatment. N6-methyladenosine (m6A)-related long noncoding RNAs (lncRNAs) have been extensively studied in recent years. However, the potential of m6A-related lncRNAs to predict recurrence in THCA is unknown.

Methods

RNA sequencing (RNA-seq) data and clinical information for THCA were downloaded from The Cancer Genome Atlas (TCGA). Differentially expressed lncRNAs (DELs) were identified using the R package DESeq2. A coexpression network based on m6A-related genes and lncRNAs was constructed. The CIBERSORT algorithm and gene set enrichment analysis (GSEA) were used for immune-infiltrating cell estimation and clustering functional enrichment analysis, respectively. A Kaplan-Meier plot was used for prognostic analysis based on m6A-associated lncRNA risk patterns. The expression of lncRNAs in recurrent and nonrecurrent THCA tissues was analyzed by real-time quantitative polymerase chain reaction (RT-qPCR).

Results

A network of m6A-related lncRNAs containing 8 lncRNAs was constructed with good predictive power for recurrence in THCA. A total of 3 clusters were obtained, and cluster 1 was most associated with THCA recurrence. We found significantly lower levels of CD8 T cells and follicular helper T cells, and significantly higher levels of dendritic cells (DCs), M2 macrophages, resting DCs, regulatory T cells, and mast cells in cluster 1 patients. Pathway analysis revealed significant enrichment in natural killer cell-mediated cytotoxicity, butyrate metabolism, and cell adhesion molecules in cluster 1. The m6A-related lncRNA risk model was effective in predicting progression-free survival (PFS) in patients with THCA recurrence. RT-qPCR analysis based on 40 THCA clinical samples from our center found the risk model to be a good predictor of recurrence in THCA patients.

Conclusions

In summary, m6A-related lncRNAs may provide a novel predictive method for prognostic relapse in THCA patients.

Molecular mechanism analysis of m6A modification-related lncRNA-miRNA-mRNA network in regulating autophagy in acute pancreatitis

This study aims to explore the molecular mechanism of N6-methyladenosine (m6A) modification-related long noncoding RNA (lncRNA)-microRNA (miRNA)-messenger RNA (mRNA) network in regulating autophagy and affecting the occurrence and development of acute pancreatitis (AP). RNA-seq datasets related to AP were obtained from Gene Expression Omnibus (GEO) database and merged after batch effect removal. lncRNAs significantly related to m6A in AP, namely candidate lncRNA, were screened by correlation analysis and differential expression analysis. In addition, candidate autophagy genes were screened through the multiple databases. Furthermore, the key pathways for autophagy to play a role in AP were determined by functional enrichment analysis. Finally, we predicted the miRNAs binding to genes and lncRNAs through TargetScan, miRDB and DIANA TOOLS databases and constructed two types of lncRNA-miRNA-mRNA regulatory networks mediated by upregulated and downregulated lncRNAs in AP. Nine lncRNAs related to m6A were differentially expressed in AP, and 21 candidate autophagy genes were obtained. Phosphoinositide 3-kinase (PI3K)-Akt signaling pathway and Forkhead box O (FoxO) signaling pathway might be the key pathways for autophagy to play a role in AP. Finally, we constructed a lncRNA-miRNA-mRNA regulatory network. An upregulated lncRNA competitively binds to 13 miRNAs to regulate 6 autophagy genes, and a lncRNA-miRNA-mRNA regulatory network in which 2 downregulated lncRNAs competitively bind to 7 miRNAs to regulate 2 autophagy genes. m6A modification-related lncRNA Pvt1, lncRNA Meg3 and lncRNA AW112010 may mediate the lncRNA-miRNA-mRNA network, thereby regulating autophagy to affect the development of AP.

Construction of a cancer-associated fibroblasts-related long non-coding RNA signature to predict prognosis and immune landscape in pancreatic adenocarcinoma

Background: Cancer-associated fibroblasts (CAFs) are an essential cell population in the pancreatic cancer tumor microenvironment and are extensively involved in drug resistance and immune evasion mechanisms. Long non-coding RNAs (lncRNAs) are involved in pancreatic cancer evolution and regulate the biological behavior mediated by CAFs. However, there is a lack of understanding of the prognostic signatures of CAFs-associated lncRNAs in pancreatic cancer patients.

Methods: Transcriptomic and clinical data for pancreatic adenocarcinoma (PAAD) and the corresponding mutation data were obtained from The Cancer Genome Atlas database. lncRNAs associated with CAFs were obtained using co-expression analysis. lncRNAs were screened by Cox regression analysis using least absolute shrinkage and selection operator (LASSO) algorithm for constructing predictive signature. According to the prognostic model, PAAD patients were divided into high-risk and low-risk groups. Kaplan-Meier analysis was used for survival validation of the model in the training and validation groups. Clinicopathological parameter correlation analysis, univariate and multivariate Cox regression, time-dependent receiver operating characteristic (ROC) curves, and nomogram were performed to evaluate the model. The gene set variation analysis (GSVA) and gene ontology (GO) analyses were used to explore differences in the biological behavior of the risk groups. Furthermore, single-sample gene set enrichment analysis (ssGSEA), tumor mutation burden (TMB), ESTIMATE algorithm, and a series of immune correlation analyses were performed to investigate the relationship between predictive signature and the tumor immune microenvironment and screen for potential responders to immune checkpoint inhibitors. Finally, drug sensitivity analyses were used to explore potentially effective drugs in high- and low-risk groups.

Results: The signature was constructed with seven CAFs-related lncRNAs (AP005233.2, AC090114.2, DCST1-AS1, AC092171.5, AC002401.4, AC025048.4, and CASC8) that independently predicted the prognosis of PAAD patients. Additionally, the high-risk group of the model had higher TMB levels than the low-risk group. Immune correlation analysis showed that most immune cells, including CD8⁺ T cells, were negatively correlated with the model risk scores. ssGSEA and ESTIMATE analyses further indicated that the low-risk group had a higher status of immune cell infiltration. Meanwhile, the mRNA of most immune checkpoint genes, including PD1 and CTLA4, were highly expressed in the low-risk group, suggesting that this population may be “hot immune tumors” and have a higher sensitivity to immune checkpoint inhibitors (ICIs). Finally, the predicted half-maximal inhibitory concentrations of some chemical and targeted drugs differ between high- and low-risk groups, providing a basis for treatment selection.

Conclusion: Our findings provide promising insights into lncRNAs associated with CAFs in PAAD and provide a personalized tool for predicting patient prognosis and immune microenvironmental landscape.

Risk Model and Immune Signature of m7G-Related lncRNA Based on Lung Adenocarcinoma

Lung cancer is a major cause of cancer-related deaths globally, with a dismal prognosis. N7-methylguanosine (m7G) is essential for the transcriptional phenotypic modification of messenger RNA (mRNA) and long noncoding RNA (lncRNA). However, research on m7G-related lncRNAs involved in lung adenocarcinoma (LUAD) regulation is still limited. Herein, we aim to establish a prognostic model of m7G-related lncRNAs and investigate their immune properties. Eight prognostic m7G-related lncRNAs were identified using univariate Cox analysis. Six m7G-related lncRNAs were identified using LASSO-Cox regression analysis to construct risk models, and all LUAD patients in The Cancer Genome Atlas (TCGA) cohort was divided into low-risk and high-risk subgroups. The accuracy of the model was verified by Kaplan-Meier analysis, time-dependent receiver operating characteristic, principal component analysis, independent prognostic analysis, nomogram, and calibration curve. Further studies were conducted on the gene set enrichment and disease ontology enrichment analyses. The gene set enrichment analysis (GSEA) revealed that the high-risk group enriched for cancer proliferation pathways, and the enrichment analysis of disease ontology (DO) revealed that lung disease was enriched, rationally explaining the superiority of the risk model. Finally, we found that the low-risk group had higher immune infiltration and checkpoint expression. It can be speculated that the low-risk group has a better effect on immunotherapy. Susceptibility to antitumor drugs in different risk subgroups was assessed, and it found that the high-risk group showed high sensitivity to first-line treatment drugs for non-small cell lung cancer. In conclusion, a risk model based on 6 m7G-related lncRNAs can not only predict the overall survival (OS) rate of LUAD patients but also guide individualized treatment for these patients.

Construction of m6A-Related lncRNA Prognostic Signature Model and Immunomodulatory Effect in Glioblastoma Multiforme

Background

Glioblastoma multiforme (GBM), the most prevalent and aggressive of primary malignant central nervous system tumors (grade IV), has a poor clinical prognosis. This study aimed to assess and predict the survival of GBM patients by establishing an m6A-related lncRNA signaling model and to validate its validity, accuracy and applicability.

Methods

RNA sequencing data and clinical data of GBM patients were obtained from TCGA data. First, m6A-associated lncRNAs were screened and lncRNAs associated with overall survival in GBM patients were obtained. Subsequently, the signal model was established using LASSO regression analysis, and its accuracy and validity are further verified. Finally, GO enrichment analysis was performed, and the influence of this signature on the immune regulation response and anticancer drug sensitivity of GBM patients was discussed.

Results

The signature constructed by four lncRNAs AC005229.3, SOX21-AS1, AL133523.1, and AC004847.1 is obtained. Furthermore, the signature proved to be effective and accurate in predicting and assessing the survival of GBM patients and could function independently of other clinical characteristics (Age, Gender and IDH1 mutation). Finally, Immunosuppression-related factors, including APC co-inhibition, T-cell co-inhibition, CCR and Check-point, were found to be significantly up-regulated in GBM patients in the high-risk group. Some chemotherapeutic drugs (Doxorubicin and Methotrexate) and targeted drugs (AZD8055, BI.2536, GW843682X and Vorinostat) were shown to have higher IC50 values in patients in the high-risk group.

Conclusion

We constructed an m6A-associated lncRNA risk model to predict the prognosis of GBM patients and provide new ideas for the treatment of GBM. Further biological experiments can be conducted on this basis to validate the clinical value of the model.