Identification of key genes and pathways for esophageal squamous cell carcinoma by bioinformatics analysis

Unveiling the potential of SLURP1 protein as a biomarker for prostate cancer screening

Background

Prostate cancer (PCa) develops slowly and lacks obvious symptoms in the early stage, which makes early screening and diagnosis difficult. Urine collection is simple and is an ideal source of biomarkers. In this study, we performed urinary proteomic studies in PCa patients to screen proteins and apply them to the non-invasive early diagnosis of PCa.

Method

Urine samples from PCa patients, benign prostatic hyperplasia (BPH) patients and normal control group were collected. Mass spectrometry was used for proteomic analysis and screening target proteins. Western blot and enzyme-linked immunosorbent assay (ELISA) were used to verify the results. Correlations with clinical indicators were explored, and receiver operating characteristic (ROC) curves were drawn to evaluate the value of target proteins in PCa.

Result

A total of 1065 proteins were identified. Urinary SLURP1 protein was significantly elevated in patients with PCa compared with normal controls and patients with BPH patients. Western blot and ELISA further verified the expression changes of SLURP1. The immunohistochemical staining results revealed a substantial increase in positive SLURP1 expression within PCa tumor tissue. Correlation analysis showed a positive correlation between the expression level of urine SLURP1 protein and serum PSA. ROC curve analysis of the SLURP1 protein in the urine of both normal individuals and PCa patients is determined to be 0.853 (95% CI=0.754 to 0.954).

Conclusion

The concentration of SLURP1 protein in urine of PCa patients is increased, which can serve as a biomarker for screening PCa.

SUMF1 overexpression promotes tumorous cell growth and migration and is correlated with the immune status of patients with glioma

BACKGROUND

Glioma is a prevalent type of malignant tumor. To date, there is a lack of literature reports that have examined the association between sulfatase modifying factor 1 (SUMF1) and glioma.

METHODS

The levels of SUMF1 were examined, and their relationships with the diagnosis, prognosis, and immune microenvironment of patients with glioma were investigated. Cox and Lasso regression analysis were employed to construct nomograms and risk models associated with SUMF1. The functions and mechanisms of SUMF1 were explored and verified using gene ontology, cell counting kit-8, wound healing, western blotting, and transwell experiments.

RESULTS

SUMF1 expression tended to increase in glioma tissues. SUMF1 overexpression was linked to the diagnosis of cancer, survival events, isocitrate dehydrogenase status, age, and histological subtype and was positively correlated with poor prognosis in patients with glioma. SUMF1 overexpression was an independent risk factor for poor prognosis. SUMF1-related nomograms and high-risk scores could predict the outcome of patients with glioma. SUMF1 co-expressed genes were involved in cytokine, T-cell activation, and lymphocyte proliferation. Inhibiting the expression of SUMF1 could deter the proliferation, migration, and invasion of glioma cells through epithelial mesenchymal transition. SUMF1 overexpression was significantly associated with the stromal score, immune cells (such as macrophages, neutrophils, activated dendritic cells), estimate score, immune score, and the expression of the programmed cell death 1, cytotoxic T-lymphocyte associated protein 4, CD79A and other immune cell marker.

CONCLUSION

SUMF1 overexpression was found to be correlated with adverse prognosis, cancer detection, and immune status in patients with glioma. Inhibiting the expression of SUMF1 was observed to deter the proliferation, migration, and invasion of cancer cells. The nomograms and risk models associated with SUMF1 could predict the prognosis of patients with glioma.

Altered genome‑wide hydroxymethylation analysis for neoadjuvant chemoradiotherapy followed by surgery in esophageal cancer

