Increased expression of TBC1D10B as a potential prognostic and immunotherapy relevant biomarker in liver hepatocellular carcinoma

The TBC1 domain family member 10B (EPI64B/TBC1D10B), a member of the RabGAP EPI64 subfamily, contains a TBC domain that confers GTPase-activating protein activity. Even though overexpression of TBC1D10B has been reported to promote tumor invasion and metastasis in gastric adenocarcinoma, the prognostic value of TBC1D10B and its correlation with DNA methylation and immune infiltration in hepatocellular carcinoma are still not known. Transcriptional expression profiles of TBC1D10B between hepatocellular carcinoma tissues and normal tissues were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus. The Clinical Proteomic Tumor Analysis Consortium and the Human Protein Atlas were used to assess the TBC1D10B protein expression. The biological functions of TBC1D10B were evaluated by the Metascape database and by Gene Set Enrichment Analysis (GSEA). Receiver operating characteristic (ROC) curve analysis was used to distinguish hepatocellular carcinoma from adjacent normal tissues. The effect of TBC1D10B on survival was estimated using the Kaplan-Meier method. DNA methylation in the TBC1D10B gene was assessed using the online MEXPRESS and MethSurv tools. The association between TBC1D10B mRNA expression and immune cell infiltration was investigated by the TIMER2 web server, tumor immune estimation resource and single-sample GSEA. This study found that TBC1D10B is highly expressed in hepatocellular carcinoma and that increased TBC1D10B mRNA expression is associated with female sex, lower Body Mass Index, high level of alpha fetal protein, and worse clinical stages. The mRNA and protein levels of TBC1D10B were verified in cells. Functional annotation indicated enrichment with negative regulation of the cell cycle, extracellular matrix, and corresponding pathways in the high-TBC1D10B phenotype. The ROC curve analysis showed that, with a cutoff level of 2.912, the accuracy, sensitive, and specificity in differentiate TBC1D10B hepatocellular carcinoma from adjacent controls were 0.931, 0.920, and 0.802, respectively. Kaplan-Meier survival analysis showed that hepatocellular carcinoma patients with high TBC1D10B had a worse prognosis than those with low TBC1D10B, especially in patients with a weight below 70 kg, height above 170 cm, and histological G2 and G3. We also found that the methylation of TBC1D10B was associated with the prognosis in patients with hepatocellular carcinoma. Moreover, correlation analysis indicated that TBC1D10B mRNA expression was positively correlated with infiltration levels of most immune cells, but negatively correlated with Th17 and cytotoxic cells infiltration. Our study indicates that increased TBC1D10B expression in hepatocellular carcinoma may play a role in tumorigenesis by regulating the cell cycle and extracellular matrix. TBC1D10B may be a novel prognostic and predictive marker and immune therapeutic target in hepatocellular carcinoma patients.