Increased expression of TNFRSF14 indicates good prognosis and inhibits bladder cancer proliferation by promoting apoptosis

DoSurvive: A webtool for investigating the prognostic power of a single or combined cancer biomarker

We present DoSurvive, a user-friendly survival analysis web tool and a cancer prognostic biomarker centered database. DoSurvive is the first database that allows users to perform multivariant survival analysis for cancers with customized gene/patient list. DoSurvive offers three survival analysis methods, Log rank test, Cox regression and accelerated failure time model (AFT), for users to analyze five types of quantitative features (mRNA, miRNA, lncRNA, protein and methylation of CpG islands) with four survival types, i.e. overall survival, disease-specific survival, disease-free interval, and progression-free interval, in 33 cancer types. Notably, the implemented AFT model provides an alternative method for genes/features which failed the proportional hazard assumption in Cox regression. With the unprecedented number of survival models implemented and high flexibility in analysis, DoSurvive is a unique platform for the identification of clinically relevant targets for cancer researcher and practitioners. DoSurvive is freely available at http://dosurvive.lab.nycu.edu.tw/.

Herpesvirus entry mediator as a potential biomarker in breast cancer compared with conventional cytotoxic T‑lymphocyte‑associated antigen 4

Breast cancer (BC) is the most common cancer in women worldwide, with 2.3 million cases recorded in 2020. Despite improvements in cancer treatment, patients with BC still succumb to the disease, due to regional and distant metastases when diagnosed at later stages. Several immune checkpoint inhibitors have been approved for BC treatment, based on their expression and role in maintaining immunosurveillance against tumors. The present study aimed to evaluate the expression of 12 immune checkpoints in patients with BC, and assess their role as diagnostic and therapeutic markers. Expression levels were measured using reverse transcription-quantitative polymerase chain reaction. Among the 12 immune markers, herpesvirus entry mediator (HVEM) was found to be significantly upregulated in patients with malignant BC compared to non-malignant controls, with a relative fold change (FC) of 1.46 and P=0.012. A similar finding was observed for cytotoxic T-lymphocyte-associated antigen 4 (CTLA4; FC=1.47 and P=0.035). In addition, receiver operating characteristic curve analysis revealed that HVEM expression allowed significant differentiation between groups, with an area under the curve of 0.74 (P=0.013). Upregulation in both HVEM and CTLA4 was revealed to be significantly associated with the human epidermal growth factor receptor-2 (HER2)-enriched phenotype (FC=3.53, P=0.009 and FC=5.98, P=0.002, respectively), while only HVEM was significantly associated with the triple-negative phenotype (FC=2.07, P=0.016). Furthermore, HVEM was significantly higher in patients with grade III tumors (FC=1.88, P=0.025) and negative vascular invasion (FC=1.67, P=0.046) compared with non-malignant controls. Serum protein levels were assessed by multiplex immunoassay, and a significant increase in HVEM was detected in patients with malignant BC compared with that in non-malignant controls (P=0.035). These data indicated that HVEM may serve as a potential biomarker and target for immunotherapy, especially for certain types of BC.

Cuproptosis correlates with immunosuppressive tumor microenvironment based on pan-cancer multiomics and single-cell sequencing analysis

Recent studies suggest that cuproptosis, a novel mode of cell death, may be associated with the development of cancer. However, no studies are showing its role in tumorigenesis, progression, and prognosis. In the present study, we comprehensively analyzed the expression difference, gene variation and methylation modification of cuproptosis-related genes (CRGs) in pan-cancer. Then, Single sample gene set enrichment analysis (ssGSEA) was used to calculate individual cuproptosis scores (CS). The association of CS with copy number variation, clinical features, immune-related genes, TMB, MSI, and tumor immune dysfunction and exclusion (TIDE) was comprehensively assessed. Single-cell transcriptome sequencing (scRNA-seq) to analyze the activation of cuproptosis in the tumor microenvironment. Immunohistochemistry (IHC) were used to validate the expression of cuproptosis hub-gene. Our study shows that CRGs were significantly expressed in a variety of tumors, and CDKN2A had the highest mutation frequency (49%) in all tumors. A significant increase in the CS was observed in most cancers and were associated with poor prognosis in the majority of tumors. CS was significantly negatively correlated with tumor microenvironment scores in more than 10 tumors and positively correlated with PD-L1 in 11 tumors, suggesting involvement in tumor immune escape. scRNA-seq suggests that CRG scores significantly increased in the cancer cells. This study opens avenues for further research on the role of cuproptosis in the occurrence and development of cancer and the development of targeted therapies based on cuproptosis.

