Identification of RNA Methylation-Related lncRNAs Signature for Predicting Hot and Cold Tumors and Prognosis in Colon Cancer

DUSP5 regulated by YTHDF1-mediated m6A modification promotes epithelial-mesenchymal transition and EGFR-TKI resistance via the TGF-β/Smad signaling pathway in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) patients have a dismal survival rate because of cancer metastasis and drug resistance. The study aims to identify the genes that concurrently modulate EMT, metastasis and EGFR-TKI resistance, and to investigate the underlying regulatory mechanisms.

METHODS

Cox regression and Kaplan-Meier analyses were applied to identify prognostic oncogenes in LUAD. Gene set enrichment analysis (GSEA) was used to indicate the biological functions of the gene. Wound-healing and Transwell assays were used to detect migratory and invasive ability. EGFR-TKI sensitivity was evaluated by assessing the proliferation, clonogenic survival and metastatic capability of cancer cells with treatment with gefitinib. Methylated RNA immunoprecipitation (MeRIP) and RNA immunoprecipitation (RIP) analyses established the level of m6A modification present on the target gene and the protein's capability to interact with RNA, respectively. Single-sample gene set enrichment (ssGSEA) algorithm used to investigate levels of immune cell infiltration.

RESULTS

Our study identified dual-specificity phosphatase 5 (DUSP5) as a novel and powerful predictor of adverse outcomes for LUAD by using public datasets. Functional enrichment analysis found that DUSP5 was positively enriched in EMT and transforming growth factor-beta (TGF-β) signaling pathway, a prevailing pathway involved in the induction of EMT. As expected, DUSP5 knockdown suppressed EMT via inhibiting the canonical TGF-β/Smad signaling pathway in in vitro experiments. Consistently, knockdown of DUSP5 was first found to inhibit migratory ability and invasiveness of LUAD cells in in vitro and prevent lung metastasis in in vivo. DUSP5 knockdown re-sensitized gefitinib-resistant LUAD cells to gefitinib, accompanying reversion of EMT progress. In LUAD tissue samples, we found 14 cytosine-phosphate-guanine (CpG) sites of DUSP5 that were negatively associated with DUSP5 gene expression. Importantly, 5'Azacytidine (AZA), an FDA-approved DNA methyltransferase inhibitor, restored DUSP5 expression. Moreover, RIP experiments confirmed that YTH N6-methyladenosine RNA binding protein 1 (YTHDF1), a m6A reader protein, could bind DUSP5 mRNA. YTHDF1 promoted DUSP5 expression and the malignant phenotype of LUAD cells. In addition, the DUSP5-derived genomic model revealed the two clusters with distinguishable immune features and tumor mutational burden (TMB).

CONCLUSIONS

Briefly, our study discovered DUSP5 which was regulated by epigenetic modification, might be a potential therapeutic target, especially in LUAD patients with acquired EGFR-TKI resistance.

LASSO regression and WGCNA-based telomerase-associated lncRNA signaling predicts clear cell renal cell carcinoma prognosis and immunotherapy response

OBJECTIVE

To investigate whether telomerase-associated lncRNA expression affects the prognosis and anti-tumor immunity of patients with renal clear cell carcinoma (ccRCC).

METHODS

A series of analyses were performed to establish a prognostic risk model and validate its accuracy. Immune-related analyses were performed to assess further the association between immune status, tumor microenvironment, and prognostic risk models.

RESULTS

Eight telomerase-associated lncRNAs associated with prognosis were identified and applied to establish a prognostic risk model. Overall survival was higher in the low-risk group.

CONCLUSION

The established prognostic risk model has a good predictive ability for the prognosis of ccRCC patients and provides a new possible therapeutic target for ccRCC.

RNA methylation-related genes INHBB and SOWAHA are associated with MSI status in colorectal cancer patients and may serve as prognostic markers for predicting immunotherapy efficacy

