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

Emerging role of RNA modification and long noncoding RNA interaction in cancer

RNA modification, especially N6-methyladenosine, 5-methylcytosine, and N7-methylguanosine methylation, participates in the occurrence and progression of cancer through multiple pathways. The function and expression of these epigenetic regulators have gradually become a hot topic in cancer research. Mutation and regulation of noncoding RNA, especially lncRNA, play a major role in cancer. Generally, lncRNAs exert tumor-suppressive or oncogenic functions and its dysregulation can promote tumor occurrence and metastasis. In this review, we summarize N6-methyladenosine, 5-methylcytosine, and N7-methylguanosine modifications in lncRNAs. Furthermore, we discuss the relationship between epigenetic RNA modification and lncRNA interaction and cancer progression in various cancers. Therefore, this review gives a comprehensive understanding of the mechanisms by which RNA modification affects the progression of various cancers by regulating lncRNAs, which may shed new light on cancer research and provide new insights into cancer therapy.

Construction and validation of a prognostic signature for mucinous colonic adenocarcinoma based on N7-methylguanosine-related long non-coding RNAs

BACKGROUND

Mucinous colonic adenocarcinoma remains a challenging disease due to its high propensity for metastasis and recurrence. N7-methylguanosine (m7G) and long non-coding RNA (lncRNA) are closely associated with the occurrence and progression of tumors. However, research on m7G-related lncRNA in mucinous colonic adenocarcinoma is lacking. Therefore, we sought to explore the prognostic impact of m7G-related lncRNAs in mucinous adenocarcinoma (MC) patients.

METHODS

In this study, Pearson analysis was used to identify m7G-related lncRNAs from transcriptome data in The Cancer Genome Atlas (TCGA). Univariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression were used to further screen m7G-related lncRNAs and incorporate them into a prognostic signature. Based on the risk model, patients were divided into low- and high-risk groups and randomly assigned to the training set and test sets in a 6:4 ratio. Kaplan-Meier, receiver operating characteristic (ROC) curve, multivariate regression, and nomogram analyses were used to confirm the accuracy of the signature. The CIBERSORT algorithm was used to calculate the degree of immune cell infiltration (ICI). Finally, the correlation of the prognostic signature with tumor mutational burden (TMB) and immunophenotype score (IPS) was evaluated.

RESULTS

A total of 432 m7G-related lncRNAs were identified by Pearson analysis. Univariate Cox regression, LASSO regression and survival analysis were performed to further select six m7G-related lncRNAs (P<0.05): AC254629.1, LINC01133, LINC01134, MHENCR, SMIM2-AS1, and XACT. Based on the risk model, heat maps, Kaplan-Meier curves, and ROC curves were constructed, and the results showed that there were significant differences in expression levels and survival status between the two risk groups. The area under the ROC curve (AUC) values for 3-, 5-, and 10-year survival in the training set were 0.944, 0.957, and 1.000, respectively. And in the test set were 0.964, 1.000, and 1.000, respectively. Subsequently, univariate and multivariate regression analyses of clinical characteristics and risk score were performed. The results of risk score were [hazard ratio (HR): 6.458, 95% confidence interval (CI): 2.708-15.403, P<0.001; HR: 7.280, 95% CI: 2.500-21.203, P<0.001], respectively. Using the risk score as an independent prognostic factor, the AUC of it over 3, 5, and 10 years was 0.911, 0.955, and 0.961, respectively. Calibration plots for the nomogram show that the model calibration line is very close to the ideal calibration line, indicating good calibration. The level of ICI was significantly different in the different risk groups. Survival analysis showed that, regardless of TMB risk, patients with MC and a high-risk score consistently had a poor overall survival (OS).

CONCLUSIONS

The m7G-related lncRNA prognostic signature has potential value for the prognosis of mucinous colonic adenocarcinoma.

Selection of M7G-related lncRNAs in kidney renal clear cell carcinoma and their putative diagnostic and prognostic role

BACKGROUND

Kidney renal clear cell carcinoma (KIRC) is a common malignant tumor of the urinary system. This study aims to develop new biomarkers for KIRC and explore the impact of biomarkers on the immunotherapeutic efficacy for KIRC, providing a theoretical basis for the treatment of KIRC patients.

METHODS

