Multi-Omics Analysis of the Tumor Microenvironment in Liver Metastasis of Colorectal Cancer Identified FJX1 as a Novel Biomarker

Predictive value of immediate early response 5 like (IER5L) in the prognosis and immune checkpoint blockade therapy of non-small cell lung cancer patients

BACKGROUND

Non-small cell lung cancer (NSCLC) is a malignancy with high mortality. Immediate early response 5 like (IER5L) has been found to associate with worse prognosis in colorectal cancer patients. However, its role in the prognosis prediction of NSCLC has remained largely unknown.

METHODS

The IER5L expression in NSCLC and normal tissues was analyzed in two public cohorts: TCGA-LUAD-LUSC and GSE159857. Additionally, functional enrichment, survival analysis, CIBERSORT and tumor mutation burden (TMB) were investigated between low- and high-IER5L level groups. The in vitro IER5L mRNA and protein levels were determined using RT-qPCR and western blot, respectively.

RESULTS

The data from TCGA-LUAD-LUSC and GSE159857 cohorts showed a high IER5L mRNA expression in NSCLC tissue samples compared to normal controls. The increased expression of IER5L in NSCLC cells were also validated by RT-qPCR and western blot analysis. Additionally, NSCLC patients with high-IER5L level had significantly worse prognosis and IER5L could be used as an independent prognostic factor for NSCLC patients. Meanwhile, patients in the high-IER5L group had higher TMB level. IER5L expression was negatively correlated with the proportion of Monocytes and T cells CD4 memory resting, and was positively related to the proportion of Tregs and M0 macrophages in tumor tissues. Besides, transcription factors TFAP4 and ZNF692 may bind to the promoter region of IER5L, and then modulate IER5L gene transcription, thereby affecting IER5L gene expression. Furthermore, GSEA results showed that IER5L gene was closely related to MAPK, PI3K-Akt, NF-kappaB signaling pathways in NSCLC.

CONCLUSION

Collectively, high IER5L expression may be a promising unfavorable prognostic biomarker and therapeutic target for NSCLC patients.

IER5L is a Prognostic Biomarker in Pan-Cancer Analysis and Correlates with Immune Infiltration and Immune Molecules in Non-Small Cell Lung Cancer

Purpose

Non-small cell lung cancer (NSCLC) accounts for the majority of lung cancer cases. Immediate early response 5 like (IER5L) plays crucial roles in progression and prognosis for several tumors, but its role in NSCLC remains unclear.

Patients and Methods

Gene expression and mutation profiles, DNA methylation data, and clinical information for cancers were downloaded from multiple databases. Relative expression, prognostic value, and correlation with disease progression of IER5L were analyzed in multiple cancers, including NSCLC. Upstream mechanisms were explored using a transcriptional network. Functional enrichment analysis, protein-protein interaction network, and gene set enrichment analysis were applied to study downstream mechanisms. Correlations of IER5L with immune infiltration, immune molecules, methylation status, and tumor mutation burden (TMB) were analyzed using R language. Finally, quantitative polymerase chain reaction (qPCR) and single-cell RNA sequencing (scRNA seq) analysis were performed to validate IER5L expression in NSCLC.

Results

Pan-cancer analysis displayed that IER5L expression was upregulated in multiple cancers and was associated with disease prognosis and progression, including NSCLC, which was validated using qPCR. scRNA seq analysis showed that multiple cells had increased IER5L expression. An EGR1-hsa-miR-8075-IER5L network was constructed for NSCLC. A total of 191 DEGs were identified between the two IER5L groups, which were significantly enriched in biological process of action potential, sodium ion transport, and regulation of membrane potential. Increased IER5L expression was primarily enriched in cell cycle, NOTCH signaling pathway, and oxidative phosphorylation pathway, and was correlated with increased regulatory T cells and neutrophils, elevated levels of immune molecules, and higher TMB.

Conclusion

Our findings show that increased IER5L expression was correlated with progression and prognosis in multiple cancers as well as with immune infiltration and immune molecules in NSCLC. Thus, IER5L is a prognostic biomarker in multiple cancers and may correlate with immunotherapeutic response in NSCLC.

Comprehensive analysis of the role of SFXN family in breast cancer