The role of RNA methylation is vital in the advancement and spread of tumors. However, its exact role in microsatellite instability in colorectal cancer (CRC) is still not fully understood. To address this gap in knowledge, this study investigated the impact of genes associated with RNA methylation on the prognosis and response to immunotherapy in individuals diagnosed with low microsatellite instability (MSI-L) or microsatellite stable (MSS) CRC. The differentially expressed genes (DEGs) in two groups of patients: those with high microsatellite instability (MSI-H) and those with MSI-L/MSS was thoroughly investigated and compared with aims of exploring the association between them and the 60 RNA methylation regulators. We employed these genes and developed an MSI-RMscore to establish a risk signature capable of forecasting patient outcomes. Furthermore, an investigation of the immunophenotypic traits was conducted encompassing patients categorized as high-risk and low-risk. By combining the MSI-RMscore and clinicopathological features, a predictive nomogram was developed, which was subsequently validated using the GEO database. Furthermore, immunohistochemistry was employed to establish the correlation between INHBB and SOWAHA and the MSI status, as well as patient prognosis. Our findings indicated that the high-risk subgroup exhibited unfavorable overall survival rates, reduced responsiveness to immune checkpoint blockers, elevated estimate scores, and increased infiltration of macrophages and fibroblasts. We also confirmed that INHBB and SOWAHA were associated with CRC patient prognosis and MSI status, as well as immunotherapy response. These findings suggest that targeting INHBB and SOWAHA could be a promising strategy to enhance patient responsiveness to immunotherapy.

Construction of RNA Methylation Modification-immune-related lncRNA Molecular Subtypes and Prognostic Scoring System in Lung Adenocarcinoma

BACKGROUND

RNA methylation modification is not only intimately interrelated with cancer development and progression but also actively influences immune cell infiltration in the tumor microenvironment (TME). RNA methylation modification genes influence the therapeutic progression of lung adenocarcinoma (LUAD), and mining RNA methylation modification prognosis-related markers in LUAD is crucial for its precise prognosis.

METHODS

RNA-Seq data and Gene sets were collected from online databases or published literature. Genomic variation analysis was conducted by the Maftools package. RNA methylation-immune-related lncRNAs were obtained by Pearson correlation analysis. Then, Consistent clustering analysis was performed to obtain RNA methylation modification- immune molecular subtypes (RMM-I Molecular subtypes) in LUAD based on selected lncRNAs. COX and random survival forest analysis were carried out to construct the RMM-I Score. The receiver operating characteristic (ROC) curve and Kaplan Meier survival analysis were used to assess survival differences. Tumor immune microenvironment was assessed through related gene signatures and CIBERSORT algorithm. In addition, drug sensitivity analysis was executed by the pRRophetic package.

RESULTS

Four RNA methylation modified-immune molecular subtypes (RMM-I1, RMM- I2, RMM-I3, RMM-I4) were presented in LUAD. Patients in RMM-I4 exhibited excellent survival advantages and immune activity. HAVCR2, CD274, and CTLA-4 expression were activated in RMM-I4, which might be heat tumors and a potential beneficial group for immunotherapy. OGFRP1, LINC01116, DLGAP1-AS2, CRNDE, LINC01137, MIR210HG, and CYP1B1-AS1 comprised the RMM-I Score. The RMM-I Score exhibited excellent accuracy in the prognostic assessment of LUAD, as patients with a low RMM- I Score exhibited remarkable survival advantage. Patients with a low RMM-I score might be more sensitive to treatment with Docetaxel, Vinorelbine, Paclitaxel, Cisplatin, and immunotherapy.

CONCLUSION

The RMM-I molecular subtype constituted the novel molecular characteristic subtype of LUAD, which complemented the existing pathological typing. More refined and accurate molecular subtypes provide help to reveal the mechanism of LUAD development. In addition, the RMM-I score offers a reliable tool for accurate prognosis of LUAD.

A methylation-related signature for predicting prognosis and sensitivity to first-line therapies in gastric cancer

Background

Methylation modification patterns play a crucial role in human cancer progression, especially in gastrointestinal cancers. We aimed to use methylation regulators to classify patients with gastric adenocarcinoma and build a model to predict prognosis, promoting the application of precision medicine.

Methods

We obtained RNA sequencing data and clinical data from The Cancer Genome Atlas (TCGA) database (n=335) and Gene Expression Omnibus (GEO) database (n=865). Unsupervised consensus clustering was used to identify subtypes of gastric adenocarcinoma. We performed functional enrichment analysis, immune infiltration analysis, drug sensitivity analysis, and molecular feature analysis to determine the clinical application for different subtypes. The univariate Cox regression analysis and the LASSO regression analysis were subsequently used to identify prognosis-related methylation regulators and construct a risk model.