Transcriptome data for KIRC was obtained from the The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Weighted gene co-expression network analysis identified KIRC-related modules of long noncoding RNAs (lncRNAs). Intersection analysis was performed differentially expressed lncRNAs between KIRC and normal control samples, and lncRNAs associated with N(7)-methylguanosine (m7G), resulting in differentially expressed m7G-associated lncRNAs in KIRC patients (DE-m7G-lncRNAs). Machine Learning was employed to select biomarkers for KIRC. The prognostic value of biomarkers and clinical features was evaluated using Kaplan-Meier (K-M) survival analysis, univariate and multivariate Cox regression analysis. A nomogram was constructed based on biomarkers and clinical features, and its efficacy was evaluated using calibration curves and decision curves. Functional enrichment analysis was performed to investigate the functional enrichment of biomarkers. Correlation analysis was conducted to explore the relationship between biomarkers and immune cell infiltration levels and common immune checkpoint in KIRC samples.

RESULTS

By intersecting 575 KIRC-related module lncRNAs, 1773 differentially expressed lncRNAs, and 62 m7G-related lncRNAs, we identified 42 DE-m7G-lncRNAs. Using XGBoost and Boruta algorithms, 8 biomarkers for KIRC were selected. Kaplan-Meier survival analysis showed significant survival differences in KIRC patients with high and low expression of the PTCSC3 and RP11-321G12.1. Univariate and multivariate Cox regression analyses showed that AP000696.2, PTCSC3 and clinical characteristics were independent prognostic factors for patients with KIRC. A nomogram based on these prognostic factors accurately predicted the prognosis of KIRC patients. The biomarkers showed associations with clinical features of KIRC patients, mainly localized in the cytoplasm and related to cytokine-mediated immune response. Furthermore, immune feature analysis demonstrated a significant decrease in immune cell infiltration levels in KIRC samples compared to normal samples, with a negative correlation observed between the biomarkers and most differentially infiltrating immune cells and common immune checkpoints.

CONCLUSION

In summary, this study discovered eight prognostic biomarkers associated with KIRC patients. These biomarkers showed significant correlations with clinical features, immune cell infiltration, and immune checkpoint expression in KIRC patients, laying a theoretical foundation for the diagnosis and treatment of KIRC.

Bioinformatics construction and experimental validation of a cuproptosis-related lncRNA prognostic model in lung adenocarcinoma for immunotherapy response prediction

Cuproptosis is a newly form of cell death. Cuproptosis related lncRNA in lung adenocarcinoma (LUAD) has also not been fully elucidated. In the present study, we aimed to construct a prognostic signature based on cuproptosis-related lncRNA in LUAD and investigate its association with immunotherapy response. The RNA-sequencing data, clinical information and simple nucleotide variation of LUAD patients were obtained from TCGA database. The LASSO Cox regression was used to construct a prognostic signature. The CIBERSORT, ESTIMATE and ssGSEA algorithms were applied to assess the association between risk score and TME. TIDE score was applied to reflect the efficiency of immunotherapy response. The influence of overexpression of lncRNA TMPO-AS1 on A549 cell was also assessed by in vitro experiments. The lncRNA prognostic signature included AL606834.1, AL138778.1, AP000302.1, AC007384.1, AL161431.1, TMPO-AS1 and KIAA1671-AS1. Low-risk group exhibited much higher immune score, stromal score and ESTIMATE score, but lower tumor purity compared with high-risk groups. Also, low-risk group was associated with a much higher score of immune cells and immune related function sets, indicating an immune activation state. Low-risk patients had relative higher TIDE score and lower TMB. External validation using IMvigor210 immunotherapy cohort demonstrated that low-risk group had a better prognosis and might more easily benefit from immunotherapy. Overexpression of lncRNA TMPO-AS1 promoted the proliferation, migration and invasion of A549 cell line. The novel cuproptosis-related lncRNA signature could predict the prognosis of LUAD patients, and helped clinicians stratify patients appropriate for immunotherapy and determine individual therapeutic strategies.

M7G-related molecular subtypes can predict the prognosis and correlate with immunotherapy and chemotherapy responses in bladder cancer patients

BACKGROUND

N7-methylguanosine (m7G) is closely associated with tumor prognosis and immune response in many cancer types. The correlation between m7G and bladder cancer (BC) needs further study. We aimed to orchestrate molecular subtypes and identify key genes for BC from the perspective of m7G.

