N6-Methyladenosine Regulators and Related LncRNAs Are Potential to be Prognostic Markers for Uveal Melanoma and Indicators of Tumor Microenvironment Remodeling

Novel methylation-related long non-coding RNA clinical outcome prediction method: the clinical phenotype and immune infiltration research in low-grade gliomas

Background

Recent studies have suggested that long non-coding RNAs (lncRNAs) may play crucial role in low-grade glioma; however, the underlying mechanisms linking them to epigenetic methylation remain unclear.

Methods

We downloaded expression level data for regulators associated with N1 methyladenosine (m1A), 5-methyladenine (m5C), and N6 methyladenosine (m6A) (M1A/M5C/M6A) methylation from the Cancer Genome Atlas-low-grade glioma (TCGA-LGG) database. We identified the expression patterns of lncRNAs, and selected methylation-related lncRNAs using Pearson correlation coefficient>0.4. Non-negative matrix dimensionality reduction was then used to determine the expression patterns of the methylation-associated lncRNAs. We constructed a weighted gene co-expression network analysis (WGCNA) network to explore the co-expression networks between the two expression patterns. Functional enrichment of the co-expression network was performed to identify biological differences between the expression patterns of different lncRNAs. We also constructed prognostic networks based on the methylation presence in lncRNAs in low-grade gliomas.

Results

We identified 44 regulators by literature review. Using a correlation coefficient greater than 0.4, we identified 2330 lncRNAs, among which 108 lncRNAs with independent prognostic values were further screened using univariate Cox regression at P< 0.05. Functional enrichment of the co-expression networks revealed that regulation of trans-synaptic signaling, modulation of chemical synaptic transmission, calmodulin binding, and SNARE binding were mostly enriched in the blue module. The calcium and CA2 signaling pathways were associated with different methylation-related long non-coding chains. Using the Least Absolute Shrinkage Selector Operator (LASSO) regression analysis, we analyzed a prognostic model containing four lncRNAs. The model's risk score was 1.12 *AC012063 + 0.74 * AC022382 + 0.32 * AL049712 + 0.16 * GSEC. Gene set variation analysis (GSVA) revealed significant differences in mismatch repair, cell cycle, WNT signaling pathway, NOTCH signaling pathway, Complement and Cascades, and cancer pathways at different GSEC expression levels. Thus, these results suggest that GSEC may be involved in the proliferation and invasion of low-grade glioma, making it a prognostic risk factor for low-grade glioma.

Conclusion

Our analysis identified methylation-related lncRNAs in low-grade gliomas, providing a foundation for further research on lncRNA methylation. We found that GSEC could serve as a candidate methylation marker and a prognostic risk factor for overall survival in low-grade glioma patients. These findings shed light on the underlying mechanisms of low-grade glioma development and may facilitate the development of new treatment strategies.

Genetics and RNA Regulation of Uveal Melanoma

Uveal melanoma (UM) is the most common intraocular malignant tumor and the most frequent melanoma not affecting the skin. While the rate of UM occurrence is relatively low, about 50% of patients develop metastasis, primarily to the liver, with lethal outcome despite medical treatment. Notwithstanding that UM etiopathogenesis is still under investigation, a set of known mutations and chromosomal aberrations are associated with its pathogenesis and have a relevant prognostic value. The most frequently mutated genes are BAP1, EIF1AX, GNA11, GNAQ, and SF3B1, with mutually exclusive mutations occurring in GNAQ and GNA11, and almost mutually exclusive ones in BAP1 and SF3B1, and BAP1 and EIF1AX. Among chromosomal aberrations, monosomy of chromosome 3 is the most frequent, followed by gain of chromosome 8q, and full or partial loss of chromosomes 1 and 6. In addition, epigenetic mechanisms regulated by non-coding RNAs (ncRNA), namely microRNAs and long non-coding RNAs, have also been investigated. Several papers investigating the role of ncRNAs in UM have reported that their dysregulated expression affects cancer-related processes in both in vitro and in vivo models. This review will summarize current findings about genetic mutations, chromosomal aberrations, and ncRNA dysregulation establishing UM biology.

Epitranscriptomic mRNA modifications governing plant stress responses: underlying mechanism and potential application

Plants inevitably encounter environmental adversities, including abiotic and biotic stresses, which significantly impede plant growth and reduce crop yield. Thus, fine-tuning the fate and function of stress-responsive RNAs is indispensable for plant survival under such adverse conditions. Recently, post-transcriptional RNA modifications have been studied as a potent route to regulate plant gene expression under stress. Among over 160 mRNA modifications identified to date, N⁶ -methyladenosine (m⁶ A) in mRNAs is notable because of its multifaceted roles in plant development and stress response. Recent transcriptome-wide mapping has revealed the distribution and patterns of m⁶ A in diverse stress-responsive mRNAs in plants, building a foundation for elucidating the molecular link between m⁶ A and stress response. Moreover, the identification and characterization of m⁶ A writers, readers and erasers in Arabidopsis and other model crops have offered insights into the biological roles of m⁶ A in plant abiotic stress responses. Here, we review the recent progress of research on mRNA modifications, particularly m⁶ A, and their dynamics, distribution, regulation and biological functions in plant stress responses. Further, we posit potential strategies for breeding stress-tolerant crops by engineering mRNA modifications and propose the future direction of research on RNA modifications to gain a much deeper understanding of plant stress biology.

Identification of N6-methylandenosine related lncRNA signatures for predicting the prognosis and therapy response in colorectal cancer patients

