Abstract

Bladder cancer (BC) is a common malignant tumor of the genitourinary system, and there are not enough tumor biomarker tests that are specific, trustworthy, and noninvasive for the diagnosis and prognosis. The purpose of this study is to investigate the clinical relevance, prognostic value, and immunological signature of Mannosidase alpha class 1B member 1(MAN1B1) expressions in BC. The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases provided the raw information that was used to analyze the expression of MAN1B1 in tumor patients. Then, a statistical study was carried out to assess the correlations of MAN1B1 expression with the clinical characteristics and the prognosis of BC. The correlation between MAN1B1 expression and tumor immune infiltration was explored via single-sample gene set enrichment analysis (ssGSEA). In human cancers, MAN1B1 expressions were shown to be generally higher in tumors than in normal specimens. We confirmed that MAN1B1 expression was distinctly increased in BC specimens compared with nontumor specimens. BC specimens with advanced T stage and M stage showed a higher level of MAN1B1. Survival analysis revealed that the overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI) of patients with high MAN1B1 expressions were distinctly worse than those with low MAN1B1 expressions. Importantly, multivariate analyses only confirmed that MAN1B1 expression was an independent prognostic factor for OS of the patients with BC. Furthermore, we observed that MAN1B1 expression level was significantly correlated with abundance of multiple immune infiltrates including Th2 cells, macrophages, Th1 cells, neutrophils, T helper cells, and NK CD56 bright cells. In conjunction with all of these findings, elevated MAN1B1 expression is associated with a poor prognosis and a higher number of immune cells in BC. MAN1B1 has the potential to act as a biomarker that can evaluate both the patient’s prognosis and the degree of immune infiltration in BC.

1. Introduction

Bladder cancer (BC) is the fourth most prevalent cancer in men and the most frequently diagnosed malignancy of the urinary system worldwide [1]. The risk of developing BC increases with age, and there has been a discernible rising trend in the overall incidence of BC over the past several years [2, 3]. Approximately 90 percent of all cases of cancer are classified as transitional cell carcinoma. The majority of bladder cancers are known as nonmuscle-invasive bladder cancers (NMIBCs), which frequently recur and progress into muscle-invasive bladder cancers [4, 5]. MIBCs account for around 30 percent of all cases of BC, and the treatment of choice for these cases is radical cystectomy in conjunction with pelvic lymph node dissection [6, 7]. In contrast to the majority of other malignancies, the diagnostic procedures, treatment options, and percentage of patients who are still alive after five years have not altered in the past three decades [8, 9]. With a recurrence rate of between 60 and 70 percent, BC continues to be a significant risk to human health all over the world. As a consequence of this, it is necessary to create methods of precise prediction in order to advance clinical diagnosis and therapy.

Mannosidase alpha class 1B member 1(MAN1B1), located on 9q34.3, encodes an enzyme belonging to the glycosyl hydrolase 47 family [10]. This enzyme has a role in N-glycan biosynthesis and belongs to the class I alpha-1,2-mannosidase family [11]. It converts Man9GlcNAc to Man8GlcNAc isomer B in a very particular manner. MAN1B1 encoded α-1,2-mannosidase α-Mannosidase is known to have an important role in both the modification of protein glycosylation and the hydrolysis of glycoprotein polysaccharides [12, 13]. It is composed of many different classes, including -Mannosidases I, -Mannosidases II, and unclassified -Mannosidase [14]. Evidence indicating alpha-mannosidase has a role in the progression of cancer has been accumulating steadily over the years [1517]. However, the potential function of MAN1B1 in tumors was rarely reported. In the current investigation, our goals were to investigate the predictive significance of MAN1B1 in BC and to investigate its connection with immune infiltration.

2. Materials and Methods

2.1. Raw Data Acquisition and Processing