The sideroflexin (SFXN) family is a group of mitochondrial membrane proteins. Although the function of the SFXN family in mitochondria has been widely recognized, the expression levels, role, and prognostic value of this family in breast cancer (BC) have not been clearly articulated and systematically analysed. In our research, SFXN1 and SFXN2 were significantly upregulated in BC versus normal samples based on Gene Expression Profiling Interactive Analysis 2 and the Human Protein Atlas databases. We found that high SFXN1 expression was significantly related to poor prognosis in BC patients and that high SFXN2 expression was significantly associated with good prognosis in BC patients. Gene Ontology analysis of the SFXN family was performed based on the STRING database to explore the potential functions of this family, including biological processes, cellular components, and molecular functions. Based on the MethSurv database, we found that two SFXN1 CpG sites (5′-UTR-S_Shelf-cg06573254 and TSS200-Island-cg17647431), two SFXN2 CpG sites (3′-UTR-Open_Sea-cg04774043 and Body-Open_Sea-cg18994254), one SFXN3 CpG site (Body-S_Shelf-cg17858697), and nine SFXN5 CpG sites (1stExon;5′-UTR-Island-cg03856450, Body-Open_Sea-cg04016113, Body-Open_Sea-cg04197631, Body-Open_Sea-cg07558704, Body-Open_Sea-cg08383863, Body-Open_Sea-cg10040131, Body-Open_Sea-cg10588340, Body-Open_Sea-cg17046766, and Body-Open_Sea-cg22830638) were significantly related to the prognosis of BC patients. According to the ENCORI database, four negative regulatory miRNAs for SFXN1 (hsa-miR-22-3p, hsa-miR-140-5p, hsa-miR-532-5p, and hsa-miR-582-3p) and four negative regulatory miRNAs for SFXN2 (hsa-miR-9-5p, hsa-miR-34a-5p, hsa-miR-532-5p, and hsa-miR-885-5p) were related to poor prognosis for BC patients. This study suggests that SFXN1 and SFXN2 are valuable biomarkers and treatment targets for patients with BC.

Identification of critical prognosis signature associated with lymph node metastasis of stomach adenocarcinomas

Lymph node metastasis (LNM) is an important factor affecting the prognosis of patients with gastric adenocarcinoma (STAD), which is the most common malignancy of the human digestive system. Current detection techniques have limited sensitivity and specificity, and there is a lack of effective biomarkers to screen for LNM. Therefore, it is critical to screen for biomarkers that predict LNM in STAD. Gene expression differential analysis (false discovery rate < 0.05, |log₂Fold change| ≥1.5) was performed on 102 LNM samples, 224 non-LNM samples, and 29 normal gastric tissue samples from The Cancer Genome Atlas (TCGA) STAD dataset, and 269 LNM-specific genes (DEGs) were obtained. Enrichment analysis showed that LNM-specific genes functioned mainly in cytokine-cytokine receptor interactions, calcium signaling, and other pathways. Ten DEGs significantly associated with overall survival in STAD patients were screened by multivariate Cox regression, and an LNM-based 10-mRNA prognostic signature was established (Logrank P < 0.0001). This 10-mRNA signature was well predicted in both the TCGA training set and the Gene Expression Omnibus validation dataset (GSE84437) and was associated with survival in patients with LNM or advanced-stage STAD. Using Kaplan-Meier survival, receiver operating characteristic curve, C-index analysis, and decision curve analysis, the 10-mRNA signature was found to be a more effective predictor of prognosis in STAD patients than the other two reported models (P < 0.0005). Protein-protein interaction network and gene set enrichment analysis of the 10-mRNA signature revealed that the signature may affect the expression of multiple biological pathways and related genes. Finally, the expression levels of prognostic genes in STAD tissues and cell lines were verified using qRT-PCR, Western blot, and the Human Protein Atlas database. Taken together, the prognostic signature constructed in this study may become an indicator for clinical prognostic assessment of LNM-STAD and provide a new strategy for future targeted therapy.

Prognostic value of SLC4A4 and its correlation with the microsatellite instability in colorectal cancer

Objective

To explore new biomarkers related to microsatellite instability in order to better predict prognosis and guide medication.

Methods

The "limma" R package was used to identify differentially expressed genes in GSE24514, and then weighted correlation network analysis was used to select key genes. Different cell types in the tumor microenvironment were identified and analyzed by single-cell sequencing, with a Lasso regression model used to screen prognostic variables. Furthermore, the correlation between microsatellite instability and potential prognostic variables was explored, as well as the expression characteristics and clinical characteristics of the prognostic variables in the TCGA, UALCAN, and HPA databases. PCR assay was used to investigate the expression of SLC4A4 in colorectal cancer cell lines. Finally, we further verified the expression of SLC4A4 by immunohistochemistry.

Results

First, 844 differentially expressed genes in GSE24514 were identified. Subsequently, weighted co-expression network analysis (WGCNA) of GSE24514 obtained all the genes significantly associated with microsatellite instability (MSI), a total of 1452. Analysis of GSE166555 single cell sequencing data set yielded 1564 differentially expressed genes. The gene sets obtained from the above three analysis processes were intersected, and 174 genes were finally obtained. The Lasso regression model revealed two potential prognostic genes, TIMP1 and SLC4A4, of which, there was a stronger correlation between microsatellite instability and SLC4A4. The mRNA and protein expression of SLC4A4 was significantly decreased in tumors, and patients with low SLC4A4 expression had a poor prognosis. In addition, SLC4A4 was specifically expressed in epithelial cells. In the microenvironment of colorectal cancer, malignant cells have a strong interaction with different stromal cells. PCR showed that SLC4A4 was significantly down-regulated in colorectal cancer cell lines Caco-2, HCT116 and HT29 compared with normal control NCM460 cell lines. Immunohistochemistry also showed low expression of SLC4A4 in colorectal cancer.

Conclusion

SLC4A4, as a tumor suppressor gene, is significantly downregulated and positively correlated with microsatellite instability, thus it may be combined with microsatellite instability to guide colorectal cancer treatment.