High STK40 Expression as an Independent Prognostic Biomarker and Correlated with Immune Infiltrates in Low-Grade Gliomas.

BACKGROUND: Expression of STK40 is observed in some cancer types, while its role in low-grade gliomas (LGG) is unclear. The present study aimed to demonstrate the relationship between STK40 and LGG based on The Cancer Genome Atlas (TCGA) database and bioinformatics analysis.
METHODS: Kruskal-Wallis test, Wilcoxon sign-rank test, and logistic regression were used to evaluate the relationship between clinicopathological features and STK40 expression. Kaplan-Meier method and Cox regression analysis were used to evaluate prognostic factors. Gene set enrichment analysis (GSEA) and immuno-infiltration analysis were used to determine the significant involvement of STK40 in function.
RESULTS: High STK40 expression in LGG was associated with WHO grade (P<0.001), IDH status (P<0.001), primary therapy outcome (P=0.027), 1p/19q codeletion (P<0.001) and histological type (P<0.001). High STK40 expression predicted a poorer overall survival (OS) (HR: 3.07; 95% CI: 2.09-4.51; P<0.001), progression-free survival (PFS) (HR:2.11; 95% CI: 1.59-2.81; P<0.001) and disease specific survival (DSS) (HR: 3.27; 95% CI: 2.17-4.92; P<0.001). STK40 expression (HR: 2.284; 95% CI: 1.125-4.637; P=0.022) was independently correlated with OS in LGG patients. GSEA demonstrated that pathways including cell cycle mitotic, neutrophil degranulation, signaling by Rho GTPases, signaling by interleukins, M phase, PI3K-Akt signaling pathway and naba secreted factors were differentially enriched in STK40 high expression phenotype. Immune infiltration analysis showed that STK40 expression was correlated with some types of immune infiltrating cells.
CONCLUSION: STK40 expression was significantly correlated with poor survival and immune infiltration in LGG, and it may be a promising prognostic biomarker in LGG.

Introduction

Gliomas are the most common intracranial tumor, accounting for approximately 70% of all primary brain tumors.1 In 2016, the WHO provided a molecular classification of gliomas and noted that molecular features have an important impact on the diagnosis and prognosis of gliomas.2 Gliomas are generally classified into low-grade gliomas (LGG) and glioblastoma multiforme (GBM).3 GBM is the most aggressive subtype of glioma, with a median overall survival (OS) of 12–16 months.4 LGG is less aggressive and the clinical outcome of LGG patients varies widely.5 In patients with LGG, gross total resection is associated with improved progression-free survival (PFS) and OS.6 The median survival time for patients with low-grade gliomas is approximately 5 to 10 years.7 Clear diagnosis and molecularly targeted therapy are important to improve the survival time and quality of life of patients with LGG. However, there are few valid and reliable biomarkers that can predict poor prognosis and direct treatment strategies in LGG. Therefore, it is of great importance to accurately predict the prognosis of LGG patients using novel and accessible biomarkers.

Serine/threonine kinase 40 (STK40), a target gene of transcriptional factor Oct4, is capable of activating the Erk/MAPK signaling pathway.8 Studies have shown that Stk40 deficiency leads to disturbances in adipogenesis, myogenesis, and erythropoiesis.9–11 STK40 could be a novel therapeutic target candidate for Triple-negative breast cancers (TNBC).12 However, the correlation between STK40 and LGG has not been studied. This study aims to explore the expression of STK40 in LGG, which may provide new directions for the development of diagnostic and therapeutic strategies for LGG.

Bioinformatics, as a science combining molecular biology and information technology, has discovered many clinically applicable tumor markers by studying the molecular mechanisms of diseases at the molecular level through data mining.13 These discoveries are important for improving the early diagnosis and prognosis of tumors and revealing the molecular mechanisms of tumor pathogenesis.14 In this study, the correlation between STK40 expression levels and clinicopathological features of LGG was explored based on the Cancer Genome Atlas (TCGA) database and LGG RNA-seq data in GTEx, comparing the difference in STK40 expression between tumor tissues and normal samples. The prognostic value of STK40 in LGG was assessed. Genomic enrichment analysis (GSEA) was performed on the high and low expression groups of STK40 to reveal its possible functions. By analyzing the correlation between STK40 expression and immune infiltration, the potential mechanisms by which STK40 regulates LGG occurrence and development were comprehensively explored and discussed.

Materials and Methods

Differential Expression of STK40

Baseline Information Sheet