The Cancer Genome Atlas (TCGA) research network has profiled and evaluated a massive collection of clinical data on > 10,000 cancer patients representing 33 distinct types of tumors. Downloading clinical data and RNA expression information from the TCGA and Genotype-Tissue Expression (GTEx) databases was accomplished through the usage of the UCSC Xena database (https://xenabrowser.net/datapages/). 33 cancer types were included: OV, PAAD, PRAD, READ, SKCM, STAD, TGCT, THCA, THYM, UCEC, UCS, ACC, BLCA, BRCA, COAD, DLBC, ESCA, GBM, HNSC, KICH, KIRC, KIRP, LAML, LGG, LIHC, LUAD, and LUSC.

2.2. Analysis of MAN1B1 Expression in Cancers

The TCGA and GTEx projects were the sources for the information regarding the abnormal expressions of MAN1B1 between cancer and normal specimens that were matched to the tumor. A tissue bank and data resource known as GTEx (https://gtexportal.org) was established by the National Institutes of Health (NIH) Common Fund. In total, 53 human normal specimens from approximately 1,000 individuals have been analyzed for RNA sequencing, genetic variation and molecular phenotypes. Plotting was done using log2 (TPM + 1) converted expression data, which was our choice for the parameter selection process.

2.3. Cox Regression Analysis and Survival Analysis

In order to investigate the connection between MAN1B1 expressions and patients’ overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI) in bladder urothelial carcinoma(BLCA) patients, a Cox regression analysis was carried out in the R environment. After using the best approach for separating patients into groups with high and low MAN1B1 expression, the Kaplan–Meier methods were applied to produce survival curves for patients. This was done after sorting patients into groups with high or low MAN1B1 expressions. An examination of the survival was carried out with the aid of survival receiver operating characteristics (ROC) and the survival package in R. The log-rank test was applied to analyze the differences between the curves, and a was considered statistically significant.

2.4. Immune Cell Infiltration

Using ssGSEA, the relative degrees of tumor infiltration by 24 different immune cell types were assessed. This allowed us to investigate the expression of genes found in published signature gene lists. The signatures that we used contained a total of 509 genes and included a wide variety of cell types involved in both adaptive and innate immune responses. The Wilcoxon rank sum test and Spearman correlation were utilized in order to investigate whether or not there is a connection between MAN1B1 and the infiltration levels of immune cells, as well as whether or not there is an association between the infiltration levels of immune cells and the various expression groups of MAN1B1.

2.5. Statistical Analysis

All statistical analyses were performed in the R package (V3.6.2). The Wilcoxon test was utilized in order to compare the differences between two groups of data that were not regularly distributed. The relationship between clinical pathologic features and MAN1B1 was analyzed with the Wilcoxon signed-rank sum test and logistic regression. ROC curves were established to evaluate the diagnostic value of MAN1B1 in BC patients. A was considered statistically significant.

3. Results

3.1. MAN1B1 Expression in Pan-Cancer

Based on TCGA datasets, we investigate the level of MAN1B1 expression in each of the many cancer types. Figure 1(a) showed that the levels of MAN1B1 expression in the tumor specimens of ACC, BLCA, BRCA, CESC, CHOL, COAD, DLBC, GBM, HNSC, KICH, KIRC, KIRP, LAML, LGG, LIHC, LUAD, LUSC, SARC, STAD, THYM, UCEC, and UCS are all higher than those in the nontumor tissues. However, the expression of MAN1B1 in the tumor tissues of SKCM and TGCT was significantly lower than that in the normal tissues (Figure 1(a)). The expression of MAN1B1 in BLCA is shown in Figure 1(b). Due to the considerable overexpression of MAN1B1 in BLCA patients, we were motivated to further investigate the diagnostic significance of this gene for BLCA patients. As shown in Figure 1(c), the results of ROC curves revealed that the AUC was 0.793 (0.6999 to 0.888, 95% CI) in a comparison between BLCA specimens and nontumor specimens. In addition, we also performed a paired t-test. A similar finding was also observed (Figures 2(a) and 2(b)). Our findings suggest that overexpression of MAN1B1 in malignancies may be a common occurrence.

3.2. Association between Clinicopathological Characteristics and MAN1B1 Expression in BLCA Patients

We evaluated the association between the expression of MAN1B1 and the clinicopathological aspects of BLCA in order to further define the relevance of MAN1B1 in BLCA. This was done in order to better understand the role that MAN1B1 plays in BLCA. On the basis of the median relative MAN1B1 expression value, all of the patients diagnosed with BLCA were separated into two groups. As shown in Table 1, we did not observe a distinct association between MAN1B1 expression and several clinical factors based on the results of the chi square test. In the T-test, we also did not observe a distinct association between MAN1B1 expression and gender and age (Figures 3(a) and 3(b)). However, we found that BLCA specimens with advanced T stage and M stage showed a higher level of MAN1B1 (Figures 3(c) and 3(d)). The expression of MAN1B1 did not change distinctly in specimens with different N stages (Figure 3(e)). Moreover, we found that the expression of MAN1B1 was distinctly increased in BLCA specimens with high grade or dead status (Figure 3(f) and 3(g)).

3.3. Associations between MAN1B1 Expression and Patient Survival

To explore the relationships between MAN1B1 expressions and survivals of BLCA patients, Kaplan–Meier methods were performed to analyze the differences in OS, DSS, and PFI. As shown in Figure 4, the OS, DSS, and PFI of patients with high MAN1B1 expressions were distinctly worse than those with low MAN1B1 expressions. In univariate analyses, we observed that MAN1B1 expression was associated with OS (Table 2), DSS (Table 3), and PFI (Table 4) of BLCA patients. However, multivariate analyses only confirmed that MAN1B1 expression (HR = 1.970, 95% CI 1.226–3.167, ) was an independent prognostic factor for OS of the patients with BLCA (Table 2). The potential prognostic values of MAN1B1 used as an independent prognostic factor for DSS and PFI of BLCA patients were not further confirmed (Tables 3 and 4).

3.4. The Correlation between MAN1B1 Expression and Immune Infiltration

The spearman correlation test was used to investigate the relationship between the expression of MAN1B1 and the amount of immune cell infiltration that was measured using ssGSEA. The abundance of acquired immunocytes, such as T helper cells and NK CD56bright cells, was found to have a negative correlation with the expression of MAN1B1, while the abundance of innate immunocytes, such as Th2 cells, Macrophages, Th1 cells, Neutrophils, NK CD56dim cells, and other such cells, had a positive correlation with the expression of MAN1B1 (Figure 5).

4. Discussion

A number of cancers, including carcinomas, share a phenomenon known as deregulation of tumor-related genes, which plays a significant part in the progression of cancer via a variety of intricate pathways [18, 19]. Alterations in the expression of a particular gene have been shown to have a strong correlation both with the development of human malignancies and their overall prognosis [20, 21]. Thus, to better understand how BC develops and progresses, finding mRNA molecular profiles related to a patient’s prognosis could reveal biological mechanisms at work. It could also lead to the discovery of new therapeutic targets for the disease.

MAN1B1 is a newly identified tumor-related gene. Previous study by Wang et al. reported that the expressions of MAN1B1 were distinctly higher in cancerous specimens than in nontumor samples. In addition, there was a correlation between increased MAN1B1 expression and a bad prognosis in patients with BC. In functional experiments, the suppression of MAN1B1 resulted in a reduction of BC cell proliferation, invasion, and migration [22]. To our best knowledge, this is the only research about the function of MAN1B1 in tumor. In this study, we first performed a pan-cancer analysis and observed that MAN1B1 expression was distinctly increased in most types of tumors, which was consistent with previous findings in BC. Our data revealed that MAN1B1 may serve as an oncogene in human tumors. Then, we further explored the relationships between MAN1B1 expressions and the prognosis of BC patients. We found that BC specimens with advanced stages showed an increased level. Interestingly, we found that the BC specimen with a dead event also showed an increased level of MAN1B1, suggesting that MAN1B1 may influence the clinical outcome of BC patients. The results of Kaplan–Meier methods confirmed that the OS, DSS, and PFI of patients with high MAN1B1 expression was significantly worse than those with low MAN1B1 expression. More importantly, multivariate analyses only confirmed that MAN1B1 expression were an independent prognostic factor for OS of the patients with BC. Our findings suggest MAN1B1 as a novel diagnostic and prognostic biomarker for BC patients.

Recent research studies have demonstrated that immune cells that infiltrate tumors, known as tumor infiltrating cells (TIICs), may be able to control the process of both the genesis and progression of tumors [23, 24]. In addition to this, TIICs have the capacity to undergo clonal expansion and are enriched preferentially in BC; a poor prognosis has been shown to correspond with an accumulation of TIICs in BC [25]. Then, our results demonstrated that MAN1B1 expression in BC was negatively associated with multiple types of immune cell infiltration. In previous studies, although it has been reported that the functions and prognostic relevance of numerous subtypes of TIICs in multiple forms of cancer are inconsistent, some of the results are still unclear. For example, infiltrating CD8+ T cells are generally considered to be tumor inhibitors that are associated with a positive prognosis in most types of cancer; however, in renal cell carcinoma and prostate cancer, CD8+ T cells are reported to be associated with poor clinical outcomes [2629]. This is because infiltrating CD8+ T cells are associated with an increased risk of death. In addition, Hald et al. observed that CD8+ T cells were a predictor of a favorable prognosis in non-small-cell lung cancer [30]; however, the findings of Tian et al., who found that CD8+ T cells serve as a predictor of a less favorable prognosis, were in direct opposition to Hald et al. [31]. In macrophages, natural killer (NK) cells, and dendritic cells, researchers have obtained results that are similarly inconsistent. Due to the specialized roles and prognostic value of TIICs, additional samples and research work that are thoroughly developed are required to confirm the prognostic significance of these cells. In this investigation, we showed that the expression of MAN1B1 was connected with the presence of immune cells in the majority of tumors, which suggests that it may have an effect on the immunological status in BC.

There were several limitations to this study. First, since all of the data used in this study were obtained by downloading them directly from public sources and analyzing them using bioinformatics approaches, the results need to be validated using additional experimental research. Secondly, in order to investigate the possible connection between MAN1B1 and the cancer-related immune microenvironment in BC, more investigations were needed.

5. Conclusion

We discovered that MAN1B1 was overexpressed in BC and that its overexpression was significantly connected to a poor prognosis. There was a possibility that MAN1B1 plays a role in the progression of tumors and metastasis. In addition, our results suggest that MAN1B1 probably plays an important part in the polarization of macrophages and the infiltration of immune cells. Thus, it has the potential to be exploited as a prognostic target for BC.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This work was supported by Jiaxing Science and Technology Program (No. 2020AD30059).