Tyrosine hydroxylase inhibits HCC progression by downregulating TGFß/Smad signaling.

The alteration of metabolic processes has been found to have significant impacts on the development of hepatocellular carcinoma (HCC). Nevertheless, the effects of dysfunction of tyrosine metabolism on the development of HCC remains to be discovered. This research demonstrated that tyrosine hydroxylase (TH), which responsible for the initial and limiting step in the bio-generation of the neuro-transmitters dopamine and adrenaline, et al. was shown to be reduced in HCC. Increased expression of TH was found facilitates the survival of HCC patients. In addition, decreased TH indicated larger tumor size, much more numbers of tumor, higher level of AFP, and the presence of cirrhosis. TH effectively impairs the growth and metastasis of HCC cells, a process dependent on the phosphorylation of serine residues (S19/S40). TH directly binds to Smad2 and hinders the cascade activation of TGFß/Smad signaling with the treatment of TGFß1. In summary, our study uncovered the non-metabolic functions of TH in the development of HCC and proposes that TH might be a promising biomarker for diagnosis as well as an innovative target for metastatic HCC.

Suppression of ITPKB degradation by Trim25 confers TMZ resistance in glioblastoma through ROS homeostasis.

Temozolomide (TMZ) represents a standard-of-care chemotherapeutic agent in glioblastoma (GBM). However, the development of drug resistance constitutes a significant hurdle in the treatment of malignant glioma. Although specific innovative approaches, such as immunotherapy, have shown favorable clinical outcomes, the inherent invasiveness of most gliomas continues to make them challenging to treat. Consequently, there is an urgent need to identify effective therapeutic targets for gliomas to overcome chemoresistance and facilitate drug development. This investigation used mass spectrometry to examine the proteomic profiles of six pairs of GBM patients who underwent standard-of-care treatment and surgery for both primary and recurrent tumors. A total of 648 proteins exhibiting significant differential expression were identified. Gene Set Enrichment Analysis (GSEA) unveiled notable alterations in pathways related to METABOLISM_OF_LIPIDS and BIOLOGICAL_OXIDATIONS between the primary and recurrent groups. Validation through glioma tissue arrays and the Xiangya cohort confirmed substantial upregulation of inositol 1,4,5-triphosphate (IP3) kinase B (ITPKB) in the recurrence group, correlating with poor survival in glioma patients. In TMZ-resistant cells, the depletion of ITPKB led to an increase in reactive oxygen species (ROS) related to NADPH oxidase (NOX) activity and restored cell sensitivity to TMZ. Mechanistically, the decreased phosphorylation of the E3 ligase Trim25 at the S100 position in recurrent GBM samples accounted for the weakened ITPKB ubiquitination. This, in turn, elevated ITPKB stability and impaired ROS production. Furthermore, ITPKB depletion or the ITPKB inhibitor GNF362 effectively overcome TMZ chemoresistance in a glioma xenograft mouse model. These findings reveal a novel mechanism underlying TMZ resistance and propose ITPKB as a promising therapeutic target for TMZ-resistant GBM.

Emerging mechanisms of the unfolded protein response in therapeutic resistance: from chemotherapy to Immunotherapy.

The accumulation of unfolded or misfolded proteins in the endoplasmic reticulum (ER) causes ER stress and activates the unfolded protein response (UPR). As an adaptive cellular response to hostile microenvironments, such as hypoxia, nutrient deprivation, oxidative stress, and chemotherapeutic drugs, the UPR is activated in diverse cancer types and functions as a dynamic tumour promoter in cancer development; this role of the UPR indicates that regulation of the UPR can be utilized as a target for tumour treatment. T-cell exhaustion mainly refers to effector T cells losing their effector functions and expressing inhibitory receptors, leading to tumour immune evasion and the loss of tumour control. Emerging evidence suggests that the UPR plays a crucial role in T-cell exhaustion, immune evasion, and resistance to immunotherapy. In this review, we summarize the molecular basis of UPR activation, the effect of the UPR on immune evasion, the emerging mechanisms of the UPR in chemotherapy and immunotherapy resistance, and agents that target the UPR for tumour therapeutics. An understanding of the role of the UPR in immune evasion and therapeutic resistance will be helpful to identify new therapeutic modalities for cancer treatment. Video Abstract.