The analysis method is described in the previous literature.15 Molecule: STK40[ENSG00000196182]. Data: RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA () LGG project. Owing to the freely accessible resource TCGA database, the study let off the institutional review board approval of the ethics committee of The Third People’s Hospital of Huai’an. The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). All data for this study were obtained from TCGA database and we did not obtain these data directly from patients or animals, nor did we intervene in these patients. Therefore, no ethical approval was necessary.

Unpaired Samples

The analysis method is described in the previous literature.15 Molecule: STK40. Data: UCSC XENA () RNAseq data in TPM format for TCGA and GTEx processed in unison by the Toil process.16 Extraction of LGG from TCGA and corresponding normal tissue data from GTEx.

ROC Analysis

The analysis method is described in the previous literature.15 Molecule: STK40. Data: UCSC XENA () RNAseq data in TPM format for TCGA and GTEx processed in unison by the Toil process.16 Extraction of LGG from TCGA and corresponding normal tissue data from GTEx.

The Relationship Between STK40 and Clinical Characteristics

Correlation of Gene Expression with Clinical Characteristics

The analysis method is described in the previous literature.15 Molecule: STK40. Clinical variables: WHO grade, 1p/19q codeletion, primary therapy outcome, IDH status, and histological type. Data: RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA LGG project.

Logistics Analysis

The analysis method is described in the previous literature.15 Dependent variable: STK40. Data: RNAseq data in level 3 HTSeq-FPKM format from the TCGA LGG project.

The Relationship Between STK40 and Prognosis

Kaplan–Meier Method

The analysis method is described in the previous literature.15,17 Molecule: STK40. Prognosis type: OS, PFS, and disease specific survival (DSS). Data: RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA LGG project.

COX Regression

The analysis method is described in the previous literature.17,18 Prognosis type: OS. Included variables: clinical characteristics and STK40. Data: RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA LGG project.

Forest Plot

The analysis method is described in the previous literature.18

Nomogram Plot

R package: rms package and survival package. Statistical method: Cox. Prognosis type: OS. Included variables: WHO grade, 1p/19q codeletion, primary therapy outcome, IDH status, age, histological type, and STK40. Data: RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA LGG project. Additional data: prognostic data from article.17

Gene Set Enrichment Analysis (GSEA)

Single Gene Differential Analysis

The analysis method is described in the previous literatures.18,19 Molecule: STK40. Data: RNAseq data in level 3 HTSeq-Counts format from the TCGA LGG project.

GSEA Analysis

The analysis method is described in the previous literatures.18,20,21

Immune Infiltration Analysis by ssGSEA

The analysis method is described in the previous literature.18,22,23 Molecule: STK40. Data: RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA LGG project.

Results

Clinical Characteristics

As shown in Table 1, a total of 528 patients were in this study. According to WHO grading, we had 224 patients of Grade II (48%) and 243 patients of Grade III (52%). According to the IDH status, we had 97 patients (WT, 18.5%) and 428 patients (Mut, 81.5%). According to the 1p/19q codeletion, we had 171 patients (codel, 32.4%), and 357 patients (non-codel, 67.6%). According to the 1p/19q codeletion, we had 110 patients (PD, 24%), 146 patients (SD, 31.9%), 64 patients (PR, 14%), and 138 patients (CR, 30.1%). According to the gender, we had 239 patients (female, 45.3%) and 289 patients (male, 54.7%). According to race, we had 487 white patients, 8 Asian patients, and 22 Black or African American patients. According to the age, we had 264 patients (≤40, 50%) and 264 patients (>40, 50%). According to the histological type, we had 195 patients (astrocytoma, 36.9%), 134 patients (oligoastrocytoma, 25.4%), and 199 patients (oligodendroglioma, 37.7%). According to the laterality, we had 256 patients (left, 48.9%), 6 patients (midline, 1.1%), and 261 patients (right, 49.9%). The age range was 32 to 53 years, with a median of 40.5 years.

Table 1

Clinical Characteristics of the LGG Patients Based on TCGA

Characteristics	Levels	Overall
n		528
WHO grade, n (%)	G2	224 (48%)
	G3	243 (52%)
IDH status, n (%)	WT	97 (18.5%)
	Mut	428 (81.5%)
1p/19q codeletion, n (%)	Codel	171 (32.4%)
	Non-codel	357 (67.6%)
Primary therapy outcome, n (%)	PD	110 (24%)
	SD	146 (31.9%)
	PR	64 (14%)
	CR	138 (30.1%)
