- Open Access
The oncogenic role of tubulin alpha-1c chain in human tumours
BMC Cancer volume 22, Article number: 498 (2022)
Tubulin alpha-1c chain (TUBA1C), a subtype of α-tubulin, has been shown to be involved in cell proliferation and cell cycle progression in several cancers and to influence cancer development and prognosis. However, a pancancer analysis of TUBA1C to reveal its immunological and prognostic roles has not been performed. In this study, we first downloaded raw data on TUBA1C expression in cancers from The Cancer Genome Atlas (TCGA) database and multiple other databases and analysed these data with R software to investigate the prognostic and immunological value of TUBA1C in cancers. Immunohistochemical analysis was performed in gliomas to further validate our findings. Overall, TUBA1C was overexpressed in most cancers, and overexpression of TUBA1C was linked to poor prognosis and higher tumour grade in patients. In addition, TUBA1C expression was associated with tumour mutation burden (TMB), microsatellite instability (MSI), the tumour microenvironment (TME) and the infiltration of immune cells. TUBA1C was also coexpressed with most immune-related genes and influenced immune-related pathways. Immunohistochemical analysis showed that TUBA1C expression was highest in glioblastoma (GBM) tissues, second highest in low-grade glioma (LGG) tissues and lowest in normal tissues. Our study indicated that TUBA1C might be a biomarker for predicting the immune status and prognosis of cancers, offering new ideas for cancer treatment.
Cancer, as a leading cause of death worldwide, is currently a major obstacle to increasing life expectancy, with its morbidity and mortality rates rising rapidly . The common therapeutic regimens for patients with advanced cancers, mainly including surgery, radiotherapy, chemotherapy and biologic therapies, are still fall short of expectations . Thus, it is absolutely imperative to identify new effective treatments to improve patient prognosis and quality of life. Cancer immunotherapy, dedicated to reactivating the antitumor immune response and blocking pathways that lead to immune escape, has received increasing attention; immune checkpoint blockade inhibitors (ICIs) have especially been noticed .
Microtubules are essential cytoskeletal components that play a key role in cell division, generation, motility and intracellular transport. Polymeric microtubules are assembled from the highly conserved subunits α- and β-tubulin . Tubulin alpha-1c chain (TUBA1C) is an isoform of α-tubulin. It has been demonstrated that TUBA1C plays a significant role in the cell cycle and immune microenvironment in lung adenocarcinoma (LUAD). Elevated expression of TUBA1C was correlated with poor outcome and with 13 tumour-infiltrating immune cells (TIICs) in LUAD . TUBA1C was found to be upregulated in hepatocellular carcinoma (HCC) and pancreatic ductal adenocarcinoma (PDAC), where it predicted poor prognosis and facilitated cell proliferation and migration [6, 7]. In addition, a previous study claimed that TUBA1C was statistically linked to the expression of RP11-480I12.5 in breast cancer (BRCA) and had prognostic value . TUBA1C has also been shown to promote aerobic glycolysis and cell growth through upregulation of YAP expression, thereby playing a role in the development of BRCA . Moreover, TUBA1C has been reported to be associated with glioma [10, 11]. These studies indicate that the TUBA1C gene has potential as a biomarker for tumour prognosis and immunotherapy outcomes. However, there has been no pancancer analysis of TUBA1C to date.
There is a complex link between malignant cancers and the tumour microenvironment (TME), which includes tumour cells, immune cells, stromal cells, endothelial cells, and cancer-associated fibroblasts . Immunotherapy focuses on recognizing and attacking cancer cells through immune cells within and outside the TME. ICIs, mainly including Programmed cell death protein 1 (PD-1)/PD-1 ligand (PD-L1) and Cytotoxic T lymphocyte antigen-4 (CTLA-4) inhibitors, act primarily by preventing the inhibition of interactions between T cells and other cells or tissues, thus allowing uninhibited activation of T cells, resulting in an antitumor effect [13, 14]. Ipilimumab, an anti-CTLA-4 antibody, was the first ICI cleared for therapeutic use by the Food and Drug Administration (FDA) in 2011 . The PD-1 blocking antibody pembrolizumab was approved by the FDA in 2014 for the treatment of patients with refractory myeloma . The use of chemotherapy in combination with ICI also holds good therapeutic promise. However, research into immunotherapy for various cancers is still in its initial stages, and further research of more general or specific immune targets is still needed to improve patient prognosis. Both tumour mutational burden (TMB) and the microsatellite instability (MSI) status can reflect the immune response to predict the outcome of cancers [17, 18].
In the current study, we performed a systematic pancancer analysis of TUBA1C through databases such as Oncomine, CCLE, HPA, and TGCA to analyse its effects on cancer prognosis, clinicopathology, the immune response, and the tumour microenvironment and validated the findings by immunohistochemistry. The potential of TUBA1C as a new immunotherapeutic target for tumour therapy was revealed.
Data download and differential analysis of TUBA1C expression in cancers
Thirty-three cancer-related RNA sequencing datasets and the clinicopathological and survival data for the corresponding patients were downloaded from the UCSC Xena website (https://xena.ucsc.edu/, derived from TCGA). Next, we analysed the mRNA expression of TUBA1C in 33 human malignancies using the online cancer microarray database Oncomine (https://www.oncomine.org/). We set the filter as we reported before . Next, we used Perl software to extract and integrate the expression information of TUBA1C in 33 cancers in TCGA (https://tcga.xenahubs.net) using the “wilcox.test” function. The R package “ggpubr” was applied to draw a box plot. In addition, the Tumor Immune Estimation Resource (TIMER) database (https://cistrome.shinyapps.io/timer/)  and Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/) database  were used for further profiling of TUBA1C expression in cancers. Moreover, we collected immunohistochemical data for cancer tissues using the HPA database (https://www.proteinatlas.org/)  and compared them with those for normal tissues to confirm the differential expression of TUBA1C in cancers at the protein level. In addition, the CGGA database (http://www.cgga.org.cn/)  was used to further investigate the expression levels and prognostic value of TUBA1C in gliomas.
Correlation analysis of TUBA1C expression with clinicopathological features and survival in human cancers
The survival information obtained from the TCGA database corresponding to each sample was used to further analyse the relationships between TUBA1C expression and clinical outcomes, including overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI), and progression-free interval (PFI). The results of survival analysis were visualized using forest plots and Kaplan–Meier (KM) curves.
Relationship between TUBA1C expression and immunity
To analyse the association between TUBA1C and TMB and MSI, we first collated the information from the TCGA database using Perl software. Next, the command “cor.test”, based on Spearman’s method, was used, and the R package “fmsb” was utilized to plot radar plots. We then calculated the immune and stromal fractions using the ESTIMATE algorithm. Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm was applied to assess tumour purity and stromal/immune cell infiltration in tumour tissues (n = 33) according to expression file . Subsequently, we evaluated the relationship between TUBA1C and TME or immune cell infiltration by applying the R packages “ggplot2”, “ggExtra” and ggpubr” (cut-off value p < 0.001). We further explored the associations between TUBA1C expression and immune cell infiltration with CIBERSORT, EPIC, TIMER, quanTIseq, xCELL and MCPcounter.
Next, CD274, CTLA4, HAVCR2, PDCD1, LAG3, SIGLEC15, PDCD1LG2, and TIGIT were chosen as immune checkpoint-related genes, and the expression data of these 8 transcripts were extracted. We used the R packages “ggplot2”, “immuneeconv” and “pheatmap” to assess immune checkpoint expression and the coexpression of TUBA1C with these immune checkpoints. The response to ICI therapy was predicted using the TIDE algorithm .
We downloaded the gene lists for immune activation, immunosuppression, chemokine receptor proteins, chemokines and Major histocompatibility complex (MHC) molecules from GSEA and processed the data in R. The results are shown in heatmaps.
Biological significance of TUBA1C in tumours
We obtained the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) gene sets  from GSEA (http://www.gsea-msigdb.org/gsea/downloads.jsp) to investigate the functional pathways of TUBA1C in tumours and then used the R packages “limma,” “org.Hs.eg.db,” “enrichplot,” and “clusterProfiler” to perform of GO functional and KEGG pathway annotation as well as pathway enrichment analysis of TUBA1C.
Immunohistochemical analysis of glioma
We collected specimens for immunohistochemistry from low-grade glioma (LGG), glioblastoma (GBM) and normal tissues from Renmin Hospital of Wuhan University. After deparaffinization, hydration and antigen repair, endogenous peroxidase activity in the sections was inhibited with 3% hydrogen peroxide for 15 min and the sections were then blocked with BSA, incubated with an anti-TUBA1C primary antibody (ab222849) overnight at 4 °C, and incubated with a secondary antibody for 50 min. Colour was developed with DAB chromogen, and the sections were restained with haematoxylin and finally subjected to dehydration and sealing. Finally, images were acquired by microscopy. We used ImageJ to quantify the immunohistochemical results.
The mRNA expression data of all genes in cancers were normalized using log2 transformation. Differences in data for each group were obtained using t tests or ANOVA. Significance was analysed with 0.05 as the threshold. Survival analysis was performed using Cox regression models, the KM method and the log-rank test. Moreover, Spearman’s or Pearson’s tests were used to identify relationships between two genes. All statistical procedures were completed in R software.
TUBA1C expression levels across cancers
We used Oncomine to assess TUBA1C gene expression in 33 cancers. It was shown that TUBA1C expression was elevated in most tumours, including bladder cancer (BLCA), breast cancer (BRCA), colorectal cancer (COAD), gastric cancer, head and neck cancer (HNSC), kidney cancer, liver cancer (LIHC), lymphoma, melanoma, myeloma, ovarian cancer (OV), pancreatic cancer (PAAD), and sarcoma (SARC). In contrast, lower TUBA1C expression was also found in BRCA and cervical cancer (CESC) (Fig. 1A). We further explored the differential expression of TUBA1C in cancers in GEPIA. Similar to the results in Oncomine, TUBA1C was upregulated in most cancers, including BLCA, BRCA, CESC, COAD, large B-cell lymphoma (DLBC), GBM, LGG, LIHC, LUAD, lung squamous cell carcinoma (LUSC), OV, PAAD, rectal cancer (READ), stomach cancer (STAD), testicular cancer (TGCT), thymoma (THYM), endometrioid cancer (UCEC) and uterine carcinosarcoma (UCS). Intriguingly, lower expression of TUBA1C was noted in acute myeloid leukaemia (LAML) (Fig. 1B).
We then analysed the data from TCGA with R software. A total of 11,057 mRNA expression profiles were obtained for 33 cancer types comprising 730 normal samples and 10,327 tumour samples. We found that TUBA1C showed higher expression in tumour samples in BLCA, BRCA, CESC, bile duct cancer (CHOL), COAD, oesophageal cancer (ESCA), GBM, HNSC, kidney chromophobe (KICH), kidney clear cell carcinoma (KIRC), kidney papillary cell carcinoma (KIRP), LIHC, LUAD, LUSC, prostate cancer (PRAD), READ, STAD and UCEC. No low expression of TUBA1C was detected in 33 cancers (Fig. 1C). Subsequently, TUBA1C was found to be upregulated in BLCA, BRCA, CESC, CHOL, COAD, ESCA, GBM, KIRC, KIRP, LIHC, LUAD, LUSC, PRAD, READ, STAD, THCA and UCEC by analysis with TIMER. In addition, the TUBA1C expression level was higher in HNSC-HPV+ tumours than in HNSC-HPV- tumours, and it was higher in SKCM tumours than in SKCM metastatic tumours (Fig. 1D). In conclusion, TUBA1C functions as a tumour promoter in most cancer types.
Immediately thereafter, we investigated the differential expression of TUBA1C at the protein level in cancerous and normal tissues in depth using Human Protein Atlas (HPA, https://www.proteinatlas.org/). The IHC results of BRCA, UCEC, renal cancer, LUAD, CESC, GBM, LIHC, PRAD and normal tissues were showed in Fig. 2A-H and these results were consisted with the expression data of TUBA1C mRNA in TCGA.
Prognostic value of TUBA1C across cancers
Subsequently, we explored the association between TUBA1C expression and the prognosis of patients across cancers. Cox regression analysis showed a significant relationship between TUBA1C expression and OS in BRAC (P = 0.031, HR = 1.288), GBM (P = 0.018, HR = 1.287), KICH (P = 0.002, HR = 7.366), KIRC (P = 0.004, HR = 1.389), KIRP (P = 0.01, HR = 2.012), LAML (P = 0.01, HR = 1.801), LGG (P < 0.001, HR = 2.373), LIHC (P < 0.001, HR = 1.675), LUAD (P = 0.002, HR = 1.373), MESO (P = 0.001, HR = 1.854), PAAD (P < 0.001, HR = 1.779), READ (P = 0.026, HR = 0.477) and SKCM (P < 0.001, HR = 1.418) (Fig. 3A). Moreover, our results suggested that high TUBA1C expression was a high-risk indicator for KIRC, KIRP, LAML, LGG, LIHC, LUAD, MESO, PAAD, SARC, and SKCM—especially KICH (HR = 7.366)—but was not a low-risk indicator for any cancer. In addition, a negative relationship between TUBA1C expression and OS was observed in 8 cancers (Fig. 3B-I): BRCA (p = 0.011), KIRC (p = 0.016), LAML (p = 0.037), LGG (p < 0.001), LIHC (p < 0.001), LUAD (p = 0.003) and SKCM (p = 0.002).
Since factors unrelated to the tumours may also contribute to death during follow-up, the relationship between TUBA1C expression and DSS in patients was subsequently analysed. It was revealed that TUBA1C expression impacted DSS in GBM (p = 0.032, HR = 1.279), KICH (p = 0.002, HR = 5.227), KIRC (p < 0.001, HR = 1.772), KIRP (p = 0.010, HR = 2.191), LGG (p < 0.001, HR = 2.498), LIHC (p < 0.001, HR = 1.538), LUAD (p = 0.034, HR =1.279), MESO (p = 0.007, HR = 1.899), PAAD (p < 0.001, HR =1.872), PRAD (p = 0.019, HR = 0.320) and SKCM (p = 0.009, HR = 1.296) (Fig. 4A). In these analyses, KICH displayed the highest HR (HR = 5.227). In addition, KM curves showed that high expression of TUBA1C was correlated with poor DSS in patients with BRCA (p = 0.044), KIRC (p < 0.001), LGG (p < 0.001), LIHC (p = 0.002), MESO (p < 0.001), PAAD (p = 0.026), SKCM (p = 0.029) and UCEC (p = 0.040) (Fig. 4B-I).
In addition, we explored the association between TUBA1C expression and the DFI in patients with different cancer types and found that increased TUBA1C expression was correlated with a decreased DFI in LGG (p < 0.001, HR = 3.229), LUAD (p = 0.007, HR = 1.445), PAAD (p = 0.005, HR = 2.281) and SARC (p = 0.026, HR = 1.314) (Fig. 5A). Moreover, in LUAD (p = 0.031), PAAD (p = 0.030) and SARC (p = 0.025), a high level of TUBA1C predicted a shorter DFI and poor prognosis in patients (Fig. 5B-D). Subsequently, we evaluated the association between TUBA1C expression and PFI. The forest plots showed that a high level of TUBA1C was correlated with a shorter PFI in COAD (p = 0.041, HR = 0.780), GBM (p = 0.007, HR = 1.893), KICH (p < 0.001, HR = 6.417), KIRC (p < 0.001, HR = 1.447), LGG (p < 0.001, HR = 1.893), LIHC (p = 0.004, HR = 1.266), LUAD (p = 0.024, HR = 1.215), MESO (p = 0.029, HR = 1.524), PAAD (p < 0.001, HR = 1.758), SARC (p = 0.004, HR = 1.288) and UCEC (p = 0.035, HR = 1.238). Among these cancers, KICH showed the highest HR (6.417, Fig. 6A). The results also indicated that high expression of TUBA1C affected the PFI unfavourably in ACC (p = 0.027) (Fig. 6B), GBM (p = 0.022) (Fig. 6C), KICH (p = 0.026) (Fig. 6D), KIRC (p = 0.005) (Fig. 6E), LGG (p < 0.001) (Fig. 6F), LUAD (p = 0.019) (Fig. 6G), MESO (p = 0.010) (Fig. 6H), PAAD (p = 0.006) (Fig. 6I) and SARC (p = 0.022) (Fig. 6K) but favourably in READ (p = 0.024) (Fig. 6J).
Relationship between TUBA1C expression and clinicopathological stage in cancers
We next assessed the relationship between TUBA1C expression and tumour stage. We observed that in BRCA (Fig. 7A, p = 0.033), HNSC (Fig. 7C, p = 0.0038), KICH (Fig. 7D, p = 0.00022), KIRC (Fig. 7E, p = 0.00025), LUAD (Fig. 7G, p = 0.034), MESO (Fig. 7H, p = 0.02) and READ (Fig. 7J, p = 0.025), significant differences in TUBA1C expression existed between stage I and stage IV tumours. In LIHC (Fig. 7F, p = 0.0067), LUAD (Fig. 7G, p = 0.0027) and MESO (Fig. 7H, p = 0.022), the expression of TUBA1C was higher in stage III than in stage I tumours. In ESCA, the TUBA1C expression in stage IV was higher than in stage III (Fig. 7B, p = 0.029). Intriguingly, in PAAD (Fig. 7I, p = 0.022), TUBA1C expression was higher in stage II tumours than in stage I tumours, but no dramatic difference existed in stage III or IV tumours. In conclusion, TUBA1C expression was correlated with tumour stage.
Relationship between TUBA1C expression and immunity
We first examined the association between TUBA1C expression level and TMB and MSI, two indicators that predict the sensitivity of cancer patients to ICIs. A positive association between TUBA1C expression and TMB was observed in 15 cancers, namely, UCS, UCEC, STAD, SKCM, SARC, PRAD, PAAD, LUAD, LGG, KIRC, KICH, COAD, CESC, BRCA and BLCA. In contrast, in THYM, a negative correlation between TUBA1C expression and TMB was observed (Fig. 8A). In ACC, UCEC, SARC and COAD, the TUBA1C expression level was positively associated with MSI, while in OV, LUSC, LUAD and LGG, a negative association existed between these factors (Fig. 8B).
Tumour development is associated not only with abnormal mutations in tumour cells but also with the composition of their microenvironment and stromal cell proportions or activation states . Therefore, the association between TUBA1C expression and the TME was investigated by assessing the relationships of stromal and immune scores to TUBA1C expression in 33 cancers using the ESTIMATE algorithm. It was shown that TUBA1C expression was positively correlated with stromal scores in GBM and LGG, while a negative correlation was observed in ESCA, STAD and TGCT (Fig. 9A). In addition, TUBA1C expression was positively related to immune scores in GBM, LGG, PCPG and THCA and negatively related to immune scores in ESCA (Fig. 9B).
TIICs, an important component of the TME, are closely associated with the development, progression and metastasis of cancer. Herein, we examined the underlying relationship between the level of infiltration of 22 immune cell types and TUBA1C expression in different cancer types. We selected the TIICs with the highest correlation with TUBA1C expression in each cancer for demonstration (Fig. 10). Our results showed that TUBA1C expression was negatively related to the infiltration of memory CD4+ T cells in COAD (Fig. 10D), STAD (Fig. 10P) and TGCT (Fig. 10Q) but positively associated in BRCA (Fig. 10B). In ESCA (Fig. 10E), TUBA1C expression was observed to be negatively correlated with infiltration of regulatory T cells. In addition, TUBA1C expression was positively related to infiltration of neutrophils in BLCA (Fig. 10A), HNSC (Fig. 10G), KIRC (Fig. 10H) and SKCM (Fig. 10O) and positively correlated with infiltration of dendritic cells in THCA (Fig. 10R). In addition, the TUBA1C expression level was correlated with infiltration of several diverse subsets of macrophages in the TMB. TUBA1C expression was positively associated with M1 macrophage infiltration in CESC (Fig. 10C), KIRP (Fig. 10I), LGG (Fig. 10J) and UCEC (Fig. 10T), with M2 macrophage infiltration in LUSC (Fig. 10M) and with M0 macrophage infiltration in SARC (Fig. 10N). In addition, TUBA1C expression was negatively related to infiltration of mast cells in LUAD (Fig. 10L) and THYM (Fig. 10S). Moreover, TUBA1C expression was positively related to infiltration of gamma delta T cells in GBM (Fig. 10F) and LIHC (Fig. 10K). We performed further analyses with other algorithms, including CIBERSORT (Fig. S1), EPIC (Fig. S2), MCPcounter (Fig. S3), quanTIseq (Fig. S4), TIMER (Fig. S5) and xCELL (Fig. S6), and also found correlations between TUBA1C and immune cell infiltration.
Next, we investigated the correlations between TUBA1C expression and the expression of immune checkpoints, including CD274, HAVCR2, CTLA4, LAG3, PDCD1LG2, PDCD1, SIGLEC15, and TIGIT, in cancers. We found that TUBA1C was correlated with immune checkpoints in all cancers except UVM, KIRC and ACC. In LUSC and ESCA, TUBA1C expression was negatively related to the expression of most immune checkpoints, while in OV, LGG, BRCA and BLCA, a positive association of TUBA1C expression with immune checkpoints was revealed (Fig. 11A), suggesting that TUBA1C has promise as a predictive factor for immunotherapeutic response in patients with these cancers. Additionally, the group with high TUBA1C responded better to ICI therapy than the group with low TUBA1C in LGG (Fig. 11B, p < 0.001), LIHC (Fig. 11C, p < 0.001), LUAD (Fig. 11D, p < 0.001) and LUSC (Fig. 11E, p = 0.011), suggesting that TUBA1C has the potential to be an immunotherapy target.
Relationship between TUBA1C expression and immune-associated genes and pathways in diverse cancers
Further gene coexpression analysis was carried out between TUBA1C and immune-associated genes, including those encoding MHC molecules, chemokines, chemokine receptors, and proteins related to immune activation and immunosuppression, in 33 cancers. The heatmaps revealed that most immune-associated genes except CCL27 were coexpressed with TUBA1C and that the main immune-associated genes were positively related to TUBA1C in LGG, LIHC, PCPG and THCA (Fig. 12). TUBA1C was coexpressed with MHC genes in several cancer types, particularly BLCA, ESCA, LGG, LICH and THCA (Fig. 12A). We also found that TUBA1C was coexpressed with immune activation genes in all cancer types and with immunosuppressive genes in most cancers except uveal melanoma (UVM) (Fig. 12B, E). TUBA1C was coexpressed with most chemokines except CCL27, and coexpression of TUBA1C with chemokines and chemokine receptors was observed in most cancers except UVM (C, D).
To explore the underlying mechanism of action of TUBA1C, we conducted GO term and KEGG pathway enrichment analyses of TUBA1C (Fig. 13). Our results suggested that TUBA1C participated in several immune-associated pathways in the GO database, such as regulation of lymphocyte activation, positive regulation of cytokine production, and T-cell activation, in PCPG. In addition, it influenced immunoglobulin in KICH and MESO (Fig. 13A). In ESCO, TUBA1C was associated with lymphocyte differentiation in GO analysis and with chemokines in KEGG analysis. In LGG, there was a tight link between TUBA1C and the GO term regulation of immune receptor processes and between TUBA1C and the KEGG pathways chemokine signaling pathway and T-cell receptor signaling pathway (Fig. 13A, B). KEGG pathway analysis also indicated that TUBA1C mediated RIG-I-like receptor signalling in OV and mediated the regulation of autophagy, antigen processing and presentation and RIG-I-like receptor signalling in LUSC (Fig. 13B). Furthermore, TUBA1C was implicated in many biological pathways, including miRNA binding and gene silencing by RNA, and in many other pathways, such as allograft rejection, drug metabolism cytochrome p40, retinol metabolism, and neuroactive ligand receptor interaction.
Expression level and prognostic value of TUBA1C in gliomas
We further investigated the expression of TUBA1C in gliomas and its impact on the prognosis of patients with gliomas. Figure 14A shows TUBA1C expression in gliomas of different histologies in the CGGA dataset. TUBA1C expression was highest in primary glioblastoma (pGBM) and recurrent pGBM (rGBM) and lowest in astrocytoma (A). In anaplastic oligodendrocytoma (AO) and anaplastic astrocytoma (AA), TUBA1C was expressed at moderate levels. Overall, TUBA1C expression increased in gliomas as WHO classification increased. In addition, TUBA1C expression was elevated in relapsed A, O, AO, and AA, similar to the results shown in Fig. 14I and Fig. 14J. Figure 14B further confirmed this phenomenon. We further found that TUBA1C expression was lower in the IDH mutant state than in the wild type (Fig. 14C), and similar results were found in WHO II, III and IV (Fig. 14D). In addition, TUBA1C expression was lower in patients with IDH mutation than in those with wild-type IDH (Fig. 14E), and similar results were found for WHO grades III and IV (Fig. 14F). In IDH wildtype LGG, the TUBA1C expression was higher than in LGG with IDH mutation and 1p/19q codeletion (Fig. 14G). Moreover, TUBA1C expression was higher in patients older than 42 years than in those aged < 42 years (Fig. 14H). Subsequently, we collected normal, LGG and GBM tissues for immunohistochemical analysis. The results showed that the IHC staining intensity was strongest in GBM (Fig. 14K, right), moderate in LGG (Fig. 14K, middle) and weakest in normal tissues (Fig. 14K, left). In primary gliomas, the probability of survival was higher in the TUBA1C low expression group than in the TUBA1C high expression group in all WHO grades (Fig. 14L, p < 0.0001) and in the WHO grade II (Fig. 14N, p = 0.032), WHO grade III (Fig. 14P, p = 0.032) and WHO grade IV (Fig. 14R, p = 0.017) subclassifications. In recurrent glioma, the probability of survival was higher in the TUBA1C low expression group in all WHO grades (Fig. 14M, p < 0.0001) and the WHO grade III (Fig. 14Q, p = 0.012) subclassification. In WHO II and IV recurrent glioma, the expression of TUBA1C was not related to the probability of survival (Fig. 14O, p = 0.98; Fig. 14S, p = 0.1).
Pancancer analysis can shed light on the expression, function and mechanism of genes in different cancers, providing new ideas and targets for improving the therapeutic efficacy in and prognosis of cancers. Microtubules, a main component of eukaryotic cells, play an important role in dynamic aggregation and depolymerization through cell division and replication , and TUBA1C is a subtype of α-tubulin. TUBA1C was reported to be related to cell proliferation and the cell cycle in diverse cancers  and was associated with the progression of several cancers , whereas no pancancer analysis has been performed for TUBA1C. Our research revealed that TUBA1C was significantly upregulated in 18 cancers, and immunohistochemistry confirmed consistent results for the protein expression levels of TUBA1C in cancers. Previous studies have shown that TUBA1C can predict poor outcome in PDAC  by modulating cell cycle signalling pathways to induce apoptosis, in LUAD  and LGG  by affecting TIICs in the tumour microenvironment, and in HCC  by cell cycle signalling pathways. Moreover, it was found that the combination of GTSE1 and TUBA1C predicted a 100% probability of triple-negative breast cancer (TNBC) in whites; NRF1, TUBA1B and BAX with EFNA4, and NRF1 and BTRC predicted a 100% probability of TNBC in blacks, confirming the potential of TUBA1C as a diagnostic marker . In our research, an association of TUBA1C expression and poor prognosis (OS, DSS, DFI and PFI) of patients was observed, particularly in LGG, LIHC, MESO and KIRC. Notably, TUBA1C was highly expressed in KIRC and is related to OS in KIRC patients, suggesting that TUBA1C has promise as a diagnostic marker for KIRC. Moreover, we found that TUBA1C was correlated with the clinicopathological stage of cancers, which helps to determine the degree of malignancy of tumours. In BRCA, HNSC, KICH, KIRC, LUAD, MESO and READ, TUBA1C was expressed at higher levels in stage IV than in stage I tumours.
In recent years, ICIs have been a hot topic in the treatment of patients with advanced cancer , but predictive biomarkers are needed. The results of studies in certain cancers indicate that TMB can predict the clinical response to ICI therapy . It has also been confirmed that patients with MSI-high colorectal cancers benefit significantly from ICI treatment compared to patients with MSI-low status . Therefore, we investigated the association between TUBA1C expression and TMB and MSI. It was shown that TUBA1C affected TMB in 16 cancers and affected MSI in 8 cancers, suggesting that the efficacy of ICI therapy can be predicted by TUBA1C expression. We can assume that in cancers where TUBA1C expression is positively correlated with TMB or MSI, a high TUBA1C level predicts a good outcome and prognosis for immunotherapy.
Tumours are strongly influenced by the surrounding normal tissue, forming a special ecological niche called the TME . The TME consists of stromal cells, fibroblasts, endothelial extracellular, innate and adaptive immune cells. Cytokines within the TME manipulate immune function, causing malfunction of the immune response and ultimately leading to tumour progression . In addition, the TME is widely accepted to be a major player in the response to immunotherapy . Therefore, it is essential to explore the relationship between TUBA1C and the TME. A higher immune score or stromal score denotes more immune or matrix components in the TME . Our results showed that TUBA1C expression was significantly negatively related to the immune components of the TME in ESCA and positively related to the immune components of GBM, LGG, PCPG and THCA. In addition, TUBA1C expression was positively related to the stromal components of the TME in GBM and LGG and negatively related to the stromal components in ESCA, STAD and TGCT.
Moreover, we investigated the association between TUBA1C expression and immune cell infiltration in the TMB. Memory T cells are involved in specific and acquired immunity and play a fundamental role in the immune system . TUBA1C expression was negatively related to memory CD4 T-cell infiltration in COAD, STAD and TGCT but positively associated in BRCA in our study. Regulatory T cells are specialized to suppress excessive immune activation and maintain immune homeostasis . TUBA1C expression was negatively correlated with regulatory T cell infiltration in ESCA. Tumours manipulate neutrophils early in their differentiation process to change the state of the tumour by creating a different phenotype and functional polarization state . TUBA1C expression was positively related to neutrophil infiltration in BLCA, HNSC, KIRC and SKCM. Macrophage infiltration in solid tumours is linked to poor prognosis and is correlated with chemoresistance in most cancers . TUBA1C expression was positively related to M1 macrophage infiltration in CESC, KIRP, LGG and UCEC, with M2 macrophage infiltration in LUSC and with M0 macrophage infiltration in SARC. Mast cells interact with the immune and nonimmune components of the TMB to regulate tissue homeostasis and immune responses and act as antitumor agents . TUBA1C expression was negatively correlated with mast cell infiltration in LUAD and THYM. Dendritic cells are natural mediators of antigen delivery, control immune tolerance and immune responses and are an important target for generating immunity against cancer . We observed a positive correlation between TUBA1C and dendritic cell infiltration in THCA.
In addition, TUBA1C was coexpressed with immune-related genes encoding chemokine proteins, immune activators, immunosuppressive factors, chemokines and MHC. Moreover, we performed mechanistic analysis of TUBA1C using the GO and KEGG pathway analyses and found that TUBA1C was associated with several immune-associated pathways, such as regulation of lymphocyte activation, positive regulation of cytokine production, T cell activation, immunoglobulin, lymphocyte differentiation, immune receptor process, RIG-I-like receptor signaling and chemokine signaling. This suggests that TUBA1C may influence cancer progression and prognosis by affecting the tumour immune response.
Moreover, we validated the expression and role of TUBA1C in gliomas. Our study showed that TUBA1C was highly expressed in glioma types with high malignancy and in recurrent gliomas. In addition, TUBA1C expression was highest in WHO IV gliomas, intermediate in WHO III gliomas and lowest in WHO II gliomas. In addition, TUBA1C expression was reduced in patients with both IDH mutation and 1p/19q codeletion. Moreover, TUBA1C overexpression was observed in older patients (age ≥ 42) and predicted a poor outcome. Immunohistochemical analysis also confirmed that TUBA1C expression was higher in high-grade gliomas than in low-grade gliomas.
In summary, the current study shows that upregulated TUBA1C is related to poor prognosis in cancer patients, is associated with a higher tumour grade, affects tumour sensitivity to and the efficacy of ICI therapy, influences immune cell infiltration in the TME, and interacts with immune-associated genes and pathways in cancers. Therefore, TUBA1C is a potential prognostic marker and may play an important role in immunotherapy.
Availability of data and materials
The datasets supporting the conclusions of this article are available on the UCSC Xena (https://xena.ucsc.edu/), Oncomine (https://www.oncomine.org/), TCGA (https://tcga.xenahubs.net), the HPA database (https://www.proteinatlas.org/) and the CGGA database (http://www.cgga.org.cn/) websites. The datasets generated and/or analyzed during the current study are also available from the corresponding authors upon reasonable request in compliance with ethical standards.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Li S, Shan X, Wang Y, Chen Q, Sun J, He Z, et al. Dimeric prodrug-based nanomedicines for cancer therapy. J Control Release. 2020;326:510–22.
Kennedy LB, Salama AKS. A review of cancer immunotherapy toxicity. CA Cancer J Clin. 2020;70(2):86–104.
Akhmanova A, Maiato H. Closing the tubulin detyrosination cycle. Science. 2017;358(6369):1381–2.
Bian T, Zheng M, Jiang D, Liu J, Sun H, Li X, et al. Prognostic biomarker TUBA1C is correlated to immune cell infiltration in the tumor microenvironment of lung adenocarcinoma. Cancer Cell Int. 2021;21(1):144.
Albahde MAH, Zhang P, Zhang Q, Li G, Wang W. Upregulated expression of TUBA1C predicts poor prognosis and promotes Oncogenesis in pancreatic ductal adenocarcinoma via regulating the cell cycle. Front Oncol. 2020;10:49.
Wang J, Chen W, Wei W, Lou J. Oncogene TUBA1C promotes migration and proliferation in hepatocellular carcinoma and predicts a poor prognosis. Oncotarget. 2017;8(56):96215–24.
Lou W, Ding B, Zhong G, Yao J, Fan W, Fu P. RP11-480I12.5-004 promotes growth and tumorigenesis of breast cancer by relieving miR-29c-3p-mediated AKT3 and CDK6 degradation. Mol Ther Nucleic Acids. 2020;21:916–31.
Wu Z, Sun S, Fan R, Wang Z. Tubulin alpha 1c promotes aerobic glycolysis and cell growth through upregulation of yes association protein expression in breast cancer. Anticancer Drugs. 2022;33(2):132–41.
Elsheikh SSM, Bakas S, Mulder NJ, Chimusa ER, Davatzikos C, Crimi A. Multi-stage association analysis of Glioblastoma gene expressions with texture and spatial patterns. Brainlesion. 2019;11383:239–50.
Zhu H, Hu X, Gu L, Jian Z, Li L, Hu S, et al. TUBA1C is a prognostic marker in low-grade Glioma and correlates with immune cell infiltration in the tumor microenvironment. Front Genet. 2021;12:759953.
Hu X, Zhu H, Zhang X, He X, Xu X. Comprehensive analysis of pan-cancer reveals potential of ASF1B as a prognostic and immunological biomarker. Cancer Med. 2021;10(19):6897–916.
Jeurling S, Cappelli LC. Treatment of immune checkpoint inhibitor-induced inflammatory arthritis. Curr Opin Rheumatol. 2020;32(3):315–20.
Zhu H, Hu X, Feng S, Jian Z, Xu X, Gu L, et al. The hypoxia-related gene COL5A1 is a prognostic and immunological biomarker for multiple human tumors. Oxidative Med Cell Longev. 2022;2022:6419695.
Ledford H. Melanoma drug wins US approval. Nature. 2011;471(7340):561.
Ribas A, Wolchok JD. Cancer immunotherapy using checkpoint blockade. Science. 2018;359(6382):1350–5.
Cristescu R, Mogg R, Ayers M, Albright A, Murphy E, Yearley J, et al. Pan-tumor genomic biomarkers for PD-1 checkpoint blockade-based immunotherapy. Science. 2018;362(6411):eaar3593.
Vilar E, Gruber SB. Microsatellite instability in colorectal cancer-the stable evidence. Nat Rev Clin Oncol. 2010;7(3):153–62.
Zhu H, Hu X, Ye Y, Jian Z, Zhong Y, Gu L, et al. Pan-cancer analysis of PIMREG as a biomarker for the prognostic and immunological role. Front Genet. 2021;12:687778.
Murray PJ. On macrophage diversity and inflammatory metabolic timers. Nat Rev Immunol. 2020;20(2):89–90.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–w102.
Thul PJ, Lindskog C. The human protein atlas: a spatial map of the human proteome. Protein Sci. 2018;27(1):233–44.
Xu J, Zhang Z, Qian M, Wang S, Qiu W, Chen Z, et al. Cullin-7 (CUL7) is overexpressed in glioma cells and promotes tumorigenesis via NF-κB activation. J Exp Clin Cancer Res. 2020;39(1):59.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Jiang P, Gu S, Pan D, Fu J, Sahu A, Hu X, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. 2018;24(10):1550–8.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–d51.
Quail DF, Joyce JA. Microenvironmental regulation of tumor progression and metastasis. Nat Med. 2013;19(11):1423–37.
Kim ND, Park ES, Kim YH, Moon SK, Lee SS, Ahn SK, et al. Structure-based virtual screening of novel tubulin inhibitors and their characterization as anti-mitotic agents. Bioorg Med Chem. 2010;18(19):7092–100.
Bates D, Eastman A. Microtubule destabilising agents: far more than just antimitotic anticancer drugs. Br J Clin Pharmacol. 2017;83(2):255–68.
Li H, Jiang X, Zhu S, Sui L. Identification of personalized dysregulated pathways in hepatocellular carcinoma. Pathol Res Pract. 2017;213(4):327–32.
Ramos J, Yoo C, Felty Q, Gong Z, Liuzzi JP, Poppiti R, et al. Sensitivity to differential NRF1 gene signatures contributes to breast cancer disparities. J Cancer Res Clin Oncol. 2020;146(11):2777–815.
Callahan MK, Postow MA, Wolchok JD. Targeting T cell co-receptors for cancer therapy. Immunity. 2016;44(5):1069–78.
Samstein RM, Lee CH, Shoushtari AN, Hellmann MD, Shen R, Janjigian YY, et al. Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nat Genet. 2019;51(2):202–6.
Lin A, Zhang J, Luo P. Crosstalk between the MSI status and tumor microenvironment in colorectal cancer. Front Immunol. 2020;11:2039.
Hessmann E, Buchholz SM, Demir IE, Singh SK, Gress TM, Ellenrieder V, et al. Microenvironmental determinants of pancreatic cancer. Physiol Rev. 2020;100(4):1707–51.
Hinshaw DC, Shevde LA. The tumor microenvironment innately modulates cancer progression. Cancer Res. 2019;79(18):4557–66.
Giannone G, Ghisoni E, Genta S, Scotto G, Tuninetti V, Turinetto M, et al. Immuno-metabolism and microenvironment in cancer: key players for immunotherapy. Int J Mol Sci. 2020;21(12):4414.
Bi KW, Wei XG, Qin XX, Li B. BTK has potential to be a prognostic factor for lung adenocarcinoma and an indicator for tumor microenvironment remodeling: a study based on TCGA data mining. Front Oncol. 2020;10:424.
Vahidi Y, Faghih Z, Talei AR, Doroudchi M, Ghaderi A. Memory CD4(+) T cell subsets in tumor draining lymph nodes of breast cancer patients: a focus on T stem cell memory cells. Cell Oncol (Dordr). 2018;41(1):1–11.
Ferreira LMR, Muller YD, Bluestone JA, Tang Q. Next-generation regulatory T cell therapy. Nat Rev Drug Discov. 2019;18(10):749–69.
Coffelt SB, Wellenstein MD, de Visser KE. Neutrophils in cancer: neutral no more. Nat Rev Cancer. 2016;16(7):431–46.
Cassetta L, Pollard JW. Targeting macrophages: therapeutic approaches in cancer. Nat Rev Drug Discov. 2018;17(12):887–904.
Rigoni A, Colombo MP, Pucillo C. Mast cells, basophils and eosinophils: from allergy to cancer. Semin Immunol. 2018;35:29–34.
Palucka K, Banchereau J. Cancer immunotherapy via dendritic cells. Nat Rev Cancer. 2012;12(4):265–77.
This study was sponsored by the National Natural Science Foundation of China (31971166).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of Renmin Hospital of Wuhan University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The correlation between TUBA1C and immune cell infiltration in CIBERSOR. Figure S2. The correlation between TUBA1C and immune cell infiltration in EPIC. Figure S3. The correlation between TUBA1C and immune cell infiltration in MCPCOUNTER. Figure S4. The correlation between TUBA1C and immune cell infiltration in QUANTISEQ. Figure S5. The correlation between TUBA1C and immune cell infiltration in TIMER. Figure S6. The correlation between TUBA1C and immune cell infiltration in XCELL.
About this article
Cite this article
Hu, X., Zhu, H., Chen, B. et al. The oncogenic role of tubulin alpha-1c chain in human tumours. BMC Cancer 22, 498 (2022). https://doi.org/10.1186/s12885-022-09595-0
- Tumour immunity
- Immune infiltration