Prognostic Impact of MITD1 and Associates With Immune Infiltration in Kidney Renal Clear Cell Carcinoma.

Kidney renal clear cell carcinoma (KIRC) is one of the most malignant diseases with poor survival rate over the world. The tumor microenvironment (TME) is highly related to the oncogenesis, development, and prognosis of KIRC. Thus, making the identification of KIRC biomarkers and immune infiltrates critically important. Microtubule Interacting and Trafficking Domain containing 1(MITD1) was reported to participate in cytokinesis of cell division. In the present study, multiple bioinformatics tools and databases were applied to investigate the expression level and clinical value of MITD1 in KIRC. We found that the expression of MITD1 was significantly increased in KIRC tissues. Further, the KIRC patients with high MITD1 levels showed a worse overall survival (OS) rate and disease free survival (DFS) rate. Otherwise, we found a significant correlation MITD1 expression and the abundance of CD8+ T cells. Functional enrichment analyses revealed that immune response and cytokine-cytokine receptor are very critical signaling pathways which associated with MITD1 in KIRC. In conclusion, our findings indicated that MITD1 may be a potential biomarker and associated with immune infiltration in KIRC.

Introduction

As one of the most common malignant cancers, renal cell carcinoma has been increasing during the last 2 decades, which accounts for 2%-3% of all new cancer patients. Among kidney cancer, kidney renal clear cell carcinoma (KIRC) is the most dominant subtype, accounting for approximately 75% of all renal cell carcinoma patients and the majority of cancer-associated deaths. It was reported that the 5-year disease specific survival rate for kidney cancer less than 10% at stage IV. Moreover, KIRC patients remain a worse prognosis compared to other kidney cancer subtypes patients. Nowadays, more and more therapies, including surgical operation, chemotherapy, radiotherapy and targeted therapies, has been used in the treatment of renal cancer, showing to be effective in the survival of patients. However, these treatments methods are still far from enough to cure the KIRC patients. Thus, it is urgent to find new and effective biomarkers for the prognosis and therapy of KIRC. For example, by comprehensive analysis of RCC clinical datasets, RGS5, TOP2A andMYC showed to play a significant role in RCC. What’s more, studies have shown that immune response is closely associated with clinical outcome in kidney cancer. Tumor-infiltrating immune cells (TIIC) play an important role in the TME to regulate tumor development and are appealing therapeutic targets. In multiple clinical and genomic investigations, KIRC has been shown to be significantly associated with immune infiltration.

MITD1 is a protein coding gene, which located in 2q11.2. And through the ESCRT-III, MITD1 is recruited to midbodies and then participate in cytokinesis abscission. ESCRT-III was shown to be disorder in several cancers, such as ovarian cancer, prostate cancer and pancreatic tumors. Moreover, cytokinesis disorders of abscission could enhance oncogene-induced mitotic stress to promote genomic instability and tumorigenesis. This may indicate that there is a certain connection between MITD1 and cytokinesis in the process of tumorigenesis. Otherwise, MITD1 has been studied to be correlated with immune infiltrates of liver cancer.

In this study, we analyzed the expression of MITD1 in KIRC by Onconmine, TIMER and the Human Protein Atlas. Then, we performed the relationship between MITD1 and clinical information by the TCGA KIRC data. Furthermore, to further investigate the biological processes and pathways associated with the regulatory network of MITD1 that could promote KIRC occurrence and development. The Gene Set Enrichment Analysis (GSEA) was performed to elevate Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Finally, the data on immune system infiltrates of KIRC from TCGA database was analyzed through the ImmuneCellAI and TIMER. The identification of prognostic MITD1 suggests the potential role of MITD1 in KIRC pathogenesis and progression.

Materials and Methods

Data Acquisition