Gender, n (%)	Female	239 (45.3%)
	Male	289 (54.7%)
Race, n (%)	Asian	8 (1.5%)
	Black or African American	22 (4.3%)
	White	487 (94.2%)
Age, n (%)	≤40	264 (50%)
	>40	264 (50%)
Histological type, n (%)	Astrocytoma	195 (36.9%)
	Oligoastrocytoma	134 (25.4%)
	Oligodendroglioma	199 (37.7%)
Laterality, n (%)	Left	256 (48.9%)
	Midline	6 (1.1%)
	Right	261 (49.9%)
Age, median (IQR)		40.5 (32, 53)

STK40 Expression is Correlated with Poor Clinical Characteristics of LGG

As shown in Figure 1, it was found that STK40 was highly expressed in LGG tissues (3.535 ± 0.023 vs 4.117 ± 0.031, P<0.001) based on 523 LGG tissues and 1152 normal brain tissues of GTEx combined TCGA database. As shown in Table 2, STK40 expression was associated with WHO grade (P<0.001), IDH status (P<0.001), 1p/19q codeletion (P<0.001), primary therapy outcome (P<0.001) and histological type (P<0.001). As shown in Figure 2 and Table 3, STK40 was significantly related to WHO grade (P<0.001), IDH status (P<0.001), primary therapy outcome (P=0.027), 1p/19q codeletion (P< 0.001) and histological type (P<0.001).

Figure 1

STK40 is significantly upregulated in LGG than normal tissues. (A) The difference expression of TLX1 in LUAD and unpaired normal lung tissues. (B) ROC curve. Significance markers: ***P<0.001.

Table 2

Correlation Between STK40 Expression and Clinicopathological Characteristics in LGG

Characteristic	Levels	Low Expression of STK40	High Expression of STK40	p	Method
n		264	264
WHO grade, n (%)	G2	147 (31.5%)	77 (16.5%)	< 0.001	Chisq.test
	G3	90 (19.3%)	153 (32.8%)
IDH status, n (%)	WT	13 (2.5%)	84 (16%)	< 0.001	Chisq.test
	Mut	249 (47.4%)	179 (34.1%)
1p/19q codeletion, n (%)	Codel	163 (30.9%)	8 (1.5%)	< 0.001	Chisq.test
	Non-codel	101 (19.1%)	256 (48.5%)
Primary therapy outcome, n (%)	PD	33 (7.2%)	77 (16.8%)	< 0.001	Chisq.test
	SD	80 (17.5%)	66 (14.4%)
	PR	34 (7.4%)	30 (6.6%)
	CR	79 (17.2%)	59 (12.9%)
Gender, n (%)	Female	125 (23.7%)	114 (21.6%)	0.382	Chisq.test
	Male	139 (26.3%)	150 (28.4%)
Race, n (%)	Asian	5 (1%)	3 (0.6%)	0.369	Fisher.test
	Black or African American	8 (1.5%)	14 (2.7%)
	White	244 (47.2%)	243 (47%)
Age, n (%)	≤40	124 (23.5%)	140 (26.5%)	0.192	Chisq.test
	>40	140 (26.5%)	124 (23.5%)
Histological type, n (%)	Astrocytoma	53 (10%)	142 (26.9%)	< 0.001	Chisq.test
	Oligoastrocytoma	59 (11.2%)	75 (14.2%)
	Oligodendroglioma	152 (28.8%)	47 (8.9%)
Laterality, n (%)	Left	125 (23.9%)	131 (25%)	0.602	Fisher.test
	Midline	2 (0.4%)	4 (0.8%)
	Right	134 (25.6%)	127 (24.3%)
Age, median (IQR)		41 (34, 53)	39 (31, 52.25)	0.151	Wilcoxon

Figure 2

Box plot evaluating STK40 expression of patients with LGG according to different clinical characteristics. (A) WHO grade. (B) 1p/19q codeletion. (C) IDH status. (D) primary therapy outcome. (E) histological type. Significance markers: **P<0.01; ***P<0.001.

Table 3

Logistic Analysis of the Association Between STK40 Expression and Clinical Characteristics

