Abstract
BACKGROUND: Cholangiocarcinoma (CHOL) is a malignant tumor that originates in the extrahepatic bile duct and can extend from the hilar region to the lower end of the common bile duct. The prognosis of CHOL patients is particularly poor; therefore, in this study, we screened mRNAs correlated with N6-methyladenosine (m(6)A) to construct a risk model for prognosis in CHOL. METHODS: The TCGA-CHOL dataset was applied to obtain and analyze the coexpression of 1281 m(6)A-related mRNAs, from which 14 were selected for further analysis through univariate proportional hazards (cox) regression analysis. Aryl hydrocarbon receptor interacting protein (AIP), CCAAT/enhancer binding protein beta (CEBPB), syndecan1 (SDC1), vacuolar protein sorting 25 homolog (VPS25) and syntaxin binding protein 2 (STXBP2) were then screened out through the least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis to develop a precise m(6)A-related mRNA prognosis risk model (MRMRPM) with an area under curve (AUC) of 0.908 and 0.923 after 1 and 2 years, respectively. We divided the samples into high-risk and low-risk groups using the m(6)A-related mRNA prognosis risk model. RESULTS: Kaplan-Meier analysis indicated poor overall survival (OS) for the high-risk group. Two Gene Expression Omnibus (GEO) datasets (GSE89748 and GSE107943) were used to validate the risk model. The results of drug sensitivity and immune cell infiltration analysis showed that the risk model could serve as a prognosis index of potential immunotherapeutic characteristics and drug sensitivity. Furthermore, the proportion of resting dendritic cells and regulatory T cells was positively associated with an increased expression of four m(6)A-related mRNAs - AIP, CEBPB, SDC1, and VPS25 - in the high-risk CHOL group. CONCLUSIONS: Our findings suggest that this model can be a prognostic indicator for CHOL patients.