E3 ligase Trim35 inhibits LSD1 demethylase activity through K63-linked ubiquitination and enhances anti-tumor immunity in NSCLC.

Targeting lysine-specific histone demethylase 1A (LSD1) can improve tumor immunogenicity of poorly immunogenic tumors, such as non-small cell lung cancer (NSCLC), with elevated T cell infiltration and sensitize tumors to anti-PD-1 therapy. However, the lack of reliable biomarkers limits utilization of LSD1 inhibitors in cancer therapy. Here, we identify an E3 ligase, Trim35, as an effective biomarker for high activity of LSD1 to predict prognosis of LSD1-targeted therapy as well as immunotherapy. Mechanistically, Trim35 represses LSD1 demethylase activity by mediating K63 ubiquitination at lysine site 422 of LSD1. Suppressed LSD1 activity facilitates ERGIC1 transcription, followed by autophagy inhibition and IFNGR1 stabilization to activate IFN-γ signaling, leading to increased MHC class I expression and immune surveillance of NSCLC cells. Furthermore, combinational use of an LSD1 inhibitor and anti-PD-1 therapy can significantly eradicate poorly immunogenic lung cancer with low Trim35. These findings strongly suggest that Trim35 is a promising biomarker for prediction of immunotherapy outcome in NSCLC.

Whole genome sequencing identifies genetic variants associated with neurogenic inflammation in rosacea.

Rosacea is a chronic inflammatory skin disorder with high incidence rate. Although genetic predisposition to rosacea is suggested by existing evidence, the genetic basis remains largely unknown. Here we present the integrated results of whole genome sequencing (WGS) in 3 large rosacea families and whole exome sequencing (WES) in 49 additional validation families. We identify single rare deleterious variants of LRRC4, SH3PXD2A and SLC26A8 in large families, respectively. The relevance of SH3PXD2A, SLC26A8 and LRR family genes in rosacea predisposition is underscored by presence of additional variants in independent families. Gene ontology analysis suggests that these genes encode proteins taking part in neural synaptic processes and cell adhesion. In vitro functional analysis shows that mutations in LRRC4, SH3PXD2A and SLC26A8 induce the production of vasoactive neuropeptides in human neural cells. In a mouse model recapitulating a recurrent Lrrc4 mutation from human patients, we find rosacea-like skin inflammation, underpinned by excessive vasoactive intestinal peptide (VIP) release by peripheral neurons. These findings strongly support familial inheritance and neurogenic inflammation in rosacea development and provide mechanistic insight into the etiopathogenesis of the condition.

Phosphorylation of PPDPF via IL6-JAK2 activates the Wnt/ß-catenin pathway in colorectal cancer.

Inflammation plays an important role in the initiation and progression of colorectal cancer (CRC) and leads to ß-catenin accumulation in colitis-related CRC. However, the mechanism remains largely unknown. Here, pancreatic progenitor cell differentiation and proliferation factor (PPDPF) is found to be upregulated in CRC and significantly correlated with tumor-node-metastasis (TNM) stages and survival time. Knockout of PPDPF in the intestinal epithelium shortens crypts, decreases the number of stem cells, and inhibits the growth of organoids and the occurrence of azoxymethane (AOM)/dextran sodium sulfate (DSS)-induced CRC. Mechanistically, PPDPF is found to interact with Casein kinase 1α (CK1α), thereby disrupting its binding to Axin, disassociating the ß-catenin destruction complex, decreasing the phosphorylation of ß-catenin, and activating the Wnt/ß-catenin pathway. Furthermore, interleukin 6 (IL6)/Janus kinase 2 (JAK2)-mediated inflammatory signals lead to phosphorylation of PPDPF at Tyr16 and Tyr17, stabilizing the protein. In summary, this study demonstrates that PPDPF is a key molecule in CRC carcinogenesis and progression that connects inflammatory signals to the Wnt/ß-catenin signaling pathway, providing a potential novel therapeutic target.

The expression profiles of signature genes from CD103+LAG3+ tumour-infiltrating lymphocyte subsets predict breast cancer survival.

