Comprehensive Analysis of Prognostic and Immune Infiltrates for E2F Transcription Factors in Human Pancreatic Adenocarcinoma

E2F transcription factors promote tumorigenicity in pancreatic ductal adenocarcinoma

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers with limited treatment options, illustrating an urgent need to identify new drugable targets in PDACs.

OBJECTIVE

Using the similarities between tumor development and normal embryonic development, which is accompanied by rapid cell expansion, we aimed to identify and characterize embryonic signaling pathways that were reinitiated during tumor formation and expansion.

METHODS AND RESULTS

Here, we report that the transcription factors E2F1 and E2F8 are potential key regulators in PDAC. E2F1 and E2F8 RNA expression is mainly localized in proliferating cells in the developing pancreas and in malignant ductal cells in PDAC. Silencing of E2F1 and E2F8 in PANC-1 pancreatic tumor cells inhibited cell proliferation and impaired cell spreading and migration. Moreover, loss of E2F1 also affected cell viability and apoptosis with E2F expression in PDAC tissues correlating with expression of apoptosis and mitosis pathway genes, suggesting that E2F factors promote cell cycle regulation and tumorigenesis in PDAC cells.

CONCLUSION

Our findings illustrate that E2F1 and E2F8 transcription factors are expressed in pancreatic progenitor and PDAC cells, where they contribute to tumor cell expansion by regulation of cell proliferation, viability, and cell migration making these genes attractive therapeutic targets and potential prognostic markers for pancreatic cancer.

High Expression of E2F4 Is an Adverse Prognostic Factor and Related to Immune Infiltration in Oral Squamous Cell Carcinoma

Background

We aimed to evaluate the prognostic value of E2F4 expression in oral squamous cell carcinoma (OSCC) and clarify its influence on immune cell infiltration and biological functions.

Methods

The Cancer Genome Atlas (TCGA) database, the STRING database, and related online tools as well as single-sample gene set enrichment analysis (ssGSEA) were used for the analyses in our study.

Results

The E2F4 expression was elevated in OSCC tumor tissue compared with paracancerous tissues. The high expression of E2F4 was closely related to the poorer overall survival, disease-free survival, and progression-free interval of OSCC. In addition, pathway enrichment analyses revealed that the top 49 genes most closely related to E2F4 were strongly associated with the cell cycle. E2F4-related proteins were closely related to the following KEGG pathways: cell cycle, cellular senescence, PI3K-Akt signaling pathway, Wnt signaling pathway, and notch signaling pathway. It was also demonstrated that the E2F4 expression was negatively correlated with immune purity and strongly related to immune cell infiltration in OSCC.

Conclusions

E2F4 can be used as a novel biomarker for the diagnosis and prognosis of OSCC.

CD155 Cooperates with PD-1/PD-L1 to Promote Proliferation of Esophageal Squamous Cancer Cells via PI3K/Akt and MAPK Signaling Pathways

BACKGROUND

Esophageal cancer is still a leading cause of death among all tumors in males, with unsatisfactory responses to novel immunotherapies such as anti-PD-1 agents. Herein, we explored the role of CD155 in esophageal squamous cell cancer (ESCA) and its underlying molecular mechanisms.

METHODS

Publicly available datasets were used for differential gene expression and immune infiltration analyses, and their correlation with patient survival. A total of 322 ESCA and 161 paracancer samples were collected and evaluated by performing immunohistochemistry and the H score was obtained by performing semiquantitative analysis. In vitro transfection of ESCA cell lines with lentivirus vectors targeting CD155 was performed to knockdown the protein. These cells were analyzed by conducting RNA sequencing, and the effects of CD155 knockdown on cell cycle and apoptosis were verified with flow cytometry and Western blotting. In addition, in vivo experiments using these engineered cell lines were performed to determine the role of CD155 in tumor formation. A small interfering RNA-mediated knockdown of Nectin3 was used to determine whether it phenocopied the profile of CD155 knockdown.

RESULTS

CD155 is highly expressed in ESCA tissues and is positively associated with PD1, PDL1, CD4, IL2RA, and S100A9 expression. Furthermore, CD155 knockdown inhibited ESCA cells' proliferation by impairing the cell cycle and inducing cell apoptosis. Bioinformatics analysis of the gene expression profile of these engineered cells showed that CD155 mainly contributed to the regulation of PI3K/Akt and MAPK signals. The downregulation of Nectin3 expression phenocopied the profile of CD155 knockdown.

DISCUSSION

CD155 may cooperate with PD-1/PD-L1 to support ESCA proliferation in ways other than regulating its underlying immune mechanisms. Indeed, CD155 downregulation can impair ESCA cell pro-cancerous behavior via the inhibition of the PI3K/Akt and MAPK signaling pathways. Moreover, Nectin3 may be a ligand of CD155 and participate in the regulation of ESCA cells' proliferation. Hence, the inhibition of CD155 may enhance the therapeutic effect of anti-PD-1 immunotherapies in ESCA.