Despite recent advances in surgical and multimodal therapies, the overall survival (OS) of advanced colorectal cancer (CRC) patients remains low. Thus, discerning sensitive prognostic biomarkers to give the optimistic treatment for CRC patients is extremely critical. N6-methyladenosine (m6A) and long noncoding RNAs (lncRNAs) play an important role in CRC progression. Nonetheless, few studies have focused on the impact of m6A-related lncRNAs on the prognosis, tumor microenvironment (TME) and treatment of CRC. In this study, 1707 m6A-related lncRNAs were identified through Pearson correlation analysis and Weighted co-expression network analysis (WGCNA) using The Cancer Genome Atlas (TCGA) cohort. Then, 28 m6A-related prognostic lncRNAs were screened by univariate Cox regression analysis, followed by identifying two clusters by consensus clustering analysis. A prognostic model consisted of 8 lncRNA signatures was constructed by the least absolute shrinkage and selection operator (LASSO). Kaplan–Meier curve analysis and a nomogram were performed to investigate the prognostic ability of this model. The risk score of prognostic model act as an independent risk factor for OS rate. Functional enrichment analysis indicated that lncRNA signatures related tumor immunity. The low-risk group characterized by increased microsatellite instability-high (MSI-H), mutation burden, and immunity activation, indicated favorable odds of OS. Moreover, the lncRNA signatures were significantly associated with the cancer stem cell (CSC) index and drug sensitivity. In addition, 3 common immune genes shared by the lncRNA signatures were screened out. We found that these immune genes were widely distributed in 2 cell types of TME. Finally, a ceRNA network was constructed to identify ZEB1-AS1 regulatory axis in CRC. We found that ZEB1-AS1 was significantly overexpressed in tumor tissues, and was related to the metastasis of EMT and the chemoresistance of 5-Fu in CRC. Therefore, our study demonstrated the important role of m6A-related lncRNAs in TME remodeling. Moreover, these results illustrated the levels of ZEB1-AS1 might be valuable for predicting the progression and prognosis of CRC, and further provided a new target for the diagnosis and treatment of CRC patients.

m6A-Related lncRNA Signature Is Involved in Immunosuppression and Predicts the Patient Prognosis of the Age-Associated Ovarian Cancer

Background

Epithelial ovarian cancers are age-associated diseases, usually diagnosed at an advanced stage. lncRNA has been discovered to interplay with N6-methyladenosine (m6A), working in tandem to promote cancer progression and worsening patient outcomes. This study is aimed at investigating the roles and mechanism of m6A-related lncRNA signature on ovarian cancers.

Methods

We retrieved TCGA and CGGA sequencing data to identify m6A-related lncRNA signature and constructed an m6A score (MS) using the LASSO algorithm. A clinical nomogram was then established to predict the overall survival of patients. Subsequently, GSEA analyses were conducted to obtain pathways involved. Expression of HLA genes, 28 tumor-infiltrating lymphocyte infiltration, and anticancer cycle were analyzed the immunological differences between high-MS and low-MS groups. Finally, immune checkpoint gene expressions and IC₅₀ of chemotherapeutic drugs were calculated, and CMap was run to identify the potential compounds and their corresponding mechanisms.

Results

We identified 16 m6A-related lncRNAs and constructed an MS model. The high-MS group showed a poor prognosis. A clinical nomogram consists of MS, and age was constructed and predicted the 1-, 3-, and 5-year survival with high accuracy. GSEA analyses presented downregulated antigen processing and presentation pathways. Immunocyte infiltrating analyses demonstrated that high-MS was associated with high infiltration of Treg cells, macrophages, and low Th1/Th2 rate. Also, high expression of immune checkpoint genes NRP1, TNFSF9, and VSIR was observed in the high-MS group. Finally, the high-MS group also predicted low IC₅₀ of vinorelbine and vorinostat.

Conclusion

This study constructed a robust prediction model for prognostic management and revealed the cross-talk between m6A and immunosuppression. Besides, the m6A lncRNA signature can predict the chemotherapeutic drug response. These will shed light on the development of novel therapeutic strategies and render survival benefits for ovarian patients.

Potential Prognostic Value of a Seven m6A-Related LncRNAs Signature and the Correlative Immune Infiltration in Colon Adenocarcinoma

Long non-coding RNAs (lncRNAs) and their N6-methyladenosine (m6A) modifications play an essential role in tumorigenesis and cancer progression. This study was designed to explore the value of m6A-related lncRNAs in prognosis and therapeutic applications of immune infiltration of colon adenocarcinoma (COAD). We downloaded the COAD gene expression and clinical data from The Cancer Genome Atlas project. By co-expression analysis, Lasso Cox regression analysis, and univariate and multivariate Cox regression, we constructed an independent prognostic signature of seven m6A-related lncRNAs. The prognostic lncRNAs were divided into two clusters by consistent clustering analysis, as well as into two groups of low–high risk based on the signature. Then we identified the relationship between the different groups with clinical features and immune cell infiltration. Cluster 2 had a higher risk score with a lower survival rate. The risk score was higher in groups with advanced clinical features, such as stage III–IV, N1-3, and M1. The expression of AC156455.1 was increased in tumor tissues and cluster 2, and the lncRNA ZEB1−AS1 was notably higher in the high-risk group. Five types of immune cells showed differences in two clusters, and most were upregulated in type 2. The expression of memory B cells was positively correlated with the risk score. The prognostic model was verified by the Gene Expression Omnibus (GEO) dataset. Besides, we found that the expression of these seven lncRNAs in tumor tissues was significantly higher than that in normal tissues, which verified the feasibility of the model. Thus, the signature of seven m6A-related lncRNAs can independently predict the prognosis of COAD. This signature is also closely associated with immune cell infiltration, and new therapeutic targets can be explored from this field.