Pyroptosis regulators exert crucial functions in prognosis, progression and immune microenvironment of pancreatic adenocarcinoma: a bioinformatic and in vitro research

Pyroptosis is an inflammatory programmed cell death, showing potentials to be a novel anti-cancer approach. However, the roles of pyroptosis-related (PR) genes (PRGs) in pancreatic adenocarcinoma (PAAD) remain elusive. In the present study, we constructed a novel PR risk signature through the lasso regression analysis. The risk signature was greatly conducive to PAAD prognostic assessment. PR risk score was identified as an independent prognostic factor and could distinguish the prognostic differences of most clinical subgroups. Meanwhile, it could improve the traditional prognostic models based on TNM-staging. Next, its prognostic value was also tested in five validation cohorts. Using CIBERSORT, ESTIMATE, and ssGSEA algorithms, the effects of PR risk signature on tumor immune microenvironment (TIM) were explored. High PR risk suppressed antitumor immune through decreasing the infiltrating levels of CD8 T and NK cells. The genomic information and histological expression of risk PRGs were uncovered by USCA and HPA databases. Somatic mutation, methylation alteration, and homozygous CNV of eight PRGs barely occurred in PAAD samples. As for therapeutic correlation, PR risk score may not predict the efficacy of PD-1/L1 inhibitors and was weakly associated with multiple drug susceptibilities. Finally, the biofunctions of toll like receptor 3 (TLR3) in pancreatic cancer (PC) cells were investigated through qPCR, MTT, colony formation, and Transwell assays. Overexpression of TLR3 could promote the proliferation, migration, and invasion of PC cells. In conclusion, PRGs play crucial roles in prognosis, progression, and immune microenvironment of PAAD. TLR3 is expected to be a promising therapeutic target.

A comprehensive multi-omics analysis reveals molecular features associated with cancer via RNA cross-talks in the Notch signaling pathway

•

Many Notch genes are identified as cancer-associated genes with an important role in tumorigenesis.

•

Dynamic expression patterns are associated with the Notch activity that are largely regulated by multiple ncRNAs.

•

Cross-talks among diverse RNAs are crucial in cancers via ceRNA network.

•

The Notch pathway shows a robust prognostic ability via integrating multi-omics features as well as their targets.

•

The Notch pathway is also correlated with immune infiltration and maybe available cancer treatment drug targets.

The Notch signaling has an important role in multiple cellular processes and is related to carcinogenic process. To understand the potential molecular features of the crucial Notch pathway, a comprehensive multi-omics analysis is performed to explore its contributions in cancer, mainly including analysis of somatic mutation landscape, pan-cancer expression, ncRNA regulation and potential prognostic power. The screened 22 Notch core genes are relative stable in DNA variation. Dynamic expression patterns are associated with the Notch activity, which are mainly regulated by multiple ncRNAs via interactions of ncRNA:mRNA and ceRNA networks. The Notch pathway shows a potential prognostic ability through integrating multi-omics features as well as their targets, and it is correlated with immune infiltration and maybe available drug targets, implying the potential role in individualized treatment. Collectively, all of these findings contribute to exploring crucial role of the key pathway in cancer pathophysiology and gaining mechanistic insights into cross-talks among RNAs and biological pathways, which indicates the possible application of the well-conserved Notch signaling pathway in precision medicine.

A combined signature of glycolysis and immune landscape predicts prognosis and therapeutic response in prostate cancer

Prostate cancer (PCa) is a common malignancy that poses a major threat to the health of men. Prostate-specific antigen (PSA) and its derivatives, as FDA-approved detection assays, are insufficient to serve as optimal markers for patient prognosis and clinical decision-making. It is widely acknowledged that aberrant glycolytic metabolism in PCa is related to tumor progression and acidifies the tumor microenvironment (TME). Considering the non-negligible impacts of glycolysis and immune functions on PCa, we developed a combined classifier in prostate cancer. The Glycolysis Score containing 19 genes and TME Score including three immune cells were created, using the univariate and multivariate Cox proportional hazards model, log-rank test, least absolute shrinkage and selection operator (LASSO) regression analysis and the bootstrap approach. Combining the glycolysis and immunological landscape, the Glycolysis-TME Classifier was then constructed. It was observed that the classifier was more accurate in predicting the prognosis of patients than the current biomarkers. Notably, there were significant differences in metabolic activity, signaling pathways, mutational landscape, immunotherapeutic response, and drug sensitivity among the Glycolysis^high/TME^low, Mixed group and Glycolysis^low/TME^high identified by this classifier. Overall, due to the significant prognostic value and potential therapeutic guidance of the Glycolysis-TME Classifier, we anticipate that this classifier will be clinically beneficial in the management of patients with PCa.

Identification of Therapeutic Targets and Prognostic Biomarkers Among Chemokine (C-C Motif) Ligands in the Liver Hepatocellular Carcinoma Microenvironment

Background: Liver hepatocellular carcinoma (LIHC) is the third leading cause of cancer-related death and the sixth most common solid tumor worldwide. In the tumor microenvironment, the cross-talk between cancer cells, immune cells, and stromal cells exerts significant effects on neoplasia and tumor development and is modulated in part by chemokines. Chemokine (C-C motif) ligands (CCL) can directly target tumor cells and stromal cells, and they have been shown to regulate tumor cell proliferation, cancer stem-like cell properties, cancer invasiveness and metastasis, which directly and indirectly affect tumor immunity and influence cancer progression, therapy and patient outcomes. However, the prognostic values of chemokines CCL in LIHC have not been clarified.

Methods: In this study, we comprehensively analyzed the relationship between transcriptional chemokines CCL and disease progression of LIHC using the ONCOMINE dataset, GEPIA, UALCAN, STRING, WebGestalt, GeneMANIA, TRRUST, DAVID 6.8, LinkedOmics, TIMER, GSCALite, and Open Targets. We validated the protein levels of chemokines CCL through western blot and immunohistochemistry.

Results: The transcriptional levels of CCL5/8/11/13/15/18/20/21/25/26/27/28 in LIHC tissues were significantly elevated while CCL2/3/4/14/23/24 were significantly reduced. A significant correlation was found between the expression of CCL14/25 and the pathological stage of LIHC patients. LIHC patients with low transcriptional levels of CCL14/21 were associated with a significantly poor prognosis. The functions of differentially expressed chemokines CCL were primarily related to the chemokine signaling pathway, cytokine–cytokine receptor interactions, and TNF-α signaling pathway. Our data suggested that RELA/REL, NFKB1, STAT1/3/6, IRF3, SPI1, and JUN were key transcription factors for chemokines CCL. We found significant correlations among the expression of chemokines CCL and the infiltration of six types of immune cells (B cells, CD8⁺ T cells, CD4⁺ T cells, macrophages, neutrophils, and dendritic cells) and immune checkpoints (PD-1. PD-L1, and CTLA-4). The western blot and immunohistochemistry results showed that protein expression levels of CCL5 and CCL20 were upregulated in LIHC. CCL5 and CCL20 were significantly correlated with the clinical outcome of patients with LIHC, and could be negatively regulated by some drugs or small molecules.

Conclusions: Our results may provide novel insights for the potential suitable targets of immunological therapy and prognostic biomarkers for LIHC.

IRF7 is a Prognostic Biomarker and Associated with Immune Infiltration in Stomach Adenocarcinoma

Background

Stomach adenocarcinoma (STAD) is one of the most prevalent malignances, ranking fifth in incidence and third in mortality among all malignances. Interferon regulatory factors (IRFs) play a vital role in immune response and tumor cellular biological process. The roles of IRFs in STAD are far from being systematically clarified.

Methods

A series of bioinformatics tools, including GEPIA, UALCAN, TIMER, Kaplan–Meier plotter and LinkedOmics, were applied to explore the expression and clinical significance of IRFs in STAD.

Results

IRF3/7 expression were upregulated in STAD in sub-group analyses based on race, gender, age, H. Pylori infection status, histological subtypes, tumor grade, individual cancer stages, and nodal metastasis status. High IRF3/7 expression were associated with poor overall survival (OS), post-progression survival (PFPS) and first progression (FP) in STAD. IRF3 and IRF7 were altered in 5% and 6% of all TCGA STAD patients. Further analysis revealed that IRF7 was significantly associated with the abundance of immune cells (B cells, Neutrophils and Dendritic cells) and the expression of most immune biomarkers. Enrichment analysis indicated that IRF7 was mainly involved in adaptive immune response, NOD-like receptor signaling pathway, Necroptosis, and Toll-like receptor signaling pathway. We also identified several IRF7-associated kinase and miRNA targets in STAD. The result of verified experiment revealed that ITF7 expression was increased in STAD tissues compared with normal tissues and prognosis analysis revealed that STAD patients with high IRF7 expression had a poor overall survival.

Conclusion

IRF7 is upregulated in STAD and associated with poor OS, PPS and FP. Moreover, IRF7 is significantly associated with the abundance of immune cells and the expression of most immune biomarkers, suggesting that IRF7 is as a prognostic biomarker and associated with immune infiltration in STAD.

In silico tools and databases for designing cancer immunotherapy