BACKGROUND: Tumour-infiltrating lymphocytes (TILs), including T and B cells, have been demonstrated to be associated with tumour progression. However, the different subpopulations of TILs and their roles in breast cancer remain poorly understood. Large-scale analysis using multiomics data could uncover potential mechanisms and provide promising biomarkers for predicting immunotherapy response. METHODS: Single-cell transcriptome data for breast cancer samples were analysed to identify unique TIL subsets. Based on the expression profiles of marker genes in these subsets, a TIL-related prognostic model was developed by univariate and multivariate Cox analyses and LASSO regression for the TCGA training cohort containing 1089 breast cancer patients. Multiplex immunohistochemistry was used to confirm the presence of TIL subsets in breast cancer samples. The model was validated with a large-scale transcriptomic dataset for 3619 breast cancer patients, including the METABRIC cohort, six chemotherapy transcriptomic cohorts, and two immunotherapy transcriptomic cohorts. RESULTS: We identified two TIL subsets with high expression of CD103 and LAG3 (CD103+LAG3+), including a CD8+ T-cell subset and a B-cell subset. Based on the expression profiles of marker genes in these two subpopulations, we further developed a CD103+LAG3+ TIL-related prognostic model (CLTRP) based on CXCL13 and BIRC3 genes for predicting the prognosis of breast cancer patients. CLTRP-low patients had a better prognosis than CLTRP-high patients. The comprehensive results showed that a low CLTRP score was associated with a high TP53 mutation rate, high infiltration of CD8 T cells, helper T cells, and CD4 T cells, high sensitivity to chemotherapeutic drugs, and a good response to immunotherapy. In contrast, a high CLTRP score was correlated with a low TP53 mutation rate, high infiltration of M0 and M2 macrophages, low sensitivity to chemotherapeutic drugs, and a poor response to immunotherapy. CONCLUSIONS: Our present study showed that the CLTRP score is a promising biomarker for distinguishing prognosis, drug sensitivity, molecular and immune characteristics, and immunotherapy outcomes in breast cancer patients. The CLTRP could serve as a valuable tool for clinical decision making regarding immunotherapy.

PHF5A regulates the expression of the DOCK5 variant to promote HNSCC progression through p38 MAPK activation.

BACKGROUND: Previously, we identified an oncogenic splicing variant of DOCK5 in head and neck squamous cell carcinoma (HNSCC); however, the mechanism for the generation of this specific DOCK5 variant remains unknown. This study aims to explore the potential spliceosome genes involved in the production of the DOCK5 variant and validate its role in regulating the progression of HNSCC. METHODS: The differentially expressed spliceosome genes involved in the DOCK5 variant were analysed in The Cancer Genome Atlas (TCGA), and the correlation between the DOCK5 variant and the potential spliceosome gene PHF5A was verified by qRT-PCR. The expression of PHF5A was detected in HNSCC cells, TCGA data and a separate primary tumour cohort. The functional role of PHF5A was examined using CCK-8, colony formation, cell scratch and Transwell invasion assays in vitro and validated in vivo in xenograft models of HNSCC. Western blot analysis was used to explore the potential mechanism of PHF5A in HNSCC. RESULTS: PHF5A was one of the top upregulated spliceosome genes in TCGA HNSCC samples with highly expressed DOCK5 variants. Knockdown or overexpression of PHF5A in HNSCC cells correspondingly altered the level of the DOCK5 variant. PHF5A was highly expressed in tumour cells and tissues and correlated with a worse prognosis of HNSCC. Loss- and gain-of-function experiments demonstrated that PHF5A could promote the proliferation, migration and invasion of HNSCC cells in vitro and in vivo. Moreover, PHF5A inhibition reversed the oncogenic effect of the DOCK5 variant in HNSCC. Western blot analysis showed that PHF5A activated the p38 MAPK pathway, and inhibition of p38 MAPK further reversed the effect of PHF5A on the proliferation, migration and invasion of HNSCC cells. CONCLUSION: PHF5A regulates the alternative splicing of DOCK5 to promote HNSCC progression through p38 MAPK activation, which provides potential therapeutic implications for HNSCC patients.

Nicotinamide adenine dinucleotide kinase promotes lymph node metastasis of NSCLC via activating ID1 expression through BMP pathway.