The role of costimulatory molecules in glioma biology and immune microenvironment

Background: Extensive research showed costimulatory molecules regulate tumor progression. Nevertheless, a small amount of literature has concentrated on the potential prognostic and therapeutic effects of costimulatory molecules in patients with glioma. Methods: The data were downloaded from The Cancer Genome Atlas (TCGA) database, Chinese Glioma Genome Atlas (CGGA) database, and Gene Expression Omnibus (GEO) database for bioinformatics analysis. R software was applied for statistical analysis. Using the FigureYa and Xiantao online tools (https://www.xiantao.love/) for mapping. Results: The Least absolute shrinkage and selection operator (LASSO) and Cox regression analysis were utilized to identify the signature consisting of five costimulatory molecules. Multivariate regression analysis revealed that the prognosis of glioma could be independently predicted by the riskscore. Furthermore, we explored clinical and genomic feature differences between the two groups. The level of tumor mutational burden (TMB) was higher in the high-risk group, while more mutation of IDH1 was observed in the low-risk group. Results of Tumor Immune Dysfunction and Exclusion (TIDE) analysis showed that high-risk patients were more prone to be responded to immunotherapy. In addition, subclass mapping analysis was performed to validate our findings and the results revealed that a significantly higher percentage of immunotherapy response rate was observed in the high-risk group. Conclusion: A novel signature with a good independent predictive capacity of prognosis was successfully identified. And our findings reveal that patients with high-risk scores were more likely to be responded to immunotherapy.

A novel chromatin regulator-related immune checkpoint related gene prognostic signature and potential candidate drugs for endometrial cancer patients

Background

Endometrial cancer (EC) is the most common gynecologic malignancy in developed countries and its prevalence is increasing. As an emerging therapy with a promising efficacy, immunotherapy has been extensively applied in the treatment of solid tumors. In addition, chromatin regulators (CRs), as essential upstream regulators of epigenetics, play a significant role in tumorigenesis and cancer development.

Methods

CRs and immune checkpoint-related genes (ICRGs) were obtained from the previous top research. The Genome Cancer Atlas (TCGA) was utilized to acquire the mRNA expression and clinical information of patients with EC. Correlation analysis was utilized for screen CRs-related ICRGs (CRRICRGs). By Cox regression and least absolute shrinkage and selection operator (LASSO) analysis, prognosis related CRRICRGs were screened out and risk model was constructed. The Kaplan–Meier curve was used to estimate the prognosis between high- and low-risk group. By comparing the IC50 value, the drugs sensitivity difference was explored. We obtained small molecule drugs for the treatment of UCEC patients based on CAMP dataset.

Results

We successfully constructed a 9 CRRICRs-based prognostic signature for patients with UCEC and found the riskscore was an independent prognostic factor. The results of functional analysis suggested that CRRICRGs may be involved in immune processes associated with cancer. Immune characteristics analysis provided further evidence that the CRRICRGs-based model was correlated with immune cells infiltration and immune checkpoint. Eight small molecule drugs that may be effective for the treatment of UCEC patients were screened. Effective drugs identified by drug sensitivity profiling in high- and low-risk groups.

Conclusion

In summary, our study provided novel insights into the function of CRRICRGs in UCEC. We also developed a reliable prognostic panel for the survival of patients with UCEC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-022-00253-w.

The Pyroptosis-Related Risk Genes APOBEC3D, TNFRSF14, and RAC2 Were Used to Evaluate Prognosis and as Tumor Suppressor Genes in Breast Cancer

Background

Pyroptosis is a type of cell death that plays an important role in predicting prognosis and immunoregulation in cancers. However, the pyroptosis-related gene signature for prognosis and immune infiltration prediction has not been studied in breast cancer (BC).

Methods

The Gene Expression Omnibus (GEO) and Cancer Genome Atlas (TCGA) databases were used to obtain the expression and clinical data of genes. 52 pyroptosis-related genes were obtained from TCGA-BC and estimated differentially expressed genes by the limma program. To categorize the molecular subtypes of pyroptosis-related genes, the ConsensusClusterPlus tool was utilized. Cox and Lasso regression analyses were used to create a signature. TCGA-BC dataset as the training set and the GSE37751 test set for risk research. Gene set enrichment analysis (GSEA) was used to conduct KEGG and GO studies of subtype groups. We also used the ssGSEA approach in the GSVA package to calculate the risk score of immune cells. Finally, pyroptosis-related genes in BC were validated using qPCR and immunohistochemical assays. Clone formation and EDU assays were used to explore the ability of signature genes to regulate the proliferation of BC cells.

Results

Based on pyroptosis-related genes, the C1 and C2 subtypes were obtained. Survival analysis results showed that the C2 group had a better prognosis. Then, a three-gene signature (APOBEC3D, TNFRSF14, and RAC2) were created by Lasso regression analysis, which had a good prediction effect in the TCGA-BC and GSE37751 datasets. Our nomogram has a fair degree of accuracy in predicting the survival rates of BC patients. The pyroptosis-related signature has a good predictive effect in evaluating the tumour microenvironment score, 28 types of immune cells and response to immune checkpoint therapy. Finally, qPCR and immunohistochemistry staining results indicated that APOBEC3D, TNFRSF14, and RAC2 expression in BC tissues was low. The results of clone formation and EdU assays showed that high expression of signature genes inhibited the proliferation ability of BC cells.

Conclusions

Based on pyroptosis-related genes (APOBEC3D, TNFRSF14, and RAC2), we built a novel prognostic molecular model for BC that might be used to assess prognostic risk and immune infiltration in BC patients. These signature genes are also tumor suppressor genes and may serve as potential targets for BC.

Single-Cell Sequencing Unveils the Heterogeneity of Nonimmune Cells in Chronic Apical Periodontitis

Chronic apical periodontitis (CAP) is a unique dynamic interaction between microbial invasions and host defense mechanisms, resulting in infiltration of immune cells, bone absorption, and periapical granuloma formation. To help to understand periapical tissue pathophysiology, we constituted a single-cell atlas for 26,737 high-quality cells from inflammatory periapical tissue and uncovered the complex cellular landscape. The eight types of cells, including nonimmune cells and immune cells, were identified in the periapical tissue of CAP. Considering the key roles of nonimmune cells in CAP, we emphasized osteo-like cells, basal/stromal cells, endothelial cells, and epithelial cells, and discovered their diversity and heterogeneity. The temporal profiling of genomic alterations from common CAP to typical periapical granuloma provided predictions for transcription factors and biological processes. Our study presented potential clues that the shift of inflammatory cytokines, chemokines, proteases, and growth factors initiated polymorphic cell differentiation, lymphangiogenesis, and angiogenesis during CAP.

A Transcriptional Signature of IL-2 Expanded Natural Killer Cells Predicts More Favorable Prognosis in Bladder Cancer

Activation of natural killer (NK) cell function is regulated by cytokines, such as IL-2, and secreted factors upregulated in the tumor microenvironment, such as platelet-derived growth factor D (PDGF-DD). In order to elucidate a clinical role for these important regulators of NK cell function in antitumor immunity, we generated transcriptional signatures representing resting, IL-2-expanded, and PDGF-DD-activated, NK cell phenotypes and established their abundance in The Cancer Genome Atlas bladder cancer (BLCA) dataset using CIBERSORT. The IL-2-expanded NK cell phenotype was the most abundant in low and high grades of BLCA tumors and was associated with improved prognosis. In contrast, PDGFD expression was associated with numerous cancer hallmark pathways in BLCA tumors compared with normal bladder tissue, and a high tumor abundance of PDGFD transcripts and the PDGF-DD-activated NK cell phenotype were associated with a poor BLCA prognosis. Finally, high tumor expression of transcripts encoding the activating NK cell receptors, KLRK1 and the CD160-TNFRSF14 receptor-ligand pair, was strongly correlated with the IL-2-expanded NK cell phenotype and improved BLCA prognosis. The transcriptional parameters we describe may be optimized to improve BLCA patient prognosis and risk stratification in the clinic and potentially provide gene targets of therapeutic significance for enhancing NK cell antitumor immunity in BLCA.

Identification of Immune-Related Genes Associated With Bladder Cancer Based on Immunological Characteristics and Their Correlation With the Prognosis

BACKGROUND

To identify the immune-related genes of bladder cancer (BLCA) based on immunological characteristics and explore their correlation with the prognosis.

METHODS

We downloaded the gene and clinical data of BLCA from the Cancer Genome Atlas (TCGA) as the training group, and obtained immune-related genes from the Immport database. We downloaded GSE31684 and GSE39281 from the Gene Expression Omnibus (GEO) as the external validation group. R (version 4.0.5) and Perl were used to analyze all data.

RESULT

Univariate Cox regression analysis and Lasso regression analysis revealed that 9 prognosis-related immunity genes (PIMGs) of differentially expressed immune genes (DEIGs) were significantly associated with the survival of BLCA patients (p < 0.01), of which 5 genes, including NPR2, PDGFRA, VIM, RBP1, RBP1 and TNC, increased the risk of the prognosis, while the rest, including CD3D, GNLY, LCK, and ZAP70, decreased the risk of the prognosis. Then, we used these genes to establish a prognostic model. We drew receiver operator characteristic (ROC) curves in the training group, and estimated the area under the curve (AUC) of 1-, 3- and 5-year survival for this model, which were 0.688, 0.719, and 0.706, respectively. The accuracy of the prognostic model was verified by the calibration chart. Combining clinical factors, we established a nomogram. The ROC curve in the external validation group showed that the nomogram had a good predictive ability for the survival rate, with a high accuracy, and the AUC values of 1-, 3-, and 5-year survival were 0.744, 0.770, and 0.782, respectively. The calibration chart indicated that the nomogram performed similarly with the ideal model.

CONCLUSION

We had identified nine genes, including PDGFRA, VIM, RBP1, RBP1, TNC, CD3D, GNLY, LCK, and ZAP70, which played important roles in the occurrence and development of BLCA. The prognostic model based on these genes had good accuracy in predicting the OS of patients and might be promising candidates of therapeutic targets. This study may provide a new insight for the diagnosis, treatment and prognosis of BLCA from the perspective of immunology. However, further experimental studies are necessary to reveal the underlying mechanisms by which these genes mediate the progression of BLCA.

Identification of a Prognosis-Related Risk Signature for Bladder Cancer to Predict Survival and Immune Landscapes

Background

Bladder cancer is the tenth most common cancer worldwide. Valuable biomarkers in the field of diagnostic bladder cancer are urgently required.

Method

Here, the gene expression matrix and clinical data were obtained from The Cancer Genome Atlas (TCGA), GSE13507, GSE32894, and Mariathasan et al. Five prognostic genes were identified by the univariate, robust, and multivariate Cox's regression and were used to develop a prognosis-related model. The Kaplan-Meier survival curves and receiver operating characteristics were used to evaluate the model's effectiveness. The potential biological functions of the selected genes were analyzed using CIBERSORT and ESTIMATE algorithms. Cancer Therapeutics Response Portal (CTRP) and PRISM datasets were used to identify drugs with high sensitivity. Subsequently, using the bladder cancer (BLCA) cell lines, the role of TNFRSF14 was determined by Western blotting, cell proliferation assay, and 5-ethynyl-20-deoxyuridine assay.

Results

GSDMB, CLEC2D, APOL2, TNFRSF14, and GBP2 were selected as prognostic genes in bladder cancer patients. The model's irreplaceable reliability was validated by the training and validation cohorts. CD8+ T cells were highly infiltrated in the high-TNFRSF14-expression group, and M2 macrophages were the opposite. Higher expression of TNFRSF14 was associated with higher expression levels of LCK, interferon, MHC-I, and MHC-II, while risk score was the opposite. Many compounds with higher sensitivity for treating bladder cancer patients in the low-TNFRSF14-expression group were identified, with obatoclax being a potential drug most likely to treat patients in the low-TNFRSF14-expression group. Finally, the proliferation of BLCA cell lines was increased in the TNFRSF14-reduced group, and the differential expression was identified. TNFRSF14 plays a role in bladder cancer progression through the Wnt/β-catenin-dependent pathway. TNFRSF14 is a potential protective biomarker involved in cell proliferation in BLCA.

Conclusion

We conducted a study to establish a 5-gene score model, providing reliable prediction for the outcome of bladder cancer patients and therapeutic drugs to individualize therapy. Our findings provide a signature that might help determine the optimal treatment for individual patients with bladder cancer.

Development of a seven-gene tumor immune microenvironment prognostic signature for high-risk grade III endometrial cancer

Uterine corpus endometrial carcinoma locally infiltrates numerous immune cells and other tumor immune microenvironment components. These cells are involved in malignant tumor growth and proliferation and the process of resistance toward immunotherapies. Here, we aimed to develop a tumor immune microenvironment-related prognostic signature for high-risk grade III endometrial carcinoma based on The Cancer Genome Atlas. The signature was systematically correlated with immune infiltration characteristics of the tumor microenvironment. The seven-gene Riskscore signature was robust and performed well in training, testing, and Gene Expression Omnibus-independent cohorts. A nomogram comprising the gene signature accurately predicted patient prognosis, with our model performing better than other endometrial cancer-related signatures. Analysis of the IMvigor210 immunotherapy cohort revealed that subgroups with a low Riskscore had a better prognosis than subgroups with a high Riskscore. Subgroups with a low Riskscore exhibited immune cell infiltration and inflammatory profiles, whereas subgroups with a high Riskscore experienced progressive disease. The receiver operating characteristic curve indicated that risk score, neoantigen, and tumor mutation burden models together accurately predicted treatment response. Taken together, we developed a tumor microenvironment-based seven-gene prognostic stratification system to predict the prognosis of patients with high-risk endometrial cancer and guide more effective immunotherapy strategies.

Investigation of IL-4, IL-10, and HVEM polymorphisms with esophageal squamous cell carcinoma: a case-control study involving 1929 participants

It is believed that an individual’s hereditary factors may be involved in the development of esophageal cancer (EC). The present study recruited 721 esophageal squamous cell carcinoma (ESCC) cases and 1208 controls and explored the roles of single nucleotide polymorphisms (SNPs) in the interleukin-4 (IL-4), IL-10, and herpesvirus entry mediator (HVEM) genes in contributing to ESCC risk. IL-4, IL-10, and HVEM SNPs were analyzed by employing an SNPscan method. After adjustment for body mass index (BMI), smoking, drinking, age and gender, we identified that the rs2070874 T>C locus in IL-4 gene decreased the risk of ESCC (CC vs. TT: P=0.008; CC vs. TT/TC: P=0.010). After a stratified analysis, we suggested that the IL-4 rs2070874 T>C variants might be a protective factor for ESCC in male, ≥63 years old, never smoking, drinking and BMI < 24 kg/m² subgroups. In addition, we identified that the rs2243263 G>C polymorphism in IL-4 gene was a risk factor for ESCC development in the BMI ≥ 24 kg/m² subgroup (GC vs. GG: P=0.030 and GC/CC vs. GG: P=0.018). We identified an association of the IL-4 rs2070874 T>C SNP with the decreased susceptibility of ESCC in stage I/II subgroup. Finally, we found an association of the IL-10 rs1800872 T>G SNP with a worse differentiation (TG vs. TT: P=0.048 and GG/TG vs. TT: P=0.032). In conclusion, the findings indicate a potential importance of IL-4 rs2070874 T>C, IL-4 rs2243263 G>C and IL-10 rs1800872 T>G SNPs in the development of ESCC.

The Analysis of PTPN6 for Bladder Cancer: An Exploratory Study Based on TCGA

PTPN6 (protein tyrosine phosphatase nonreceptor type 6), a tyrosine phosphatase, is known to be signaling molecules that regulate a variety of cellular processes including cell growth, differentiation, mitotic cycle, and oncogenic transformation. Previous studies have demonstrated that PTPN6 expression is relatively elevated in several malignancies. However, the role of PTPN6 in bladder cancer (BC) remains unclear. The purpose of this study was to explore the prognostic value of PTPN6 in BC. RNA-seq data from The Cancer Genome Atlas (TCGA) was used to identify the expression level of PTPN6 in BC. The relationship between clinical pathologic features and PTPN6 were analyzed with the Wilcoxon signed-rank test. The prognostic and predictive value of PTPN6 was evaluated by survival analysis and nomogram. Gene Set Enrichment Analysis (GSEA) was conducted to explore the potential molecular mechanisms of PTPN6 in BC. Finally, Tumor Immune Estimation Resource (TIMER) was applied to investigate the relationship between PTPN6 and immune cell infiltration in the tumor microenvironment. Results indicated that PTPN6 was overexpressed in BC tissues compared with normal bladder tissues and was significantly correlated with grade, stage, T, and N. Survival analysis showed that low expression of PTPN6 was significantly related to the poor overall survival (OS) in BC patients. Coexpression analysis showed that PTPN6 and TNFRSF14 (Tumor necrosis factor receptor superfamily member 14) have a close correlation in BC. GSEA showed that multiple cancer-associated signaling pathways are differentially enriched in the PTPN6 high expression phenotype. Moreover, the expression level of PTPN6 was positively associated with the infiltration of B cells, CD4+T cells, dendritic cells, and neutrophils and negatively associated with CD8+ T cells and macrophages in BC. In conclusion, we identified that PTPN6 may be a novel prognostic biomarker in BC based on the TCGA database. Further clinical trials are needed to confirm our observations and mechanisms underlying the prognostic value of PTPN6 in BC also deserve further experimental exploration.

HVEM has a broader expression than PD-L1 and constitutes a negative prognostic marker and potential treatment target for melanoma

HVEM (Herpes Virus Entry Mediator) engagement of BTLA (B and T Lymphocyte Attenuator) triggers inhibitory signals in T cells and could play a role in evading antitumor immunity. Here, HVEM expression levels in melanoma metastases were analyzed by immunohistochemistry, correlated with overall survival (OS) in 116 patients, and validated by TCGA transcriptomic data. Coincident expression of HVEM and its ligand BTLA was studied in tumor cells and tumor-infiltrating lymphocytes (TILs) by flow cytometry (n = 21) and immunofluorescence (n = 5). Candidate genes controlling HVEM expression in melanoma were defined by bioinformatics studies and validated by siRNA gene silencing. We found that in patients with AJCC stage III and IV melanoma, OS was poorer in those with high HVEM expression on melanoma cells, than in those with a low expression, by immunohistochemistry (p = .0160) or TCGA transcriptomics (p = .0282). We showed a coincident expression of HVEM at the surface of melanoma cells and of BTLA on TILs. HVEM was more widely expressed than PD-L1 in melanoma cells. From a mechanistic perspective, in contrast to PDL1, HVEM expression did not correlate with an IFNγ signature but with an aggressive gene signature. Interestingly, this signature contained MITF, a key player in melanoma biology, whose expression correlated strongly with HVEM. Finally, siRNA gene silencing validated MITF control of HVEM expression. In conclusion, HVEM expression seems to be a prognosis marker and targeting this axis by checkpoint-inhibitors may be of interest in metastatic melanoma.