Integrated analysis of N6-methyladenosine- and 5-methylcytosine-related long non-coding RNAs for predicting prognosis in cervical cancer

N6-methyladenosine (m6A) and 5-methylcytosine (m5C) play a role in modifying long non-coding RNAs (lncRNAs) implicated in tumorigenesis and progression. This study was performed to evaluate prognostic value of m6A- and m5C-related lncRNAs and develop an efficient model for prognosis prediction in cervical cancer (CC). Using gene expression data of TCGA set, we identified m6A- and m5C-related lncRNAs. Consensus Clustering Analysis was performed for samples subtyping based on survival-related lncRNAs, followed by analyzing tumor infiltrating immune cells (TIICs). Optimal signature lncRNAs were obtained using lasso Cox regression analysis for constructing a prognostic model and a nomogram to predict prognosis. We built a co-expression network of 23 m6A-related genes, 15 m5C-related genes, and 62 lncRNAs. Based on 9 m6A- and m5C-related lncRNAs significantly associated with overall survival (OS) time, two molecular subtypes were obtained, which had significantly different OS time and fractions of TIICs. A prognostic model based on six m6A- and m5C-related signature lncRNAs was constructed, which could dichotomize patients into two risk subgroups with significantly different OS time. Prognostic power of the model was successfully validated in an independent dataset. We subsequently constructed a nomogram which could accurately predict survival probabilities. Drug sensitivity analysis found preferred chemotherapeutic agents for high and low-risk patients, respectively. Our study reveals that m6A- and m5C-related lncRNAs are associated with prognosis and immune microenvironment of CC. The m6A- and m5C-related six-lncRNA signature may be a useful tool for survival stratification in CC and open new avenues for individualized therapies.