Immunotherapy is a rapidly growing therapy for cancer which have numerous benefits over conventional treatments like surgery, chemotherapy, and radiation. Overall survival of cancer patients has improved significantly due to the use of immunotherapy. It acts as a novel pillar for treating different malignancies from their primary to the metastatic stage. Recent preferments in high-throughput sequencing and computational immunology leads to the development of targeted immunotherapy for precision oncology. In the last few decades, several computational methods and resources have been developed for designing immunotherapy against cancer. In this review, we have summarized cancer-associated genomic, transcriptomic, and mutation profile repositories. We have also enlisted in silico methods for the prediction of vaccine candidates, HLA binders, cytokines inducing peptides, and potential neoepitopes. Of note, we have incorporated the most important bioinformatics pipelines and resources for the designing of cancer immunotherapy. Moreover, to facilitate the scientific community, we have developed a web portal entitled ImmCancer (https://webs.iiitd.edu.in/raghava/immcancer/), comprises cancer immunotherapy tools and repositories.

Multiomics profiling of the expression and prognosis of MCMs in endometrial carcinoma

Minichromosome maintenance (MCM) family members are a group of genes involved in regulating DNA replication and cell division and have been identified as oncogenes in various cancer types. Several experimental studies have suggested that MCMs are dysregulated in endometrial carcinoma (EC). However, the expression pattern, clinical value and functions of different MCMs have yet to be analyzed systematically and comprehensively. We analyzed expression, survival rate, DNA alteration, PPT network, GGI network, functional enrichment cancer hallmarks and drug sensitivity of MCMs in patients with EC based on diverse datasets, including Oncomine, GEPIA, Kaplan–Meier Plotter, HPA, Sangerbox and GSCALite databases. The results indicated that most MCM members were increased in EC and showed a prognostic value in survival analysis, which were considerately well in terms of PFS and OS prognostic prediction. Importantly, functional enrichment, PPI network and GGI network suggested that MCMs interact with proteins related to DNA replication and cell division, which may be the mechanism of MCM promote EC progression. Further data mining illustrated that MCMs have broad DNA hypomethylation levels and high levels of copy number aberrations in tumor tissue samples, which may be the mechanism causing the high expression level of MCMs. Moreover, MCM2 can activate or suppress diverse cancer-related pathways and is implicated in EC drug sensitivity. Taking together, our findings illustrate the expression pattern, clinical value and function of MCMs in EC and imply that MCMs are potential targets for precision therapy and new biomarkers for the prognosis of patients with EC.

Microarray analysis to explore the effect of CXCL12 isoforms in a pancreatic pre-tumor cell model

CXCL12 expression was significantly lower in tumor samples than in corresponding normal samples. CXCL12 expression was significantly positively related to the infiltration levels of T cells, dendritic cells (DCs), immature DCs, cytotoxic cells, Tfh cells, mast cells, B cells, Th1 cells, natural killer (NK) cells, pDCs, neutrophils, and T helper cells (Spearman correlation coefficient > 0.5, P < 0.001) and negatively correlated with the infiltration level of NK CD56bright cells. In addition, pancreatic hTERT-HPNE cells treated with three diverse CXCL12 isoforms exhibited changes mainly in the regulation of the epithelial-mesenchymal transition activation pathway.

Comprehensive Analysis Identified ETV7 as a Potential Prognostic Biomarker in Bladder Cancer

Background

The tumor microenvironment (TME) plays a crucial role in the initiation and progression of cancer. Bladder cancer (BLCA) is a malignant tumor of the genitourinary system. Its heterogeneity results in significant differences in the prognosis of patients. To date, this is still a huge challenge for clinical treatment. In recent years, more and more evidence showed that dysregulation of transcription factors (TFs) plays an important role in tumor progression, invasion, and metastasis. Unfortunately, the role of TFs on the tumor microenvironment in bladder cancer is unclear.

Methods

The original data of BLCA and corresponding adjacent tissues were obtained from The Cancer Genome Atlas (TCGA) database. TFs were downloaded from the Animal Transcription Factor DataBase (Animal TFDB). Intersection analysis was used to obtain TFs that were differentially expressed between tumor and adjacent tissues. Gene Set Cancer Analysis (GSCALite) and CIBERSORT software were used to reveal the key differentially expressed TFs (DE-TFs). Subsequently, UALCAN and Human Protein Atlas (HPA) databases were used to disclose the expression of key DE-TFs in BLCA. The K-M curve divulged the relationship between the key DE-TFs and the patient's overall survival (OS), and the univariate and multivariate Cox regression analyses were conducted to explore independent prognostic factors. The cluster profiler package and Gene Set Enrichment Analysis (GSEA) were used for functional enrichment of genes related to the key DE-TFs. Finally, CIBERSORT software analyzed the immune landscape of BLCA.

Results

We obtained a total of 117 BLCA-related DE-TFs. Among them, ETV7 was identified as the key DE-TFs due to its association with the autophagy activation pathway and various immune cells in cancer. Online databases of UALCAN and HPA indicated that ETV7 was overexpressed in tumors and negatively correlated with tumor severity. The K-M curve showed that the OS of patients with high expression of ETV7 was poor, which indicated that it was an independent prognostic factor. Functional enrichment of 87 DEGs between ETV7-high and -low expression groups indicated that it was closely related to the immune response and the functions of a variety of immune cells. Finally, CIBERSORT results proved that the high and low expression of ETV7 also caused significant differences in the tumor immune microenvironment of patients.

Conclusion

Overall, we proved that the transcription factor ETV7 was a novel prognostic factor, which may improve the individualized outcome prediction in BLCA by regulating the tumor immune microenvironment.

Multi-Omics Reveal the Immunological Role and the Theragnostic Value of miR-216a/GDF15 Axis in Human Colon Adenocarcinoma

Colon adenocarcinoma (COAD) is the most common type of gastrointestinal cancer and is still the third leading cause of cancer-related mortality worldwide. Accurate screening tools for early diagnosis and prediction of prognosis and precision treatment strategies are urgently required to accommodate the unmet medical needs of COAD management. We herein aimed to explore the significance of the microRNA (miR)-216a/growth differentiation factor 15 (GDF15) axis in terms of clinical value, tumor immunity, and potential mechanisms in COAD by using multi-omic analysis. The gene expression levels of miR-216a and GDF15 showed an increase in the COAD group compared to those of the normal group. The expression of miR-216a presented a negative correlation with GDF15 in COAD tumor tissue. The use of an in vitro luciferase reporter assay and bioinformatic prediction revealed that miR-216a-3p acted toward translational inhibition on GDF15 by targeting its 3′untranslated region (UTR) site. High miR-216a expression was associated with decreased overall survival (OS), while the high expression of GDF15 was associated with increased OS. Enriched type 1 T-helper (Th1), enriched regulatory T (Treg), enriched eosinophils, and decreased nature killer T-cells (NKTs) in COAD tumor tissue may play counteracting factors on the tumor-regulatory effects of miR-216a and GDF15. In addition, high GDF15 expression had associations with suppressed immunoinhibitory genes and negative correlations with the infiltration of macrophages and endothelial cells. The enrichment analysis revealed that GDF15 and its co-expression network may be implicated in mitochondrial organization, apoptosis signaling, and endoplasmic reticulum (ER) stress response. The Genomics of Drug Sensitivity in Cancer (GDSC) and Cancer Therapeutics Response Portal (CTRP) analysis identified that Gemcitabine acted as a precision treatment for COAD when GDF15 expression was low. This study supports the miR-216a/GDF15 axis as a diagnostic/prognostic panel for COAD, identifies Th1, Treg, eosinophils, and NKTs as counteracting factors, indicates potential relationships underlying immunomodulation, mitochondrial organization, apoptotic signaling, and ER stress and unveil Gemcitabine as a potential drug for the development of treatment strategy when combined with targeting GDF15.

Identification of Novel Biomarkers for Predicting Prognosis and Immunotherapy Response in Head and Neck Squamous Cell Carcinoma Based on ceRNA Network and Immune Infiltration Analysis

Objectives

Patients with head and neck squamous cell carcinoma (HNSCC) have poor prognosis and show poor responses to immune checkpoint (IC) inhibitor (ICI) therapy. Competing endogenous RNA (ceRNA) networks, tumor-infiltrating immune cells (TIICs), and ICIs may influence tumor prognosis and response rates to ICI therapy. This study is aimed at identifying prognostic and IC-related biomarkers and key TIIC signatures to improve prognosis and ICI therapy response in HNSCC patients.

Methods and Results

Ninety-five long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and 1746 mRNAs were identified using three independent methods. We constructed a ceRNA network and estimated the proportions of 22 immune cell subtypes. Ten ceRNAs were related to prognosis according to Kaplan–Meier analysis. Two risk signatures based, respectively, on nine ceRNAs (ANLN, CFL2, ITGA5, KDELC1, KIF23, NFIA, PTX3, RELT, and TMC7) and three immune cell types (naïve B cells, neutrophils, and regulatory T cells) via univariate Cox regression, least absolute shrinkage and selection operator, and multivariate Cox regression analyses could accurately and independently predict the prognosis of HNSCC patients. Key mRNAs in the ceRNA network were significantly correlated with naïve B cells and regulatory T cells and with stage, grade, and immune and molecular subtype. Eight IC genes exhibited higher expression in tumor tissues and were correlated with eight key mRNAs in the ceRNA network in HNSCC patients with different HPV statuses according to coexpression and TIMER 2.0 analyses. Most drugs were effective in association with expression of these key signatures (ANLN, CFL2, ITGA5, KIF23, NFIA, PTX3, RELT, and TMC7) based on GSCALite analysis. The prognostic value of key biomarkers and associations between key ceRNAs and IC genes were validated using online databases. Eight key ceRNAs were confirmed to predict response to ICI in other cancers based on TIDE analysis.

Conclusions

We constructed two risk signatures to accurately predict prognosis in HNSCC. Key IC-related signatures may be associated with response to ICI therapy. Combinations of ICIs with inhibitors of eight key mRNAs may improve survival outcomes of HNSCC patients.

Identification of the Nerve-Cancer Cross-Talk-Related Prognostic Gene Model in Head and Neck Squamous Cell Carcinoma

The incidence of head and neck squamous cell carcinoma (HNSC) is increasing year by year. The nerve is an important component of the tumor microenvironment, which has a wide range of cross-talk with tumor cells and immune cells, especially in highly innervated organs, such as head and neck cancer and pancreatic cancer. However, the role of cancer-nerve cross-talk-related genes (NCCGs) in HNSC is unclear. In our study, we constructed a prognostic model based on genes with prognostic value in NCCGs. We used Pearson’s correlation to analyze the relationship between NCCGs and immune infiltration, microsatellite instability, tumor mutation burden, drug sensitivity, and clinical stage. We used single-cell sequencing data to analyze the expression of genes associated with stage in different cells and explored the possible pathways affected by these genes via gene set enrichment analysis. In the TCGA-HNSC cohort, a total of 23 genes were up- or downregulated compared with normal tissues. GO and KEGG pathway analysis suggested that NCCGs are mainly concentrated in membrane potential regulation, chemical synapse, axon formation, and neuroreceptor-ligand interaction. Ten genes were identified as prognosis genes by Kaplan-Meier plotter and used as candidate genes for LASSO regression. We constructed a seven-gene prognostic model (NTRK1, L1CAM, GRIN3A, CHRNA5, CHRNA6, CHRNB4, CHRND). The model could effectively predict the 1-, 3-, and 5-year survival rates in the TCGA-HNSC cohort, and the effectiveness of the model was verified by external test data. The genes included in the model were significantly correlated with immune infiltration, microsatellite instability, tumor mutation burden, drug sensitivity, and clinical stage. Single-cell sequencing data of HNSC showed that CHRNB4 was mainly expressed in tumor cells, and multiple metabolic pathways were enriched in high CHRNB4 expression tumor cells. In summary, we used comprehensive bioinformatics analysis to construct a prognostic gene model and revealed the potential of NCCGs as therapeutic targets and prognostic biomarkers in HNSC.

Development and Clinical Validation of Novel 8-Gene Prognostic Signature Associated With the Proportion of Regulatory T Cells by Weighted Gene Co-Expression Network Analysis in Uterine Corpus Endometrial Carcinoma

Background

Uterine corpus endometrial carcinoma (UCEC) is a gynecological malignant tumor with low survival rate and poor prognosis. The traditional clinicopathological staging is insufficient to estimate the prognosis of UCEC. It is necessary to select a more effective prognostic signature of UCEC to predict the prognosis and immunotherapy effect of UCEC.

Methods

CIBERSORT and weighted correlation network analysis (WGCNA) algorithms were combined to screen modules related to regulatory T (Treg) cells. Subsequently, univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses were used to identify the genes in key modules. The difference in overall survival (OS) between high- and low-risk patients was analyzed by Kaplan-Meier analysis. The Tregs-related risk signature (TRRS) was screened by uni- and multivariate Cox analyses. Afterward, we analyzed the expression difference of TRRS and verified its ability to predict the prognosis of UCEC and the effect of immunotherapy.

Results

Red module has the highest correlation with Tregs among all clustered modules. Pathways enrichment indicated that the related processes of UCEC were primarily associated to the immune system. Eight genes (ZSWIM1, NPRL3, GOLGA7, ST6GALNAC4, CDC16, ITPK1, PCSK4, and CORO1B) were selected to construct TRRS. We found that this TRRS is a significantly independent prognostic factor of UCEC. Low-risk patients have higher overall survival than high-risk patients. The immune status of different groups was different, and tumor-related pathways were enriched in patients with higher risk score. Low-risk patients are more likely take higher tumor mutation burden (TMB). Meanwhile, they are more sensitive to chemotherapy than patients with high-risk score, which indicated a superior prognosis. Immune checkpoints such as PD-1, CTLA4, PD-L1, and PD-L2 all had a higher expression level in low-risk group. TRRS expression really has a relevance with the sensitivity of UCEC patients to chemotherapeutic drugs.

Conclusion

We developed and validated a TRRS to estimate the prognosis and reflect the immune status of UCEC, which could accurately assess the prognosis of patients with UCEC and supply personalized treatments for them.

A Pan-Cancer Analysis of the Oncogenic and Immunogenic Role of m6Am Methyltransferase PCIF1

Background

Phosphorylated CTD-interacting factor 1 (PCIF1) is identified as the only known methyltransferase of N6,2′-O-dimethyladenosine (m6Am) in mRNA. However, its oncogenic and immunogenic role in cancer research is at an initial stage.

Methods

Herein, we carried out a pan-cancer analysis of PCIF1, with a series of datasets (e.g., TIMER2.0, GEPIA2, cBioPortal).

Results

PCIF1 expression was higher in most cancers than normal tissues and was discrepant across pathological stages. Highly expressed PCIF1 was positively correlated with overall survival (OS) or disease-free survival (DFS) of some tumors. PCIF1 expression had a positive correlation with CD4⁺ T-cell infiltration in kidney renal clear cell carcinoma (KIRC), CD8⁺ T cells, macrophages, and B cells in thyroid carcinoma (THCA), and immune checkpoint genes (ICGs) in LIHC but a negative correlation with CD4⁺ T cells, neutrophils, myeloid dendritic cells, and ICGs in THCA. It also affected tumor mutational burden (TMB) and microsatellite instability (MSI) of most tumors.

Conclusion

PCIF1 expression was correlated with cancer prognosis and immune infiltration, suggesting it to be a potential target for cancer therapy.

Hypomethylation-driven AKT Serine/Threonine Kinase 3 promotes testicular germ cell tumors proliferation and negatively correlates to immune infiltration

AKT Serine/Threonine Kinase 3 (AKT3) has been reported to play an important role in different tumors. However, its clinical value, biological function, and molecular mechanism in testicular germ cell tumors (TGCT) remains unclear. In the current study, we applied the Gene Set Cancer Analysis (GSCA), UCSC XENA, Gene Expression Omnibus (GEO), the Human Protein Atlas (HPA), LinkedOmics, DiseaseMeth version 2.0, TISIDB, and other databases for TGCT data mining. Then, we investigated AKT3’s mechanism of action and clinical survival significance via bioinformatics followed by in vitro experiments. We found that AKT3 was upregulated and had frequent copy number amplifications in TGCT, which were associated with poor survival outcomes of patients. On the other hand, mutations that led to AKT3 loss-of-function were correlated to a better prognosis in patients. Moreover, AKT3 silencing significantly inhibited the proliferation, DNA synthesis and colony formation of NCCIT cells (a TGCT cell line). AKT3 might participate in TGCT progression through multiple signaling pathways, such as ErbB, oxidative phosphorylation, and affecting tumor immune infiltration. Also, the upregulation of AKT3 mRNA expression might be driven by the hypomethylation of its promoter region. Overall, AKT3 is a potential TGCT oncogene and can be further used as a therapeutic target.

The role of AKR1 family in tamoxifen resistant invasive lobular breast cancer based on data mining

Background

Tamoxifen (TAM) resistance to invasive lobular cell carcinoma is a challenge for breast cancer treatment. This study explored the role of Aldo-keto reductase family 1 (AKR1) family in tamoxifen-resistant aggressive lobular breast cancer based on data mining.

Methods

TAM-resistant invasive lobular breast cancer gene chip was downloaded from the Gene Expression Omnibus (GEO) database (accession-numbered as GSE96670). The online analytical tool GEO2R was used to screen for differentially expressed genes in TAM-resistant invasive lobular breast cancer cells and TAM-sensitive counterparts. A protein-protein interaction (PPI) networks were constructed using the STRING online platform and the Cytoscape software. GeneMANIA and GSCALite online tools were used to reveal the potential role of these hub genes in breast cancer progression and TAM resistance development. And the used the GSE67916 microarray data set to verify the differentially expression of these hub genes in breast cancer. The protein expression levels of AKR1C1, AKR1C2 and AKR1C3 in TAM-sensitive and resistant breast cancer cells were compared. The TAM sensitivity of breast cancer cells with or without AKR1C1, AKR1C2 or AKR1C3 gene manipulation was evaluated by cell viability assay.

Results

A total of 184 differentially expressed genes were screened. Compared with TAM sensitive breast cancer cells, 162 were up-regulated and 22 were down-regulated. The study identified several hub genes in the PPI network that may be involved in the development of TAM resistance of breast cancer, including signal transducer and activator of transcription 1 (STAT1), estrogen receptor alpha (ESR1), fibronectin1 (FN1), cytochrome P4501B1 (CYP1B1), AKR1C1, AKR1C2, AKR1C3 and uridine diphosphate glucuronosyltransferase (UGT) 1A family genes (UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A7, UGT1A8, UGT1A9, UGT1A10). Compared with TAM-sensitive counterparts, the expression levels of AKR1C1, AKR1C2, and AKR1C3 were up-regulated in TAM-resistant breast cancer cells.

Conclusions

Overexpression of each of these three genes significantly increased the resistance of breast cancer cells to TAM treatment, while their knockdown showed opposite effects, indicating that they are potential therapeutic target for the treatment of TAM-resistant breast cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-09040-8.

Omics- and Pharmacogenomic Evidence for the Prognostic, Regulatory, and Immune-Related Roles of PBK in a Pan-Cancer Cohort

PDZ-binding kinase (PBK) is known to regulate tumor progression in some cancer types. However, its relationship to immune cell infiltration and prognosis in different cancers is unclear. This was investigated in the present study by analyzing data from TCGA, GEO, GETx, TIMER, CPTAC, GEPIA2, cBioPortal, GSCALite, PROGNOSCAN, PharmacoDB, STRING, and ENCORI databases. PBK was overexpressed in most tumors including adenocortical carcinoma (hazard ratio [HR] = 2.178, p < 0.001), kidney renal clear cell carcinoma (KIRC; HR = 1.907, p < 0.001), kidney renal papillary cell carcinoma (HR = 3.024, p < 0.001), and lung adenocarcinoma (HR = 1.255, p < 0.001), in which it was associated with poor overall survival and advanced pathologic stage. PBK methylation level was a prognostic marker in thyroid carcinoma (THCA). PBK expression was positively correlated with the levels of BIRC5, CCNB1, CDC20, CDK1, DLGAP5, MAD2L1, MELK, PLK1, TOP2A, and TTK in 32 tumor types; and with the levels of the transcription factors E2F1 and MYC, which regulate apoptosis, the cell cycle, cell proliferation and invasion, tumorigenesis, and metastasis. It was also negatively regulated by the microRNAs hsa-miR-101-5p, hsa-miR-145-5p, and hsa-miR-5694. PBK expression in KIRC, liver hepatocellular carcinoma, THCA, and thymoma was positively correlated with the infiltration of immune cells including B cells, CD4+T cells, CD8⁺ T cells, macrophages, monocytes, and neutrophils. The results of the functional enrichment analysis suggested that PBK and related genes contribute to tumor development via cell cycle regulation. We also identified 20 drugs that potentially inhibit PBK expression. Thus, PBK is associated with survival outcome in a variety of cancers and may promote tumor development and progression by increasing immune cell infiltration into the tumor microenvironment. These findings indicate that PBK is a potential therapeutic target and has prognostic value in cancer treatment.

SKA3 Serves as a Biomarker for Poor Prognosis in Kidney Renal Papillary Cell Carcinoma

Background

There is a surprising paucity of studies investigating the potential mechanism of SKA3 in the progression and prognosis of kidney renal papillary cell carcinoma (KIRP).

Methods

We used TCGA and other databases to analyze the expression, clinical value, and potential mechanisms of SKA3 in KIRP patients. We also explored therapeutic agents for KIRP through GSCALite.

Results

SKA3 mRNA expression was significantly upregulated and the area under the curve was 0.792 (95% CI 0.727–0.856). Increased SKA3 expression was related to shorter overall survival, disease-specific survival and progression-free survival. Hub genes in protein–protein interactions were CDK1, CDC20, CCNB1, CCNA2, BUB1, AURKB, BUB1B, PLK1, CCNB2, and MAD2L1, which were differentially expressed and also associated with KIRP prognosis. Gene-set enrichment analysis indicated that E2F targets, epithelial–mesenchymal transition, glycolysis, the WNT signaling pathway, and other pathways were highly enriched upon SKA3 upregulation. Gene-set variation analysis of SKA3 and its ten hub genes showed that the significant correlation of cancer-related pathways included the cell cycle, DNA damage, hormone androgen receptor, hormone estrogen receptor, PI3K/Akt, and Ras/MAPK. In addition, we found that MEK inhibitors, ie, trametinib, selumetinib, PD0325901, and RDEA119, may be feasible targeting agents for KIRP patients.

Conclusion

SKA3 might contribute to poor prognosis of KIRP through cell cycle, DNA damage, hormone androgen receptor, hormone estrogen receptor, PI3K/Akt, and RAS/MAPK. SKA3 potentially serves as a prognostic biomarker and target for KIRP.

A new prognostic risk model based on autophagy-related genes in kidney renal clear cell carcinoma

This study aimed to explore the potential role of autophagy-related genes in kidney renal clear cell carcinoma (KIRC) and develop a new prognostic-related risk model. In our research, we used multiple bioinformatics methods to perform a pan-cancer analysis of the CNV, SNV, mRNA expression, and overall survival of autophagy-related genes, and displayed the results in the form of heat maps. We then performed cluster analysis and LASSO regression analysis on these autophagy-related genes in KIRC. In the cluster analysis, we successfully divided patients with KIRC into five clusters and found that there was a clear correlation between the classification and two clinicopathological features: tumor, and stage. In LASSO regression analysis, we used 13 genes to create a new prognostic-related risk model in KIRC. The model showed that the survival rate of patients with KIRC in the high-risk group was significantly lower than that in the low-risk group, and that there was a correlation between this grouping and the patients' metastasis, tumor, stage, grade, and fustat. The results of the ROC curve suggested that this model has good prediction accuracy. The results of multivariate Cox analysis show that the risk score of this model can be used as an independent risk factor for patients with KIRC. In summary, we believe that this research provides valuable data supporting future clinical treatment and scientific research.

The Integrative Analysis of Thrombospondin Family Genes in Pan-Cancer Reveals that THBS2 Facilitates Gastrointestinal Cancer Metastasis

Recent cancer studies have found that the thrombospondin (THBS) family, including THBS1, THBS2, THBS3, THBS4, and THBS5, play vital roles in the development and progression of human cancers. However, their relationships with tumor stage, prognosis, and tumor immunity in pan-cancer have not been systematically reported. In the present study, we employed versatile public databases to assess the expression and mutations of different THBSs in pan-cancer and performed functional experiments to analyze the roles of THBS2 in gastrointestinal cancer metastasis. Our findings indicate that THBS genes are frequently mutated in various cancers and the dysregulation of THBS family members is associated with the progression of some cancers such as gastric cancer, colon cancer, and lung cancer. Further analyses indicate that THBS genes are associated with cancer hallmarks such as cell cycle and epithelial-mesenchymal transition (EMT). Importantly, thrombospondins, especially THBS1 and THBS2, are correlated with the immune cell infiltration level in gastrointestinal cancers. Our experiments further verified that THBS2 participates in tumor metastasis by enhancing EMT. Therefore, the overall analyses reveal that THBSs might offer us potential chances for tumor diagnosis and therapy.

Identification and Validation of Afatinib Potential Drug Resistance Gene BIRC5 in Non-Small Cell Lung Cancer

Introduction

Resistance to second-generation epidermal growth factor receptor-tyrosine kinase inhibitor (EGFR-TKI), afatinib, is the most significant challenge in the clinical management of non-small cell lung cancer (NSCLC), and the underlying mechanisms remain unclear.

Methods

Genomic signatures that may confer afatinib resistance in NSCLC were identified via data mining of public databases and integrative bioinformatic analyses. Furthermore, acquired afatinib-resistant lung adenocarcinoma cell lines (HCC827 AR) were established by long-term exposure under afatinib in vitro for stepwise escalation. The expression of baculovirus IAP repeat protein 5 (BIRC5) was detected by western blot, and cellular viability of HCC827 AR was determined by CCK8.

Results

Through integrative bioinformatic analyses of public datasets, overexpression of baculovirus IAP repeat protein 5 (BIRC5) was identified in both afatinib-resistant NSCLC cells and tissues, and BIRC5 overexpression was positively correlated with lymph node metastasis as well as pathological stage in NSCLC. Furthermore, NSCLC patients with BIRC5 overexpression showed poor survival outcomes. Immune infiltration analysis suggested that BIRC5 expression was significantly inversely correlated with tumor-infiltrating cell numbers and immune biomarker expression in NSCLC. The functions of genes co-expressed with BIRC5 were mainly enriched in cell cycle mitotic phase transition, double-strand break repair, and negative regulation of the cell cycle process signaling pathway. In addition, overexpression of BIRC5 protein was detected in afatinib-resistant cells by western blot, while BIRC5-expressing cells treated with BIRC5 inhibitor, YM155, were sensitive to afatinib.

Conclusions

In this study, we showed that overexpression of BIRC5 resulted in resistance to afatinib in NSCLC and BIRC5-specific inhibitors may overcome the resistant phenotype, indicating that dysregulation of the apoptotic cell death pathway may be the key mechanism underlying TKI resistance in the development of NSCLC.

Prognostic Values and Clinical Significance of S100 Family Member's Individualized mRNA Expression in Pancreatic Adenocarcinoma

Objective: Pancreatic adenocarcinoma (PAAD) is a common malignant tumor worldwide. S100 family (S100s) is wildly involved in regulating the occurrence, development, invasion, metastasis, apoptosis, and drug resistance of many malignant tumors. However, the expression pattern, prognostic value, and oncological role of individual S100s members in PAAD need to be elucidated.

Methods: The transcriptional expression levels of S100s were analyzed through the Oncomine and GEPIA, respectively. The protein levels of S100s members in PAAD were studied by Human Protein Atlas. The correlation between S100 mRNA expression and overall survival and tumor stage in PAAD patients was studied by GEPIA. The transcriptional expression correlation and gene mutation rate of S100s members in PAAD patients were explored by cBioPortal. The co-expression networks of S100s are identified using STRING and Gene MANIA to predict their potential functions. The correlation of S100s expression and tumor-infiltrating immune cells was tested by TIMER. Pathway activity and drug target analyzed by GSCALite.

Results: 13 S100s members were upregulated in PAAD tissues. 15 S100s members were associated with TP53 mutation. Expression levels of S100A3/A5/A6/A10/A11/A14/A16/B/P/Z were significantly correlated with the pathological stage. Prognosis analysis demonstrated that PAAD patients with low mRNA levels of S100A1/B/Z or high levels of S100A2/A3/A5/A10/A11/A14/A16 had a poor prognosis. Immuno-infiltration analysis showed that the mRNA levels of S100A10/A11/A14/A16 were correlated with the infiltration degree of macrophages in PAAD. Drug sensitivity analysis showed that PAAD expressing high levels of S100A2/A6/A10/A11/A13/A14/A16 maybe resistant to small molecule drugs.

Conclusion: This study identifies the clinical significance and biological functions of the S100s in PAAD, which may provide novel insights for the selection of prognostic biomarkers.

Prognostic Value of E2F Transcription Factor Expression in Pancreatic Adenocarcinoma

Background

Pancreatic adenocarcinoma (PAAD) is one of the deadliest types of cancer. In the early stages, patients often have atypical symptoms, making diagnosis difficult. The prognosis of diagnosed patients is very poor and treating PAAD is challenging. Therefore, determining reliable risk factors related to PAAD development is critical for improving patient prognosis. E2F family transcription factors (TFs) are essential regulators of DNA synthesis and cell cycle progression in eukaryotic cells, and they have been identified as prognostic biomarkers associated with multiple cancer types. However, further research is necessary to establish the prognostic relevance of these TFs in PAAD patients.

Material/Methods

We assessed PAAD patient transcriptional and outcome data using the TIMER, ONCOMINE, STRING, GEPIA, cBioPortal, Kaplan-Meier Plotter, GSCALite, and starBase databases.

Results

PAAD tumor tissues exhibited increased expression of E2F1/3/5/7/8 relative to that in normal tissues, while the expression of E2F2/3/6/8 was associated with a more advanced tumor stage. Survival analyses indicated that PAAD patients expressing higher levels of E2F1/2/3/7/8 exhibited shorter overall survival (OS) and disease-free survival (DFS) than patients expressing lower levels of these TFs. In addition, E2F4 and E2F6 overexpression was associated with poorer DFS and OS, respectively. We also found that the expression of E2Fs was significantly correlated with immune infiltrates, including CD8+ T cells, CD4+ T cells, B cells, dendritic cells, neutrophils, and macrophages.

Conclusions

Our study may provide new insights into the optimal choice of immunotherapy and promising novel targets for therapeutic intervention in PAAD patients.

The Prognostic and Immunological Value of Guanylate-Binding Proteins in Lower-Grade Glioma: Potential Markers or Not?

Seven guanylate-binding proteins (GBPs, GBP1–7), identified as a subfamily of interferon-γ-induced guanosine triphosphate hydrolases (GTPases), has been reported to be closely associated with tumor progression, metastasis, and prognosis of cancer patients in recent years. However, the expression patterns, prognostic value, immune infiltration relevance, and biological functions of GBPs in lower-grade glioma (LGG) remain elusive. In this study, by analysis and verification through multiple public data platforms, we found that GBP1, 2, 3, 4 were significantly upregulated in LGG tissues vs normal brain tissue. Analysis based on the Cox proportional hazard ratio and Kaplan–Meier plots demonstrated that the high expressions of GBP 1, 2, 3, 4 were significantly correlated with the poor prognosis of LGG patients. Correlation analysis of clinical parameters of LGG patients indicated that the expressions of GBP 1, 2, 3, 4 were significantly associated with the histological subtype and tumor histological grade of LGG. Furthermore, the correlation analysis of immune infiltration showed that the expressions of GBP1, 2, 3, 4 were significantly and positively correlated with the level of tumor immune-infiltrating cells. In particular, GBP1, 2, 3, 4 expressions were strongly correlated with the infiltration levels of monocyte, TAM, and M1/M2 macrophage, revealing their potential to regulate the polarity of macrophages. Finally, we used the GSEA method to explore the signaling pathways potentially regulated by GBP1, 2, 3, 4 and found that they were all closely associated with immune-related signaling pathways. Collectively, these findings suggested that GBP1, 2, 3, 4 were potent biomarkers to determine the prognosis and immune cell infiltration of LGG patients.

Overexpression of NREP Promotes Migration and Invasion in Gastric Cancer Through Facilitating Epithelial-Mesenchymal Transition

The identification of biomarkers and effective therapeutic targets for gastric cancer (GC), the most common cause of cancer-related deaths around the world, is currently a major focus area in research. Here, we examined the utility of Neuronal Regeneration Related Protein (NREP) as a prognostic biomarker and therapeutic target for GC. We assessed the clinical relevance, function, and molecular role of NREP in GC using bioinformatics analysis and experimental validation. Our results showed that in GC, NREP overexpression was significantly associated with a poor prognosis. Our findings also suggested that NREP may be involved in the activation of cancer-associated fibroblasts and the epithelial-mesenchymal transition (EMT), with transforming growth factor β1 mediating both processes. In addition, NREP expression showed a positive correlation with the abundance of M2 macrophages, which are potent immunosuppressors. Together, these results indicate that NREP is overexpressed in GC and affects GC prognosis. Thus, NREP could be a prognostic biomarker and therapeutic target for GC.

Comprehensive bioinformatics analyses of APOBECs family and identification of APOBEC3D as the unfavorable prognostic biomarker in clear cell renal cell carcinoma

Purpose: At present, how early screening for ccRCC is still a thorny issue for urologists. Probing the mechanisms underlying the development of ccRCC and finding relevant prognostic biomarkers remains crucial. Therefore, we systematically analyzed the APOBEC family in this study and identified APOBEC3D as a prognostic biomarker. Methods: In this study, based on the TCGA database, we systematically assessed the expression and prognosis of the APOBEC family and analyzed potential bioinformatic pathways. We then constructed nomograms to predict the prognosis of ccRCC patients better. Afterward, we further focused on APOBEC3D in our data on ccRCC specimens. The APOBEC3D should be extensively studied in ccRCC in the future. Results: The results showed that the APOBEC family showed the most significant changes in expression in ccRCC. The pathway enrichment analysis showed that APOBEC3 family members mainly regulated cytidine and cytosine-related processes. Subsequently, the Cox regression was used to construct prognostic signature, and validated in ICGC and GEO databases. Next, a nomogram was created integrating clinical parameters showing good predictive performance. Finally, we screened for APOBEC3D and found in our clinical sample that patients with high expression of APOBEC3D had a worse prognosis. Conclusion: Based on these results, APOBEC family members play important roles in the development of ccRCC, and APOBEC3D could serve as the biomarker for predicting patient prognosis.

Hypomethylation and downregulation of miR-23b-3p are associated with upregulated PLAU: a diagnostic and prognostic biomarker in head and neck squamous cell carcinoma

Background

DNA methylation and miRNA-target genes play an important part in the early development of various tumors and have been studied as tumor biomarkers. Although previous studies have reported a cluster of molecular events (such as aberrant alterations of genomics and epigenetics), little is known of the potential biomarkers for early diagnosis and prognostic evaluation in head and neck squamous cell carcinoma (HNSCC).

Methods

Multiple bioinformatics tools based on The Cancer Genome Atlas (TCGA) database and clinical samples were applied to evaluate the beneficial biomarkers in HNSCC. We focused on the role of plasminogen activator urokinase (PLAU), including diagnostic and prognostic significance, gene expression analysis, aberrant DNA methylation characteristics, interaction of miRNAs and associated signaling pathways.

Results

We found that PLAU was aberrantly upregulated in HNSCC, regardless of the mRNA or protein level. The results of receiver operating characteristic (ROC) curve and Cox regression analysis revealed that PLAU was a diagnostic and independent prognostic factor for patients with HNSCC. Hypomethylation of PLAU was closely related to poor survival in HNSCC. Additionally, miR-23b-3p was predicted to target PLAU and was significantly downregulated in HNSCC tissues. Therefore, our findings suggested that PLAU functioned as a promoter in the pathological process of HNSCC. DNA hypomethylation and downregulation of miR-23b-3p were associated with PLAU overexpression. Finally, our findings provided evidence of a significant interaction between PLAU-target and miRNAs-target pathways, indicating that miR-23b-3p suppresses malignant properties of HNSCC by targeting PLAU via Ras/MAPK and Akt/mTOR signaling pathways.

Conclusions

PLAU is overexpressed and may serve as an independent diagnostic and prognostic biomarker in HNSCC. Hypomethylation and downregulation of miR-23b-3p might account for the oncogenic role of PLAU in HNSCC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02251-w.

Transcriptomic-Based Identification of the Immuno-Oncogenic Signature of Cholangiocarcinoma for HLC-018 Multi-Target Therapy Exploration

Cholangiocarcinomas (CHOLs), hepatobiliary malignancies, are characterized by high genetic heterogeneity, a rich tumor microenvironment, therapeutic resistance, difficulty diagnosing, and poor prognoses. Current knowledge of genetic alterations and known molecular markers for CHOL is insufficient, necessitating the need for further evaluation of the genome and RNA expression data in order to identify potential therapeutic targets, clarify the roles of these targets in the tumor microenvironment, and explore novel therapeutic drugs against the identified targets. Consequently, in our attempt to explore novel genetic markers associated with the carcinogenesis of CHOL, five genes (SNX15, ATP2A1, PDCD10, BET1, and HMGA2), collectively termed CHOL-hub genes, were identified via integration of differentially expressed genes (DEGs) from relatively large numbers of samples from CHOL GEO datasets. We further explored the biological functions of the CHOL-hub genes and found significant enrichment in several biological process and pathways associated with stem cell angiogenesis, cell proliferation, and cancer development, while the interaction network revealed high genetic interactions with a number of onco-functional genes. In addition, we established associations between the CHOL-hub genes and tumor progression, metastasis, tumor immune and immunosuppressive cell infiltration, dysfunctional T-cell phenotypes, poor prognoses, and therapeutic resistance in CHOL. Thus, we proposed that targeting CHOL-hub genes could be an ideal therapeutic approach for treating CHOLs, and we explored the potential of HLC-018, a novel benzamide-linked small molecule, using molecular docking of ligand-receptor interactions. To our delight, HLC-018 was well accommodated with high binding affinities to binding pockets of CHOL-hub genes; more importantly, we found specific interactions of HLC-018 with the conserved sequence of the AT-hook DNA-binding motif of HMGA2. Altogether, our study provides insights into the immune-oncogenic phenotypes of CHOL and provides valuable information for our ongoing experimental validation.

Comprehensive landscape of the renin-angiotensin system in Pan-cancer: a potential downstream mediated mechanism of SARS-CoV-2

Background: SARS-CoV-2, the cause of the worldwide COVID-19 pandemic, utilizes the mechanism of binding to ACE2 (a crucial component of the renin-angiotensin system [RAS]), subsequently mediating a secondary imbalance of the RAS family and leading to severe injury to the host. However, very few studies have been conducted to reveal the mechanism behind the effect of SARS-CoV-2 on tumors.

Methods: Demographic data extracted from 33 cancer types and over 10,000 samples were employed to determine the comprehensive landscape of the RAS. Expression distribution, pretranscriptional and posttranscriptional regulation and posttranslational modifications (PTMs) as well as genomic alterations, DNA methylation and m6A modification were analyzed in both tissue and cell lines. The clinical phenotype, prognostic value and significance of the RAS during immune infiltration were identified.

Results: Low expression of AGTR1 was common in tumors compared to normal tissues, while very low expression of AGTR2 and MAS1 was detected in both tissues and cell lines. Differential expression patterns of ACE in ovarian serous cystadenocarcinoma (OV) and kidney renal clear cell carcinoma (KIRC) were correlated with ubiquitin modification involving E3 ligases. Genomic alterations of the RAS family were infrequent across TCGA pan-cancer program, and ACE had the highest alteration frequency compared with other members. Low expression of AGTR1 may result from hypermethylation in the promoter. Downregulation of RAS family was linked to higher clinical stage and worse survival (as measured by disease-specific survival [DSS], overall survival [OS] or progression-free interval [PFI]), especially for ACE2 and AGTR1 in KIRC. ACE-AGTR1, a classical axis of the RAS family related to immune infiltration, was positively correlated with M2-type macrophages, cancer-associated fibroblasts (CAFs) and immune checkpoint genes in most cancers.

Conclusion: ACE, ACE2, AGT and AGTR1 were differentially expressed in 33 types of cancers. PTM of RAS family was found to rely on ubiquitination. ACE2 and AGTR1 might serve as independent prognostic factors for LGG and KIRC. SARS-CoV-2 might modify the tumor microenvironment by regulating the RAS family, thus affecting the biological processes of cancer.

The Expression Patterns and Prognostic Value of the Proteasome Activator Subunit Gene Family in Gastric Cancer Based on Integrated Analysis

Increasing evidence supports that proteasome activator subunit (PSME) genes play an indispensable role in multiple tumors. The diverse expression patterns, prognostic value, underlying mechanism, and the role in the immunotherapy of PSME genes in gastric cancer (GC) have yet to be fully elucidated. We systematically demonstrated the functions of these genes in GC using various large databases, unbiased in silico approaches, and experimental validation. We found that the median expression levels of all PSME genes were significantly higher in GC tissues than in normal tissues. Our findings showed that up-regulated PSME1 and PSME2 expression significantly correlated with favorable overall survival, post-progression survival, and first progression survival in GC patients. The expression of PSME1 and PSME2 was positively correlated with the infiltration of most immune cells and the activation of anti-cancer immunity cycle steps. Moreover, GC patients with high PSME1 and PSME2 expression have higher immunophenoscore and tumor mutational burden. In addition, a receiver operating characteristic analysis suggested that PSME3 and PSME4 had high diagnostic performance for distinguishing GC patients from healthy individuals. Moreover, our further analysis indicated that PSME genes exert an essential role in GC, and the present study indicated that PSME1 and PSME2 may be potential prognostic markers for enhancing survival and prognostic accuracy in GC patients and may even act as potential biomarkers for GC patients indicating a response to immunotherapy. PSME3 may serve as an oncogene in tumorigenesis and may be a promising therapeutic target for GC. PSME4 had excellent diagnostic performance and could serve as a good diagnostic indicator for GC.

Identification and Validation of Potential Candidate Genes of Colorectal Cancer in Response to Fusobacterium nucleatum Infection

Objective: Recent investigations revealed the relationship between Fusobacterium nucleatum (Fn) infection and colorectal cancer (CRC). However, how the host genes changes contribute to CRC in response to Fn infection remains largely unknown.

Materials and methods: In the present study, we aimed to comprehensively analyze microarray data obtained from a Caco-2 infection cell model using integrated bioinformatics analysis and further identify and validate potential candidate genes in Fn-infected Caco-2 cells and CRC specimens.

Results: We identified 10 hub genes potentially involved in Fn induced tumor initiation and progression. Furthermore, we demonstrated that the expression of centrosomal protein of 55 kDa (CEP55) is significantly higher in Fn-infected Caco-2 cells. Knocking down of CEP55 could arrest the cell cycle progression and induce apoptosis in Fn-infected Caco-2 cells. The expression of CEP55 was positively correlated with the Fn amount in Fn-infected CRC patients, and these patients with high CEP55expression had an obviously poorer differentiation, worse metastasis and decreased cumulative survival rate.

Conclusion: CEP55 plays an important role in Fn-infected colon cancer cell growth and cell cycle progression and could be used as a new diagnostic and prognostic biomarker for Fn-infected CRC.

Dysregulation of ECRG4 is associated with malignant properties and of prognostic importance in human gastric cancer

BACKGROUND

We have previously characterized esophageal carcinoma-related gene 4 (ECRG4) as a novel tumor suppressor gene, which is frequently inactivated in nasopharyngeal carcinoma and breast cancer. Nevertheless, the expression status and prognostic significance of ECRG4 maintain elusive in human gastric cancer. Herein, we examined ECRG4 expression profile in gastric cancer and assessed its association with clinicopathological characteristics and patient survival.

METHODS

Online data mining, real-time RT-PCR and immunohistochemistry were employed to determined ECRG4 expression at transcriptional and protein levels in tumors vs. noncancerous tissues. Statistical analyses including the Kaplan-Meier survival analysis and the Cox hazard model were utilized to detect the impact on clinical outcome. Moreover, ECRG4 expression was silenced in gastric cancer SGC7901 cells, and cell proliferation, colony formation and invasion assays were carried out.

RESULTS

ECRG4 mRNA and protein levels were obviously downregulated in cancer tissues than noncancerous tissues. Statistical analyses demonstrated that low ECRG4 expression was found in 34.5% (58/168) of primary gastric cancer tissues, which was associated with higher histological grade (P= 0.018), lymph node metastasis (P= 0.011), invasive depth (P= 0.020), advanced tumor stage (P= 0.002) and poor overall survival (P< 0.001). Multivariate analysis showed ECRG4 expression is an independent prognostic predictor (P< 0.001). Silencing ECRG4 expression promoted gastric cancer cell growth and invasion. Western blot analysis revealed the anti-metastatic functions of ECRG4 by downregulating of E-cadherin and α-Catenin, as well as upregulating N-cadherin and Vimentin.

CONCLUSIONS

Our observations reveal that ECRG4 expression is involved in gastric cancer pathogenesis and progression, and may serve as a candidate prognostic biomarker for this disease.

Identification of potential target genes of non-small cell lung cancer in response to resveratrol treatment by bioinformatics analysis

Non-small cell lung cancer (NSCLC) is the most common type in lung cancer in the world, and it severely threatens the life of patients. Resveratrol has been reported to inhibit cancer. However, mechanisms of resveratrol inhibiting NSCLC were unclear. The aim of this study was to identify differentially expressed genes (DEGs) of NSCLC treated with resveratrol and reveal the potential targets of resveratrol in NSCLC. We obtained mRNA expression profiles of two datasets from the National Center for Biotechnology Information Gene Expression Omnibus (NCBI-GEO) and 271 DEGs were selected for further analysis. Data from STRING shown that 177 nodes and 342 edges were in the protein-protein interaction (PPI) network, and 10 hub genes (ANPEP, CD69, ITGAL, PECAM1, PTPRC, CD34, ITGA1, CCL2, SOX2, and EGFR) were identified by Cytoscape plus-in cytoHubba. Survival analysis revealed that NSCLC patients showing low expression of PECAM1, ANPEP, CD69, ITGAL, and PTPRC were associated with worse overall survival (OS) (P < 0.05), and high expression of SOX2 and EGFR was associated with worse OS for NSCLC patients (P < 0.05). Overall, we identified ANPEP, CD69, ITGAL, and PTPRC as potential candidate genes which were main effects of resveratrol on the treatment of NSCLC. ANPEP, ITGAL, CD69, and PTPRC are all clusters of differentiation (CD) antigens, might be the targets of resveratrol. The bioinformatic results suggested that the inhibitory effect of resveratrol on lung cancer may be related to the immune signaling pathway. Further studies are needed to validate these findings and to explore their functional mechanisms.

Downregulated Expression of CLEC9A as Novel Biomarkers for Lung Adenocarcinoma

Purpose

Abnormal CLEC9A expression is concerned with carcinogenesis. However, the role of CLEC9A in lung adenocarcinoma (LUAD) remains unknown. The goal of this study was to reveal the role of CLEC9A in LUAD based on bioinformatics and cellular functional experiments.

Materials and methods

Data available from The Cancer Genome Atlas (TCGA) were employed to study CLEC9A expression and mutations in LUAD. Expression and alterations of CLEC9A were analyzed using UALCAN and cBioPortal, respectively. Kaplan-Meier analysis was used to analyze the effect of CLEC9A on the survival of LUAD. Protein-protein interaction (PPI) network was built using GeneMANIA analysis. The similar genes of CLEC9A were obtained using GEPIA analysis, while co-expression genes correlated with CLEC9A were identified using LinkedOmics analysis. The effects of CLEC9A expression on immune cell infiltration was assessed. The effect of CLEC9A on the proliferation, apoptosis, cell cycle distribution, and invasion of human LUAD cells was detected in the LUAD cell line.

Results

CLEC9A was downregulated and the CLEC9A gene was often altered in LUAD. The survival of LUAD patients was correlated with the expression level of CLEC9A. The similar genes of CLEC9A were linked to functional networks involving positive regulation of interleukin-12 production, plasma membrane and CD40 receptor binding, primary immunodeficiency, intestinal immune network for IgA production, and cell adhesion molecules pathways. Cell cycle, apoptosis, EMT, and RAS/MAPK were significantly enriched pathways in positive and negative correlation genes with CLEC9A. A difference in the immune infiltration level of immune cell between the high and low CLEC9A expression groups was observed. Somatic cell copy number alternations (CNAs) of the CLEC9A, including arm-level gain and arm-level deletion, observably changed the infiltration levels of B cells, CD4+ T cells, macrophages, and neutrophils in LUAD. Except for LAG3, the expression of CD274, CTLA4, PDCD1, and TIGIT was positively correlated with the expression level of CLEC9A. After transfection, overexpression and knockdown of CLEC9A could affect the proliferation, apoptosis, cell cycle distribution, and invasion of LUAD cells.

Conclusion

CLEC9A is associated with prognosis and tumor immune microenvironment of LUAD, suggesting that CLEC9A may be considered as a novel biomarker for LUAD.

Identification of the Role and Clinical Prognostic Value of Target Genes of m6A RNA Methylation Regulators in Glioma

m6A RNA methylation regulators can regulate the growth, progression, and invasion of glioma cells by regulating their target genes, which provides a reliable support for the m6A regulator–target axes as the novel therapeutic targets and clinical prognostic signature in glioma. This study aimed to explore the role and prognostic value of m6A RNA methylation regulators and their targets. Expression profiles and clinicopathological data were obtained from the Chinese Glioma Genome Atlas (CGGA), The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and Clinical Proteome Tumor Analysis Consortium (CPTAC) datasets. Differential expression and correlation analyses were performed between normal and glioma tissues at mRNA and protein levels. Univariate Cox regression, survival, and Lasso Cox regression analyses were conducted to identify and establish the prognostic gene signature. Kaplan–Meier curve, multivariate Cox regression analysis, and ROC were utilized to evaluate the prognostic capacity of the prognostic gene signature. The correlation analysis, systematic bioinformatics analysis, and cell experiment were performed to further understand the potential underlying molecular mechanisms and drug sensitivity. Our results suggested that IGF2BP2, KIAA1429, METTL16, and METTL3, as well as 208 targets are involved in the occurrence of glioma, GBM, and LGG. YTHDF1 and 78 targets involved the occurrence of glioma and GBM, not LGG, among which 181 genes were associated with overall survival. From other findings and our cell experiment results, we demonstrated that METTL3 can activate Notch pathway and facilitate glioma occurrence through regulating its direct targets NOTCH3, DLL3, and HES1, and Notch pathway genes may serve as the potential treatment targets for glioma. Our study established and validated a seven-gene signature comprising METTL3, COL18A1, NASP, PHLPP2, TIMP1, U2AF2, and VEGFA, with a good capability for predicting glioma survival, which may guide therapeutic customization and clinical decision-making. These genes were identified to influence 81 anticancer drug responses, which further contributes to the early phase clinical trials of drug development.

MXD3 as an onco-immunological biomarker encompassing the tumor microenvironment, disease staging, prognoses, and therapeutic responses in multiple cancer types

MAX dimerization (MXD) protein 3 (MXD3) is a member of the MXD family of basic-helix-loop-helix-leucine-zipper (bHLHZ) transcription factors that plays pivotal roles in cell cycle progression and cell proliferation. However, there is insufficient scientific evidence on the pathogenic roles of MXD3 in various cancers and whether MXD3 plays a role in the immuno-oncology context of the tumor microenvironment, pathogenesis, prognosis, and therapeutic response of different tumors through certain common molecular mechanisms; thus, we saw a need to conduct the present in silico pan-cancer study. Using various computational tools, we interrogated the role of MXD3 in tumor immune infiltration, immune evasion, tumor progression, therapy response, and prognosis of cohorts from various cancer types. Our results indicated that MXD3 was aberrantly expressed in almost all The Cancer Genome Atlas (TCGA) cancer types and subtypes and was associated with the tumor stage, metastasis, and worse prognoses of various cohorts. Our results also suggested that MXD3 is associated with tumor immune evasion via different mechanisms involving T-cell exclusion in different cancer types and by tumor infiltration of immune cells in thymoma (THYM), liver hepatocellular carcinoma (LIHC), and head and neck squamous cell carcinoma (HNSC). Methylation of MXD3 was inversely associated with messenger (m)RNA expression levels and mediated dysfunctional T-cell phenotypes and worse prognoses of cohorts from different cancer types. Finally, we found that genetic alterations and oncogenic features of MXD3 were concomitantly associated with deregulation of the DBN1, RAB24, SLC34A1, PRELID1, LMAN2, F12, GRK6, RGS14, PRR7, and PFN3 genes and were connected to phospholipid transport and ion homeostasis. Our results also suggested that MXD3 expression is associated with immune or chemotherapeutic outcomes in various cancers. In addition, higher MXD3 expression levels were associated with decreased sensitivity of cancer cell lines to several mitogen-activated protein kinase kinase (MEK) inhibitors but led to increased activities of other kinase inhibitors, including Akt inhibitors. Interestingly, MXD3 exhibited higher predictive power for response outcomes and overall survival of immune checkpoint blockade sub-cohorts than three of seven standardized biomarkers. Altogether, our study strongly suggests that MXD3 is an immune-oncogenic molecule and could serve as a biomarker for cancer detection, prognosis, therapeutic design, and follow-up.

Low Levels of TRIM28-Interacting KRAB-ZNF Genes Associate with Cancer Stemness and Predict Poor Prognosis of Kidney Renal Clear Cell Carcinoma Patients

Simple Summary

This is the first report investigating the involvement of TRIM28-interacting KRAB-ZNFs in kidney cancer progression. We demonstrate a significant negative association between KRAB-ZNFs and cancer stemness followed by an attenuated immune-suppressive response and reveal the prognostic role for several KRAB-ZNFs. Our findings may help better understand the molecular basis of kidney cancer and ultimately pave the way to more appropriate prognostic tools and novel therapeutic strategies directly eradicating the dedifferentiated compartment of the tumor.

Abstract

Krüppel-associated box zinc finger (KRAB-ZNF) proteins are known to regulate diverse biological processes, such as embryonic development, tissue-specific gene expression, and cancer progression. However, their involvement in the regulation of cancer stemness-like phenotype acquisition and maintenance is scarcely explored across solid tumor types, and to date, there are no data for kidney renal clear cell cancer (KIRC). We have harnessed The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database transcriptomic data and used several bioinformatic tools (i.e., GEPIA2, GSCALite, TISIDB, GSEA, CIBERSORT) to verify the relation between the expression and genomic alterations in KRAB-ZNFs and kidney cancer, focusing primarily on tumor dedifferentiation status and antitumor immune response. Our results demonstrate a significant negative correlation between KRAB-ZNFs and kidney cancer dedifferentiation status followed by an attenuated immune-suppressive response. The transcriptomic profiles of high KRAB-ZNF-expressing kidney tumors are significantly enriched with stem cell markers and show a depletion of several inflammatory pathways known for favoring cancer stemness. Moreover, we show for the first time the prognostic role for several KRAB-ZNFs in kidney cancer. Our results provide new insight into the role of selected KRAB-ZNF proteins in kidney cancer development. We believe that our findings may help better understand the molecular basis of KIRC.

ASF1B Promotes Oncogenesis in Lung Adenocarcinoma and Other Cancer Types

Anti-silencing function 1B histone chaperone (ASF1B) is known to be an important modulator of oncogenic processes, yet its role in lung adenocarcinoma (LUAD) remains to be defined. In this study, an integrated assessment of The Cancer Genome Atlas (TCGA) and genotype-tissue expression (GTEx) datasets revealed the overexpression of ASF1B in all analyzed cancer types other than LAML. Genetic, epigenetic, microsatellite instability (MSI), and tumor mutational burden (TMB) analysis showed that ASF1B was regulated by single or multiple factors. Kaplan-Meier survival curves suggested that elevated ASF1B expression was associated with better or worse survival in a cancer type-dependent manner. The CIBERSORT algorithm was used to evaluate immune microenvironment composition, and distinct correlations between ASF1B expression and immune cell infiltration were evident when comparing tumor and normal tissue samples. Gene set enrichment analysis (GSEA) indicated that ASF1B was associated with proliferation- and immunity-related pathways. Knocking down ASF1B impaired the proliferation, affected cell cycle distribution, and induced cell apoptosis in LUAD cell lines. In contrast, ASF1B overexpression had no impact on the malignant characteristics of LUAD cells. At the mechanistic level, ASF1B served as an indirect regulator of DNA Polymerase Epsilon 3, Accessory Subunit (POLE3), CDC28 protein kinase regulatory subunit 1(CKS1B), Dihydrofolate reductase (DHFR), as established through proteomic profiling and Immunoprecipitation-Mass Spectrometry (IP-MS) analyses. Overall, these data suggested that ASF1B serves as a tumor promoter and potential target for cancer therapy and provided us with clues to better understand the importance of ASF1B in many types of cancer.

Upregulation of the APOBEC3 Family Is Associated with a Poor Prognosis and Influences Treatment Response to Raf Inhibitors in Low Grade Glioma

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) has been identified as a group of enzymes that catalyze cytosine deamination in single-stranded (ss) DNA to form uracil, causing somatic mutations in some cancers. We analyzed the APOBEC3 family in 33 TCGA cancer types and the results indicated that APOBEC3s are upregulated in multiple cancers and strongly correlate with prognosis, particularly in low grade glioma (LGG). Then we constructed a prognostic model based on family expression in LGG where the APOBEC3 family signature is an accurate predictive model (AUC of 0.85). Gene mutation, copy number variation (CNV), and a differential gene expression (DEG) analysis were performed in different risk groups, and the weighted gene co-expression network analysis (WGCNA) was employed to clarify the role of various members in LGG; CIBERSORT algorithm was deployed to evaluate the landscape of LGG immune infiltration. We found that upregulation of the APOBEC3 family expression can strengthen Ras/MAPK signaling pathway, promote tumor progression, and ultimately reduce the treatment benefits of Raf inhibitors. Moreover, the APOBEC3 family was shown to enhance the immune response mediated by myeloid cells and interferon gamma, as well as PD-L1 and PD-L2 expression, implying that they have immunotherapy potential. Therefore, the APOBEC3 signature enables an efficient assessment of LGG patient survival outcomes and expansion of clinical benefits by selecting appropriate individualized treatment strategies.

