The DNA methylation profile of non-coding RNAs improves prognosis prediction for pancreatic adenocarcinoma

Jie Zhang,Keqing Shi,Weiguo Huang,Wanqing Weng,Zhongjing Zhang,Yangyang Guo,Tuo Deng,Yukai Xiang,Xiaofeng Ni,Bicheng Chen,Mengtao Zhou

DOI: https://doi.org/10.1186/s12935-019-0828-8

2019-04-23

Abstract:Background: Compelling lines of evidence indicate that DNA methylation of non-coding RNAs (ncRNAs) plays critical roles in various tumour progression. In addition, the differential methylation of ncRNAs can predict prognosis of patients. However, little is known about the clear relationship between DNA methylation profile of ncRNAs and the prognosis of pancreatic adenocarcinoma (PAC) patients. Methods: The data of DNA methylation, RNA-seq, miRNA-seq and clinical features of PAC patients were collected from TCGA database. The DNA methylation profile was obtained using the Infinium HumanMethylation450 BeadChip array. LASSO regression was performed to construct two methylation-based classifiers. The risk score of methylation-based classifiers was calculated for each patient, and the accuracy of the classifiers in predicting overall survival (OS) was examined by ROC curve analysis. In addition, Cox regression models were utilized to assess whether clinical variables and the classifiers were independent prognostic factors for OS. The targets of miRNA and the genes co-expressed with lncRNA were identified with DIANA microT-CDS and the Multi-Experiment Matrix (MEM), respectively. Moreover, DAVID Bioinformatics Resources were applied to analyse the functional enrichment of these targets and co-expressed genes. Results: A total of 4004 CpG sites of miRNA and 11,259 CpG sites of lncRNA were screened. Among these CpG sites, 8 CpG sites of miRNA and 7 CpG sites of lncRNA were found with regression coefficients. By multiplying the sum of methylation degrees of the selected CpGs with these coefficients, two methylation-based classifiers were constructed. The classifiers have shown good performance in predicting the survival rate of PAC patients at varying follow-up times. Interestingly, both of these two classifiers were predominant and independent factors for OS. Furthermore, functional enrichment analysis demonstrated that aberrantly methylated miRNAs and lncRNAs are related to calcium ion transmembrane transport and MAPK, Ras and calcium signalling pathways. Conclusion: In the present study, we identified two methylation-based classifiers of ncRNA associated with OS in PAC patients through a comprehensive analysis of miRNA and lncRNA profiles. We are the first group to demonstrate a relationship between the aberrant DNA methylation of ncRNAs and the prognosis of PAC, and this relationship would contribute to individualized PAC therapy.