Skip to main content

Identification of HM13 as a prognostic indicator and a predictive biomarker for immunotherapy in hepatocellular carcinoma



Histocompatibility minor 13 (HM13) is a signal sequence stubbed intramembrane cleavage catalytic protein that is essential for cell signaling, intracellular communication, and cancer. However, the expression of HM13 and its prognostic value, association with tumor-infiltrating immune cells (TIICs) in the microenvironment, and potential to predict immunotherapeutic response in HCC are unknown.


The HM13 expression, clinicopathology analysis, and its influence on survival were analyzed in multiple public databases and further verified in collected HCC and normal tissues by qRT-PCR and immunohistochemistry staining assay (IHC). Furthermore, the lentivirus vector encoding HM13-shRNA to manipulate HM13 expression was selected to investigate whether HM13 could influence the malignant growth and metastasis potential of HCC cells. Finally, significant impacts of HM13 on the HCC tumor microenvironment (TME) and reaction to immune checkpoint inhibitors were analyzed.


Upregulated HM13 was substantially correlated with poor prognosis in patients with HCC, and could facilitate the proliferation and migratory potential of HCC cells. Additionally, patients with high HM13 expression might be more sensitive to immunotherapy.


HM13 might be a prognostic biomarker and potential molecular therapeutic target for HCC.

Peer Review reports


Hepatocellular carcinoma (HCC), which accounts for almost 90% of all primary liver cancers, is one of the most frequent tumors and the leading cause of cancer-related deaths in the world [1]. The incidence and development of HBV or HCV infection, alcohol abuse, HIV co-infection, cirrhosis, and aflatoxin exposure may vary depending on the prevalence of risk factors such as HBV or HCV infection, alcohol abuse, HIV co-infection, cirrhosis, and aflatoxin exposure [2,3,4]. Furthermore, a growing body of evidence links diabetes, metabolic syndrome, and obesity to an increased risk of HCC in patients with non-alcoholic fatty liver disease (NAFLD) [5,6,7,8], implying that people who follow healthy lifestyle habits, such as a regular exercise routine and avoiding metabolic syndrome, may have a lower risk of developing HCC [9]. Most patients with early and intermediate stage HCC are known to receive therapeutic options such as liver resection, liver transplantation, radiofrequency ablation, and embolization [10]. However, due to a lack of biomarkers for early detection and the insidious nature of HCC, more than half of patients with HCC are already in an advanced stage by the time they are detected, at which point they are often unsuitable for curative therapy and the prognosis is still poor [11, 12]. Furthermore, many patients continue to experience tumor recurrence months or years later [13]. Although liver biopsy is the most reliable means of diagnosing HCC, it is an intrusive procedure with risks [14]. As a result, novel approaches to prognostic biomarkers and therapeutic targets are still needed to ameliorate dismal prognosis.

Histocompatibility small 13 (HM13) encodes a protein that localizes to the endoplasmic reticulum and catalyzes its intramembrane proteolysis following cleavage of some signal peptides from preproteins, which is required for the generation of lymphocyte surface (HLA-E) epitopes derived from MHC class I signal peptides [15, 16]. Growing evidence supports the association between HM13 expression, tumor-infiltrating immune cells (TIICs) [16, 17], and cancer [18, 19]. However, the function of HM13 in HCC has never been reported before. This study aimed to investigate the relationship between HM13 expression profiles, clinicopathological features, immune infiltration, and immunotherapy responses in HCC patients and to explore the role of HM13 in the development of HCC, which may contribute to our understanding of the mechanisms of HCC.


Bioinformatics analysis by public databases