METHODS

RNA-seq and clinical data of BC patients were extracted from TCGA and GSE13507 datasets. The patients were subtyped by "ConsensusClusterPlus" and "limma." The clusters were validated by the Kaplan‒Meier curves, univariable and multivariate Cox regression models, the concordance index, and calibration curves. The immunotherapy response was evaluated by immune checkpoints, immune infiltration, TIDE score, and IMvigor210 cohort. Genomics of Drug Sensitivity in Cancer was utilized to predict the chemotherapy response between the clusters.

RESULTS

The m7G-related cluster was ultimately established by EIF4G1, NUDT11, NUDT10, and CCNB1. The independent prognostic value of the m7G-related cluster was validated by the TCGA and GSE13507 datasets. The cluster was involved in immune-associated pathways, such as neutrophil degranulation, antigen processing cross-presentation, and signaling by interleukins pathways. Meanwhile, cluster 2 was positively correlated with many immune checkpoints, such as CD274, CTLA4, HAVCR2, LAG3, PDCD1, and PDCD1LG2. The cluster 2 was significantly correlated with a higher TIDE score than the cluster 1. Furthermore, in the IMvigor210 cohort, patients in the cluster 1 had a higher response rate than those in the cluster 2. Patients in the cluster 2 were sensitive to many chemotherapies.

CONCLUSIONS

We successfully determined molecular subtypes and identified key genes for BC from the perspective of m7G, thereby providing a roadmap for the evolution of immunotherapy and precision medicine.

Analysis of m7G methylation modification patterns and pulmonary vascular immune microenvironment in pulmonary arterial hypertension

Background

M⁷G methylation modification plays an important role in cardiovascular disease development. Dysregulation of the immune microenvironment is closely related to the pathogenesis of PAH. However, it is unclear whether m⁷G methylation is involved in the progress of PAH by affecting the immune microenvironment.

Methods

The gene expression profile of PAH was obtained from the GEO database, and the m⁷G regulatory factors were analyzed for differences. Machine learning algorithms were used to screen characteristic genes, including the least absolute shrinkage and selection operator, random forest, and support vector machine recursive feature elimination analysis. Constructed a nomogram model, and receiver operating characteristic was used to evaluate the diagnosis of disease characteristic genes value. Next, we used an unsupervised clustering method to perform consistent clustering analysis on m⁷G differential genes. Used the ssGSEA algorithm to estimate the relationship between the m⁷G regulator in PAH and immune cell infiltration and analyze the correlation with disease-characteristic genes. Finally, the listed drugs were evaluated through the screened signature genes.

Results

We identified 15 kinds of m⁷G differential genes. CYFIP1, EIF4E, and IFIT5 were identified as signature genes by the machine learning algorithm. Meanwhile, two m⁷G molecular subtypes were identified by consensus clustering (cluster A/B). In addition, immune cell infiltration analysis showed that activated CD4 T cells, regulatory T cells, and type 2 T helper cells were upregulated in m⁷G cluster B, CD56 dim natural killer cells, MDSC, and monocyte were upregulated in the m⁷G cluster A. It might be helpful to select Calpain inhibitor I and Everolimus for the treatment of PAH.

Conclusion

Our study identified CYFIP1, EIF4E, and IFIT5 as novel diagnostic biomarkers in PAH. Furthermore, their association with immune cell infiltration may facilitate the development of immune therapy in PAH.

A Novel m7G-Related Gene Signature Predicts the Prognosis of Colon Cancer

Colon cancer (CC), one of the most common malignancies worldwide, lacks an effective prognostic prediction biomarker. N7-methylguanosine (m7G) methylation is a common RNA modification type and has been proven to influence tumorigenesis. However, the correlation between m7G-related genes and CC remains unclear. The gene expression levels and clinical information of CC patients were downloaded from public databases. Twenty-nine m7G-related genes were obtained from the published literature. Via unsupervised clustering based on the expression levels of m7G-related genes, CC patients were divided into three m7G clusters. Based on differentially expressed genes (DEGs) from the above three groups, CC patients were further divided into three gene clusters. The m7G score, a prognostic model, was established using principal component analysis (PCA) based on 15 prognosis-associated m7G genes. KM curve analysis demonstrated that the overall survival rate was remarkably higher in the high-m7G score group, which was much more significant in advanced CC patients as confirmed by subgroup analysis. Correlation analysis indicated that the m7G score was associated with tumor mutational burden (TMB), PD-L1 expression, immune infiltration, and drug sensitivity. The expression level of prognosis-related m7G genes was further confirmed in human CC cell lines and samples. This study established an m7G gene-based prognostic model (m7G score), which demonstrated the important roles of m7G-related genes during CC initiation and progression. The m7G score could be a practical biomarker to predict immunotherapy response and prognosis in CC patients.

