Construction of N-7 methylguanine-related mRNA prognostic model in uterine corpus endometrial carcinoma based on multi-omics data and immune-related analysis

Abstract

AbstractN-7 methylguanine (m7G) is one of the most common RNA base modifications in post-transcriptional regulation, which participates in multiple processes such as transcription, mRNA splicing and translation during the mRNA life cycle. However, its expression and prognostic value in uterine corpus endometrial carcinoma (UCEC) have not been systematically studied. In this paper, the data such as gene expression profiles, clinical data of UCEC patients, somatic mutations and copy number variants (CNVs) are obtained from the cancer genome atlas (TCGA) and UCSC Xena. By analyzing the expression differences of m7G-related mRNA in UCEC and plotting the correlation network maps, a risk score model composed of four m7G-related mRNAs (NSUN2, NUDT3, LARP1 and NCBP3) is constructed using least absolute shrinkage and selection operator (LASSO), univariate and multivariate Cox regression in order to identify prognosis and immune response. The correlation of clinical prognosis is analyzed between the m7G-related mRNA and UCEC via Kaplan–Meier method, receiver operating characteristic (ROC) curve, principal component analysis (PCA), t-SNE, decision curve analysis (DCA) curve and nomogram etc. It is concluded that the high risk is significantly correlated with (P < 0.001) the poorer overall survival (OS) in patients with UCEC. It is one of the independent risk factors affecting the OS. Differentially expressed genes are identified by R software in the high and low risk groups. The functional analysis and pathway enrichment analysis have been performed. Single sample gene set enrichment analysis (ssGSEA), immune checkpoints, m6A-related genes, tumor mutation burden (TMB), stem cell correlation, tumor immune dysfunction and rejection (TIDE) scores and drug sensitivity are also used to study the risk model. In addition, we have obtained 3 genotypes based on consensus clustering, which are significantly related to (P < 0.001) the OS and progression-free survival (PFS). The deconvolution algorithm (CIBERSORT) is applied to calculate the proportion of 22 tumor infiltrating immune cells (TIC) in UCEC patients and the estimation algorithm (ESTIMATE) is applied to work out the number of immune and matrix components. In summary, m7G-related mRNA may become a potential biomarker for UCEC prognosis, which may promote UCEC occurrence and development by regulating cell cycles and immune cell infiltration. It is expected to become a potential therapeutic target of UECE.