Results

Through unsupervised consensus clustering, patients were divided into two subtypes (cluster A and cluster B) with different clinical outcomes. Cluster B included patients with a better prognosis outcome and who were more likely to respond to immunotherapy. We then successfully built a predictive model and found five methylation-related genes (CHAF1A, CPNE8, PHLDA3, SPARC, and EHF) potentially significant to the prognosis of patients. The 1-, 3-, and 5-year areas under the curve of the risk model were 0.712, 0.696, and 0.759, respectively. The risk score was an independent prognostic factor and had the highest concordance index among common clinical indicators. Meanwhile, the tumor microenvironment, sensitivity of chemotherapeutic drugs, molecular features, and oncogenic dedifferentiation differed significantly across the risk groups and subtypes.

Conclusions

We classified patients with gastric adenocarcinoma based on methylation regulators, which has positive implications for first-line clinical treatment. The prognostic model could predict the prognosis of patients and help to promote the development of precision medicine.

A novel m6A/m5C/m1A score signature to evaluate prognosis and its immunotherapy value in colon cancer patients

BACKGROUND

Colon cancer features strong heterogeneity and invasiveness, with high incidence and mortality rates. Recently, RNA modifications involving m6A, m5C, and m1A play a vital part in tumorigenesis and immune cell infiltration. However, integrated analysis among various RNA modifications in colon cancer has not been performed.

METHODS

RNA-seq profiling, clinical data and mutation data were obtained from The Cancer Genome Atlas and Gene Expression Omnibus. We first explored the mutation status and expression levels of m6A/m5C/m1A regulators in colon cancer. Then, different m6A/m5C/m1A clusters and gene clusters were identified by consensus clustering analysis. We further constructed and validated a scoring system, which could be utilized to accurately assess the risk of individuals and guide personalized immunotherapy. Finally, m6A/m5C/m1A regulators were validated by immunohistochemical staining and RT-qPCR.

RESULTS

In our study, three m6A/m5C/m1A clusters and gene clusters were identified. Most importantly, we constructed a m6A/m5C/m1A scoring system to assess the clinical risk of the individuals. Besides, the prognostic value of the score was validated with three independent cohorts. Moreover, the level of the immunophenoscore of the low m6A/m5C/m1A score group increased significantly with CTLA-4/PD-1 immunotherapy. Finally, we validated that the mRNA and protein expression of VIRMA and DNMT3B increased in colon cancer tissues.

CONCLUSIONS

We constructed and validated a stable and powerful m6A/m5C/m1A score signature to assess the survival outcomes and immune infiltration characteristics of colon cancer patients, which further guides optimization of personalized treatment, making it valuable for clinical translation and implementation.

Identification of N1 methyladenosine-related biomarker predicting overall survival outcomes and experimental verification in ovarian cancer

AIM

This study aimed to construct a N1-methyladenosine (m1A)-related biomarker model for predicting the prognosis of ovarian cancer (OVCA).

METHODS

OVCA samples were clustered into two subtypes using the Non-Negative Matrix Factorization (NMF) algorithm, including TCGA (n = 374) as the training set and GSE26712 (n = 185) as the external validation set. Hub genes, which were screened to construct a risk model, and nomogram to predict the overall survival of OVCA were explored and validated through various bioinformatic analysis and quantitative real-time PCR.

RESULTS

Following bootstrap correction, the C-index of nomogram was 0.62515, showing reliable performance. The functions of DEGs in the high- and low-risk groups were mainly enriched in immune response, immune regulation, and immune-related diseases. The immune cells relevant to the expression of hub genes were explored, for example, Natural Killer (NK) cells, T cells, activated dendritic cells (aDC).

CONCLUSIONS

AADAC, CD38, CACNA1C, and ATP1A3 might be used as m1A-related biomarkers for OVCA, and the nomogram labeled with m1A for the first time had excellent performance for predicting overall survival in OVCA.

Identification and validation of a novel prognostic model for gastric cancer based on m7G-related genes

Background

The role of N7-methyladenosine (m7G)-related genes in the progression and prognosis of gastric cancer (GC) remains unclear. This study aimed to explore prognostic biomarkers for GC based on m7G methylation regulators and to construct a prognostic risk model.

Methods