TBC1D3 family is a prognostic biomarker and correlates with immune infiltration in kidney renal clear cell carcinoma

The TBC1D3 family is overexpressed in many cancers, including kidney renal clear cell carcinoma (KIRC), which is associated with tumor-infiltrating lymphocytes. However, the expression and prognosis of TBC1D3 family and tumor-infiltrating lymphocytes in KIRC remain unknown. In the present study, we systematically explored and validated the expression and prognostic value of TBC1D3 family expression in KIRC using multiple public databases. In addition, the function of the TBC1D3 family members and the correlations between TBC1D3 family expression and KIRC immune infiltration levels were investigated. We found that TBC1D3 family members were rarely mutated (less than 5 frequencies). TBC1D3 family was overexpressed in KIRC; high expression of the TBC1D3 family members was correlated with poor prognosis. In addition, TBC1D3D may positively regulate proliferation, and overexpression of TBC1D3 promoted clear cell renal cell carcinoma proliferation in vitro. In terms of immune infiltrating levels, TBC1D3 family expression was positively associated with CD4⁺ T cells infiltrating levels. These findings suggest that the TBC1D3 family expression is correlated with prognosis and immune infiltrating levels. Therefore, the TBC1D3 family can be used as a biomarker for KIRC and a prognostic biomarker for determining the prognosis and immune infiltration levels in KIRC.