Esophageal cancer has high incidence rate in China. Neoadjuvant chemoradiotherapy (nCRT) has become the standard treatment for esophageal squamous cell carcinoma (ESCC). However, there are few reliable epigenetic parameters for patients with ESCC undergoing neoadjuvant therapy. Genomic extract from tumor tissue was amplified and sequenced using the Illumina HiSeq4000 to quantify genes associated methylation or hydromethylation in 12 patients with ESCC undergoing nCRT. The genome-wide hydroxymethylation were analyzed by methylated and hydroxymethylated DNA immunoprecipitation sequencing by MACS2 software and UCSC RefSeq database. Abnormal DNA methylation was statistically different between nCRT-well (showed a pathological complete response to nCRT) and nCRT-poor (showed incomplete pathological response to nCRT) patients. Levels of ten-eleven translocation 1, 2 and 3 mRNA and protein were higher in tumor tissue in nCRT-well group patients than in nCRT-poor group patients. Illumina HiSeq 4000 sequencing identified 2925 hypo-differentially hydroxymethylated region (DhMRs) and 292 hyper-DhMRs in promoter between nCRT-well and nCRT-poor patients. Biological processes associated with hyper-DhMRs included 'snRNA processing', 'hormone-mediated signaling pathway' and 'cellular response'. Metabolic processes were associated with hypo-DhMRs. These data may explain the functional response to nCRT in patients with abnormal promoter of methylation gene-associated mRNA expression. The present results implied that hyper-DhMRs and hypo-DhMRs affect molecular pathways, such as hippo and Notch signaling pathways, highlighting epigenetic modifications associated with clinical response to nCRT in patients with esophageal cancer.

Overexpression of cyclin-dependent kinase 1 in esophageal squamous cell carcinoma and its clinical significance

Cyclin‐dependent kinase 1 (CDK1) plays a significant role in certain malignancies. However, it remains unclear whether CDK1 plays a role in esophageal squamous cell carcinoma (ESCC). The aim of this study was to analyze the expression and clinical value of CDK1 in ESCC. CDK1 protein in 151 ESCC tissues and 138 normal esophageal tissues was detected by immunohistochemistry. RNA‐seq of eight pairs of ESCC and adjacent esophageal specimens was performed to evaluate the levels of CDK1 mRNA. Microarray and external RNA‐seq data from 664 cases of ESCC and 1733 cases of control tissues were used to verify the difference in CDK1 expression between the two groups. A comprehensive analysis of all data was performed to evaluate the difference in CDK1 between ESCC tissues and control tissues. Further, functional enrichment analyses were performed based on differentially expressed genes (DEGs) of ESCC and co‐expressed genes (CEGs) of CDK1. In addition, a lncRNA‐miRNA‐CDK1 network was constructed. The expression of CDK1 protein was obviously increased in ESCC tissues (3.540 ± 2.923 vs. 1.040 ± 1.632, P < 0.001). RNA‐seq indicated that the mRNA level of CDK1 was also highly expressed in ESCC tissues (5.261 ± 0.703 vs. 2.229 ± 1.161, P < 0.0001). Comprehensive analysis revealed consistent up‐regulation of CDK1 (SMD = 1.41; 95% CI 1.00–1.83). Further, functional enrichment analyses revealed that the functions of these genes were mainly concentrated in the cell cycle. A triple regulatory network of PVT1‐hsa‐miR‐145‐5p/hsa‐miR‐30c‐5p‐CDK1 was constructed using in silico analysis. In summary, overexpression of CDK1 is closely related to ESCC tumorigenesis.

Oncogenic enhancers drive esophageal squamous cell carcinogenesis and metastasis

The role of cis-elements and their aberrations remains unclear in esophageal squamous cell carcinoma (ESCC, further abbreviated EC). Here we survey 28 H3K27ac-marked active enhancer profiles and 50 transcriptomes in primary EC, metastatic lymph node cancer (LNC), and adjacent normal (Nor) esophageal tissues. Thousands of gained or lost enhancers and hundreds of altered putative super-enhancers are identified in EC and LNC samples respectively relative to Nor, with a large number of common gained or lost enhancers. Moreover, these differential enhancers contribute to the transcriptomic aberrations in ECs and LNCs. We also reveal putative driver onco-transcription factors, depletion of which diminishes cell proliferation and migration. The administration of chemical inhibitors to suppress the predicted targets of gained super-enhances reveals HSP90AA1 and PDE4B as potential therapeutic targets for ESCC. Thus, our epigenomic profiling reveals a compendium of reprogrammed cis-regulatory elements during ESCC carcinogenesis and metastasis for uncovering promising targets for cancer treatment.