Characteristics	Total (N)	Odds Ratio (OR)	P value
WHO grade (G3 vs G2)	467	3.245 (2.228–4.760)	<0.001
1p/19q codeletion (non-codel vs codel)	528	51.644 (26.025–117.948)	<0.001
Primary therapy outcome (CR vs PD&SD&PR)	458	0.635 (0.423–0.948)	0.027
IDH status (Mut vs WT)	525	0.111 (0.058–0.199)	<0.001
Gender (Male vs Female)	528	1.183 (0.840–1.669)	0.336
Race (White vs Asian & Black or African American)	517	0.762 (0.355–1.596)	0.473
Age (>40 vs ≤40)	528	0.784 (0.557–1.104)	0.164
Laterality (Midline & Right vs Left)	523	0.919 (0.652–1.295)	0.630
Histological type (Oligodendroglioma vs Astrocytoma & Oligoastrocytoma)	528	0.160 (0.106–0.236)	<0.001

Role of STK40 in LGG Patient Survival

As shown in Figure 3, STK40 was positively correlated with poor OS (HR: 3.07; 95% CI: 2.09–4.51; P<0.001), PFS (HR:2.11; 95% CI: 1.59–2.81; P<0.001), and Disease specific survival (DSS) (HR: 3.27; 95% CI: 2.17–4.92; P<0.001) of LGG patients. As shown in Table 4, high STK40 expression levels were associated with worse OS (HR: 3.07; 95% CI: 2.09–4.51; P<0.001), WHO grade (HR: 3.059; 95% CI: 2.046–4.573; P<0.001), 1p/19q codeletion (HR: 2.493; 95% CI: 1.590–3.910; P<0.001*), primary therapy outcome (HR: 0.203; 95% CI: 0.094–0.437; P<0.001), IDH status (HR: 0.186; 95% CI: 0.130–0.265; P<0.001), age (HR: 2.889; 95% CI: 2.009–4.155; P<0.001), histological type (HR: 0.673; 95% CI: 0.471–0.961; P=0.029). As shown in Table 4 and Figure 4, STK40 (HR: 2.284; 95% CI: 1.125–4.637; P=0.022), WHO grade (HR: 1.735; 95% CI: 1.073–2.808; P=0.025), primary therapy outcome (HR: 0.276; 95% CI: 0.127–0.601; P=0.001), IDH status (HR: 0.353; 95% CI: 0.215–0.580; P<0.001), age (HR: 3.315; 95% CI: 2.092–5.254; P<0.001) were independently correlated with OS in multivariate analysis. The above data indicated STK40 is a prognostic factor and increased STK40 level is associated with poor OS. A nomogram was constructed to predict the 1-, 3-, and 5-year survival probability of LGG patients by combining the expression level of STK40 with clinical variables (Figure 5).

Figure 3

High expression of STK40 is associated with poor OS and DSS in patients with LGG. (A) The Kaplan–Meier curves and number at risk of OS in LGG. (B) The Kaplan–Meier curves and number at risk of PFS in LGG. (C) The Kaplan–Meier curves and number at risk of DSS in LGG.

Table 4

Univariate and Multivariate Cox Regression Analyses of Clinical Characteristics Associated with OS

Characteristics	Total (N)	Univariate Analysis		Multivariate Analysis
		Hazard Ratio (95% CI)	P value	Hazard Ratio (95% CI)	P value
WHO grade (G3 vs G2)	466	3.059 (2.046–4.573)	<0.001	1.735 (1.073–2.808)	0.025
1p/19q codeletion (non-codel vs codel)	527	2.493 (1.590–3.910)	<0.001	0.862 (0.377–1.975)	0.726
Primary therapy outcome (CR vs PD&SD&PR)	457	0.203 (0.094–0.437)	<0.001	0.276 (0.127–0.601)	0.001
IDH status (Mut vs WT)	524	0.186 (0.130–0.265)	<0.001	0.353 (0.215–0.580)	<0.001
Gender (Male vs Female)	527	1.124 (0.800–1.580)	0.499
Race (White vs Asian & Black or African American)	516	0.849 (0.395–1.823)	0.675
Age (>40 vs ≤40)	527	2.889 (2.009–4.155)	<0.001	3.315 (2.092–5.254)	<0.001
Histological type (Oligodendroglioma vs Astrocytoma & Oligoastrocytoma)	527	0.673 (0.471–0.961)	0.029	0.762 (0.457–1.272)	0.299
Laterality (Midline & Right vs Left)	522	0.776 (0.548–1.099)	0.154
STK40 (High vs Low)	527	3.071 (2.089–4.515)	<0.001	2.284 (1.125–4.637)	0.022

Figure 4

Forest plot of the multivariate Cox regression analysis for LGG.