A New Survival Model Based on ADAMTSs for Prognostic Prediction in Clear Cell Renal Cell Carcinoma

The main purpose of this study was to explore the genetic variation, gene expression, and clinical significance of ADAMTSs (a disintegrin and metalloprotease domains with thrombospondin motifs) across cancer types. Analysis of data from the TCGA (The Cancer Genome Atlas) database showed that the ADAMTSs have extensive CNV (copy number variation) and SNV (single nucleotide variation) across cancer types. Compared with normal tissues, the methylation of ADAMTSs in cancer tissues is also significantly different, which affects the expression of ADAMTS gene and the prognosis of cancer patients. Through gene expression analysis, we found that ADAMTS family has significant changes in gene expression across cancer types and is closely related to the prognosis of carcinoma, especially in ccRCC (clear cell renal cell carcinoma). LASSO regression analysis was used to establish a prognostic model based on the ADAMTSs to judge the prognosis of patients with ccRCC. Multiple Cox regression analysis suggested that age, grade, stage, and risk score of the prognostic model of ccRCC were independent prognostic factors in patients with renal clear cell carcinoma. These findings indicate that the ADAMTSs-based survival model can accurately predict the prognosis of patients with ccRCC and suggest that ADAMTSs are a potential prognostic biomarker and therapeutic target in ccRCC.