Dual immunological and proliferative regulation of immune checkpoint FGL1 in lung adenocarcinoma: The pivotal role of the YY1–FGL1–MYH9 axis

Rational

Lung cancer is the most common tumor worldwide, with the highest mortality rate and second highest incidence. Immunotherapy is one of the most important treatments for lung adenocarcinoma (LUAD); however, it has relatively low response rate and high incidence of adverse events. Herein, we explored the therapeutic potential of fibrinogen-like protein 1 (FGL1) for LUAD.

Methods

Data from GEPIA and ACLBI databases were assessed to explore gene–gene correlations and tumor immune infiltration patterns. A total of 200 patients with LUAD were recruited. FGL1 levels in the serum and cellular supernatant were determined by enzyme-linked immunosorbent assay. In vitro and in vivo experiments were performed to assess the effect FGL1 on the proliferation of LUAD cells. Cocultures were performed to explore the effect of FGL1 knockdown in lung cancer cells on T cells, concerning cytokine secretion and viability. PROMO and hTFtarget databases were used for transcription factor prediction. Quantitative polymerase chain reaction (qPCR), chromatin immunoprecipitation, and dual luciferase reporter assays were performed to validate the identified transcription factor of FGL1. Immunoprecipitation, mass spectrometry and gene ontology analysis were performed to explore the downstream partners of FGL1.

Results

FGL1 expression in LUAD was positively associated with PDL1, but not for PD1 expression. Moreover, FGL1 was positively associated with the CD3D expression and negatively associated with FOXP3, S100A9, and TPSB2 within the tumor site. FGL1 promotes the secretion of interleukin-2 by T cells in vitro, simultaneously inducing their apoptosis. Indeed, YY1 is the upstream molecule of FGL1 was found to be transcriptionally regulated by YY1 and to directly by to MYH9 to promote the proliferation of LUAD cells in vitro and in vivo.

Conclusions

FGL1 is involved in the immunological and proliferative regulation of LUAD cells by controlling the secretion of important immune-related cytokines via the YY1–FGL1–MYH9 axis. Hence, targeting FGL1 in LUAD may pave the way for the development of new immunotherapies for tackling this malignancy.

Nucleophosmin 1 is a prognostic marker of gastrointestinal cancer and is associated with m6A and cuproptosis

Background: NPM1 is highly expressed in a variety of solid tumors and promotes tumor development. However, there are few comprehensive studies on NPM1 analysis in gastrointestinal cancer.

Methods: We used bioinformatics tools to study the expression difference of NPM1 between gastrointestinal cancer and control group, and analyzed the relationship between its expression level and the diagnosis, prognosis, functional signaling pathway, immune infiltration, m6A and cuproptosis related genes of gastrointestinal cancer. At the same time, the expression difference of NPM1 between esophageal carcinoma (ESCA) samples and control samples was verified by in vitro experiments.

Results: NPM1 was overexpressed in gastrointestinal cancer. In vitro experiments confirmed that the expression of NPM1 in ESCA samples was higher than that in normal samples. The expression of NPM1 has high accuracy in predicting the outcome of gastrointestinal cancer. The expression of NPM1 is closely related to the prognosis of multiple gastrointestinal cancers. Go and KEGG enrichment analysis showed that NPM1 co-expressed genes involved in a variety of biological functions. NPM1 expression is potentially associated with a variety of immune cell infiltration, m6A and cuproptosis related genes in gastrointestinal cancers.

Conclusion: NPM1 can be used as a diagnostic and prognostic marker of gastrointestinal cancer, which is related to the immune cell infiltration and the regulation of m6A and cuproptosis.

Comprehensive Analysis of the Immune Cell Infiltration Landscape and Immune-Related Methylation in Retinoblastoma

Retinoblastoma is a common pediatric intraocular cancer, originating from cone precursors. The development of immunotherapies can help eradicate the tumor without vision loss, which would largely improve the quality of life of patients with retinoblastoma. Investigation of the tumor immune microenvironment provides knowledge for developing novel immunotherapies in cancer. However, the immune cell infiltrative landscape of retinoblastoma is unknown. Here, we compared the relative expression of immune gene signatures among 59 patients with retinoblastoma. The patients were divided into two subgroups according to the 28 types of immune cell infiltration (ICI) scores. We found that a subgroup with high ICI scores had increased expression levels of late cone markers, while the other subgroup exhibited larger tumor size and metastasis propensity. Furthermore, hypermethylated genes in the high-ICI subgroup were associated with immune regulation in the tumor microenvironment, suggesting that DNA methylation may play a vital regulatory role in retinoblastoma immunity. Our study provides a comprehensive framework for the systemic analysis of the influences of epigenetic events on the tumor immune microenvironment. We anticipate that our assay can not only provide insights into tumor immune regulation but also open up the perspectives for the identification of novel immunotherapy targets for retinoblastoma.

