Abstract
BACKGROUND: Neoadjuvant chemoradiotherapy (nCRT) followed by surgery benefits survival for patients with esophageal squamous cell carcinomas (ESCC) compared with surgery alone, but the clinical outcomes of nCRT are heterogeneous. This study aimed to elucidate transcriptional factor (TF) regulation network and competitive endogenous RNA (ceRNA) network determining response of ESCC to nCRT. MATERIALS AND METHODS: RNA microarray data of GSE59974 and GSE45670 were analyzed to investigate the significant changes of lincRNAs, miRNAs, mRNAs in responders and non-responders of nCRT in ESCC. Functional and enrichment analyses were conducted by clusterProfiler. The target lincRNAs and mRNAs of miRNAs were predicted by miRWalk. The ceRNA and TF regulatory networks were constructed using Cytoscape. RESULTS: Differentially expressed genes between responders and non-responders mainly enriched in biological process including Wnt signaling pathway and regulation of cell development and morphogenesis involved in differentiation. Besides, these genes showed enrichment in molecular function of glycosaminoglycan binding, metalloendopeptidase inhibitor and growth factor activity. KEGG analysis enriched these genes in pathways of neurotrophin signaling pathway, cell adhesion molecules and Wnt signaling pathway. We also constructed ceRNA network and TF network regulating response of ESCC to nCRT. Core regulatory miRNAs were miR-520a, miR-548am, miR-3184, miR-548d, miR-4725, miR-148a, miR-4659a and key regulatory TFs included MBNL1, SLC26A3, BMP4, ZIC1 and ANKRD7. CONCLUSION: We identified significantly altered lincRNAs, miRNAs and mRNAs involved in the nCRT response of ESCC. In addition, the ceRNA regulatory network of lincRNA-miRNA-mRNA and TF regulatory network were constructed, which would elucidate novel molecular mechanisms determining nCRT response of ESCC, thus providing promising clues for clinical therapy.