Umbrella Review on Associations Between Single Nucleotide Polymorphisms and Lung Cancer Risk

The aim is to comprehensively and accurately assess potential relationships between single nucleotide polymorphisms (SNP) and lung cancer (LC) risk by summarizing the evidence in systematic reviews and meta-analyses. This umbrella review was registered with the PROSPERO international prospective register of systematic reviews under registration number CRD42020204685. The PubMed, Web of Science, and Embase databases were searched to identify eligible systematic reviews and meta-analyses from inception to August 14, 2020. The evaluation of cumulative evidence was conducted for associations with nominally statistical significance based on the Venice criteria and false positive report probability (FPRP). This umbrella review finally included 120 articles of a total of 190 SNP. The median number of studies and sample size included in the meta-analyses were five (range, 3–52) and 4 389 (range, 354–256 490), respectively. A total of 85 SNP (in 218 genetic models) were nominally statistically associated with LC risk. Based on the Venice criteria and FPRP, 13 SNP (in 22 genetic models), 47 SNP (in 99 genetic models), and 55 SNP (in 94 genetic models) had strong, moderate, and weak cumulative evidence of associations with LC risk, respectively. In conclusion, this umbrella review indicated that only 13 SNP (of 11 genes and one miRNA) were strongly correlated to LC risk. These findings can serve as a general and helpful reference for further genetic studies.