A pan-cancer study of HM13 expression was performed using the Gene Expression Profiling Interactive Analysis (GEPIA) database [20]. The GEO, TCGA, and ICGC databases were used to investigate the relationship between HM13 expression and clinicopathologic features in samples with and without HCC. The human protein atlas database (HPA) was used to compare the protein expression of HM13 in normal and HCC tissues ( The researchers utilized R software to do a Cox regression analysis to see if HM13 was an independent predictive factor in the TCGA study. The potential of HM13 in predicting immunotherapeutic response was investigated using the IMvigor210 cohort.

Clinical data

Fifty-nine HCC samples and 50 normal tissues were selected from patients who underwent liver biopsy. Before surgery, none of the patients had had radiation or chemotherapy, and those with missing data were eliminated. Two prominent pathologists evaluated and co-diagnosed all pathological data.

Quantitative real-time PCR, immunohistochemistry staining, and western-blotting assay

As previously described [21], we used quantitative real-time PCR (qRT-PCR, F603101, Sangon Biotech, China) and immunohistochemical staining test (IHC) to compare the mRNA and protein expression levels of the HM13 and CHOP genes in HCC and normal tissues. The following are the primer sequences: HM13, Forward: GATCCGCATAACGGCAGTG, Reverse GAGCCCCAAGAGTGTGCAG, CHOP, Forward: GGAAACAGAGTGGTCATTCCC, Reverse CTGCTTGAGCCGTTCATTCTC, and β-ACTIN, Forward: CGTGGGCCGCCCTAGGCACCA, Reverse TTGGCTTAGGGTTCAGGGGGG. The sections were subjected to IHC after being heated in citrate buffer for antigen retrieval, quenching endogenous peroxidase with 2% hydrogen peroxide, and blocking with 5% bovine serum albumin (BSA, A602449, Sangon Biotech, China). After that, the slices were incubated at 4 °C for the duration of the night with an anti-HM13 (ab247061, Abcam) antibody and an isotype control. DAB chromogen (A690009, Sangon Biotech, China) was used to stain the sections progressively after washing with PBST (Phosphate Buffer Solution with Tween-20). Two competent pathologists performed the IHC findings assessment using the single-blind and unified criteria procedures. To distinguish between low/loss (4) and high (> 4) expression of HM13 in HCC and normal tissues, a final score of the sum of the extent of expression score (no positive cells = 0, 10% = 1, 10–50% = 2, positive staining of > 50% = 3) and intensity score (negative = 0, weak = 1, moderate = 2, strong = 3) was used. Finally, the western-blotting (WB) assay was used to examine the knockdown effect of HM13 protein in MHCC97L cells. RIPA buffer (20–188, Sigma) was applied to the whole-cell lysates of MHCC97L cells. A 12% SDS-PAGE gel was used to separate 30 g of protein, which was then transferred to a PVDF membrane. Membranes were incubated with anti-HM13 (1:1000, ab247061, Abcam) and anti-β-ACTIN (1:2000, ab8227, Abcam) antibodies at 4 °C overnight. Blocking was performed with a 5% blocking buffer. Following that, membranes were treated for 1 hour with the secondary antibody Anti-HRP Rabbit (1:10000, 207674, Abcam). To produce protein bands, an ECL reagent (P0018S, Beyotime, China) was utilized. Using Image Lab software, the relative expression of HM13 was normalized to β-ACTIN. The experiments covered in our study were all set up with 2 additional independent controls.

Cell lines and cell culture

The Shanghai Institute of Cells (Chinese Academy of Science, Shanghai, China) provided cell lines (HepG2, Hep3B, HuH7, and MHCC97L). Cells were grown in prescribed DMEM medium (E600008, Sangon Biotech, China) with 10% fetal bovine serum (FBS, E600001, Sangon Biotech, China) at 37 °C with 5% CO2 in 100% humidity.

Plasmids construction and cell infection

The Lentiviral vector encoding HM13’s Lenti-shRNA vector system (PCDH-GFP) (HM13-shRNA sequence: CCGGGTGCCTGAAACAATCACCAGCCACTCGAGTGGCTGGTGATTGTTTCAGGCACTTTTT) were chosen. Control vectors were shRNA controls with a non-targeting sequence. The HM13-shRNA was selected for creating stable knockdown (KD) cell lines using lipofectamine™ 3000 transfection reagent (L3000001, Invitrogen, USA) according to the manufacturer’s instructions, and MHCC97L cells were grown in DMEM (E600008, Sangon Biotech, China) supplemented with 10% FBS (E600001, Sangon Biotech, China) for cell infection.

Cell proliferation assay

The cell counting kit-8 assay was used to determine the vitality of the cells (CCK-8, E606335, Sangon Biotech, China). Cells were seeded at a density of 1 × 105 cells per well in 96-well cell culture clusters and cultivated for 1 day, 2 days, and 3 days, respectively. After culture, 10 μL CCK-8 solution was applied to each well, and the absorbance was measured with a microplate reader at 450 nm within 4 hours.

Wound-healing assay

In a serum-free medium, cells were grown to about 80% confluence for 12 hours before being scratched with a sterile pipette tip. After scratching, cells were grown for the following 24 hours. The migrative ability was quantified and analyzed by measuring the scratch area which was not covered by cells at 100 x magnification. Three independent duplicates were needed to ensure the accuracy of this assay and the wound closure rate was calculated as [1 - (wound area / original wound area)] from photographs acquired at 24 h.

Transwell assay

The 24-Transwell chambers were used in the invasion experiment (3244, Corning Inc., Corning, USA). After coating transwell filters with Matrigel (BD Biosciences, USA), cultured cells were resuspended at a density of 1 × 104 cells per mL in 200 μL serum-free DMEM (E600008, Sangon Biotech, China), applied to the Matrigel membrane, and plated in the upper chamber while the bottom chamber was filled with 500 μL DMEM supplemented with 10% FBS (E600001, Sangon Biotech, China). Cells adhering to the bottom chamber were stained with 0.1% crystal violet in PBS for 15 minutes and counted under microscopy at 200 x magnification after being cultured at 37 °C for 48 hours.

Immune infiltrates analysis and signaling pathway analysis

To evaluate the potential relationships between HM13 expression profiles and TIICs in the HCC microenvironment, the TCGA database was used to gauge the proportions of TIICs via CIBERSORT [22] and the Timer database [23]. Subsequently, to explore the signaling pathways involved in HM13 mediated effects, Gene set enrichment analysis (GSEA) was performed to distinguish the altogether cautioned Hallmark items with FDR < 0.05.

Genetic alterations and TMB analysis

The somatic variants data of 350 HCC patients were downloaded from TCGA to analyze the difference of genetic alterations between the high- and low-HM13 expression subgroups with R package “maftools”, and the tumor mutation burden (TMB) of each patient was subsequently assessed.

HM13 methylation level analysis

The data on HM13 copy number variation (CNV) and methylation level was obtained from the cBioPortal web platform (, and a comparison of the varying HM13 gene expressions in HM13 copy number variation groups and the correlation between HM13 methylation level and HM13 gene expression were performed. TCGA data was utilized to determine changes in the promoter methylation level of HM13 between HCC and normal tissues using the UALCAN web tool ( The MethSurv online tool ( was utilized to investigate the HM13 methylation level’s predictive value in the TCGA-LIHC cohort.

Statistical analysis

R software was used to conduct statistical analysis (Version 4.0.2). When applicable, categorical data were compared using the Pearson chi-square test or Fisher exact test, and quantitative data between the two groups were compared using the t-test. The influence of HM13 expression on the OS rate of patients with HCC was assessed using Kaplan-Meier survival analysis with the log-rank test. The link between HM13 expression and overall survival, as well as other clinical variables, was investigated using a Cox proportional hazards regression model. One-way analysis of variance (ANOVA) was performed to analyze several groups. For HM13 expression, a P value of 0.05 was considered statistically significant. A flow chart was shown in Fig. 1.

Fig. 1
figure 1

A flow chart


Pan-cancer analysis of HM13 expression in TCGA database

The results of the pan-cancer analysis performed in GEPIA showed that HM13 was significantly increased in various tumor tissues, including BLCA, BRCA, COAD, DLBC, GBM, KAML, LGG, LIHC, OV, READ, SKCM, STAD, TGCT, THYM, UCEC, and UCS, when compared to normal tissues (Fig. 2A, p < 0.01), however, HM13 was exclusively linked to poor overall survival and disease-free survival in patients with LGG (Fig. 2B, p < 0.05) and LIHC (Fig. 2C, p < 0.05). Considering that clinical samples of LGG were difficult to collect, and we had collected clinical samples of HCC and normal tissues, subsequently, we focused on the role of HM13 in HCC.

Fig. 2
figure 2

Pan-cancer analysis of HM13 expression in TCGA. A HM13 expression in 33 types of cancers. High HM13 expression was related to overall survival and disease-free survival in patients with LGG (B) and LIHC (C)

HM13 expression is significantly increased in HCC

The results analyzed in multiple databases including GSE25097, GSE62232, GSE76427, GSE101685, GSE36376, TCGA, ICGC, and collected clinical samples showed that HM13 expression level in HCC tissues was significantly higher than that in normal liver tissues (Fig. 3A, p < 0.05), and differential protein expression of HM13 between normal and HCC tissues analyzed by HPA (Fig. 3B) and IHC (Fig. 3C) were consistent with the above results. Among these clinical samples, we detected 3 cases of HM13High and 47 cases of HM13Low/loss in normal tissues, 26 cases of HM13High and 33 cases of HM13Low/loss in tumor tissues (p < 0.05), respectively. In summary, HM13 expression levels, both mRNA and protein expression, were significantly increased in HCC.

Fig. 3
figure 3

HM13 was overexpressed in HCC. A HM13 mRNA expression was analyzed in multiple databases. B HM13 protein expression was analyzed in HPA databases. C HM13 protein expression was analyzed by IHC

HM13 expression increased with the progression of HCC

The relationship between HM13 expression and tumor progression was explored and the results revealed that both HM13 mRNA (Fig. 4A, p < 0.05) and protein expression (Fig. 4B, p < 0.05) significantly increased in later stages of HCC compared with that in early stages. Among the 59 tumor tissues, we detected 12 cases of HM13High and 30 cases of HM13Low/loss in early stages tissues, and 14 cases of HM13High and 3 cases of HM13Low/loss in later stages tissues (p < 0.05), respectively.

Fig. 4
figure 4

HM13 expression increased with the progression of HCC. A HM13 mRNA expression in early and later stages of HCC. B HM13 protein expression in early and later stages of HCC

Moreover, HM13 expression was also related to grade and relapse in HCC patients (Table 1).

Table 1 Relationship between HM13 expression and clinical pathological features

Increased HM13 expression indicated a poor survival outcome of HCC

Kaplan-Meier analysis was performed to explore the relationship between HM13 expression and overall survival outcome and the results revealed that patients with lower HM13 expression had a better overall survival rate (Fig. 5A, p < 0.05). Furthermore, we found that HM13 expression was an independent prognostic factor for HCC patients by Cox proportional regression model in TCGA (Fig. 5B, p < 0.05), ICGC (Fig. 5C, p < 0.05), and clinical sample cohort (Fig. 5D, p < 0.05).

Fig. 5
figure 5

Overexpressed HM13 indicated a poor prognosis for HCC. A Kaplan–Meier plots of HM13 expression in three HCC datasets. B-D Univariate and multivariate Cox regression analysis revealed that HM13 expression was an independent prognostic factor in three HCC databases

HM13 expression was involved in immune infiltration in HCC

The TIMER database was used to assess the correlation between HM13 expression and TIICs levels using the Pearson method and the result showed that HM13 was significantly related to infiltration levels of B cell, CD4+ T cell, CD8+ T cell, macrophage, neutrophil, and dendritic cell (Fig. 6A). Furthermore, difference analysis of TIICs infiltration levels between high- and low-HM13 expression subgroups was assessed by Wilcoxon signed-rank test based on the CIBERSORT algorithm and the result demonstrated that HCC patients with higher HM13 expression had modestly increased ratios of macrophages M0 and resting dendritic cells, while patients with lower HM13 expression had significantly elevated ratios of resting NK cell, macrophages. M2 and resting mast cells (Fig. 6B).

Fig. 6
figure 6

HM13 expression was involved in immune infiltration in HCC. TIICs analysis by Timer database (A) and CIBERSORT (B)

Relationship between HM13 and immune checkpoint genes and the potential of HM13 in predicting immunotherapeutic response in HCC

The results of genetic alterations analysis indicated that the mutation rates of the top 10 most significantly mutated genes were remarkably different in the two subgroups (Fig. 7A). The most significantly mutated gene was TP53 in the HM13-high subgroup, while CTNNB1 was the most significantly mutated gene in the HM13-low subgroup. Subsequently, the TMB of each patient was assessed. HM13 expression was closely related to TMB, and patients with higher HM13 expression had modestly increased ratios of TMB (Fig. 7B, p < 0.05). In the following, we found the expression of PD1 and CTLA4 were significantly related to the expression of HM13 (Fig. 7C, p < 0.001), and patients with higher HM13 expression had significantly increased PD1 and CTLA4 (Fig. 7D, p < 0.001), indicating that immune checkpoint inhibitors (ICIs) treatment was more effective for patients with higher HM13 expression. Finally, the IMvigor210 cohort was used to analyze the potential of HM13 in predicting immunotherapeutic response and we found that patients with high HM13 expression might be more sensitive to immunotherapy (Fig. 7E, p < 0.05).

Fig. 7
figure 7

Relationship between HM13 and immune checkpoint genes and the potential of HM13 in predicting immunotherapeutic response in HCC. A Oncoplots of top 10 mutated genes in high- and low-HM13 expression subgroups. B Correlation and difference analysis of HM13 expression and TMB. C-D Correlation and difference analysis of HM13 and PD1 and CTLA4 expression. E HM13 expression in different anti-PD1 and antiCTLA4 responses of patients from the IMvigor210 cohort

HM13 promoted the proliferation and migratory potential of HCC cells

The HM13 expression level was assessed in four HCC cell lines including HepG2, Hep3B, HuH7, and MHCC97L (Fig. 8A, p < 0.05). Next, we selected the lentivirus vector encoding HM13-shRNA to manipulate HM13 expression in MHCC97L cell lines. HM13 expressions were confirmed by qRT-PCR (Fig. 8B, p < 0.001). To investigate the influence of HM13 on the malignant growth and metastasis potential, cell proliferation assay, wound-healing assay, and transwell assay were performed. For cell proliferation ability, the results of CCK-8 showed that the proliferation rate of down-regulated HM13 was decreased (Fig. 8C, p < 0.05). For the wound-healing assay, we found that HM13 knockdown remarkably inhibited wound closure (Fig. 8D, p < 0.01). Meanwhile, in transwell chambers with Matrigel assay, HM13 knockdown dramatically reduced cell invasion (Fig. 8E, p < 0.001).

Fig. 8
figure 8

HM13 promoted the proliferation and migratory potential of HCC cells. A HM13 expression in four HCC cell lines. B Stable cell lines for knockdown of HM13. C CCK-8 assay. D Wound-healing assay. E Transwell assay

Signaling pathway analysis of HM13-mediated effects

According to the results of the GSEA, the expression of HM13 was significantly correlated with HALLMARK DNA REPAIR and HALLMARK UNFOLDED PROTEIN RESPONSE (Fig. 9A), suggesting that HM13 may be associated with endoplasmic reticulum stress and unresponsive protein response (UPR). Excitingly, consistent with the above analysis, we used qRT-PCR detection to find that HM13 was significantly correlated with the expression of URP-related pathway molecule CHOP (Fig. 9B, p < 0.001), and CHOP was also highly expressed in HCC tissue samples with high HM13 expression (Fig. 9C, p < 0.05). In addition, the expression of CHOP was also significantly reduced after HM13 knockdown (Fig. 9D, p < 0.01). All these suggest that HM13 may play an important regulatory role in the activation of PERK/CHOP pathway.

Fig. 9
figure 9

Signaling pathway analysis of HM13-mediated effects A GSEA analysis. B Correlation analysis. C Difference analysis. D CHOP expression was significantly reduced after HM13 knockdown

HM13 Hypomethylation in HCC

Patients with HM13 CNV gain had a greater level of HM13 expression (Fig. 10A, p < 0.05), according to the cBioPortal, however HM13 methylation was not connected to HM13 gene expression (Fig. 10B, p > 0.05). The promoter methylation of HM13 in HCC tissues was substantially greater than in surrounding normal tissues (Fig. 10C, p < 0.05), according to the UALCAN. Patients with decreased HM13 methylation had a better overall survival (Fig. 10D, p < 0.01), according to the MethSurv online tool.

Fig. 10
figure 10

HM13 Hypomethylation in HCC. A Expression levels in different CNV of HM13. B Correlation analysis between HM13 methylation and its expression. C The promoter methylation of HM13 in HCC. D Kaplan-Meier survival of the promoter methylation of HM13


With the characteristics of high malignancy and rapid progress, more than half of patients with HCC were diagnosed at a late stage with a worse clinical outcome and would die within 3–6 months after diagnosis, which mainly resulted from tumor location, tumor cell invasion, and metastasis, tumor recurrence, and body characteristics, although significant advances have been achieved in the medical fields of anticancer therapeutics [1, 24,25,26]. Various genetic variations have a role in the start and progression of HCC. Identifying genes that are dysregulated during tumor growth can therefore aid in improving prognosis and therapy options. HM13, an intramembrane cleavage catalytic protein of signal sequence stubs, was proved to play an important role in the maturation of tail-anchored protein [16] and dislocation of TGF-β1from the endoplasmic reticulum [15] and be necessary for the removal of the signal peptide that remains attached to the hepatitis C virus core protein after the initial proteolytic processing of the polyprotein [27]. Meanwhile, HM13 was identified to increase the levels of cytokines in EGFRvIII secretion profiles and facilitate tumor progression and development of Glioblastoma in vivo and vitro while EGFRvIII was the most common mutation [18, 28, 29]. Interestingly, the focus of our research was on elucidating the mechanism by which down-regulated HM13 might impact the malignancy of HCC cells.

As can be observed from our findings, HM13 was considerably elevated and connected to poor prognosis and survival rate in patients with HCC, regardless of whether the data came from public databases or clinical samples obtained, which was consistent with the findings in GBM patients [18]. It is worth noting that in normal tissues, we only detected 3 cases of high expression of HM13, and the rest were extremely low, and even no expression of HM13 was detected. To better highlight the high expression level of HM13 in HCC, we selected normal tissues that do not express HM13 for display. In addition, we reconfirmed our results with public databases including Kaplan-Meier Plotter and the Human Protein Atlas database. We found that patients with increased HM13 expression had poor overall survival, whether by Kaplan-Meier Plotter or Kaplan-Meier Plotter (Supplementary Fig. S1). Subsequently, we explored the role of HM13 methylation in HCC. According to the cBioPortal, we found that although HM13 methylation was not associated with HM13 gene expression, patients with increased HM13 CNV had higher levels of HM13 expression. According to the UALCAN, the promoter methylation of HM13 was significantly higher in HCC tissues than in surrounding normal tissues. According to the MethSurv online tool, patients with reduced HM13 methylation had better overall survival. These indicate that HM13 methylation is also involved in the occurrence and progression of HCC, which deserves more in-depth research in future work. Finally, we aimed to investigate whether HM13 could influence the malignant growth and metastasis potential of HCC cells after selecting the lentivirus vector encoding HM13-shRNA to manipulate HM13 expression. The results obtained indicated that down-regulated HM13 could inhibit the proliferative and metastatic potential of HCC cells, further supporting the notion of the enhancement role played by HM13 in HCC.

Importantly, our study expounded on the relationship between HM13 expression and TME in HCC. TIICs have been reported to play a key regulatory role in tumorigenesis and the progression of HCC [30,31,32]. We found a significant difference in the proportion of some TIICs between the HM13-high and HM13-low subgroups in TCGA. Furthermore, we found that the TMB, expression of PD1 and CTLA4 were significantly and positively correlated with the expression of HM13, and their levels were significantly higher in the HM13-high subgroup than in the HM13-low subgroup. We also found that patients with high HM13 expression might be more sensitive to immunotherapy in the IMvigor210 cohort. The above data suggest that HM13 expression can accurately predict the immune profile of HCC, and patients with high HM13 expression may be more responsive to ICI therapy.

Last but not least, there are certain limits to our research. First, after screening public databases, the expression levels and survival prediction values of HM13 were validated in a clinical sample cohort; however, the sample size was small, and additional verification in randomized controlled trials was required. Furthermore, in public databases, the expression of PD1 and CTLA4 was utilized to predict patient response to ICI treatment; thus, HCC patients treated with ICIs were required to evaluate the true predictive usefulness of HM13.


In conclusion, our study revealed for the first time that upregulation of HM13 was significantly associated with poor survival outcomes in HCC patients and promoted cell proliferation and migration capacity. In addition, HM13 may play an important role in the regulation of TIICs and immune checkpoint gene expression in the immune microenvironment of HCC. HM13 has the potential to be a predictive biomarker as well as a molecular therapeutic target for HCC.

Availability of data and materials

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



Histocompatibility minor 13


Tumor-infiltrating immune cells


Immunohistochemistry staining assay


Tumor microenvironment


Hepatocellular carcinoma


Non-alcoholic fatty liver disease


Gene Expression Profiling Interactive Analysis


Human protein atlas


Tumor mutation burden


Immune checkpoint inhibitors


Fetal bovine serum


Cell counting kit-8 assay


Copy number variation


  1. Harkus U, Wankell M, Palamuthusingam P, McFarlane C, Hebbard L. Immune checkpoint inhibitors in HCC: cellular, molecular and systemic data. Semin Cancer Biol. 2022.

  2. Shen H, Wang Z, Ren S, Wang W, Duan L, Zhu D, et al. Prognostic biomarker MITD1 and its correlation with immune infiltrates in hepatocellular carcinoma (HCC). Int Immunopharmacol. 2020;81:106222.

    Article  CAS  Google Scholar 

  3. Bravi F, Tavani A, Bosetti C, Boffetta P, La Vecchia C. Coffee and the risk of hepatocellular carcinoma and chronic liver disease: a systematic review and meta-analysis of prospective studies. Eur J Cancer Prev. 2017;26(5):368–77.

    Article  CAS  Google Scholar 

  4. Ioannou GN, Bryson CL, Weiss NS, Miller R, Scott JD, Boyko EJ. The prevalence of cirrhosis and hepatocellular carcinoma in patients with human immunodeficiency virus infection. Hepatology. 2013;57(1):249–57.

    Article  Google Scholar 

  5. Schlesinger S, Aleksandrova K, Pischon T, Jenab M, Fedirko V, Trepo E, et al. Diabetes mellitus, insulin treatment, diabetes duration, and risk of biliary tract cancer and hepatocellular carcinoma in a European cohort. Ann Oncol. 2013;24(9):2449–55.

    Article  CAS  Google Scholar 

  6. Tsilidis KK, Kasimis JC, Lopez DS, Ntzani EE, Ioannidis JP. Type 2 diabetes and cancer: umbrella review of meta-analyses of observational studies. BMJ. 2015;350:g7607.

    Article  Google Scholar 

  7. Calle EE, Rodriguez C, Walker-Thurmond K, Thun MJ. Overweight, obesity, and mortality from cancer in a prospectively studied cohort of U.S. adults. N Engl J Med. 2003;348(17):1625–38.

    Article  Google Scholar 

  8. Foerster F, Gairing SJ, Müller L, Galle PR. NAFLD-driven HCC: safety and efficacy of current and emerging treatment options. J Hepatol. 2022;76(2):446–57.

    Article  CAS  Google Scholar 

  9. Li WQ, Park Y, McGlynn KA, Hollenbeck AR, Taylor PR, Goldstein AM, et al. Index-based dietary patterns and risk of incident hepatocellular carcinoma and mortality from chronic liver disease in a prospective study. Hepatology. 2014;60(2):588–97.

    Article  CAS  Google Scholar 

  10. Forner A, Reig M, Bruix J. Hepatocellular carcinoma. Lancet. 2018;391(10127):1301–14.

    Article  Google Scholar 

  11. Hollebecque A, Malka D, Ferté C, Ducreux M, Boige V. Systemic treatment of advanced hepatocellular carcinoma: from disillusions to new horizons. Eur J Cancer. 2015;51(3):327–39.

    Article  Google Scholar 

  12. Yun UJ, Bae SJ, Song YR, Kim YW. A Critical YAP in Malignancy of HCC Is Regulated by Evodiamine. Int J Mol Sci. 2022;23(3).

  13. Fujiwara N, Friedman SL, Goossens N, Hoshida Y. Risk factors and prevention of hepatocellular carcinoma in the era of precision medicine. J Hepatol. 2018;68(3):526–49.

    Article  Google Scholar 

  14. Yin L, He N, Chen C, Zhang N, Lin Y, Xia Q. Identification of novel blood-based HCC-specific diagnostic biomarkers for human hepatocellular carcinoma. Artif Cells Nanomed Biotechnol. 2019;47(1):1908–16.

    Article  CAS  Google Scholar 

  15. Loureiro J, Lilley BN, Spooner E, Noriega V, Tortorella D, Ploegh HL. Signal peptide peptidase is required for dislocation from the endoplasmic reticulum. Nature. 2006;441(7095):894–7.

    Article  CAS  Google Scholar 

  16. Lemberg MK, Bland FA, Weihofen A, Braud VM, Martoglio B. Intramembrane proteolysis of signal peptides: an essential step in the generation of HLA-E epitopes. J Immunol. 2001;167(11):6441–6.

    Article  CAS  Google Scholar 

  17. Fang H, Sheng S, Chen B, Wang J, Mao D, Han Y, et al. A Pan-Cancer analysis of the oncogenic role of cell division cycle-associated protein 4 (CDCA4) in human tumors. Front Immunol. 2022;13:826337.

    Article  CAS  Google Scholar 

  18. Wei JW, Cai JQ, Fang C, Tan YL, Huang K, Yang C, et al. Signal peptide peptidase, encoded by HM13, contributes to tumor progression by affecting EGFRvIII secretion profiles in glioblastoma. CNS Neurosci Ther. 2017;23(3):257–65.

    Article  CAS  Google Scholar 

  19. van Es JH, van Gijn ME, Riccio O, van den Born M, Vooijs M, Begthel H, et al. Notch/gamma-secretase inhibition turns proliferative cells in intestinal crypts and adenomas into goblet cells. Nature. 2005;435(7044):959–63.

    Article  Google Scholar 

  20. 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.

    Article  CAS  Google Scholar 

  21. Su L, Zhang G, Kong X. Prognostic significance of pregnancy zone protein and its correlation with immune infiltrates in hepatocellular carcinoma. Cancer Manag Res. 2020;12:9883–91.

    Article  CAS  Google Scholar 

  22. Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med. 2015;21(8):938–45.

    Article  CAS  Google Scholar 

  23. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.

    Article  CAS  Google Scholar 

  24. Li Y, Bai W, Zhang J. MiR-200c-5p suppresses proliferation and metastasis of human hepatocellular carcinoma (HCC) via suppressing MAD2L1. Biomed Pharmacother. 2017;92:1038–44.

    Article  CAS  Google Scholar 

  25. Bruix J, Reig M, Sherman M. Evidence-based diagnosis, staging, and treatment of patients with hepatocellular carcinoma. Gastroenterology. 2016;150(4):835–53.

    Article  Google Scholar 

  26. Llovet JM, Villanueva A, Lachenmayer A, Finn RS. Advances in targeted therapies for hepatocellular carcinoma in the genomic era. Nat Rev Clin Oncol. 2015;12(8):436.

    Article  Google Scholar 

  27. McLauchlan J, Lemberg MK, Hope G, Martoglio B. Intramembrane proteolysis promotes trafficking of hepatitis C virus core protein to lipid droplets. EMBO J. 2002;21(15):3980–8.

    Article  CAS  Google Scholar 

  28. Zhou J, Cheng T, Li X, Hu J, Li E, Ding M, et al. Epigenetic imprinting alterations as effective diagnostic biomarkers for early-stage lung cancer and small pulmonary nodules. Clin Epigenetics. 2021;13(1):220.

    Article  CAS  Google Scholar 

  29. Ji P, Wang H, Cheng Y, Liang S. Prognostic prediction and gene regulation network of EIF2S2 in hepatocellular carcinoma based on data mining. J Gastrointest Oncol. 2021;12(6):3061–78.

    Article  Google Scholar 

  30. Ringelhan M, Pfister D, O'Connor T, Pikarsky E, Heikenwalder M. The immunology of hepatocellular carcinoma. Nat Immunol. 2018;19(3):222–32.

    Article  CAS  Google Scholar 

  31. Wei J, Liu Y, Zhao C. Integrated analysis of FAM57A expression and its potential roles in hepatocellular carcinoma. Front Oncol. 2021;11:719973.

    Article  Google Scholar 

  32. Zhou P, Lu Y, Zhang Y, Wang L. Construction of an immune-related six-lncRNA signature to predict the outcomes, immune cell infiltration, and immunotherapy response in patients with hepatocellular carcinoma. Front Oncol. 2021;11:661758.

    Article  Google Scholar 

Download references


Not applicable.


No funding.

Author information

Authors and Affiliations



LHC designed the study and critically reviewed the manuscript. ZGH downloaded and analyzed the data. LXP wrote the manuscript. YQK completed cellular and animal experiments. The final manuscript has been approved by all authors.

Corresponding author

Correspondence to Hongchun Liu.

Ethics declarations

Ethics approval and consent to participate

This study was supported by the Ethics Committees of Zhengzhou University. Written informed consent was obtained from all patients. All methods were performed following the relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

All authors declare no conflict of interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Additional file 2: Fig. S1.

Increased HM13 expression had poor overall survival, whether in Kaplan-Meier Plotter (A) or HPA (B) database.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, G., Lv, X., Yang, Q. et al. Identification of HM13 as a prognostic indicator and a predictive biomarker for immunotherapy in hepatocellular carcinoma. BMC Cancer 22, 888 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: