HAUS Augmin-Like Complex Subunit 1 Influences Tumour Microenvironment and Prognostic Outcomes in Glioma

Background The expression of HAUS Augmin-like complex subunit 1 (HAUS1), a protein-coding gene, is low in normal samples among various cancers with pan-cancer analysis. The depletion of HAUS1 in cells decreases the G2/M cell compartment and induces apoptosis. However, the detailed expression pattern of HAUS1 and its correlation with immune infiltration in glioma (LGG and GBM) (LGG: low-grade glioma; GBM: glioblastoma) remain unknown. Therefore, in this study, we examined the role and prognostic value of HAUS1 in glioma. Methods Transcriptional expression data of HAUS1 were collected from the CGGA and TCGA databases. The Kaplan–Meier analysis, univariate and multivariate Cox analyses, and receiver operating characteristic (ROC) curves were used to analyse the clinical significance of HAUS1 in glioma. The STRING database was used to analyse protein-protein interactions (PPI), and the “ClusterProfiler” package was used for functional enrichment analysis to examine the possible biological roles of HAUS1. In addition, the HAUS1 promoter methylation modification was analysed using MEXPRESS, and the association between HAUS1 expression and tumour-infiltrating immune cells was investigated using CIBERSORT. Results Based on the data retrieved from TCGA (703 samples) and CGGA (1018 samples), an elevated expression of HAUS1 was observed in glioma samples, which was associated with poorer survival of patients, unfavourable clinical characteristics, 1p/19q codeletion status, WHO grade, and IDH mutation status. Furthermore, multivariate and univariate Cox analyses revealed that HAUS1 was an independent predictor of glioma. HAUS1 expression level was associated with several tumour-infiltrating immune cells, such as Th2 cells, macrophages, and activated dendritic cells. The outcomes of ROC curve analysis showed that HAUS1 was good to prognosticate immune infiltrating levels in glioma with a higher area under the curve (AUC) value (AUC = 0.974). Conclusions HAUS1 was upregulated and served as a biomarker for poor prognosis in patients with glioma. High HAUS1 expression was associated with several tumour-infiltrating immune cells such as Th2 cells, macrophages, and activated dendritic cells, which had high infiltration levels. Therefore, these findings suggest that HAUS1 is a potential biomarker for predicting the prognosis of patients with glioma and plays a pivotal role in immune infiltration in glioma.


Introduction
Glioma is a life-threatening highly malignant primary brain tumour with a high mortality rate [1,2]. However, owing to diffusive and invasive growth, the treatment of glioma often leads to a dismal prognosis [3,4]. Moreover, the overall survival of patients with glioma is poor, with the 5-year survival rate of adult patients being <5% [5,6]. Chemotherapy resistance of glioma may be the main cause of poor prognosis [7]. Nowadays, lacking of corresponding molecular markers and understanding of the mechanism of glioma pathogenesis bring about a huge challenge to the early diagnosis of glioma [8]. Although there are some predictive molecular markers such as 1p/19q codeletion and IDH mutations, the high complexity of the molecular mechanism of glioma urgently requires more prognostic markers for glioma [9]. erefore, it is imperative to discover new therapeutic strategies to enhance the overall survival of patients with glioma. e correlation between immunotherapy and the tumour microenvironment (TME) has been receiving increasing attention [10]. TME plays an important role in chemotherapy, radiotherapy, and immunotherapy [11], and the efficacy of immunotherapy can be determined by TME, which differs among organs [12]. In addition, immunotherapy has been manifested to gradually become an integral component of cancer therapy [13]. Studies have revealed that analysing the characteristics of TME, especially immune cell infiltration, can improve the understanding of the likelihood of immunotherapeutic response or survival [14]. In recent years, owing to the emergence of novel immunotherapy strategies involving PD-1 and PD-L1 and the clinical efficacy of PD-1 and PD-L1 blockade in various solid cancers, blocking of the PD-1 and PD-L1 pathways has emerged as an effective strategy for cancer treatment [15,16]. PD-L1 has been recently identified as a frequently used biomarker for anti-PD-1-based immunotherapy [17]. Moreover, clinical trials have revealed that integrating PD-1/PD-L1 with CTLA4 blockade treatment can provide improved therapeutic benefits compared to single blockade treatment [18]. However, only a small proportion of patients can benefit from anti-CTLA4 or anti-PD-1/PD-L1 immunotherapy [19]. erefore, identifying novel immune-related biomarkers is necessary to improve the efficacy of immunotherapy.
HAUS augmin-like complex subunit 1 (HAUS1) is a protein-coding gene that belongs to the HAUS protein family that contains HAUS1-8 genes. HAUS1 encodes 1 of the 8 subunits of the 390 kDa human augmin complex, which participates in the generation of microtubules inside the mitotic spindle [20,21]. In addition, HAUS1 may play a major role in centrosome integrity and completion of cytokinesis [22]. Microtubules are widely involved in intracellular life activities, including cell morphology maintenance, cell transport, cell migration, cell division, and cell signal transduction. erefore, as one of the augmin complexes involved in microtubule generation, HAUS1 can participate in multiple human cell life activities. However, the role of HAUS1 in tumours remains unclear, especially, its relationship with immune infiltration and TME.
In this study, we investigated the prognostic significance of HAUS1 mRNA expression and methylation in patients with glioma using data from the Chinese Glioma Genome Atlas (CGGA) and e Cancer Genome Atlas (TCGA) databases. In addition, we evaluated the enrichment functions of HAUS1 using GSEA and investigated the correlation between HAUS1 expression and infiltration levels of different tumour-infiltrating lymphocytes (TILs).
e clinical parameters of the patients with complete clinical data in TCGA and CGGA databases are listed in Table 1, respectively. In addition, pan-cancer analyses of HAUS1 were conducted using TCGA dataset and GEPIA (http://gepia.cancer-pku. cn/). All downloaded samples were divided into the low-and high-expression groups based on the median HUAS1 expression. Because all data were downloaded from public open-access databases (TCGA and CGGA), it was not necessary to obtain ethics approval or informed consent. Figure 1 shows the workflow of this study.

Differential Expression
Analysis. GEPIA (http://gepia. cancer-pku.cn/) is an online database for gene expression profiling and interactive analyses of tumour and normal samples; it comprises 9736 tumour samples and 8587 normal samples [24,25]. In this study, we used GEPIA to examine the expression profile of HAUS1 in glioma and conduct a pan-cancer analysis to examine HAUS1 expression level in 33 cancers.

Protein-Protein Interaction Analysis and Functional Enrichment Analysis.
e STRING database (http://stringdb.org/) is the largest database of protein-protein interactions for predicting functional relationships between proteins [26,27]. In this study, we used the STRING database to screen for co-expression genes and build PPI networks with an interaction score of >0.4. Subsequently, we used the "ClusterProfiler" and "ggplot2" packages for Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis to examine the enrichment function of these genes [28].

Methylation
Analysis of HAUS1. MEXPRESS (https:// mexpress.be/) is an easy-to-use online tool for integration and visualisation of clinical expression and DNA methylation data of cancer and normal tissues from TCGA database [29,30]. In this study, we used the MEXPRESS online tool to estimate and visualise DNA methylation of HAUS1 in TCGA dataset.

Analysis of Differentially Expressed Genes (DEGs).
Data retrieved from TCGA database were divided into the low-and high-expression groups based on the median HAUS1 expression level. Differentially expressed genes (DEGs) between the two groups were screened using the limma, heatmap, and ggplot2 packages in R. ese DEGs were visualised on a volcano map, and the top 15 genes were visualised on heat maps. e cut-off values for screening DEGs were set as |log 2 fold change| > 1.5 and adjusted P value <0.001.
2.6. Functional Enrichment Analysis. GSEA was performed to investigate the biological pathways that differed substantially between the high-and low-HAUS1-expression groups. e "c2.cp.kegg.v7.4.symbols.gmt" and "c5.all.v7.4.symbols.gmt" gene sets were used to evaluate biological processes related to HAUS1. e analysis was performed using 1000 permutations, and P-value <0.05 and (FDR) q value <0.05 were set as the cut-off criteria. e CIBERSORT algorithm was used to examine the association between the expression of HAUS1 and the infiltration of 22 types of TILs in glioma. CIBERSORT is an online tool that uses a deconvolution algorithm to sensitively and specifically calculate the proportion of TILs [31,32]. We divided CGGA samples into the low-and high-expression groups according to the median expression level of HAUS1.

Quantitative RT-PCR.
Total RNA was isolated from the paraneoplastic and tumour tissues of patients with glioma of different grades using the TRIzol reagent (Sigma-Aldrich, St. Louis, MO, USA). RNA from each sample (2 μg) was reverse transcribed to cDNA, which was used as a template in a 20-μL reaction mixture (10 μL of PCR mixture, 0.5 μL of forward and reverse primers each, 2 μL of cDNA template, and an appropriate volume of water). Subsequently, quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed using FastStart Universal SYBR ® Green Master (Roche, Germany) on an ABI QuantStudio5 Q5 realtime PCR System ( ermo Fisher Scientific, USA). e PCR conditions were maintained as follows: initial DNA denaturation at 95°C for 30 s, followed by 45 cycles at 94°C for 15 s, 56°C for 30 seconds, and 72°C for 20 seconds. Each sample was examined in triplicate. e threshold cycle (CT) readings were recorded and normalised to glyceraldehyde 3phosphate dehydrogenase (GAPDH) levels in all samples using the 2 −ΔΔCT method. e mRNA expression levels of samples were compared to those of para-cancerous tissue controls. e sequences of primer pairs for target genes are shown in Table 2. 2.9. Western Blot. Paraneoplastic tissue and tumour tissue samples from patients with glioma of different grades were lysed in RIPA buffer (Solarbio, Beijing, China) supplemented with protease and phosphatase inhibitors for 30 min, followed by denaturation at 95°C for 10 min. e protein samples were then separated via SDS-PAGE and transferred onto polyvinylidene fluoride (PVDF) membranes. Subsequently, the membranes were blocked with 5% skim milk powder solution for 1 h and incubated overnight with primary antibodies, including anti-HAUS1 antibody (1 : 500, Abclonal, A17250) and anti-GAPDH antibody (1 : 5000, Abcam, ab8245). e following day, the membranes were incubated with secondary antibodies (1 : 5000) for 2 h at room temperature. ereafter, an ECL kit (Billerica Millipore, MA, USA) was used for observation, and protein bands were visualised on a ChemiDoc detection system (Bio-Rad, Hercules, CA, USA) and quantified using the ImageJ software (National Institutes of Health, USA).

Statistical
Analysis. Statistical analyses of TCGA and CGGA data were conducted using various R packages (version 3.6.1). e relationship between HAUS1 and different clinical characteristics was evaluated using logistic regression. In addition, multivariate and univariate Cox analyses were performed to analyse the correlation between different clinical parameters and HAUS1 expression. e results of Cox analyses are presented in Table 3. e Kaplan-Meier curve revealing progression-free interval (PFI), disease-specific survival (DSS), and overall survival (OS) was constructed using the survival package in R, and time-dependent receiver operating characteristic (ROC) curves were constructed using the survival ROC package. Statistical significance was set as P < 0.05.

Pan-cancer Analysis of the Expression Pattern and
Prognostic Significance of HAUS1. Based on the data retrieved from TCGA and GTEx databases, we discovered that HAUS1 exhibited significantly elevated expression in a majority of human cancers, except for ACC, ESCA, KICH, MESO, PCPG, LAML, OV, PRAS, SARC, and UVM (Supplementary Figure 1A). Subsequently, we used GEPIA2  Figure 1B) and found that HAUS1 expression was significantly high in DLBC, GBM, LGG, LIHC, PAAD, SKCM, TGCT, and THYM. erefore, these results revealed that HAUS1 expression was significantly higher in DLBC, GBM, LGG, LIHC, PAAD, SKCM, TGCT, and THYM. Owing to the overexpression of HAUS1, we further evaluated its prognostic significance in these cancers. Because GBM and LGG are types of brain tumours, their data were combined for analysis. We analysed the association between OS and HAUS1 expression in DLBC, GBM, LGG, LIHC, PAAD, SKCM, TGCT, and THYM using the GEPIA2 database. As shown in Supplementary Figures 1C-F and 2A-B and Figure 2(b), HAUS1 was significantly associated  with OS in glioma (LGG and GBM) and LIHC. However, HAUS1 expression was higher and more significantly associated with OS in glioma (LGG and GBM) than in LIHC.
As K-M depends on the cut-off value, and the efficiency of the test is not as high as Cox, we further investigated the prognostic value of HAUS1 in pan cancer by using the tool of Sangerbox (Supplementary Figure 2C). And the results showed that the expression level of HAUS1 plays a valuable role in glioma. e high expression level and positive prognostic value of HAUS1 in glioma (LGG and GBM) suggested that patients with LGG and GBM might be eligible candidates for anti-HAUS1 immunotherapy. erefore, in this study, we mainly discussed the underlying mechanisms and biological functions of HAUS1 in patients with glioma.

Expression Level and Prognostic Value of HAUS1 in
Glioma Samples from TCGA and CGGA. We compared the mRNA level of HAUS1 between normal and glioma tissue samples obtained from GTEx and TCGA databases. In Table 1, we can see that HAUS1 expression level in glioma patients has a prognostic value in WHO grade (P ≤ 0.01), IDH mutation status (P ≤ 0.01), and 1p/19q codeletion status (P ≤ 0.01) based on the CGGA and TCGA database. As shown in Figure 2(a), HAUS1 expression was higher in glioma tissues than in normal tissues (P < 0.001). Furthermore, to assess whether HAUS1 could serve as a prognostic marker for glioma, we performed Kaplan-Meier survival analysis to compare the PFI, DSS, and OS of patients with glioma between the low-and high-HAUS1-expression groups. In TCGA dataset, patients in the high-HAUS1expression group had a worse OS (HR � 3.58, P < 0.001) (Figure 2(b)). Similarly, in the CGGA dataset, reduced expression levels of HAUS1 were correlated with a better OS (P < 0.001) (Figure 2(c)). In addition, elevated HAUS1 expression level was correlated with poorer DSS (HR � 3.79, P < 0.001) and PFI (HR � 2.72, P < 0.001) (Figures 2(d) and 2(e)).
Furthermore, time-dependent ROC curves were constructed to investigate the prognostic value of HAUS1 in predicting 1-, 3-, and 5-year OS, DSS, and PFI based on TCGA data. As demonstrated in Figure 2(f ), the AUC values for predicting the 1-, 3-, and 5-year OS of patients with glioma were 0.751, 0.815, and 0.756, respectively. As shown in Figure 2(g), the AUC values for predicting the 1-, 3-, and 5-year DSS of patients with glioma were 0.750, 0.807, and 0.771, respectively. In addition, the AUC values for predicting the 1-, 3-, and 5-year PFI of patients with glioma were 0.741, 0.718, and 0.705, respectively ( Figure 2(h)). Furthermore, based on the CGGA data, the AUC values for predicting 1-, 3-, and 5-year OS were 0.646, 0.722, and 0.744, respectively ( Figure 2(i)). erefore, these results revealed that HAUS1 is a moderately sensitive index for predicting the prognosis of glioma and can serve as an effective biomarker for glioma.

Relationship between HAUS1 Expression and Different
Clinical Characteristics. A total of 325 and 693 glioma samples from the CGGA database and 749 glioma samples from TCGA database were used to assess the relationship between HAUS1 overexpression and malignant behaviour of glioma. First, we investigated the differential expression of   0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Time ( HAUS1 in different subgroups stratified based on 1p/19q codeletion status, IDH mutation status, and WHO grade. e expression level of HAUS1 in Dataset 1 (ID: mRNA-seq_325) revealed an increasing trend from LGG (Grade II) to HGG (Grade IV) (Figure 4(a)). In Dataset 2 (ID: mRNAseq_693), the expression level of HAUS1 considerably increased from Grade II to Grade IV (Figure 4(b)). Similarly, in TCGA dataset, the expression level of HAUS1 showed an increasing trend from Grade II to Grade IV (Figure 4(c)). ese results suggested that advanced WHO grades were correlated with elevated expression levels of HAUS1. By analysing the relationship between HAUS1 and methylation, we found that HAUS1 was expression higher in higher MGMT methylation group (Figure 4(d)).
Furthermore, we investigated the relationship between mRNA expression of HAUS1 and WHO grades using western blot and qRT-PCR (Figures 4(e) and 4(f )). e results revealed that HAUS1 had a higher mRNA expression level in glioma tissues than in normal tissues. Moreover, the mRNA expression of HAUS1 increased with the increasing WHO grades. A high tumour grade often predicted a worse prognosis, and the mRNA expression of HAUS1 increased with the malignant progression of glioma.
ese results demonstrated that HAUS1 overexpression was related to malignant clinicopathological characteristics of glioma.    Elevated HAUS1 expression (red and green, Figure 6(a)) in the IDH1 mutation group indicated a dismal survival outcome, implying that HAUS1 was a significant prognostic factor for patients with glioma with the corresponding IDH1 genotypes (P < 0.001). Furthermore, we examined the relationship between HAUS1 expression and chemotherapy or radiotherapy. e findings revealed that lower expression levels of HAUS1 without chemotherapy indicated optimal survival outcomes (purple, Figure 6(b)), whereas higher expression levels of HAUS1 were correlated with the worst survival outcomes (red, Figure 6(b)). erefore, these results suggested that improved survival outcomes were correlated with low-HAUS1 expression without chemotherapy. As shown in Figure 6(c), we could obviously see that low-expression level of HAUS1 with or without radiotherapy revealed better survival outcomes compared to the highexpression level of HAUS1. Finally, we examined the relationship between HAUS1 expression and 1p/19q codeletion status and found that the poorest survival outcomes were correlated with elevated HAUS1 expression without 1p/19q codeletion (green, Figure 6(d)).

PPI Networks and Functional
Annotations. e STRING database was used to establish PPI networks, and GO and KEGG were used for functional enrichment analyses. As shown in the PPI network in Figure 7(a), 10 co-expressed genes were closely correlated with HAUS1. As demonstrated in Figure 7(b), the biological mechanisms related to HAUS1 were associated with microtubules, spindle, and spindle organisation. e functional annotations revealed that these genes were associated with tubulin binding, microtubule binding, and microtubule minus-end binding. e relationship between HAUS1 expression and co-expressed genes in glioma based on TCGA data is demonstrated in Figures 7(c)-7(l).

Functional Enrichment Analysis of HAUS1.
After analysing the correlation between HAUS1 and its coexpressed genes, we performed functional enrichment analysis in the low-and high-HAUS1-expression groups. To assess the underlying mechanisms of HAUS1 in promoting tumour progression, we screened for DEGs in the low-and high-HAUS1-expression groups. A total of 1457 DEGs were identified; of which, 1203 were upregulated and 254 were downregulated. A heat map and volcano plot were used to demonstrate the expression of these DEGs (Figures 8(a) and 8(b)). Furthermore, GSEA was performed to investigate key pathways associated with (e) Representative western blot images of HAUS1 and GAPDH from three groups. GAPDH was used as a protein-loading control. Relative protein levels of HAUS1 in the three groups are expressed as mean ± S.E.M of three independent experiments ( * P < 0.05, * * P < 0.01). (f ) qRT-PCR results showing significant differences in mRNA expression of HAUS1; data are expressed as fold change compared to the NORM group and as mean ± SEM of three independent experiments ( * * P < 0.01, * * * * P < 0.0001). 8 Journal of Oncology HAUS1. As shown in Table 4, P < 0.05 and FDR < 0.05 were used as the cut-off criteria. KEGG analysis revealed five categories that were positively associated with elevated expression of HAUS1: cell cycle, P53 signalling pathway, nucleotide sugar, and amino sugar metabolism, pyrimidine metabolism, and degradation of other glycans. In addition, KEGG analysis revealed four negatively correlated categories: long-term depression; neuroactive ligand-receptor interaction; long-term potentiation and alanine, aspartate, and glutamate metabolism (Figure 9(a)). GO analysis revealed five categories that were positively correlated with high expression of HAUS1: modulation of haematopoietic stem cell differentiation, modulation of stem cell differentiation, G1-S phase transition in the cell cycle, and transduction of innate immune response-activating signal. In addition, GO analysis revealed five negatively correlated categories: neurotransmitter receptor activity, glutamate secretion, glutamate metabolism, glutamate receptor signalling pathway, and modulation of glutamate receptor signalling pathway (Figure 9(b)). ese results showed that pathways regulating cell and body metabolism and amino acid catabolism, which are important for patients with glioma, were strongly associated with HAUS1 expression.

Relationship between HAUS1 Expression and DNA Methylation and Replication.
To further examine the underlying mechanisms of HAUS1 upregulation in glioma samples, we analysed the association between HAUS1 expression and its methylation. We found a methylation site (cg06121461) in the DNA sequence of HAUS1 that was negatively correlated with its expression in GBM samples (Supplementary Figure 3A). In LGG samples, we found six methylation sites (cg18069568, cg20479805, cg26626598, cg21686188, cg26713775, and cg04376617) that were negatively correlated with HAUS1 expression (Supplementary Figure 3B). Furthermore, we performed differential expression analysis of four DNA methyltransferases (DNMT1, DNMT2, DNMT3A, and DNMT3B) and assessed their correlation with HAUS1 expression by using Sangerbox tool. As demonstrated in Supplementary Figure   LGG samples, whereas, it was positively correlated with only DNMT2 (R � 0.22, P � 0.014) in GBM samples. In addition, we evaluated the association between HAUS1 expression and DNA replication. Mismatch repair (MMR) is a wellrecognised intracellular mechanism. When a critical gene function in MMR is lost, it leads to errors in DNA replication that cannot be repaired, resulting in the production of a large number of somatic mutations. In this study, we used TCGA expression profile data to assess the association between HAUS1 expression and mutations of five MMR genes, namely, EPCAM, PMS2, MLH1, MSH2, and MSH6. In GBM samples, HAUS1 expression was significantly positively correlated with MLH1 and MSH2. In LGG samples, HAUS1 expression was positively and strongly correlated with PMS2, MLH1, MSH2, and MSH6 and negatively correlated with EPCAM (Supplementary Figure 3D). ese results indicated that HAUS1 could mediate tumorigenesis and progression of glioma.

Relationship between HAUS1 Expression and Tumour
Immune Microenvironment in Glioma. TIL, a vital component of TME, plays a pivotal role in modulating tumour growth and progression. We used CIBERSORT to extract and process the signature gene expression profiles of immune cells to create a systematic representation of immune cell patterns. A total of 22 types of TILs were screened using TCGA and CGGA samples, and their relationship with HAUS1 expression in patients with glioma was further assessed (Supplementary Figure 4A-B). In addition, we investigated whether the tumour-immune microenvironment was different between patients with elevated HAUS1 expression and those with low-HAUS1 expression. e 689 glioma samples from TCGA and 749 glioma samples from CGGA were divided into the high-and low-expression groups based on the median HAUS1 expression. e infiltration levels of M0 macrophages, M2 macrophages, and gammadelta T cells were higher in the high-expression group than in the low-expression group, whereas, those of monocytes were lower in TCGA and CGGA datasets (Figures 10(a) and 10(b)). In addition, we investigated the relationship among the 22 types of TILs (Figures 11(a) and 11(b)). e heat map revealed a weak-to-moderate correlation among the ratios of various types of immune cells. HAUS1 was negatively correlated with monocytes; therefore, we subsequently examined the relationship between HAUS1 expression and marker genes of monocytes (Figure 11(c)). e findings revealed that HAUS1 was negatively associated with HIVEP2 (R � −0.43) and MBP (R � −0.29).
In addition, as previous studies have confirmed that the infiltration of immune cells in TME is closely correlated with checkpoint blockade therapy [33], we investigated the association between the expression of HAUS1 and distinct immune checkpoint molecules (Figure 12). e results showed that most immune checkpoints examined were positively correlated with HAUS1 expression, such as CD276 (R � 0.

Discussion
Gliomas are highly heterogeneous tumours, including lowgrade glioma (LGG) to high-grade glioma (HGG), and are the most common primary central nervous system (CVS) tumours in humans [34,35]. Based on the classification proposed by the World Health Organization (WHO), glioma can be classified as astrocytoma, oligodendroglioma, or mixed oligoastrocytoma [34]. Although multiple clinical treatments are available for glioma, such as radiotherapy, surgery, oral medication, and chemotherapy, the prognosis of patients remains unsatisfactory [36]. Moreover, currently available therapeutic interventions are known to exhibit severe short-term and long-term side effects, such as neurocognitive deficits, sterility, endocrinopathies, and postoperative mutism. Furthermore, patients with recurring illnesses have a notably unsatisfactory prognosis after first therapy, with a median survival of <6 months [37]. erefore, developing more efficient treatment methods for glioma is necessary. Previous studies have shown that glioma cells may act against various constituents in their microenvironment, resulting in the establishment of an immunosuppressive microenvironment and the progression of glioma [38]. In this context, modulating the glioma microenvironment while targeting glioma cells with key biomarkers simultaneously appears to be a feasible therapeutic approach to treat patients with glioma [39].      IKZF1   MARCKSL1  MMP15  ASGR1  PQBP1  MBP  PNPLA6  ATP6V1B2  HIVEP2  TMBIM6  DAZAP2  CFP  TEX264  HAUS1  ACTG1  CFL1  UPK3A  ANXA5  CTBS  CD1D  ASGR2  In this study, we discovered that the mRNA expression of HAUS1 was higher in glioma tissues than in normal tissues and increased with tumour grade. However, little information is available regarding the association between HAUS1 expression and glioma risk and prognosis. e pan-cancer analysis in this study revealed that HAUS1 expression was high in most cancers but only had a high statistical significance in the glioma tissue samples. We found that HAUS1 overexpression was associated with PRS type, histological characteristics, age, grade, 1p/19q codeletion status, chemotherapy status, and IDH mutation status. erefore, we speculated that HAUS1 overexpression might be correlated with the malignant behaviour of glioma. Furthermore, multivariate and univariate Cox analyses confirmed that HAUS1 is an independent prognostic factor for glioma.
e Kaplan-Meier curves demonstrated that the elevated HAUS1 expression was associated with a shorter OS, PFI, and DSS, and ROC curves demonstrated that HAUS1 might serve as an effective diagnostic biomarker for the differentiation of glioma from healthy tissues. Given that HAUS1 has a good prognostic value, we built a nomogram that integrated HAUS1 expression with various clinical characteristics, and the results revealed that HAUS1 could accurately predict the 1-, 3-, and 5-year survival of patients with glioma. e presence or absence of the 1p/19q codeletion and IDH mutation were proposed as a novel glioma classification approach by the WHO [40]. erefore, in this study, we investigated the correlation of HAUS1 with 1p/19q codeletion status and IDH mutation status. Moreover, the correlation of HAUS1 expression with chemotherapy and radiotherapy was also investigated. Functional enrichment analysis revealed that elevated HAUS1 expression was associated with cell cycle, P53 signalling pathway, and regulation of stem cell differentiation, which was consistent with the  TNFSF9   CD200   TNFSF18   TNFRSF8   VSIR   TMIGD2   VTCN1   CD160   ADORA2A   LAG3   TNFRSF25   ICOSLG  HAUS1  CD276  NRP1  TNFRSF4  CD244  TNFRSF18  TNFSF15  CD200R1  CD28  TNFSF4  CD70  TNFRSF9  IDO1  CD40  LGALS9  LAIR1  CD86  HAVCR2  CD80  PDCD1  CD44  CD274  PDCD1LG2  TNFSF14  TNFRSF14  CD48  BTLA  CTLA4  CD27  CD40LG  ICOS  IDO2  KIR3DL1  TIGIT  HHLA2  BTNL2  TNFSF9  CD200  TNFSF18  TNFRSF8  VSIR  TMIGD2  VTCN1  CD160  ADORA2A  LAG3  ese findings suggested that HAUS1 served as a prognostic biomarker for glioma.
HAUS1 is a protein-coding gene that is involved in microtubule nucleation, which is important in assembling mitotic spindle, and regulates NuMA, a crucial protein in spindle orientation [41]. In addition, HAUS1 plays a significant role in the maintenance of centrosome integrity and completion of cytokinesis [22]. In this study, we found that HAUS1 plays a crucial role in cell division and functional activity of cells. However, high HAUS1 expression was correlated with a poor OS in the patients with glioma. erefore, we speculated that high HAUS1 expression stimulated tumour cell division and survival, thus exacerbating the malignant behaviour of glioma.
Furthermore, we investigated the correlation of HAUS1 expression with the TME and immunotherapy. HAUS1 overexpression was correlated with high infiltration of gammadelta T cells, M0 macrophages, and M2 macrophages but with low infiltration of monocytes. Monocytes are major cells of the innate immune system and play a crucial role in adaptive immune response [42,43]. In addition, they play a significant role in TME and are correlated with cancer progression, immune escape and initiation, and regulation of immune and inflammatory responses [44,45]. Previous studies have reported that monocytes are used as immune targets in arterial hypertension and act as pivotal immune cells in sepsis [46,47]. In glioma, the most infiltrating immune cells are monocytes, along with microglia, accounting for approximately one-third of all the immune cells in glioma tissues [48]. erefore, we investigated the relationship between HAUS1 and glioma. However, this study has a few limitations. First, we did not further explore and verify the functional role of HAUS1 in the immune microenvironment of glioma. However, this subject is rather novel and warrants further investigation. Second, no clinical samples were used to verify the results. However, we analysed and verified the results based on a large amount of data from public databases and in vivo, which may provide guidance and a reference for future studies.