A Comprehensive Multiomics Analysis Identified Ubiquilin 4 as a Promising Prognostic Biomarker of Immune-Related Therapy in Pan-Cancer

Recently, it was reported that ubiquilin 4 (UBQLN4) alteration was associated with genomic instability in some cancers. However, whether UBQLN4 is a valuable biomarker for the prognosis of immunotherapy in pan-cancer was not identified. We evaluated the biologic and oncologic significance of UBQLN4 in pan-cancer at multiomics level, such as expression, mutation, copy number variation (CNV), methylation, and N⁶-methyladenosine (m⁶A) methylation. These omics data were obtained from several public databases, including Oncomine, The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), the Genotype-Tissue Expression (GTEx), the Human Protein Atlas (HPA), Gene Set Cancer Analysis (GSCA), m⁶A-Atlas, CancerSEA, and RNAactDrug. We found that UBQLN4 mRNA and protein were overexpressed in most cancer types, and the expression, mutation, CNV, and methylation of UBQLN4 were associated with the prognosis of some cancers. Mechanistically, UBQLN4 was involved in angiogenesis, DNA damage, apoptosis, and the pathway of PI3K/AKT and TSC/mTOR. Moreover, UBQLN4 mRNA was significantly correlated with immune checkpoints, tumor mutational burden (TMB), microsatellite instability (MSI), and mismatch repair (MMR). And, the correlation among UBQLN4 mRNA, CNV, and methylation and immune microenvironment was also identified. Furthermore, UBQLN4 was associated with the sensitivity of chemotherapy and targeted drugs at multiomics level. In conclusion, UBQLN4 was a promising prognostic biomarker of immune-related therapy in pan-cancer.

