Untargeted metabolomics analysis of esophageal squamous cell cancer progression

90% of esophageal cancer are esophageal squamous cell carcinoma (ESCC) and ESCC has a very poor prognosis and high mortality. Nevertheless, the key metabolic pathways associated with ESCC progression haven't been revealed yet. Metabolomics has become a new platform for biomarker discovery over recent years. We aim to elucidate dominantly metabolic pathway in all ESCC tumor/node/metastasis (TNM) stages and adjacent cancerous tissues. We collected 60 postoperative esophageal tissues and 15 normal tissues adjacent to the tumor, then performed Liquid Chromatography with tandem mass spectrometry (LC-MS/MS) analyses. The metabolites data was analyzed with metabolites differential and correlational expression heatmap according to stage I vs. con., stage I vs. stage II, stage II vs. stage III, and stage III vs. stage IV respectively. Metabolic pathways were acquired by Kyoto Encyclopedia of Genes and Genomes. (KEGG) pathway database. The metabolic pathway related genes were obtained via Gene Set Enrichment Analysis (GSEA). mRNA expression of ESCC metabolic pathway genes was detected by two public datasets: gene expression data series (GSE)23400 and The Cancer Genome Atlas (TCGA). Receiver operating characteristic curve (ROC) analysis is applied to metabolic pathway genes. 712 metabolites were identified in total. Glycerophospholipid metabolism was significantly distinct in ESCC progression. 16 genes of 77 genes of glycerophospholipid metabolism mRNA expression has differential significance between ESCC and normal controls. Phosphatidylserine synthase 1 (PTDSS1) and Lysophosphatidylcholine Acyltransferase1 (LPCAT1) had a good diagnostic value with Area under the ROC Curve (AUC) > 0.9 using ROC analysis. In this study, we identified glycerophospholipid metabolism was associated with the ESCC tumorigenesis and progression. Glycerophospholipid metabolism could be a potential therapeutic target of ESCC progression.