RNA sequencing profiles with corresponding clinicopathological information associated with GC of which the histological type was stomach adenocarcinoma (STAD) were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), respectively. A total of 29 m7G regulators were extracted from previous studies. According to the expression similarity of m7G regulators, the GC samples obtained from TCGA were further classified into 2 clusters demonstrating different overall survival (OS) rates and genetic heterogeneity, and the differentially expressed genes (DEGs) between these 2 clusters were defined as m7G-related genes. Univariate regression analysis and regression analysis were then used to obtain the prognostic m7G-related genes. The samples in TCGA and Genotype-Tissue Expression (GTEx) were used to verify the differential expression and prognostic value of these m7G-related genes contained in the prognostic model. Subsequently, the risk score was combined with other prognostic factors to develop a nomogram. The predictive ability of the nomogram was evaluated by the standard receiver operating characteristic (ROC) curve. Gene set enrichment analysis (GSEA) was used to identify activation pathways in both groups. Finally, the association between the prognostic model and the immune characteristics of GC were appraised.

Results

A prognostic model consisting of 11 m7G-related genes was constructed. GC patients in the high-risk group were shown to have a poor prognosis and this result was further demonstrated in each group. The risk model can be applied for patients with different clinical features. The results of GSEA showed that cell adhesion, cell junction, and focal adhesion were highly enriched in the high-risk group. In addition, we found that the expression of programmed cell death ligand 1 (PD-L1) was significantly elevated in the low-risk group, whereas programmed cell death ligand 2 (PD-L2) and tumor necrosis factor receptor superfamily member 4 (TNFRSF4) were overexpressed in the high-risk group.

Conclusions

We successfully built and verified a m7G relevant prognostic model for predicting prognosis and providing a new train of thought for improving the treatment of GC.

Identification and development of a novel risk model based on cuproptosis-associated RNA methylation regulators for predicting prognosis and characterizing immune status in hepatocellular carcinoma

BACKGROUND

Cuproptosis, a novel cell death caused by excess copper, is quite obscure in hepatocellular carcinoma (HCC) and needs more investigation.

METHODS

RNA-seq and clinical data of HCC patients TCGA database were analyzed to establish a predictive model through LASSO Cox regression analysis. External dataset ICGC was used for the verification. GSEA and CIBERSORT were applied to investigate the molecular mechanisms and immune microenvironment of HCC. Cuproptosis induced by elesclomol was confirmed via various in vitro experiments. The expression of prognostic genes was verified in HCC tissues using qRT-PCR analysis.

RESULTS

Initially, 18 cuproptosis-associated RNA methylation regulators (CARMRs) were selected for prognostic analysis. A nine-gene signature was created by applying the LASSO Cox regression method. Survival and ROC assays were carried out to validate the model using TCGA and ICGC database. Moreover, there exhibited obvious differences in drug sensitivity in terms of common drugs. A higher tumor mutation burden was shown in the high-risk group. Additionally, significant discrepancies were found between the two groups in metabolic pathways and RNA processing via GSEA analysis. Meanwhile, CIBERSORT analysis indicated different infiltrating levels of various immune cells between the two groups. Elesclomol treatment caused a unique form of programmed cell death accompanied by loss of lipoylated mitochondrial proteins and Fe-S cluster protein. The results of qRT-PCR indicated that most prognostic genes were differentially expressed in the HCC tissues.

CONCLUSION

Overall, our predictive signature displayed potential value in the prediction of overall survival of HCC patients and might provide valuable clues for personalized therapies.

Identification and validation of RNA methylation-related alternative splicing gene signature for low-grade glioma to predict survival and immune landscapes

BACKGROUND

Low-grade glioma (LGG) is a crucial pathological type of glioma. The present study aimed to explore multiple RNA methylation regulator-related AS events and investigate their prognostic values in LGG.

METHODS

The prognostic model for low-grade glioma was established using the LASSO regression analysis. To validate prognostic value, we performed Kaplan-Maier survival analysis, ROC curves and nomograms. The ESTIMATE algorithm, the CIBERSORT algorithm and the ssGSEA algorithm were utilized to explore the role of the immune microenvironment in LGG. Subsequently, we then used GO, KEGG and GSEA enrichment analysis to explore the functional roles of these genes. In addition, we employed the GDSC database to screen potential chemotherapeutic agents.

RESULTS