Metastasis is a significant cause of high mortality in lung cancer. Lymph node (LN) metastasis is the most common metastatic pathway in non-small cell lung cancer and the most crucial factor affecting the prognosis of NSCLC. Nevertheless, the molecular mechanism underlying metastasis is unknown. We demonstrated that higher NADK expression suggests worsened survival prognosis, and NADK expression positively correlates with the lymph node metastasis rate and TNM and AJCC stages in NSCLC patients. Moreover, patients with LN metastasis show higher NADK expression than those without LN metastasis. NADK can promote NSCLC progression by enhancing the migration, invasion, lymph node metastasis and growth of NSCLC cells. Mechanistically, NADK inhibits the ubiquitination and degradation of BMPR1A by interacting with Smurf1, further activating the BMPs signalling pathway and promoting ID1 transcription. In conclusion, NADK may be a potential diagnostic indicator and a novel therapeutic target for metastatic NSCLC.

SLC3A2 promotes tumor-associated macrophage polarization through metabolic reprogramming in lung cancer.

Tumor-associated macrophages (TAMs) are one of the most abundant immunosuppressive cells in the tumor microenvironment and possess crucial functions in facilitating tumor progression. Emerging evidence indicates that altered metabolic properties in cancer cells support the tumorigenic functions of TAMs. However, the mechanisms and mediators the underly the cross-talk between cancer cells and TAMs remain largely unknown. In the present study, we revealed that high solute carrier family 3 member 2 (SLC3A2) expression in lung cancer patients was associated with TAMs and poor prognosis. Knockdown of SLC3A2 in lung adenocarcinoma cells impaired M2 polarization of macrophages in a coculture system. Using metabolome analysis, we identified that SLC3A2 knockdown altered the metabolism of lung cancer cells and changed multiple metabolites, including arachidonic acid, in the tumor microenvironment. More importantly, we showed that arachidonic acid was responsible for SLC3A2-mediated macrophage polarization in the tumor microenvironment to differentiate into M2 type both in vitro and in vivo. Our data illustrate previously undescribed mechanisms responsible for TAM polarization and suggest that SLC3A2 acts as a metabolic switch on lung adenocarcinoma cells to induce macrophage phenotypic reprogramming through arachidonic acid.

The Association of R-Loop Binding Proteins Subtypes with CIN Implicates Therapeutic Strategies in Colorectal Cancer.

Chromosomal instability (CIN) covers approximately 65 to 70% of colorectal cancer patients and plays an essential role in cancer progression. However, the molecular features and therapeutic strategies related to those patients are still controversial. R-loop binding proteins (RLBPs) exert significant roles in transcription and replication. Here, integrative colorectal cancer proteogenomic analysis identified two RLBPs subtypes correlated with distinct prognoses. Cluster I (CI), represented by high expression of RLBPs, was associated with the CIN phenotype. While Cluster II (CII) with the worst prognosis and low expression of RLBPs was composed of a high percentage of patients with mucinous adenocarcinoma or right-sided colon cancer. The molecular feature analysis revealed that the active RNA processing, ribosome synthesis, and aberrant DNA damage repair were shown in CI, a high inflammatory signaling pathway, and lymphocyte infiltration was enriched in CII. In addition, we revealed 42 tumor-associated RLBPs proteins. The CI with high expression of tumor-associated proteins was sensitive to drugs targeting genome integrity and EGFR in both cell and organoid models. Thus, our study unveils a significant molecular association of the CIN phenotype with RLBPs, and also provides a powerful resource for further functional exploration of RLBPs in cancer progression and therapeutic application.

Girdin Promotes Tumorigenesis and Chemoresistance in Lung Adenocarcinoma by Interacting with PKM2.

Girdin, an Akt substrate, has been reported to promote tumorigenesis in various tumors. However, the role of Girdin in a spontaneous tumor model has not yet been explored. Here, we studied the role of Girdin in lung adenocarcinoma (LUAD) using the autochthonous mouse model and found that Girdin led to LUAD progression and chemoresistance by enhancing the Warburg effect. Mechanistically, Girdin interacted with pyruvate kinase M2 (PKM2), which played a vital role in aerobic glycolysis. Furthermore, Girdin impaired Platelet Derived Growth Factor Receptor Beta (PDGFRß) degradation, which in turn, promoted PKM2 tyrosine residue 105 (Y105) phosphorylation and inhibited PKM2 activity, subsequently promoting aerobic glycolysis in cancer cells. Taken together, our study demonstrates that Girdin is a crucial regulator of tumor growth and may be a potential therapeutic target for overcoming the resistance of LUAD cells to chemotherapy.

Drebrin promotes lung adenocarcinoma cell migration through inducing integrin ß1 endocytosis.