The RNA-seq data of KIRC and normal tissue were obtained from the The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/). Four kidney cancer datasets were obtained from Oncomine (http://www.oncomine.org) database, a web-based data mining platform with RNA-seq and microarray data of human tumors. Some samples with incomplete clinical information were deleted when KIRC samples were analyzed by comprehensive bioinformatics. The RNA-seqdata and clinical information of 534 patients were downloaded from TCGA database. Subsequently, by chi-square test, the statistical difference of these characteristics with MITD1 expression was analyzed.

TIMER Database Analysis

TIMER (http://cistrome.dfci.harvard.edu/TIMER/) is a comprehensive database for systematic analysis of immune infiltrates with different kinds of cancer. It used a previously published statistical method to analyze RNA-seq expression profile data to detect the infiltration of immune cells in tumor tissues. The current release of TIMER incorporated 10,009 samples across 23 cancer types from TCGA dataset to estimate the abundance of immune infiltrates.

Assessment of Prognosis Clinical Characteristics and MITD1 in KIRC Cohort

Based on the expression level of MITD1, we divided the 534 KIRC samples into 2 groups: high MITD1 expression group and low MITD1 expression group. Then we performed Kaplan-Meier analysis on 534 KIRC patients through SPSS software (version 24.0, Chicago, IL, United States). A univariate Cox proportional hazard regression analysis was used to evaluate the association between overall survival (OS) time and clinical characteristics in the trial cohort. P-value < 0.05 was considered to have statistically significant prognostic value.

The Human Protein Atlas

The Human Protein Atlas database (HPA) (https://www.proteinatlas.org) was used to analyze protein expression of MITD1 between normal and KIRC tissues. HPA portal is an open database with millions of high-resolution images showing the spatial distribution of proteins in 44 different normal human tissues and 20 different cancer types.

Gene Set Enrichment Analysis

Gene Set Enrichment Analysis (GSEA) is a computational method used to determine whether a set of priori defined genes show statistically significant, concordant differences between 2 biological states. GSEA preliminarily classified the genes according to their correlation with the MITD1 expression. This computational method illustrates the significant difference in survival observed between the high and low MITD1 groups. For each analysis, we ran 1000 repetitions of gene set permutations. In addition, nominal P values and standardized enrichment scores (NES) were used to classify enrichment pathways in each phenotype.

Correlation Analysis Between MITD1 and Immune Infiltration in KIRC

The composition and abundance of immune cells in the tumor microenvironment have a great influence on tumor progression and the effect of immunotherapy. GraphPad Prism software (version 8.4.2, https://www.graphpad.com/updates/) was used to analyze the correlation between MITD1 expression and TIICs, the cutoff was P-value < 0.01. The Immune Cell Abundance Identifier (ImmuCellAI) detected the abundance of immune cells, and the differences in the infiltration immune cells of 534 KIRC were analyzed. ImmuneCellAI is a tool for estimating the abundance of 24 immune cells from gene expression datasets, including RNA-Seq and microarray data, of which 24 immune cells are composed of 6 immune cell types and 18 T-cell subtypes (http://bioinfo.life.hust.edu.cn/web/ImmuCellAI/).

Statistical Analysis

All statistical analyses were enforced with corresponding packages of R (Version 3.6.3). Kaplan-Meier survival analysis and chi square test was conducted using SPSS (version 24.0). To calculate the 95% CI and HR, univariate Cox regression analyses were used to evaluate MITD1 expression between clinical factors and OS. For all analysis, the threshold P-value < 0.05 was considered significant statistically.

Results

MITD1 Was Upregulated in KIRC

The mRNA-seq data from TCGA datasets were analyzed by the TIMER tools, which showed the expression levels of MITD1 in different kinds of tumors. MITD1 with a higher expression level existed in most cancers, including BLCA, ESCA, COAD, GBM, HNSC, KIRP, LIHC, LUAD, LUSC, PRAD, STAD, READ as well as KIRC (Figure 1A). Then, the expression level of MITD1 in tumors and normal tissues of kidney were analyzed by Oncomine database. As shown in Figure 1B, higher expression level of MITD1 was presented in kidney cancers which included 4 studies. Four studies on MITD1 mRNA levels in Renal cancer vs normal kidney tissue was shown in Figure 1C to F. Further, immunohistochemistry (IHC) results available from the HPA database showed that MITD1 has a higher expression in KIRC tumor tissues compared to nontumor tissues (Figure 2).

Figure 1.

MITD1 mRNA expression levels in normal and cancer tissues. A, Human MITD1 expression levels in different tumor types from TCGA database were determined by TIMER database (*P < 0.05, **P < 0.01, ***P < 0.001). B, MITD1 in datasets of different cancers compared with normal tissues in the oncomine database. Red color indicates upregulated expression in cancer, and blue color indicates downregulated expression in cancer. C-F, Meta-analysis of the 4 datasets from 4 studies on MITD1 mRNA levels in renal cancer vs normal kidney tissue searched by oncomine database.

Figure 2.

MITD1 protein expression levels in normal and KIRC tissues. The expression of MITD1 protein was visualized using immunohistochemistry via the human protein atlas database.

Connection Between the Expression of MITD1 and Clinicopathology

Then, univariate analyses were used to evaluate whether the MITD1 expression level and various clinic pathological characteristics were independent prognostic parameters of KIRC patient overall survival. We found that high expression of the MITD1 gene, age, hemoglobin, laterality, neoplasm, pathologic M, pathologic N, pathologic T, platelet, tumor stage, intermediate dimension and longest dimension were independent prognostic factors in patients with KIRC (Figure 3A). Next, the potentially prognostic significance of MITD1 expression levels was investigated by TCGA KIRC data. The patients’ clinical classification information is shown in Table 1, as a result, high MITD1 expression found to be closely correlated with pathologic T stage (P = 0), platelet qualitative (P = 0.009), live status (P = 0) and tumor stage (P = 0.001). These suggested that MITD1 could be an independent prognostic indicator for KIRC. The distribution of MITD1 expression and survival status of KIRC patients were shown in Figure 3B and C. Kaplan-Meier survival analysis of the KIRC dataset showed that MITD1 expression was significantly with OS (P = 0.001) and DFS (P = 0.001). Low MITD1 expression was associated with better prognosis, and high expression of MITD1 was associated with poor prognosis in KIRC patients (Figure 3D and E). As the time dependent ROC curve showed that the 5-year and 10-year AUC of MITD1 expression were 0.599 and 0.722 for predicting survival, respectively (Figure 3F), which indicated that MITD1 has great prognostic value.

Figure 3.

Univariate Cox analysis and survival curve of differential M1TD1 expression. A, Univariate Cox analysis of MITD1 expression and other clinicopathological variables. B-C, MITD1 expression distribution and survival status. D, Overall survival curve of differential MITD1 expression. E, Disease free survival curve of differential MITD1expression. F, ROC curves of MITD1 of 5 year and 10 year.

Table 1.

MITD1 Expression Level and Clinical Pathological Factors in 534 KIRC Patient.a

Characteristics	Case (No.)	No. of patients (%)	Expression level		P value
			Low (No. cases)	High (No. cases)
Age	534				0.193
<60		247 (46.3)	131	116
≥60		287 (53.7)	136	151
Hemoglobin_result	454				0.191
Normal		187 (41.2)	102	85
Non-normal		267 (58.8)	129	138
Laterality	533				0.83
Left		250 (46.9)	124	126
Right		283 (53.1)	143	140
Neoplasm_histologic_grade	526				0.387
G1		14 (2.7)	7	7
G2		231 (43.9)	122	109
G3		206 (39.2)	101	105
G4		75 (14.3)	31	44
Pathologic_M	502				0.081
M0		424 (84.5)	225	199
M1		78 (15.5)	33	45
Pathologic_N	256				0.771
N0		240 (93.8)	114	126
N1		16 (6.3)	7	9
Pathologic_T	534				0
T1		274 (51.3)	162	112
T2		70 (13.1)	25	45
T3		179 (33.5)	75	104
T4		11 (2.1)	5	6
Platelet_qualitative_result	444				0.009
Low		47 (10.6)	20	27
Normal		360 (81.1)	195	165
Elevated		37 (8.3)	11	26
Serum_calcium_result	364				0.957
Normal		151 (41.5)	74	77
Non-normal		213 (58.5)	105	108
Gender.demographic	534				0.586
female		186 (34.8)	96	90
male		348 (65.2)	171	177
status	534				0
Alive		361 (67.6)	200	161
Dead		173 (32.4)	67	106
Tumor_stage	531				0.001
Stage I		268 (50.5)	158	110
Stage II		58 (10.9)	23	35
Stage III		123 (23.2)	50	73
Stage IV		82 (15.4)	35	47
Intermediate_dimension	495				0.4
Low		290 (58.6)	154	136
High		205 (41.4)	101	104
Tumor dimension	495				0.059
Low		267 (53.9)	148	119
High		228 (46.1)	107	121

aThe boldface values indicate significance with P value < 0.05.

Biological Function Analysis of M1TD1

In order to investigate the potential biological functions of MITD1 in KIRC, the GO term and KEGG pathway enrichment analyses were used to address our needs. We used the GSEA to reveal the significantly enriched (FDR < 0.05, P-value < 0.05) GO terms and KEGG pathways with high MITD1 expression. The top 5 GO terms positively correlated with high expression of MITD1were innate immune response, cytokine-mediated signaling pathway, lymphocyte activation, regulation of leukocyte activation and T cell activation. The top 5 5 negatively correlated categories: regulation of polysaccharide biosynthetic process, regulation of smooth muscle cell differentiation, retina morphogenesis in camera-type eye, cornification and positive regulation of smooth muscle cell migration. As shown in Table 2 and Figure 4A, the biological processes (BPs) strongly associated with high expression of MITD1 were immune response and cell life. And the most significantly positively changed 5 pathways were Ribosome, Viral protein interaction with cytokine and cytokine receptor, Primary immunodeficiency, Cytokine-cytokine receptor interaction and Autoimmune thyroid disease. Also, the top 5 negative correlated pathways were Aldosterone-regulated sodium reabsorption Citrate cycle (TCA cycle), Olfactory transduction, Regulation of lipolysis in adipocytes, Gastric acid secretion, as shown in Table 3 and Figure 4B. These results indicated that cytokine receptor and citrate cycle are very critical pathways which associated with MITD1 in KIRC.

Table 2.

Top 5 Significantly Correlated GO Term With MITD1 Expression Based on P-Value.

Regulation	GO term	NES	P-value
Positive	Innate immune response	1.416384989	0.001124859
	Cytokine-mediated signaling pathway	1.435599139	0.001128668
	Lymphocyte activation	1.562070926	0.001148106
	Regulation of leukocyte activation	1.465387692	0.001157407
	T cell activation	1.695728707	0.001187648
Negative	Regulation of polysaccharide biosynthetic process	−2.17911054	0.002304147
	Regulation of smooth muscle cell differentiation	−1.87121347	0.002304147
	Retina morphogenesis in camera-type eye	−1.88561074	0.002304147
	Cornification	−1.82220123	0.002304147
	Positive regulation of smooth muscle cell migration	−2.01859177	0.002336449

Figure 4.

Enrichment plots from gene set enrichment analysis. A, GSEA results showing differential enrichment of genes in GO biological process term with high MITD1 expression. B, GSEA results showing differential enrichment of genes in KEGG pathway with high MITD1 expression.

Table 3.

Top 5 Significantly Correlated KEGG Pathway With MITD1 Expression Based on P-Value.

Regulation	KEGG pathway	NES	P-value
Positive	Ribosome	2.111089495	0.000140095
	Viral protein interaction with cytokine and cytokine receptor	2.123368383	0.000152975
	Primary immunodeficiency	2.055747293	0.000166667
	Cytokine-cytokine receptor interaction	1.704887989	0.000268132
	Autoimmune thyroid disease	2.040880331	0.000332557
Negative	Aldosterone-regulated sodium reabsorption	−2.17485499	0.000249377
	Citrate cycle (TCA cycle)	−1.99797266	0.000249875
	Olfactory transduction	−2.063311	0.000249875
	Regulation of lipolysis in adipocytes	−2.07478781	0.000273373
	Gastric acid secretion	−2.11046891	0.000277855

Immune Infiltration of MITD1 in KIRC

In order to investigate the possible effect of MITD1 on different immune cell types in the KIRC microenvironment, the ImmuCellAI database was utilized to calculate the ingredients of 24 immune cells in 534 KIRC sample. As was shown in Figure 5A, the immune cells varied significantly among different samples. Therefore, we compared the level of immune cells in the high expression MITD1 group and low expression group in the above study (Figure 5B). Most of the immune cells were significantly correlated with MITD1 expression. We found that CD8 navie, Th17, NKT, Monocyte and Neutrophil cells were negative correlation with MITD1 expression. Cytotoxic, Th1, Thf, effector memory, NK, Gamma delta and CD8 T cells were positively correlated with MITD1 expression. Further, the TIMER web tool was performed to investigate the association between MITD1 and tumor immune infiltrating levels (Figure 6A). Results indicated that MITD1 was not significantly correlated with tumor purity (r = −0.036, P > 0.05), CD4+ T cells (r = 0.024, P > 0.05), DC cells (r = 0.103, P > 0.05), B cell (r = 0.118, P > 0.05) and Macrophage (r = −0.098, P > 0.05) but positively correlated with infiltrating levels of CD8+ T cells (r = 0.258, P = 1.82e−8) and neutrophil (r = 0.255, P = 2.88e−8). All the analyses above may indicate that MITD1 expression is correlated with immune infiltration level in KIRC, especially by CD8+ T cells.

Figure 5.

The correlation of MITD1 and immune infiltration in KIRC. (A) Immune cells estimated by ImmuCellAI in KIRC; (B) the comparison of immune cells with high and low expression of MITD1.

Figure 6.

Correlations between MITD1 expression and immune infiltration levels in KIRC by TIMER.

Discussion

KIRC is the main kind of renal cell carcinoma, which takes up for 60%-85% in all renal cell carcinoma cases. Despite the development in diagnosis, surgery, drug therapy, radiotherapy and chemotherapy, the clinical outcome of KIRC still remains unsatisfactory. Although many studies have reported possible diagnostic or therapeutic targets for renal cancer, but remain mostly theoretical. As a well-known heterogeneous disease, KIRC still lacks specific diagnostic biomarkers for individualized therapy, especially in immunotherapy.

MITD1 protein is considered to associate with the midbody. A previous study reported that ESCRT-III subunits are very important for the recruitment of MITD1 to the midbody. Also, ESCRT-III complex helps MITD1 to participate in the abscission of cytokinesis. ESCRT-III is in synergy with oncogene-induced mitotic stress, the abnormal regulation of abscission may lead to oncogenesis and genetic instability. During the process, MITD1 could coordinate ESCRT-III activity during abscission with early events in the final phase of cell division. The genomic and clinical datasets on the public database have been used to mined to find potential novel biomarkers and underlying mechanism of kidney cancer. In this study, we integrated a variety of bioinformatics approaches to prove that MITD1 could be a prognostic biomarker for KIRC. The TCGA RNA-seq data of the KIRC cases were used to analyze the prognostic value of MITD1. All the relationship of MITD1 with various tumor characteristics and immune cell responses support that high expression of MITD1 showed a poor survival of KIRC. Through our analyses, we found that high MITD1 expression of KIRC patients are more likely to present a more advanced pathologic T stage, platelet qualitative, live status and tumor stage than those patients with low expression of MITD1. We speculated that dysregulated expression of MITD1 could play an important role of tumorigenesis and tumor immunology in KIRC progression. Coincidentally, the founding was consistent with the expression model of MITD1 in liver cancer.

The functional enrichment analyses of MITD1 co-expressed related genes by GSEA were investigated to further figure out the function and associated signaling pathways in KIRC. Our work suggested that higher expression of MITD1 in KIRC patients could mainly alter immune response and cytokine-cytokine receptor interaction. Tumor immune infiltration by different kinds of immune cells is an active immune response against various tumors. Cytokine-cytokine receptor interaction is usually enriched during tumorigenesis, such as liver cancer and colorectal cancer. All these above suggested that MITD1 may be a potential biomarker of prognosis and immune therapeutic target of KIRC.

Over the past few years, the medical treatment of RCC has transited from nonspecific immune approaches to novel immunotherapy agents. Targeted immunotherapy is considered to be a promising new strategy for hereditary and sporadic disease. Subsequently, ImmuCellAI tools and TIMER database were used to reveal the relationship between MITD1 expression and immune cells infiltration levels in KIRC. By comprehensive analyses, we found that the relationship betweenMITD1 and CD8+ T cells is the strongest. CD8+ T cells are the preferred immune cells for targeting cancer cell therapy. Macrophage type 2 (M2) cells, Cancer-associated fibroblasts (CAFs), and regulatory T cells (Tregs) could form immune barriers against CD8+ T cell-mediated anti-tumor immune responses. Thus, CD8+ T cells need to initiate and activate effector CTLs in a process known as the tumor immune cycle to produce a durable and effective anti-tumor immune response. Also, a study has shown that infiltrating CD 8+ T cells was correlated with kidney injury in patients with anti-glomerular basement membrane disease. In a word, these results indicated that MITD1 could regulate and recruit immune infiltrating cells in KIRC. However, in order to more accurately understand the relationship between MITD1 and CD8 + T cells in vivo and in vitro, controlled and multicenter clinical trials are needed.

Taken together, we found that MITD1 is highly expressed in KIRC through integrated bioinformatics analyses. Otherwise, increased MITD1 expression indicates poor prognosis and increased immune infiltration levels of CD8+ T cells. Therefore, MITD1 may be a potential prognostic biomarker for KIRC and play a vital role in immune cell infiltration.