Figure 5

Nomogram for predicting the probability of patients with 1-, 3- and 5-year overall survival.

STK40-Related Pathways Based on GSEA

There were 111 data sets which showed significantly differential enrichment in STK40 high expression phenotype, and we selected the top 9 data sets with high values of normalized enrichment score (NES), in Table 5 and Figure 6, including cell cycle mitotic, neutrophil degranulation, signaling by Rho GTPases, signaling by interleukins, M phase, Rho GTPase effectors, extracellular matrix organization, PI3K-Akt signaling pathway and naba secreted factors.

Table 5

Gene Sets Enriched in the High STK40 Expression Phenotype

Description	NES	P Adjust	FDR
REACTOME_CELL_CYCLE_MITOTIC	1.725	0.015	0.011
REACTOME_NEUTROPHIL_DEGRANULATION	1.949	0.015	0.011
REACTOME_SIGNALING_BY_RHO_GTPASES	1.548	0.015	0.011
REACTOME_SIGNALING_BY_INTERLEUKINS	2.249	0.015	0.011
REACTOME_M_PHASE	1.660	0.015	0.011
REACTOME_RHO_GTPASE_EFFECTORS	1.711	0.015	0.011
REACTOME_EXTRACELLULAR_MATRIX_ORGANIZATION	2.238	0.015	0.011
WP_PI3KAKT_SIGNALING_PATHWAY	1.637	0.015	0.011
NABA_SECRETED_FACTORS	1.748	0.015	0.011

Figure 6

Enrichment plots from gene set enrichment analysis (GSEA). (A) neutrophil degranulation, (B) signaling by interleukins, (C) GPCR-ligand binding, (D) G alpha (I) signaling events, (E) VEGFAVEGFR-2 signaling pathway, (F) Class A 1 Rhodopsin-Like Receptors, (G) naba secreted factors, (H) PI3K-Akt signaling pathway, and (I) Focal Adhesion-PI3K-Akt-mTOR-signaling pathway, were differentially enriched in STK40-related OC.

The Correlation Between STK40 Expression and Immune Infiltration

The correlation between expression of STK40 and immune infiltration by ssGSEA with Spearman r in Figure 7 and Table 6 showed that STK40 expression was negatively correlated with infiltration levels of (P=0.001), NK CD56bright cells (P<0.001), pDC (Plasmacytoid DC) (P<0.001) and Treg (P<0.001), and positively correlated with that of aDC (activated DC) (P<0.001), B cells (P=0.002), cytotoxic cells (P<0.001), eosinophils (P<0.001), iDC (immature DC) (P<0.0018), macrophages (P<0.001), neutrophils (P<0.001), NK CD56dim cells (P<0.001), NK cells (P<0.001), T cells (P<0.001), T helper cells (P<0.001), Tem (T effector memory) (P=0.044), Tgd (T gamma delta) (P<0.001), Th1 cells (P=0.031), Th17 cells (P<0.001) and Th2 cells (P<0.001).

Figure 7

Correlation of STK40 expression levels with immune infiltration (Spearman r).

Table 6

The Correlation Between STK40 Expression and Immune Cells Detected by Spearman Correlation Method

Gene Name	Cell Type	Correlation Coefficient (Spearman)	P value (Spearman)
STK40	aDC	0.357	<0.001
STK40	B cells	0.136	0.002
STK40	CD8 T cells	−0.047	0.284
STK40	Cytotoxic cells	0.301	<0.001
STK40	DC	−0.174	<0.001
STK40	Eosinophils	0.517	<0.001
STK40	iDC	0.345	<0.0018
STK40	Macrophages	0.57	<0.001
STK40	Mast cells	−0.044	0.316
STK40	Neutrophils	0.475	<0.001
STK40	NK CD56bright cells	−0.305	<0.001
STK40	NK CD56dim cells	0.179	<0.001
STK40	NK cells	0.346	<0.001
STK40	pDC	−0.357	<0.001
STK40	T cells	0.311	<0.001
STK40	T helper cells	0.315	<0.001
STK40	Tcm	−0.001	0.981
STK40	Tem	0.088	0.044
STK40	TFH	−0.033	0.455
STK40	Tgd	0.228	<0.001
STK40	Th1 cells	0.094	0.031
STK40	Th17 cells	0.282	<0.001
STK40	Th2 cells	0.313	<0.001
STK40	Treg	−0.298	<0.001

Discussion