Eight RNA methylation related AS events were involved in construct a survival and prognosis model, which had good ability of independent prediction for patients with LGG. Patients in the high-risk group had shorter life expectancy and higher mortality, while patients in the low-risk group had a better prognosis. We constructed a nomogram which showed an excellent predictive performance for individual OS. The risk score exhibited a close correlation with some immune cells and expression of immune checkpoints. Patients in high-risk group were characterized by immunosuppressive microenvironment and poor response to immunotherapy, and were sensitive to more chemotherapeutic drugs. Pathway and functional enrichment analyses further confirmed that significant differences existed in immune landscape between the two subgroups.

CONCLUSION

The prognostic RNA methylation-related alternative splicing signature constructed could constitute a promising prognostic biomarker, which could serve to optimize treatment regimens.

A telomerase regulation-related lncRNA signature predicts prognosis and immunotherapy response for gastric cancer

BACKGROUND

Telomeres are involved in the development and progression of gastric cancer (GC). However, the association of telomerase regulation-related lncRNAs with prognosis and immunotherapy responsiveness in gastric cancer is unclear.

METHODS

This study systematically evaluated the relationship between lncRNAs co-expressed with 67 telomerase regulatory genes and gastric cancer prognosis. The risk scores of the samples were calculated based on telomerase regulation-related lncRNAs with prognostic value, and the samples were classified into high-/low-risk groups. The prognostic value of risk groups was then evaluated, a GC prognostic prediction model based on risk groups and clinical characteristics was established, and the prediction accuracy of the model was clarified by receiving operating characteristic (ROC) curves and calibration curves. Finally, the value of risk grouping in GC immunotherapy sensitivity was predicted by comparing MSI status and tumor mutation load between the high- and low-risk groups.

RESULTS

We identified 13 lncRNAs with prognostic value co-expressed with telomerase regulatory genes and observed that the prognosis of the low-risk group was significantly better than that of the high-risk group. Meanwhile, a GC overall survival (OS) prediction model based on risk grouping and clinical characteristics was developed, and ROC curves and calibration curves confirmed the good predictive ability of the model. In addition, the low-risk group exhibited a higher tumor mutation load and MSI-H, suggesting a possible benefit of immunotherapy.

CONCLUSION

We found that telomerase regulation-related lncRNAs have prognostic value in GC patients and contribute to the exploration of more effective immunotherapeutic strategies.

m6A/ m1A /m5C/m7G-related methylation modification patterns and immune characterization in prostate cancer

Methylation has a close relationship with immune reactions, metastasis, and cancer cell growth. Additionally, RNA methylation-related proteins have emerged as potential cancer therapeutic targets. The connection between the tumor microenvironment (TME) and methylation-related genes (MRGs) remains unclear. We explored the expression patterns of the MRGs in the genome and transcriptional fields of 796 prostate cancer (PCa) samples using two separate data sets. We identified a relationship between patient clinicopathological characteristics, prognosis, TME cell infiltrating qualities, and different MRG changes, as well as the identification of two distinct molecular groupings. Then, we formed an MRGs model to predict overall survival (OS), and we tested the accuracy of the model in patients with PCa. In addition, we developed a very accurate nomogram to improve the MRG model’s clinical applicability. The low-risk group had fewer tumor mutational burden (TMB), greater tumor immune dysfunction and exclusion (TIDE) ratings, fewer mutant genes, and better OS prospects. We discuss how MGRs may affect the prognosis, clinically important traits, TME, and immunotherapy responsiveness in PCa. In order to get a better understanding of MRGs in PCa, we could further explore the prognosis and create more effective immunotherapy regimens to open new avenues.

RNA Modification in Inflammatory Bowel Diseases

Inflammatory bowel disease (IBD) is a chronic inflammatory disorder characterized by damage to the intestinal mucosa, which is caused by a combination of factors. These include genetic and epigenetic alterations, environmental influence, microorganism interactions, and immune conditions. Some populations with IBD show a cancer-prone phenotype. Recent studies have provided insight into the involvement of RNA modifications in the specific pathogenesis of IBD through regulation of RNA biology in epithelial and immune cells. Studies of several RNA modification-targeting reagents have shown preferable outcomes in patients with colitis. Here, we note a new awareness of RNA modification in the targeting of IBD and related diseases, which will contribute to early diagnosis, disease monitoring, and possible control by innovative therapeutic approaches.