Lung adenocarcinoma (LUAD) is the most common type of lung cancers, which remains the leading cause of cancer-related death worldwide. Drebrin can promote cell migration and invasion with poor prognosis, but its roes in LUAD tumor progression remains unknown. We showed that the expression of Drebrin was upregulated in clinical LUAD samples. A Kaplan-Meier survival analysis showed that a high expression of Drebrin predicated poor prognosis in LUAD. In vitro, Drebrin promoted anchorage-independent growth and migration of LUAD cells. Drebrin interacted with dynamin through CT domain, and served as an adaptor to promote LUAD cell migration through inducing integrin ß1 endocytosis. Thus, this study demonstrated the critical role of Drebrin in LUAD and associated mechanism.

Asparagine synthetase regulates lung-cancer metastasis by stabilizing the ß-catenin complex and modulating mitochondrial response.

The availability of asparagine is the limitation of cell growth and metastasis. Asparagine synthetase (ASNS) was an essential enzyme for endogenous asparagine products. In our study, ASNS-induced asparagine products were essential to maintain tumor growth and colony formations in vitro. But mutated ASNS which defected endogenous asparagine products still upregulated cell invasiveness, which indicated that ASNS promoted invasiveness by alternative pathways. Mechanically, ASNS modulated Wnt signal transduction by promoting GSK3ß phosphorylation on ser9 and stabilizing the ß-catenin complex, as result, ASNS could promote more ß-catenin translocation into nucleus independent of endogenous asparagine. At the same time, ASNS modulated mitochondrial response to Wnt stimuli with increased mitochondrial potential and membrane fusion. In summary, ASNS promoted metastasis depending on Wnt pathway and mitochondrial functions even without endogenous asparagine products.

Engineering an enthesis-like graft for rotator cuff repair: An approach to fabricate highly biomimetic scaffold capable of zone-specifically releasing stem cell differentiation inducers.

Rotator cuff (RC) attaches to humerus across a triphasic yet continuous tissue zones (bone-fibrocartilage-tendon), termed "enthesis". Regrettably, rapid and functional enthesis regeneration is challenging after RC tear. The existing grafts bioengineered for RC repair are insufficient, as they were engineered by a scaffold that did not mimic normal enthesis in morphology, composition, and tensile property, meanwhile cannot simultaneously stimulate the formation of bone-fibrocartilage-tendon tissues. Herein, an optimized decellularization approach based on a vacuum aspiration device (VAD) was developed to fabricate a book-shaped decellularized enthesis matrix (O-BDEM). Then, three recombinant growth factors (CBP-GFs) capable of binding collagen were synthesized by fusing a collagen-binding peptide (CBP) into the N-terminal of BMP-2, TGF-ß3, or GDF-7, and zone-specifically tethered to the collagen of O-BDEM to fabricate a novel scaffold (CBP-GFs/O-BDEM) satisfying the above-mentioned requirements. After ensuring the low immunogenicity of CBP-GFs/O-BDEM by a novel single-cell mass cytometry in a mouse model, we interleaved urine-derived stem cell-sheets into this CBP-GFs/O-BDEM to bioengineer an enthesis-like graft. Its high-performance on regenerating enthesis was determined in a canine model. These findings indicate this CBP-GFs/O-BDEM may be an excellent scaffold for constructing enthesis-like graft to patch large/massive RC tears, and provide breakthroughs in fabricating graded interfacial tissue.

Fatty acid oxidation fuels glioblastoma radioresistance with CD47-mediated immune evasion.

Glioblastoma multiforme (GBM) remains the top challenge to radiotherapy with only 25% one-year survival after diagnosis. Here, we reveal that co-enhancement of mitochondrial fatty acid oxidation (FAO) enzymes (CPT1A, CPT2 and ACAD9) and immune checkpoint CD47 is dominant in recurrent GBM patients with poor prognosis. A glycolysis-to-FAO metabolic rewiring is associated with CD47 anti-phagocytosis in radioresistant GBM cells and regrown GBM after radiation in syngeneic mice. Inhibition of FAO by CPT1 inhibitor etomoxir or CRISPR-generated CPT1A-/-, CPT2-/-, ACAD9-/- cells demonstrate that FAO-derived acetyl-CoA upregulates CD47 transcription via NF-κB/RelA acetylation. Blocking FAO impairs tumor growth and reduces CD47 anti-phagocytosis. Etomoxir combined with anti-CD47 antibody synergizes radiation control of regrown tumors with boosted macrophage phagocytosis. These results demonstrate that enhanced fat acid metabolism promotes aggressive growth of GBM with CD47-mediated immune evasion. The FAO-CD47 axis may be targeted to improve GBM control by eliminating the radioresistant phagocytosis-proofing tumor cells in GBM radioimmunotherapy.

