Abstract
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. 1 Among kidney cancer, kidney renal clear cell carcinoma (KIRC) is the most dominant subtype, accounting for approximately 75% of all renal cell carcinoma patients 2,3 and the majority of cancer-associated deaths. 4 It was reported that the 5-year disease specific survival rate for kidney cancer less than 10% at stage IV. 5 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. 6,7 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. 8,9 What’s more, studies have shown that immune response is closely associated with clinical outcome in kidney cancer. 10,11 Tumor-infiltrating immune cells (TIIC) play an important role in the TME to regulate tumor development and are appealing therapeutic targets. 10 In multiple clinical and genomic investigations, KIRC has been shown to be significantly associated with immune infiltration. 12,13
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. 14 ESCRT-III was shown to be disorder in several cancers, such as ovarian cancer, 15 prostate cancer 16 and pancreatic tumors. 17 Moreover, cytokinesis disorders of abscission could enhance oncogene-induced mitotic stress to promote genomic instability and tumorigenesis. 18 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. 19
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, 20 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. 21 It used a previously published statistical method to analyze RNA-seq expression profile data to detect the infiltration of immune cells in tumor tissues. 22 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.
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. 23,24
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.
25
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
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
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
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).

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 (*

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 (

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.
MITD1 Expression Level and Clinical Pathological Factors in 534 KIRC Patient.a
aThe boldface values indicate significance with
Top 5 Significantly Correlated GO Term With MITD1 Expression Based on
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,

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.
Top 5 Significantly Correlated KEGG Pathway With MITD1 Expression Based on
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,

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.

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. 28 Despite the development in diagnosis, surgery, drug therapy, radiotherapy and chemotherapy, the clinical outcome of KIRC still remains unsatisfactory. 29 -31 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. 14 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. 18 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. 8,32 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. 19
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. 33 Cytokine-cytokine receptor interaction is usually enriched during tumorigenesis, such as liver cancer 34 and colorectal cancer. 35 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. 36 Targeted immunotherapy is considered to be a promising new strategy for hereditary and sporadic disease. 37 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. 38 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. 38 Also, a study has shown that infiltrating CD 8+ T cells was correlated with kidney injury in patients with anti-glomerular basement membrane disease. 39 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.
Footnotes
Authors’ Note
Chujie Chen and Yiyu Sheng conceived, designed, and performed the project. All data will be provided by reasonable request from the corresponding author. This study did not require an ethical board approval because it was a study based on bioinformatics.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