TLR3 Serves as a Prognostic Biomarker and Associates with Immune Infiltration in the Renal Clear Cell Carcinoma Microenvironment

Background

Clear cell renal cancer (KIRC) is one of the most common cancers globally, with a poor prognosis. TLRs play a vital role in anticancer immunity and the regulation of the biological progress of tumour cells. However, the precise role of TLRs in KIRC is still ambiguous.

Methods

Various bioinformatics analysis and clinical validation of tissues were performed to evaluate the prognostic value of TLRs and their correlation with immune infiltration in KIRC.

Results

The expression of TLR2/3/7/8 was increased at both mRNA and protein levels in KIRC. TLRs in KIRC were involved in the activation of apoptosis, EMT, RAS/MAPK, and RTK pathways, as well as the inhibition of the cell cycle and the hormone AR pathway. Drug sensitivity analysis revealed that high expression of TLR3 and low expression of TLR7/9/10 were resistant to most of the small molecules or drugs from CTRP. Enrichment analyses showed that TLRs were mainly involved in innate immune response, toll-like receptor signalling pathway, NF-kappa B signalling pathway, and TNF signalling pathway. Furthermore, a high-level TLR3 expression was associated with a favourable prognosis in KIRC. Validation research further confirmed that TLR3 expression was increased in KIRC tissues, and high TLR3 levels were associated with poor overall survival. Moreover, TLR3 in KIRC showed a positive association with an abundance of immune cells, including B-cells, CD4+ T-cells, CD8+ T-cells, macrophage, neutrophils, and dendritic cells, and the expression of the immune biomarker sets. Several TLR3-associated kinase, miRNA, or transcription factor targets were also identified in KIRC.

Conclusion