Pantothenate Kinase 1 Inhibits the Progression of Hepatocellular Carcinoma by Negatively Regulating Wnt/ß-catenin Signaling.

Hyperactivation of Wnt/ß-catenin signaling has been reported in hepatocellular carcinoma (HCC). However, the mechanisms underlying the hyperactivation of Wnt/ß-catenin signaling are incompletely understood. In this study, Pantothenate kinase 1 (PANK1) is shown to be a negative regulator of Wnt/ß-catenin signaling. Downregulation of PANK1 in HCC correlates with clinical features. Knockdown of PANK1 promotes the proliferation, growth and invasion of HCC cells, while overexpression of PANK1 inhibits the proliferation, growth, invasion and tumorigenicity of HCC cells. Mechanistically, PANK1 binds to CK1α, exerts protein kinase activity and cooperates with CK1α to phosphorylate N-terminal serine and threonine residues in ß-catenin both in vitro and in vivo. Additionally, the expression levels of PANK1 and ß-catenin can be used to predict the prognosis of HCC. Collectively, the results of this study highlight the crucial roles of PANK1 protein kinase activity in inhibiting Wnt/ß-catenin signaling, suggesting that PANK1 is a potential therapeutic target for HCC.

SIK2 maintains breast cancer stemness by phosphorylating LRP6 and activating Wnt/ß-catenin signaling.

Breast cancer stem cells (BCSCs) are the main drivers of recurrence and metastasis. However, commonly used drugs rarely target BCSCs. Via screenings, we found that Salt-inducible kinase 2 (SIK2) participated in breast cancer (BC) stemness maintenance and zebrafish embryos development. SIK2 was upregulated in recurrence samples. Knockdown of SIK2 expression reduced the proportion of BCSCs and the tumor initiation of BC cells. Mechanistically, SIK2, phosphorylated by CK1α, directly phosphorylated LRP6 in a SIK2 kinase activity-dependent manner, leading to Wnt/ß-catenin signaling pathway activation. ARN-3236 and HG-9-91-01, inhibitors of SIK2, inhibited LRP6 phosphorylation and ß-catenin accumulation and disturbed stemness maintenance. In addition, the SIK2-activated Wnt/ß-catenin signaling led to induction of IDH1 expression, causing metabolic reprogramming in BC cells. These findings demonstrate a novel mechanism whereby Wnt/ß-catenin signaling pathway is regulated by different kinases in response to metabolic requirement of CSCs, and suggest that SIK2 inhibition may potentially be a strategy for eliminating BCSCs.

Natural killer cell-related gene signature predicts malignancy of glioma and the survival of patients.

BACKGROUND: Natural killer (NK) cells-based therapies are one of the most promising strategies against cancer. The aim of this study is to investigate the natural killer cell related genes and its prognostic value in glioma. METHODS: The Chinese Glioma Genome Atlas (CGGA) was used to develop the natural killer cell-related signature. Risk score was built by multivariate Cox proportional hazards model. A cohort of 326 glioma samples with whole transcriptome expression data from the CGGA database was included for discovery. The Cancer Genome Atlas (TCGA) datasets was used for validation. GO and KEGG were used to reveal the biological process and function associated with the natural killer cell-related signature. We also collected the clinical pathological features of patients with gliomas to analyze the association with tumor malignancy and patients' survival. RESULTS: We screened for NK-related genes to build a prognostic signature, and identified the risk score based on the signature. We found that NK-related risk score was independent of various clinical factors. Nature-killer cell gene expression is correlated with clinicopathological features of gliomas. Innovatively, we demonstrated the tight relation between the risk score and immune checkpoints, and found NK-related risk score combined with PD1/PDL1 patients could predict the patient outcome. CONCLUSION: Natural killer cell-related gene signature can predict malignancy of glioma and the survival of patients, these results might provide new view for the research of glioma malignancy and individual immunotherapy.