Novel Insights into the Prognosis and Immunological Value of the SLC35A (Solute Carrier 35A) Family Genes in Human Breast Cancer

According to statistics 2020, female breast cancer (BRCA) became the most commonly diagnosed malignancy worldwide. Prognosis of BRCA patients is still poor, especially in population with advanced or metastatic. Particular functions of each members of the solute carrier 35A (SLC35A) gene family in human BRCA are still unknown regardless of awareness that they play critical roles in tumorigenesis and progression. Using integrated bioinformatics analyses to identify therapeutic targets for specific cancers based on transcriptomics, proteomics, and high-throughput sequencing, we obtained new information and a better understanding of potential underlying molecular mechanisms. Leveraging BRCA dataset that belongs to The Cancer Genome Atlas (TCGA), which were employed to clarify SLC35A gene expression levels. Then we used a bioinformatics approach to investigate biological processes connected to SLC35A family genes in BRCA development. Beside that, the Kaplan-Meier estimator was leveraged to explore predictive values of SLC35A family genes in BCRA patients. Among individuals of this family gene, expression levels of SLC35A2 were substantially related to poor prognostic values, result from a hazard ratio of 1.3 (with 95 percent confidence interval (95% CI: 1.18-1.44), the p for trend (ptrend) is 3.1 × 10-7). Furthermore, a functional enrichment analysis showed that SLC35A2 was correlated with hypoxia-inducible factor 1A (HIF1A), heat shock protein (HSP), E2 transcription factor (E2F), DNA damage, and cell cycle-related signaling. Infiltration levels observed in specific types of immune cell, especially the cluster of differentiation found on macrophages and neutrophils, were positively linked with SLC35A2 expression in multiple BRCA subclasses (luminal A, luminal B, basal, and human epidermal growth factor receptor 2). Collectively, SLC35A2 expression was associated with a lower recurrence-free survival rate, suggesting that it could be used as a biomarker in treating BRCA.

Comprehensive Analysis of Hexokinase 2 Immune Infiltrates and m6A Related Genes in Human Esophageal Carcinoma

Background: Hexokinase 2 not only plays a role in physiological function of human normal tissues and organs, but also plays a vital role in the process of glycolysis of tumor cells. However, there are few comprehensive studies on HK2 in esophageal carcinoma (ESCA) needs further study.

Methods: Oncomine, Tumor Immune Estimation Resource (TIMER), The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database were used to analyze the expression differences of HK2 in Pan-cancer and ESCA cohort, and to analyze the correlation between HK2 expression level and clinicopathological features of TCGA ESCA samples. GO/KEGG, GGI, and PPI analysis of HK2 was performed using R software, LinkedOmics, GeneMANIA and STRING online tools. The correlation between HK2 and ESCA immune infiltration was analyzed TIMER and TCGA ESCA cohort. The correlation between HK2 expression level and m6A modification of ESCA was analyzed by utilizing TCGA ESCA cohort.

Results: HK2 is highly expressed in a variety of tumors, and its high expression level in ESCA is closely related to the weight, cancer stages, tumor histology and tumor grade of ESCA. The analysis results of GO/KEGG showed that HK2 was closely related to cell adhesion molecule binding, cell-cell junction, ameboidal-type cell migration, insulin signaling pathway, hif-1 signaling pathway, and insulin resistance. GGI showed that HK2 associated genes were mainly involved in the glycolytic pathway. PPI showed that HK2 was closely related to HK1, GPI, and HK3, all of which played an important role in tumor proliferation. The analysis results of TIMER and TCGA ESCA cohort indicated that the HK2 expression level was related to the infiltration of various immune cells. TCGA ESCA cohort analyze indicated that the HK2 expression level was correlated with m6A modification genes.

Conclusion: HK2 is associated with tumor immune infiltration and m6A modification of ESCA, and can be used as a potential biological target for diagnosis and therapy of ESCA.

Identification of CDK2-Related Immune Forecast Model and ceRNA in Lung Adenocarcinoma, a Pan-Cancer Analysis

Background

Tumor microenvironment (TME) plays important roles in different cancers. Our study aimed to identify molecules with significant prognostic values and construct a relevant Nomogram, immune model, competing endogenous RNA (ceRNA) in lung adenocarcinoma (LUAD).

Methods

