Correspondence: Yurong Sun (supiaoy66@163.com) Background: Esophageal cancer is one of the most poorly diagnosed and fatal cancers in the world. Although a series of studies on esophageal cancer have been reported, the molecular pathogenesis of the disease is still elusive. Aim: To investigate the molecular process of esophageal cancer comprehensively and deeply. Methods: Differential expression analysis was performed to identify differentially expressed genes (DEGs) in different stages of esophageal cancer. Then exacting gene interaction modules and hub genes were identified in module interaction network. Further, though survival analysis, methylation analysis, pivot analysis, and enrichment analysis, some important molecules and related function or pathway were identified to elucidate potential mechanism in esophageal cancer. Results: A total of 7457 DEGs and 14 gene interaction modules were identified. These module genes were significantly involved in the positive regulation of protein transport, gastric acid secretion, insulin-like growth factor receptor binding and other biological processes (BPs), as well as p53 signaling pathway, ERBB signaling pathway and epidermal growth factor receptor (EGFR) signaling pathway. Then, transcription factors (TFs) (including HIF1A) and ncRNAs (including CRNDE and hsa-mir-330-3p) significantly regulate dysfunction modules were identified. Further, survival analysis showed that GNGT2 was closely related to survival of esophageal cancer. And DEGs with strong methylation regulation ability were identified, including SST and SH3GL2. Conclusion: These works not only help us to reveal the potential regulatory factors in the development of disease, but also deepen our understanding of its deterioration mechanism.