With the development of molecular detection and sequencing technologies, the distinction between low-grade and high-grade gliomas has shifted from histology to molecular markers.24 KLF11 expression is an independent prognostic factor for poor survival in glioma patients.25 IRF7 could be a potential biomarker for the prognosis of LGG patients.26 APOB methylation may serve as a novel and promising prognostic biomarker for LGG patients.3 KIF4A, 9, 18A, 23 contribute to improved treatment and prognosis of LGG.27 The IRX1 gene exhibits oncogenic effects in glioma disease and may be clinically important as a future therapeutic target.28 SEMA3B is a potential target for the targeted treatment of gliomas.29 SASH1 expression was positively associated with better postoperative survival in glioma patients.30 Therefore, it is important to study mRNAs as new LGG biomarkers and therapeutic targets in the future.

STK40 is overexpressed in TNBC.12 In this study, High STK40 expression in LGG was associated with WHO grade (P< 0.001), IDH status (P< 0.001), primary therapy outcome (P=0.027), 1p/19q codeletion (P< 0.001) and histological type (P< 0.001). LGG expressed more STK40 than adjacent brain tissue, especially in patients with high grade (G3). These phenomena suggested that STK40 may be involved in tumorigenesis and promote proliferation. Primary therapy outcome (PD&SD&PR) may be due to a large tumor burden due to high STK40. ELF3 expression was high in LGG patients with WT IDH status and non-codel 1p/19q codeletion, the reasons for which need to be further investigated. High STK40 expression predicted a poorer OS (HR: 3.07; 95% CI: 2.09–4.51; P<0.001), PFS (HR:2.11; 95% CI: 1.59–2.81; P<0.001) and DSS (HR: 3.27; 95% CI: 2.17–4.92; P<0.001). And STK40 expression (HR: 2.284; 95% CI: 1.125–4.637; P=0.022) was independently correlated with OS in LGG patients.

The miR-31 target Stk40 was de-inhibited, thereby inhibiting the STK40-NF-κΒ-controlled inflammatory pathway, with the result that cell proliferation was reduced and apoptosis was activated in esophageal squamous cell carcinoma (ESCC).31 STK40 is a pseudokinase that binds the E3 biotin ligase COP1.32 STK40 is one of two genes specifically down-regulated shortly after LRP5 depletion.12 In this study, based on GESA, STK40 was related to pathways including cell cycle mitotic, neutrophil degranulation, signaling by Rho GTPases, signaling by Interleukins, M phase, Rho GTPase effectors, extracellular matrix organization, PI3K-Akt signaling pathway and naba secreted factors.

Another important aspect of this study was to investigate the relationship between STK40 expression and the level of diverse immune infiltration in LGG. Immune infiltration with LGG is currently a hot topic and knowledge of immune infiltrating cells is beneficial to the development of immunotherapy for LGG. The results showed modest correlations between STK40 expression and infiltration levels of DC, NK CD56bright cells, pDC and Treg, aDC, B cells, cytotoxic cells, eosinophils, iDC, macrophages, neutrophils, NK CD56dim cells, NK cells, T cells, T helper cells, Tem, Tgd, Th1 cells, Th17 cells and Th2 cells in LGG. These correlations may suggest potential mechanisms by which STK40 inhibits the function of DC, NK CD56bright cells, pDC, and Treg, and subsequently promotes the function of aDC, B cells, cytotoxic cells, eosinophils, iDC, macrophages, neutrophils, NK CD56dim cells, NK cells, T cells, T helper cells, Tem, Tgd, Th1 cells, Th17 cells and Th2 cells.

Despite some limitations, this is the first work to explore the relationship between STK40 and LGG. This study was performed based on RNA sequencing from the TCGA database. We were unable to illustrate the expression of STK40 at the protein level or to clearly evaluate the direct mechanism of STK40 involvement in LGG. However, the results of bioinformatics analysis need to be demonstrated by some experiments. Further studies are needed regarding the direct mechanism of STK40-mediated LGG.

Conclusion

STK40 was highly expressed in LGG relative to normal tissue and related to poor OS, PFS, and DSS. STK40 might participate in the development of LGG by pathways including cell cycle mitotic, neutrophil degranulation, signaling by Rho GTPases, signaling by interleukins, M phase, Rho GTPase effectors, extracellular matrix organization, PI3K-Akt signaling pathway and naba secreted factors and immune infiltrating cells. This study partly revealed the role of STK40 in LGG, providing a potential biomarker for the diagnosis and prognosis of LGG.