Time series expression pattern of key genes reveals the molecular process of esophageal cancer

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.

Crosstalk between MMP-13, CD44, and TWIST1 and its role in regulation of EMT in patients with esophageal squamous cell carcinoma

Matrix metalloproteinases (MMPs) play key roles in epithelial-mesenchymal transition (EMT) for the development of cancer cell invasion and metastasis. MMP-13 is an extracellular matrix (ECM)-degrading enzyme that plays crucial roles in angiogenesis, cell cycle regulation, niche maintenance, and transforming squamous epithelial cells in various tissues. CD44, a transmembrane glycoprotein expressed on esophageal tumor cells, is required for EMT induction and invasion in esophageal squamous cell carcinoma (ESCC). The transcription factor TWIST1, as EMT and stemness marker, regulates MMPs expression and is identified as the downstream target of CD44. In this study, we aimed to investigate the probable interplay between the expression of key genes contributing to ESCC development, including MMP-13, TWIST1, and CD44 with clinical features for introducing novel diagnostic and therapeutic targets in the disease. The gene expression profiling of MMP-13, TWIST1, and CD44 was performed using quantitative real-time PCR in tumor tissues from 50 ESCC patients compared to corresponding margin non-tumoral tissues. Significant overexpression of MMP-13, CD44S, CD44V3, CD44V6, and TWIST1 were observed in 74%, 36%, 44%, 44%, and 52% of ESCC tumor samples, respectively. Overexpression of MMP-13 was associated with stage of tumor progression, metastasis, and tumor location (P < 0.05). There was a significant correlation between TWIST1 overexpression and grade (P < 0.05). Furthermore, overexpression of CD44 variants was associated with stage of tumor progression, grade, tumor invasion, and location (P < 0.05). The results indicated the significant correlation between concomitant expression of MMP-13/TWIST1, TWIST1/CD44, and CD44/MMP-13 with each other in a variety of clinicopathological traits, including depth of tumor invasion, tumor location, stage of tumor, and lymph node involvement in ESCC tissue samples (P < 0.05). Collectively, our results indicate that the TWIST1-CD44-MMP-13 axis is involved in tumor aggressiveness, proposing these genes as regulators of EMT, diagnostic markers, and therapeutic targets in ESCC.

Cellular Processes Involved in Jurkat Cells Exposed to Nanosecond Pulsed Electric Field

Recently, nanosecond pulsed electric field (nsPEF) has been considered as a new tool for tumor therapy, but its molecular mechanism of function remains to be fully elucidated. Here, we explored the cellular processes of Jurkat cells exposed to nanosecond pulsed electric field. Differentially expressed genes (DEGs) were acquired from the GEO2R, followed by analysis with a series of bioinformatics tools. Subsequently, 3D protein models of hub genes were modeled by Modeller 9.21 and Rosetta 3.9. Then, a 100 ns molecular dynamics simulation for each hub protein was performed with GROMACS 2018.2. Finally, three kinds of nsPEF voltages (0.01, 0.05, and 0.5 mV/mm) were used to simulate the molecular dynamics of hub proteins for 100 ns. A total of 1769 DEGs and eight hub genes were obtained. Molecular dynamic analysis, including root mean square deviation (RMSD), root mean square fluctuation (RMSF), and the Rg, demonstrated that the 3D structure of hub proteins was built, and the structural characteristics of hub proteins under different nsPEFs were acquired. In conclusion, we explored the effect of nsPEF on Jurkat cell signaling pathway from the perspective of molecular informatics, which will be helpful in understanding the complex effects of nsPEF on acute T-cell leukemia Jurkat cells.