Expression and survival analysis of the STAT gene family in diffuse gliomas using integrated bioinformatics

STAT4 facilitates PD-L1 level via IL-12R/JAK2/STAT3 axis and predicts immunotherapy response in breast cancer

Signal transducer and activator of transcription 4 (STAT4) is a critical transcription factor for T helper cell differentiation and tumor cells. Although its prognostic role and gene function have been reported in several carcinomas, the role of STAT4 in vitro and in vivo in breast cancer remains poorly understood. The effect of STAT4 in immunotherapy is also unclear. Therefore, we integrated bulk transcriptomics, experiments, and single-cell transcriptomics to systematically analyze its function in prognosis and signaling pathway. Several clinical breast cancer cohorts confirmed STAT4 as a T-cell relevant prognostic biomarker. Overexpressed STAT4 increased programmed cell death ligand 1 (PD-L1) and major histocompatibility complex class II levels in breast cancer cells. In molecular mechanism, transcriptional synergy between STAT4 and STAT3 transactivated interleukin (IL)-12R and involved a positive feedback loop: STAT4/IL-12R/JAK2-STAT3-STAT4, which contributed to the upregulation of PD-L1 expression. The above signaling axis was defined as the STAT4-related pathway and its score was used to predict T-cell expansion and anti-PD1 treatment response. These findings highlight a novel molecular mechanism indirectly regulating PD-L1 through the STAT4-related pathway: IL-12R/JAK2-STAT3-STAT4/PD-L1, and it has potential application in predicting anti-PD-1 immunotherapy response, which may pave the way for stratified immunotherapy in breast cancer.

A novel defined risk signature of interferon response genes predicts the prognosis and correlates with immune infiltration in glioblastoma

BACKGROUND

Interferons (IFNs) have been implemented as anti-tumor immunity agents in clinical trials of glioma, but only a subset of glioblastoma (GBM) patients profits from it. The predictive role of IFNs stimulated genes in GBM needs further exploration to investigate the clinical role of IFNs.

METHODS

This study screened 526 GBM patients from three independent cohorts. The transcriptome data with matching clinical information were analyzed using R. Immunohistochemical staining data from the Human Protein Atlas and DNA methylation data from MethSurv were used for validation in protein and methylation level respectively.

RESULTS

We checked the survival effect of all 491 IFNs response genes, and found 54 genes characterized with significant hazard ratio in overall survival (OS). By protein-protein interaction analysis, 10 hub genes were selected out for subsequent study. And based on the expression of these 10 genes, GBM patients could be divided into two subgroups with significant difference in OS. Furthermore, the least absolute shrinkage and selection operator cox regression model was utilized to construct a multigene risk signature, including STAT3, STAT2 and SOCS3, which could serve as an independent prognostic predictor for GBM. The risk model was validated in two independent GBM cohorts. The GBM patients with high risk scores mainly concentrated in the GBM Mesenchymal subtype. The higher risk group was enriched in hypoxia, angiogenesis, EMT, glycolysis and immune pathways, and had increased Macrophage M2 infiltration and high expression of immune checkpoint CD274 (namely PD-L1).

CONCLUSIONS

Our findings revealed the three-gene risk model could be an independent prognostic predictor for GBM, and they were crucial participants in immunosuppressive microenvironment of GBM.

The Emerging Role of MTHFD Family Genes in Regulating the Tumor Immunity of Oral Squamous Cell Carcinoma

Objective

To investigate the function and regulatory mechanisms of methylenetetrahydrofolate dehydrogenase (MTHFD) family genes in oral squamous cell carcinoma (OSCC), especially focus on their regulating role in tumor immunity.

Methods

The publicly available data from the TCGA database were used to investigate the expression pattern and regulatory role of MTHFD family genes in OSCC. More importantly, the involvement of MTHFD family genes in tumor immunity was investigated in terms of immune and stromal cell infiltration in tumor microenvironment, tumor-infiltrating immune cells, and immunomodulatory genes (e.g., immunoinhibitory genes and immunostimulatory genes). Statistical analysis was performed using R software packages and public web servers.

Results

MTHFD family genes were considerably upregulated in OSCC as compared with normal oral tissue. Patients with high MTHFD2 expression presented worse survival outcomes than those with low MTHFD2 expression. Functional enrichment analysis showed that the top 100 positively and negatively correlated genes of the MTHFD family genes were significantly enriched in several KEGG pathways, including cell cycle, spliceosome, DNA replication, and Th17 cell differentiation. As a result of tumor immunity analysis, MTHFD2L expression was found to be negatively related to the Estimate-Stromal-Immune score in OSCC; however, there was no statistical significance between the Estimate-Stromal-Immune score and MTHFD1, MTHFD1L, or MTHFD2 in OSCC. Additionally, MTHFD family genes were found to be significantly positively correlated with tumor-infiltrating immune cells, including Treg and Th17 cells. Moreover, MTHFD family genes were significantly correlated with several immune inhibitory genes such as CD274 and CTLA4 and several immune-stimulatory genes such as CXCL12, CXCR4, and TMIGD2.

Conclusion

Given the expression pattern, prognostic value, biological functions, and involvement in tumor immunity, MTHFD family genes could serve as potential therapeutic biomarkers in targeting tumor immunity in oral cancer.

Implications of Human Antimicrobial Peptide Defensin Beta-1 in Clinical Oral Squamous Cell Carcinoma Patients via an Integrated Bioinformatics Approach

BACKGROUND

The human antimicrobial peptide defensin beta 1 (DEFB1) has been found to play antimicrobial and anti-inflammatory roles in oral diseases; however, its tumor-regulating role in oral squamous cell carcinoma (OSCC) has not yet been researched by using an integrative bioinformatics approach.

OBJECTIVE

To investigate the regulating mechanisms of the DEFB1 gene in OSCC in terms of its expression patterns, prognostic values, biological functions, and implication for tumor immunity.

METHODS

The DEFB1 gene expression pattern and regulatory involvement in OSCC were investigated using publically accessible data from TCGA database. R software tools and public web servers were utilized to conduct statistical analysis of data from cancer and noncancerous samples.

RESULTS

DEFB1 was found to be significantly downregulated in OSCC tumor samples compared with healthy control oral samples. The DEFB1 gene was found associated with the prognostic outcomes of OSCC, and its upregulation represented better survival outcome. Gene set enrichment analysis (GSEA) results showed that DEFB1-significantly correlated genes were mainly enriched in four signaling pathways mediating the antitumor role of DEFB1 in OSCC, including extracellular matrix-related pathway, RTK/PI3K/AKT/mTOR pathway, keratinization, and cytokine-related pathway. The gene-gene interaction network showed that DEFB1 was closely correlated with several genes, for example, CCR6 (C-C motif chemokine receptor 6), CXCL1 (C-X-C motif chemokine ligand 1), MAP4K2 (mitogen-activated protein kinase kinase kinase kinase 2), PTGER3 (prostaglandin E receptor 3), and MMP7 (matrix metallopeptidase 7). Moreover, DEFB1 was found to be involved in the tumor immunity of OSCC by regulating the function of tumor macrophage cells, mast cells, T cells, and NK cells.

CONCLUSIONS

Given the dysregulation, prognostic value, and tumor progression-related biological pathway alteration, indicating the tumor immune-modulatory role of DEFB1 in OSCC, the DEFB1 gene should be regarded as a potential therapeutic target for treating oral cancer.