Comprehensive pan-cancer analysis of N7-methylguanosine regulators: Expression features and potential implications in prognosis and immunotherapy

Although immunotherapy has made great strides in cancer therapy, its effectiveness varies widely among individual patients as well as tumor types, and there is an urgent need to develop biomarkers for effectively assessing immunotherapy response. In recent years, RNA methylation regulators have demonstrated to be novel potential biomarkers for prognosis as well as immunotherapy of cancers, such as N6-methyladenine (m6A) and 5-methylcytosine (m5C). N7-methylguanosine (m7G) is a prevalent RNA modification in eukaryotes, but the relationship between m7G regulators and prognosis as well as tumor immune microenvironment is still unclear. In this study, a pan-cancer analysis of 26 m7G regulators across 17 cancer types was conducted based on the bioinformatics approach. On the one hand, a comprehensive analysis of expression features, genetic variations and epigenetic regulation of m7G regulators was carried out, and we found that the expression tendency of m7G regulators were different among tumors and their aberrant expression in cancers could be affected by single nucleotide variation (SNV), copy number variation (CNV), DNA methylation and microRNA (miRNA) separately or simultaneously. On the other hand, the m7Gscore was modeled based on single sample gene set enrichment analysis (ssGSEA) for evaluating the relationships between m7G regulators and cancer clinical features, hallmark pathways, tumor immune microenvironment, immunotherapy response as well as pharmacotherapy sensitivity, and we illustrated that the m7Gscore exhibited tight correlations with prognosis, several immune features, immunotherapy response and drug sensitivity in most cancers. In conclusion, our pan-cancer analysis revealed that m7G regulators may exert critical roles in the tumor progression and immune microenvironment, and have the potential as biomarkers for predicting prognosis, immunotherapy response as well as candidate drug compounds for cancer patients.

N7-methylguanosine-related lncRNAs: Distinction between hot and cold tumors and construction of predictive models in colon adenocarcinoma

Colon adenocarcinoma (COAD) is a prevalent malignant tumor that severely threatens human health across the globe. Immunotherapy is an essential need for patients with COAD. N7-methylguanosine (m7G) has been associated with human diseases, and non-coding RNAs (lncRNAs) regulate various tumor-related biological processes. Nonetheless, the m7G-related lncRNAs involved in COAD regulation are limited. This study aims to construct the clustering features and prognostic model of m7G-related lncRNAs in COAD. First, The Cancer Genome Atlas (TCGA) database was used to identify m7G-related differentially expressed lncRNAs (DELs), based on which COAD cases could be classified into two subtypes. Subsequently, univariate Cox analysis was used to identify 9 prognostic m7G-related lncRNAs. Further, Five candidates were screened by LASSO-Cox regression to develop new models. The patients were divided into high-risk and low-risk groups based on the median risk score. Consequently, the Kaplan-Meier survival curve demonstrated a statistically significant overall survival (OS) between the high- and low-risk groups (P&lt;0.001). Multivariate Cox regression analysis revealed that risk score is an independent prognostic factor in COAD patients (P&lt;0.001). This confirms the clinical applicability of the model. Additionally, we performed Gene Set Enrichment Analysis (GSEA), which uncovered the biological and functional differences between risk subgroups, i.e., enrichment of immune-related diseases in the high-risk group and enrichment of metabolic-related pathways in the low-risk group. In a drug sensitivity analysis, high-risk group were more sensitive to some chemotherapeutics and targeted drugs than low-risk group. Eventually, the stability of the model was confirmed by qRT-PCR. Our study unraveled the features of different immune states of COAD and established a prognostic model, including five m7G-related lncRNAs for COAD patients. These results will bolster clinical treatment and survival prediction of COAD.