Our results indicate that TLR3 serves as a prognostic biomarker and associated with immune infiltration in KIRC. This work lays a foundation for further studies on the role of TLR3 in the carcinogenesis and progression of KIRC.

A Novel Intercellular Communication-Associated Gene Signature for Prognostic Prediction and Clinical Value in Patients With Lung Adenocarcinoma

Background

Lung cancer remains the leading cause of cancer death globally, with lung adenocarcinoma (LUAD) being its most prevalent subtype. This study aimed to identify the key intercellular communication-associated genes (ICAGs) in LUAD.

Methods

Eight publicly available datasets were downloaded from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. The prognosis-related ICAGs were identified and a risk score was developed by using survival analysis. Machine learning models were trained to predict LUAD recurrence based on the selected ICAGs and clinical information. Comprehensive analyses on ICAGs and tumor microenvironment were performed. A single-cell RNA-sequencing dataset was assessed to further elucidate aberrant changes in intercellular communication.

Results

Eight ICAGs with prognostic potential were identified in the present study, and a risk score was derived accordingly. The best machine-learning model to predict relapse was developed based on clinical information and the expression levels of these eight ICAGs. This model achieved a remarkable area under receiver operator characteristic curves of 0.841. Patients were divided into high- and low-risk groups according to their risk scores. DNA replication and cell cycle were significantly enriched by the differentially expressed genes between the high- and the low-risk groups. Infiltrating immune cells, immune functions were significantly related to ICAGs expressions and risk scores. Additionally, the changes of intercellular communication were modeled by analyzing the single-cell sequencing dataset.

Conclusion

The present study identified eight key ICAGs in LUAD, which could contribute to patient stratification and act as novel therapeutic targets.

Comprehensive Analysis of RUNX and TGF-β Mediated Regulation of Immune Cell Infiltration in Breast Cancer

Runt-related transcription factors (RUNXs) can serve as both transcription activators and repressors during biological development, including immune cell maturation. RUNX factors have both tumor-promoting and tumor-suppressive roles in carcinogenesis. Immune cell infiltration and the tumor immune microenvironment have been found to be key regulators in breast cancer progression, treatment response, and patient outcome. However, the relationship between the RUNX family and immune cell infiltration in breast cancer remains unclear. We performed a comprehensive analysis to reveal the role of RUNX factors in breast cancer. Analysis of patient data in the Oncomine database showed that the transcriptional levels of RUNX proteins in breast cancer were elevated. Kaplan–Meier plotter (KM plotter) analysis showed that breast cancer patients with higher expression of RUNX proteins had better survival outcomes. Through analysis of the UALCAN database, we found that the transcriptional levels of RUNX factors were significantly correlated with some breast cancer patient characteristics. cBio Cancer Genomics Portal (cBioPortal) analysis showed the proportions of different RUNX genomic alterations in various subclasses of breast cancer. We also performed gene ontology (GO) and pathway analyses for the significantly differentially expressed genes that were correlated with RUNX factors in breast cancer. TIMER database analysis showed that immune cell infiltration in breast cancer could be affected by the transcriptional level, mutation, and gene copy number of RUNX proteins. Using the Gene Set Cancer Analysis (GSCA) database, we analyzed the effects of RUNX gene methylation on the level of immune cell infiltration in breast cancer. We found that the methylation level changes of RUNX2 and RUNX3 had opposite effects on immune cell infiltration in breast cancer. We also analyzed the relationship between the methylation level of RUNX genes and the TGF-β signaling pathway using the TISIDB database. The results showed that the methylation levels of RUNX1 and RUNX3 were correlated with the expression of TGF-β1. In summary, our analysis found that the RUNX family members can influence the infiltration of various immune cells in breast cancer depending on their expression level, mutation, gene copy number, and methylation. The RUNX family is an important regulator of immune cell infiltration in breast cancer and may serve as a potential prognostic biomarker.

A pan-cancer analysis of the expression of gasdermin genes in tumors and their relationship with the immune microenvironment

BACKGROUND

Gasdermins (GSDMs) are a class of proteins related to pyrolysis and in humans, consist of GSDMA, GSDMB, GSDMC, GSDMD, DFNA5, and DFNB59. The inflammatory factors and cell contents released during pyrolysis can recruit immune cells and change the microenvironment. However, to date, there is a paucity of studies examining the relationship between GSDMs and the immune microenvironment in tumors. Therefore, this current report analyzed the expression of GSDM genes in tumors and their relationship with the immune microenvironment.

METHODS

Apply GSCALite and GEPIA2 online analysis tools to analyze the gene expression levels and the Single nucleotide variant (SNV), copy number variation (CNV), and methylation characteristics of GSDM genes respectively. Use R software or TISIDB online analysis tool to carry out the correlation analysis required in the article. Furthermore, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted to examine the role of these GSDM genes in various cancers.

RESULTS

The results demonstrated that CNV can cause an increase in GSDM gene expression, and methylation can inhibit GSDM gene expression. The elevated expression of GSDMA, GSDMB, GSDMC, GSDMD, and DFNA5 in some or most tumors was often accompanied by elevated immune scores, increased immune cell infiltration, and high expression of major histocompatibility complex (MHC) molecules, chemokines and their receptors, and immune checkpoint-related genes. However, DFNB59 was often negatively correlated with these indicators in tumors. GSDMD was the most highly expressed GSDM protein in various normal tissues and tumors, and showed the strongest correlation with immune microenvironment-related genes. Moreover, the methylation of GSDMD was accompanied by low immune cell infiltration, low expression of MHC molecule-related genes, low expression of chemokines and receptor-related genes, and low expression of immune checkpoint-related genes.

CONCLUSIONS

Therefore, the expression of GSDM-related genes is associated with the tumor immune microenvironment. The GSDM genes, especially GSDMD, may be used as therapeutic targets to predict or change the tumor microenvironment and as biomarkers to predict the therapeutic efficacy of immune checkpoint inhibitors.

Glutamine Metabolism Regulators Associated with Cancer Development and the Tumor Microenvironment: A Pan-Cancer Multi-Omics Analysis

BACKGROUND

In recent years, metabolic reprogramming has been identified as a hallmark of cancer. Accumulating evidence suggests that glutamine metabolism plays a crucial role in oncogenesis and the tumor microenvironment. In this study, we aimed to perform a systematic and comprehensive analysis of six key metabolic node genes involved in the dynamic regulation of glutamine metabolism (referred to as GLNM regulators) across 33 types of cancer.

METHODS

We analyzed the gene expression, epigenetic regulation, and genomic alterations of six key GLNM regulators, including SLC1A5, SLC7A5, SLC3A2, SLC7A11, GLS, and GLS2, in pan-cancer using several open-source platforms and databases. Additionally, we investigated the impacts of these gene expression changes on clinical outcomes, drug sensitivity, and the tumor microenvironment. We also attempted to investigate the upstream microRNA-mRNA molecular networks and the downstream signaling pathways involved in order to uncover the potential molecular mechanisms behind metabolic reprogramming.

RESULTS

We found that the expression levels of GLNM regulators varied across cancer types and were related to several genomic and immunological characteristics. While the immune scores were generally lower in the tumors with higher gene expression, the types of immune cell infiltration showed significantly different correlations among cancer types, dividing them into two clusters. Furthermore, we showed that elevated GLNM regulators expression was associated with poor overall survival in the majority of cancer types. Lastly, the expression of GLNM regulators was significantly associated with PD-L1 expression and drug sensitivity.

CONCLUSIONS

The elevated expression of GLNM regulators was associated with poorer cancer prognoses and a cold tumor microenvironment, providing novel insights into cancer treatment and possibly offering alternative options for the treatment of clinically refractory cancers.

Pan-Cancer Analysis of Immune Complement Signature C3/C5/C3AR1/C5AR1 in Association with Tumor Immune Evasion and Therapy Resistance

Despite the advances in our understanding of the genetic and immunological basis of cancer, cancer remains a major public health burden with an ever-increasing incidence rate globally. Nevertheless, increasing evidence suggests that the components of the complement system could regulate the tumor microenvironment (TME) to promote cancer progression, recurrence, and metastasis. In the present study, we used an integrative multi-omics analysis of clinical data to explore the relationships between the expression levels of and genetic and epigenetic alterations in C3, C5, C3AR1, and C5AR1 and tumor immune evasion, therapy response, and patient prognosis in various cancer types. We found that the complements C3, C5, C3AR1, and C5AR1 have deregulated expression in human malignancies and are associated with activation of immune-related oncogenic processes and poor prognosis of cancer patients. Furthermore, we found that the increased expression levels of C3, C5, C3AR1, and C5AR1 were primarily predicted by copy number variation and gene methylation and were associated with dysfunctional T-cell phenotypes. Single nucleotide variation in the gene signature co-occurred with multiple oncogenic mutations and is associated with the progression of onco-immune-related diseases. Further correlation analysis revealed that C3, C5, C3AR1, and C5AR1 were associated with tumor immune evasion via dysfunctional T-cell phenotypes with a lesser contribution of T-cell exclusion. Lastly, we also demonstrated that the expression levels of C3, C5, C3AR1, and C5AR1 were associated with context-dependent chemotherapy, lymphocyte-mediated tumor killing, and immunotherapy outcomes in different cancer types. In conclusion, the complement components C3, C5, C3AR1, and C5AR1 serve as attractive targets for strategizing cancer immunotherapy and response follow-up.

Identification and Characterization of Alcohol-related Hepatocellular Carcinoma Prognostic Subtypes based on an Integrative N6-methyladenosine methylation Model

Background: Alcohol consumption increases the risk of hepatocellular carcinoma (HCC), and associated with a high mortality rate and poor prognosis. N6-methyladenosine (m6A) methylations play key roles in tumorigenesis and progression. However, our current knowledge about m6A in alcohol-related HCC (A-HCC) remains elucidated. Herein, the authors construct an integrative m6A model based on A-HCC subtyping and mechanism exploration workflow. Methods: Based on the m6A expressions of A-HCC and in vivo experiment, different prognosis risk A-HCC subtypes are identified. Meanwhile, multiple interdependent indicators of prognosis including patient survival rate, clinical pathological prognosis and immunotherapy sensitivity. Results: The m6A model includes LRPPRC, YTHDF2, KIAA14219, and RBM15B, classified A-HCC patients into high/low-risk subtypes. The high-risk subtype compared to the low-risk subtype showed phenotypic malignancy, poor prognosis, immunosuppression, and activation of tumorigenesis and proliferation-related pathways, including the E2F target, DNA repair, and mTORC1 signalling pathways. The expression of Immunosuppressive cytokines DNMT1/EZH2 was up-regulated in A-HCC patients, and teniposide may be a potential therapeutic drug for A-HCC. Conclusion: Our model redefined A-HCC prognosis risk, identified potential m6As linking tumour progress and immune regulations and selected possible therapy target, thus promoting understanding and clinical applications about A-HCC.