“GEO2R,” “limma” R packages were used to identify all differentially expressed mRNAs from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. Genes with P-value <0.01, LogFC>2 or <-2 were included for further analyses. The function analysis of 250 overlapping mRNAs was shown by DAVID and Metascape software. By UALCAN, Oncomine and R packages, we explored the expression levels, survival analyses of CDK2 in 33 cancers. “Survival,” “survminer,” “rms” R packages were used to construct a Nomogram model of age, gender, stage, T, M, N. Univariate and multivariate Cox regression were used to establish prognosis-related immune forecast model in LUAD. CeRNA network was constructed by various online databases. The Genomics of Drug Sensitivity in Cancer (GDSC) database was used to explore correlations between CDK2 expression and IC50 of anti-tumor drugs.

Results

A total of 250 differentially expressed genes (DEGs) were identified to participate in many cancer-related pathways, such as activation of immune response, cell adhesion, migration, P13K-AKT signaling pathway. The target molecule CDK2 had prognostic value for the survival of patients in LUAD (P = 5.8e-15). Through Oncomine, TIMER, UALCAN, PrognoScan databases, the expression level of CDK2 in LUAD was higher than normal tissues. Pan-cancer analysis revealed that the expression, stage and survival of CDK2 in 33 cancers, which were statistically significant. Through TISIDB database, we selected 13 immunodepressants, 21 immunostimulants associated with CDK2 and explored 48 genes related to these 34 immunomodulators in cBioProtal database (P < 0.05). Gene Set Enrichment Analysis (GSEA) and Metascape indicated that 49 mRNAs were involved in PUJANA ATM PCC NETWORK (ES = 0.557, P = 0, FDR = 0), SIGNAL TRANSDUCTION (ES = –0.459, P = 0, FDR = 0), immune system process, cell proliferation. Forest map and Nomogram model showed the prognosis of patients with LUAD (Log-Rank = 1.399e-08, Concordance Index = 0.7). Cox regression showed that four mRNAs (SIT1, SNAI3, ASB2, and CDK2) were used to construct the forecast model to predict the prognosis of patients (P < 0.05). LUAD patients were divided into two different risk groups (low and high) had a statistical significance (P = 6.223e-04). By “survival ROC” R package, the total risk score of this prognostic model was AUC = 0.729 (SIT1 = 0.484, SNAI3 = 0.485, ASB2 = 0.267, CDK2 = 0.579). CytoHubba selected ceRNA mechanism medicated by potential biomarkers, 6 lncRNAs-7miRNAs-CDK2. The expression of CDK2 was associated with IC50 of 89 antitumor drugs, and we showed the top 20 drugs with P < 0.05.

Conclusion

In conclusion, our study identified CDK2 related immune forecast model, Nomogram model, forest map, ceRNA network, IC50 of anti-tumor drugs, to predict the prognosis and guide targeted therapy for LUAD patients.

Prognostic Values of E2F1/2 Transcriptional Expressions in Chromophobe Renal Cell Carcinoma Patients: Evidence from Bioinformatics Analysis

Background

Numerous studies on the E2F transcription factors have led to increasing insights that E2Fs could be an important driver of the formation and progression of many human cancers. Little is known about the function of distinct E2Fs in chromophobe renal cell carcinoma (chRCC).

Methods

We utilized the UALCAN, GEPIA, Cancer Genome Atlas (TCGA) database, cBioPortal, Metascape, STRING, Cytoscape, GeneMANIA, TIMER, TISIDB, GSCALite, and MEXPRESS databases to investigate the transcription level, genetic alteration, methylation, and biological function of E2Fs in chRCC patients, and its association with the occurrence, progress, prognosis, and immune cell infiltration in patients with chRCC.

Results

We found that E2F1/2/4/7/8 were more expressed in chRCC tissues than in normal tissues, while the expression of E2F5/6 was lower in the former than in the latter, and the expression levels of E2F1/2/4/5/6//7/8 were also associated with the histological parameters of chRCC, including T-stage and N-stage. Higher expression of E2F1/2/7/8 was found to be significantly correlated with worse overall survival (OS) in chRCC patients. Cox regression and time-dependent ROC analysis further suggested that E2F1/2 could be the potential independent biomarkers for chRCC prognosis. Besides, a moderate mutation rate of E2Fs (34%) was noticed in chRCC, and the genetic mutations in E2Fs were associated with poor survival of chRCC patients. We noticed that the expression of E2Fs was statistically correlated with the immune cell infiltration in chRCC. Moreover, we also found that the expression of E2F1 was significantly correlated with tumor-infiltrating lymphocytes and immunomodulators, E2F7 expression was associated with MHC molecules, and the expression of E2F1/8 was correlated to their methylation levels.

Conclusion

Our results provide novel insights for selecting the prognostic biomarkers for chRCC and suggest that E2F1/2 could act as potential prognostic biomarkers for the survival of chRCC patients. However, more in-depth experiments are required to identify the underlying mechanisms and verify the clinical value of E2F1/2 in the prognosis of chRCC.