Abstract:BACKGROUND Esophageal squamous cell carcinoma (ESCC) is one of the most common malignancies worldwide, and its development comprises a multistep process from intraepithelial neoplasia (IN) to carcinoma (CA). However, the critical regulators and underlying molecular mechanisms remain largely unknown. AIM To explore the genes and infiltrating immune cells in the microenvironment that are associated with the multistage progression of ESCC to facilitate diagnosis and early intervention. METHODS A mouse model mimicking the multistage development of ESCC was established by providing warter containing 4-nitroquinoline 1-oxide (4NQO) to C57BL/6 mice. Moreover, we established a control group without 4NQO treatment of mice. Then, transcriptome sequencing was performed for esophageal tissues from patients with different pathological statuses, including low-grade IN (LGIN), high-grade IN (HGIN), and CA, and controlled normal tissue (NOR) samples. Differentially expressed genes (DEGs) were identified in the LGIN, HGIN, and CA groups, and the biological functions of the DEGs were analyzed via Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses. The CIBERSORT algorithm was used to detect the pattern of immune cell infiltration. Immunohistochemistry (IHC) was also conducted to validate our results. Finally, the Luminex multiplex cytokine analysis was utilized to measure the serum cytokine levels in the mice. RESULTS Compared with those in the NOR group, a total of 681541, and 840 DEGs were obtained in the LGIN, HGIN, and CA groups, respectively. Using the intersection of the three sets of DEGs, we identified 86 genes as key genes involved in the development of ESCC. Enrichment analysis revealed that these genes were enriched mainly in the keratinization, epidermal cell differentiation, and interleukin (IL)-17 signaling pathways. CIBERSORT analysis revealed that, compared with those in the NOR group, M0 and M1 macrophages in the 4NQO group showed stronger infiltration, which was validated by IHC. Serum cytokine analysis revealed that, compared with those in the NOR group, IL-1β and IL-6 were upregulated, while IL-10 was downregulated in the LGIN, HGIN, and CA groups. Moreover, the expression of the representative key genes, such as S100a8 and Krt6b, was verified in external human samples, and the results of immunohistochemical staining were consistent with the findings in mice. CONCLUSION We identified a set of key genes represented by S100a8 and Krt6b and investigated their potential biological functions. In addition, we found that macrophage infiltration and abnormal alterations in the levels of inflammation-associated cytokines, such as IL-1β, IL-6, and IL-10, in the peripheral blood may be closely associated with the development of ESCC.

Comprehensive multi-omics analysis identified core molecular processes in esophageal cancer and revealed GNGT2 as a potential prognostic marker

Identification of Significant Biomarkers and Pathways Associated with Gastric Carcinogenesis by Whole Genome-Wide Expression Profiling Analysis

Identification of Key Genes and Evaluation of Clinical Outcomes in Lung Squamous Cell Carcinoma Using Integrated Bioinformatics Analysis

Identification of Key Genes and Signaling Pathways Associated with the Progression of Gastric Cancer

Comprehensive Analysis and Identification of Key Driver Genes for Distinguishing Between Esophageal Adenocarcinoma and Squamous Cell Carcinoma

Integrative genomics analysis of hub genes and their relationship with prognosis and signaling pathways in esophageal squamous cell carcinoma

Prognostic Biomarkers Identification in Esophageal Cancer Based on WGCNA and Single-Cell Analysis

A Genome‐wide Analysis of Long Noncoding RNA Profile Identifies Differentially Expressed Lncrnas Associated with Esophageal Cancer

Identification of serum prognostic marker miRNAs and construction of microRNA-mRNA networks of esophageal cancer

Transcriptome sequencing reveals novel biomarkers and immune cell infiltration in esophageal tumorigenesis

Identification of crucial genes correlated with esophageal cancer by integrated high-throughput data analysis

Identification of Key Genes Associated with Esophageal Adenocarcinoma Based on Bioinformatics Analysis

Microarray analyses reveal genes related to progression and prognosis of esophageal squamous cell carcinoma.

Identification and Validation of Potential Pathogenic Genes and Prognostic Markers in ESCC by Integrated Bioinformatics Analysis

Landscape of Expression Profiles in Esophageal Carcinoma by the Cancer Genome Atlas Data.

Identification of Key Modules and Hub Genes Involved in Esophageal Squamous Cell Carcinoma Tumorigenesis Using WCGNA

Integrated Bioinformatics Analysis to Identify Key Genes Related to the Prognosis of Esophageal Squamous Cell Carcinoma

Identification of Potential Biomarkers Associated with Prognosis in Gastric Cancer via Bioinformatics Analysis

Construction and investigation of β3GNT2-associated regulatory network in esophageal carcinoma

Analysis of Significant Genes and Pathways in Esophageal Cancer Based on Gene Expression Omnibus Database

Diagnostic model of combined ceRNA and DNA methylation related genes in esophageal carcinoma