SAMD13 serves as a useful prognostic biomarker for hepatocellular carcinoma
European Journal of Medical Research volume 28, Article number: 514 (2023)
Hepatocellular carcinoma (HCC) is the most common form of liver cancer and the 5-year relative overall survival (OS) rate is less than 20%. Since there are no specific symptoms, most patients with HCC are diagnosed in an advanced stage with poor prognosis. Therefore, identifying novel prognostic biomarkers to improve the survival of patients with HCC is urgently needed. In the present study, we attempted to identify SAMD13 (Sterile Alpha Motif Domain-Containing Protein 13) as a novel biomarker associated with the prognosis of HCC using various bioinformatics tools. SAMD13 was found to be highly expressed pan-cancer; however, the SAMD13 expression was significantly correlated with the worst prognosis in HCC. Clinicopathological analysis revealed that SAMD13 upregulation was significantly associated with advanced HCC stage and high-grade tumor type. Simultaneously, high SAMD13 expression resulted in association with various immune markers in the immune cell subsets by TIMER databases and efficacy of immunotherapy. Methylation analysis showed SAMD13 was remarkably associated with prognosis. Furthermore, a six-hub gene signature associated with poor prognosis was correlated with the cell cycle, transcription, and epigenetic regulation and this analysis may support the connection between SAMD13 expression and drug-resistance. Our study illustrated the characteristics of SAMD13 role in patients with HCC using various bioinformatics tools and highlights its potential role as a therapeutic target and promising biomarker for prognosis in HCC.
Hepatocellular carcinoma (HCC) has one of the highest morbidity and mortality globally and the incidence of HCC is significantly increased annually. HCC ranked second and third in major causes of cancer-related deaths due to the complexity, heterogeneity, and high recurrence following surgical resection [1, 2]. In patients with early stage HCC, surgical or percutaneous resection is recommended as a first-line treatment option; however, the recurrence rate is around 50% 5-year post-surgery [3, 4]. The prognostic markers play an important role in monitoring tumor recurrence or metastasis and judging prognosis . Therefore, identification of specific prognostic biomarker of HCC helps to provide a great opportunity to improve the prognosis in patients.
SAMD13 gene is cytogenetically located on the short (p) arm of chromosome 1, at position 31.1 (1p31.1) and one of the sterile alpha motif domain-containing (SAMD)-proteins. The SAMD is the most common protein interaction module which has been implicated in diverse biological roles for cellular processes, such as binding to self- or other SAMD, non-SAMD, RNA, DNA, and even lipids [6,7,8,9,10,11,12]. A previous study revealed that SAMD13 was downregulated on the micropapillary area in invasive micropapillary carcinoma (IMPC) . This study suggested that SAMD13 expression is relevant to the pathological features and immune microenvironment in IMPC; however, the expression profile and biological role of SAMD13 in hepatocellular carcinoma remain totally unknown.
This study aimed to explore the expression profile and prognostic effects of SAMD13 using pan-cancer database and various bioinformatics in patients with HCC. Furthermore, we examined the relationship between SAMD13 expression and clinicopathological criteria of The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA–LIHC) using UALCAN. Subsequently, we performed methylation profiling analyses to predict the overall survival rate of patients with HCC. The correlation between SAMD13 and tumor-infiltrating immune cells in the tumor microenvironment of HCC was also determined. Networks and functional enrichment analyses were further performed to evaluate whether SAMD13 expression is associated with transcriptional and epigenetic regulation and its related genes in clinical features. Finally, based on the functional enrichment analyses, the importance of SAMD13 gene for drug resistance and its value as targets for novel drug development using GEO data set was validated. Taken together, our results provided a better understanding of clinical significance of SAMD13 in HCC, and prognostic assessment implication in patients with HCC.
Materials and methods
Data acquisition and processing
UALCAN database is an open-access web for analyzing cancer OMICS data based on tumors and normal samples from TCGA, MET500, CPTAC, and CBTTC . Here, the UALCAN was used to compare the expression levels of SAMD13 in pan-cancer, including HCC and clinical data, as well as the associations between them. GEPIA and GEPIA2 is a web tool for gene expression profiling and interactive analyses based on TCGA and The Genotype-Tissue Expression (GTEx) data . This database was used to access overall survival (OS) and disease-free survival (DFS) based on the expression of SAMD13, and the further association between this gene and the expression of immune-related marker genes was verified.
The TCGA–LIHC (HCC) data set was used as the training set, while three other independent data sets (GSE22058, GSE25097, and GSE45436) and Gene Expression database of Normal and Tumor tissues 2 (GENT2) (which is web-based tool from public gene expression data sets) were used as validation sets. Gene expression profiling data sets, including GSE22058 (97 non-tumor and 100 tumor), GSE25097 (6 normal, 40 cirrhosis, and 268 tumor), and GSE45436 (39 non-tumor and 95 tumor), and chemo-resistant gene profiling data sets, including GSE54175 (cisplatin- or doxorubicin-resistant subline), GSE93595 (anti-angiogenic therapy (JNJ-28841072)-resistant subline), GSE121153 (sorafenib-resistant subline), GSE125180 (doxorubicin-resistant subline), GSE109211 (clinical parameters from responders and non-responders to sorafenib in patients with HCC) were downloaded from the gene expression omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/). cBioportal for Cancer Genomics database which analyzes multi-omics data from The Cancer Genome Atlas was used for analyzing genetic alteration for SAMD13 gene in a liver study (TCGA, PanCancer Atlas) . The KM–plotter database was used for survival curves in various cancers . Hazard ratios (HR) and p values (from the log-rank test) were calculated online.
Characterization of tumor microenvironment in HCC
Human Liver Browser and Single-cell Atlas in Liver Cancer (scAtlasLC) were used to characterize SAMD13 expression of single-cell transcriptomic profiles in HCC [18, 19]. The Tumor Immune Estimation Resource database (TIMER 2.0) used analysis and visualization of immune infiltrates from TCGA in TIMER Database . In addition to exploring the associations of HCC-infiltrating immune cells with SAMD13 gene expression, the degree of immune infiltration of B cell, CD8 + T cell, CD4 + T cell, macrophages, neutrophil, and dendritic cell (DC) was analyzed in TIMER database. The strength of correlations evaluated by the purity-adjusted partial Spearman’s rho value and estimated statistical significance in TIMER 2.0.
Comprehensive analysis of SAMD13 methylation
MethSurv database could provide the initial assessment of DNA methylation biomarkers using TCGA . In this study, single CpG methylation of SAMD13 and Heatmap analysis were performed in patients with HCC. Next, Shiny Methylation Analysis Resource Tool (SMART)  was used to analyze differential methylation by each SAMD13 probe and Spearman’s correlation between methylation level and mRNA level. The CpG-aggregated methylation values were measured by mean (β values). To characterize the methylation patterns, the significant CpGs were classified according to their functional roles in genomic locations, such as promoters within 1,500 bps of a transcription start site (TSS) (TSS1500); within 200 bps of a TSS (TSS200); 5’ untranslated regions (5'UTR); first exon (1stExon); body (non-promoter); 3'UTR (non-promoter).
Functional and pathway enrichment analysis
The 22 genes with the strongest correlation with SAMD13 were selected using Pathway Commons  and this allowed the creation of a protein–protein interaction (PPI) network for SAMD13 genes as well as binding and targeted genes. To determine the functional and pathway meaning, cellular component (CC), molecular function (MF), biological process (BP), and reactome pathway (RP) for 22 selected genes forming a cluster with SAMD13 was performed using geneontology.org . Functional enrichment analyses of Gene Ontology (GO) terms were performed using Fisher's Exact test and adjusted by a False Discovery Rate (FDR) correction for multiple testing in pantherdb.org (http://pantherdb.org/).
SAMD13 Immunohistochemistry staining
Cell lines and cell culture
Huh7, HepG2, SK-Hep1 cells were purchased from the American Type Culture Collection. SNU475 and SNU449 cells were Korean Cell Line Bank (Seoul, Republic of Korea). Huh7, HepG2, SK-Hep1 Cells were maintained in Dulbecco’s modified Eagle’s medium containing 10% FBS (GIBCO, Grand Island, NY, USA).SNU475 and SNU449 cells were maintained in RPMI1640 medium containing 10% FBS in a humidified incubator with 5% CO2 at 37 °C. HepaRG cells were cultured in William’s E medium supplemented with 10% of FBS, 5 µg/mL insulin, 2‐mM Glutamax, 1% penicillin–streptomycin, and 50 µM hydrocortisone hemisuccinate (Sigma-Aldrich).
Real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR) analysis
Total RNA was isolated using a TRIzol reagent-based kit (Intron Biotech, Seongnam-Si, Republic of Korea). Reverse transcription was performed using the SuperScript IV First-Strand Synthesis System for RT-PCR (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s protocol. cDNA was amplified with specific primers and SYBR Premix Ex Taq (Takara Bio, Otsu, Shiga, Japan and Agilent Technologies, Santa Clara, CA, USA). The SAMD13 mRNA levels were normalized to the amount of GAPDH mRNA. qPCR was performed according to the manufacturer’s instructions (Applied Biosystems, Foster City, CA, USA, and Agilent Technologies). Relative quantification of gene expression was performed using the 2-ΔΔCT method. The primer pair used for SAMD13 (forward: 5’-CAAGGAAAATGGCTCTGTCGGTG-3’; reverse: 5’-AGCTTGCTCCTCAAATCCCACG-3’) and GAPDH (forward: 5’-AACGGGAAGCTTGTCATCAATGGAAA-3’; reverse: 5’- GCATCAGCAGAGGGGGCAGAG-3’).
The statistical analysis was calculated automatically based on the online database above. Student's t test implemented by GraphPad Prism (Ver7).Correlations were analyzed by Spearman and Pearson’s correlation. p < 0.05 was considered statistically significant.
Impact of expression and prognosis of SAMD13 in various human cancers
To explore the expression pattern of SAMD13between tumor and normal tissues in various types of cancer, we examined the SAMD13 differences using the UALCAN database. The comparison of expression level between each type of normal and cancer cells revealed that SAMD13 increased in breast invasive carcinoma, cholangiocarcinoma, glioblastoma multiforme, kidney renal papillary cell carcinoma, liver hepatocellular carcinoma, lung adenocarcinoma, prostate adenocarcinoma, and stomach adenocarcinoma, while decreased in colon adenocarcinoma, head and neck squamous cell carcinoma, kidney chromophobe, rectum adenocarcinoma, thyroid carcinoma, and uterine corpus endometrial carcinoma (Fig. 1A). To validate these findings, we further examined the expression levels of SAMD13 using the GENT2 portal. As shown in Additional file 1: Figure S1, SAMD13 was upregulated in blood, brain, liver, and bone cancers, while downregulated in colon, breast, kidney, head and neck, tongue, thyroid, adrenal gland, pharynx, testis, and endometrial cancers. To assess the prognostic significance of SAMD13 expression in various types of cancers, log-rank test was performed using pan-cancer RNA-seq in KM-plotter. Among these 21 types of cancer, four were significantly associated with the prognosis of all cancers (Fig. 1B). High mRNA expression level of SAMD13 were observed to be significantly associated with worse overall survival in lung adenocarcinoma (HR, 1.51; 95% CI 1.15–2.07; p = 0.004) and hepatocellular carcinoma (HR, 2.15; 95% CI 1.15–3.06; p < 0.001), whereas thymoma (HR, 0.23; 95% CI 0.05–1.12; p = 0.048) and kidney renal clear cell carcinoma (HR, 0.011; 95% CI 0.49–0.91; p = 0.011) were significantly associated better prognosis. In addition, these observations were validated by GEPIA2 in pan-cancer database. As shown in Fig. 1C, high SAMD13 expression is associated with a poorer overall survival LIHC and LUAD, while low expression of SAMD13 indicated better prognosis in THYM. When analyzing disease-free survival, the high expression of SAMD13 is associated with a worse prognosis in LGG, LICH, and LUAD, while the low expression of SAMD13 indicated better prognosis in COAD, GBM, KIRC, and THCA.
Validation and survival analysis of SAMD13 in HCC
Since SAMD13 gene expression was not only upregulated in various carcinomas but also linked to worst prognosis in HCC, we further performed validation on additional three independent GEO data sets. As expected, these GEO data sets also revealed the significant augmentation of SAMD13 mRNA expression in the tumor group compared to their normal/non-tumor group. Interestingly, increased mRNA expression of SAMD13 was also revealed in the cirrhosis group at GSE25097 (Fig. 2A). Next, this trend was evaluated at the protein level between normal liver and HCC tissues, and as a result, the immunohistochemical data show that SAMD13 expression was higher than that compared to normal tissue (Fig. 2B).To further validate the differential expression of SAMD13, mRNA expression of SAMD13 was conducted on 1 immortalized hepatic cell line (HepaRG), 2 well-differentiated HCC cell lines (Huh7 andHepG2), and 3 poorly differentiated HCC cell lines (SK-Hep-1, SNU475, and SNU449). As shown Fig. 3C, the SAMD13 expression were increased in poorly differentiated HCC cells relative to well-differentiated HCC cells suggesting that SAMD13 may be acted HCC tumors to more aggressive tumor phenotypes. To better understand the relationship between SAMD13 expression and prognosis in HCC, we investigated the correlation between SAMD13 expression and overall survival (OS), disease specific survival (DSS), disease free interval (DFI), and progression free interval (PFI). As shown in Fig. 2D, high expression of SAMD13 was significantly associated with shorter OS (log-rank test = 11.6; p = 0.0006605), DSS (log-rank test = 7.32; p = 0.006818), and PFI (log-rank test = 6.127; p = 0.01331); meanwhile, DFI was correlated with shorter tendency (log-rank test = 3.343; p = 0.0675).
Clinicopathological features and SAMD13 expression
We then determined a relationship between clinicopathological features and SAMD13 expression using TCGA data through UALCAN in HCC, and summarized in Table1. According to tumor stage, SAMD13 expression were significantly upregulated in HCC group classified as stage I to stage III compared to the corresponding normal group. Based on tumor grade, a significant increase in SAMD13 levels was observed in the HCC group in all tumor grades. In addition, upregulation of SAMD13 expression was observed in patients with HCC and fibrolamellarcarcinoma types. In addition, SAMD13 expression was statistically significant in all other categories of patients with HCC, including gender, age, and race. These results suggest that SAMD13 was significantly correlated with clinical parameters in patients with HCC.
Association of SAMD13 with expression profiling of tumor microenvironment (TME) in HCC
To further assess the effect of SAMD13 on the TME, we investigated the association between SAMD13 and expression profiling of TME in HCC using the Human Liver Browser and scAtlasLC data sets. It was revealed that SAMD13 was mainly expressed in scar-associated macrophages (SAMs), T cells, tumor-associated macrophages (TAMs), tumor-associated endothelial cells (TECs), and carcinoma cells, while moderately expressed in hepatocytes and cholangiocytes (Fig. 3A, B). Then, we further analyzed between SAMD13 and immune infiltrates in HCC using TIMER database. The results showed that the expression of SAMD13 was significantly correlated with B cell (r = 0.159, p = 3.19e − 03),CD8 + T cells (r = 0.106, p = 5.03e-02), CD4 + T cells (r = 0.282, p = 1.02e − 07), macrophage (r = 0.291, p = 4.32e-08), neutrophils (r = 0.148, p = 5.93e − 03), and dendritic cells (r = 0.249, p = 3.32e-06)in HCC (Fig. 3C). To study the relationship between SAMD13 and infiltrated immune cells, TIMER and GEPIA database were use. Immune gene markers of B cell, T cell (general), CD8 + cell, CD4 + cell, tumor-associated macrophage (TAM), M1 and M2 macrophage, neutrophil, natural killer cell (NK) cell, DC, Th1/2/fh/17cell, Treg, and T cell exhaustion in HCC were assessed. The results revealed that SAMD13 was a significant positive relationship including all immune infiltrates, including B cell (CD17,CD79A), T cell general (CD3D, CD3E, CD2), CD8 + T cell (CD8A, CD8B), CD4 + T cell (CD4), TAM (CD68, IL10), M1 macrophage (IRF5, PTGS2), M2 macrophage (CD163, VSIG4, MS4A4A), neutrophil (ITGAM,CCR7), DC (HLA-DPB1, HLA-DRA, HLA-DPA1, CD1C, NRP1, ITGAX), Th1 (STAT4, STAT1, IFNG, TNF), Th2 (GATA3, STAT6, STAT5A), Th17 (STAT3), Treg (CCR8, TGFB1, STAT5B), and T cell exhaustion (PDCD1, CTLA4, HAVCR2). However, NK cell and Tfh was not significantly associated between SAMD13 and immune infiltrates (Table 2). In addition, the comprehensive prognostic value analysis was performed on SAMD13 expression and immune cell infiltration in HCC. The results of Fig. 3D showed that low expression of SAMD13 along with low immune infiltration of neutrophil, macrophage, macrophage M0, macrophage M2, and myeloid-derived suppressor cells (MDSC) had better prognosis than the high expression of SAMD13 group in HCC. However, there were no significantly changed inCD8 + T cells, CD4 + T cells, B cell, DC, and NK cells (data not shown). These results indicated that SAMD13 is not only closely related to immune cell infiltration, but also correlated with prognosis in HCC.
Association between SAMD13 profiles and clinical characteristics of patients with HCC
To prove into the role of SAMD13 methylation in HCC, we studied the correlation between SAMD13 methylation state and its expression levels. Based on SMART database, SAMD13 methylation level is significantly lower in HCC tissues compared to normal liver tissues (Fig. 4A) and specific methylation sites were presented in heatmaps using MethSurv database (Fig. 4B). According to the heatmap, there were 15 CpG sites in SAMD13 and we found that average methylation of all CpG sitesof SAMD13 (Overall; Aggregation, p = 2.4e-05), including N_Shore (Aggregation, p = 0.00048), S_Shore(Aggregation, p = 0.058), S_Shelf (Aggregation, p = 0.039), and Open_Sea (Aggregation, p = 6.3e-14) was significantly lower in tumor tissues than in the normal counterpart (Fig. 4C). At the same time, we performed survival analyses using methylation profiling of SAMD13 CpGs in HCC. The results indicated that CpGs of SAMD13 gene containing cg15103960 (HR, 0.71; 95% CI, 0.51-1; p = 0.052), cg02041547 (HR, 0.46; 95% CI, 0.32-0.64; p < 0.001), and cg23086720 (HR, 0.51; 95% CI, 0.36-0.72; p < 0.001) was associated with better prognosis, while cg23694882 (HR, 1.52; 95% CI, 1.08-2.14; p = 0.016), cg15089272 (HR, 1.5; 95% CI, 1.07-2.12; p = 0.019), and cg23925111 (HR, 1.62; 95% CI, 1.15-2.27; p = 0.006) showed opposite (Table 3). We also found that cg20426713, cg11555919, cg23694882, cg13305246, cg22529952, and cg22529952 probes were significantly lower in tumor tissues than in the normal counterpart, while cg23086720 probe showed higher in tumor compared to normal control (Additional file 1: Figure S2). Given the close relationship between the methylation of SAMD13 and SAMD13 expression, we found a positive correlation between the methylation level and the expression of SAMD13 (Overall; Aggregation, p = 1.1e-05), including N_Shore (Aggregation, p = 0.012), S_Shore (Aggregation, p = 0.00023), and Open_Sea (Aggregation, p = 2.8e-09) (Fig. 4D and Additional file 1: Figure S3). Moreover, we examined the copy number alterations and revealed that SAMD13 mRNA had higher copy number gain compare to diploid or shallow deletion (Fig. 5A). Interestingly, later tumor stages including T3, T3A, T3B, and T4 showed a relatively high tendency for deletion or gain of copy number alterations, although absolute copy number alterations were decreased (Fig. 5B). Since the relationship between copy number amplification and individual DNA methylation patterns affects gene regulation, we next explored the relationship between methylation value and copy number amplification in SAMD13. Interestingly, we found that the average methylation value of CpGs of SAMD13 was significantly negative correlation with copy number amplification in HCC (Aggregation, p = 2.3e-05) (Fig. 5C). In methylation value on clinical relevance, we found that cg23694882 was much hyper-methylated in stage-dependent manner (p = 0.0022), while cg23086720 showed hypo-methylated in stage dependent manner (p = 0.00073). These results indicated that methylation value and copy number alterations of SAMD13 gene plays a crucial role in the prognosis that might reflect the complexity/heterogeneity with HCC patients in SAMD13.
Networks analyses and functional enrichment analyses of SAMD13 in HCC
To further investigate the biological role in SAMD13, genes were found to be interacted with SAMD13 based on Pathway Commons and STRING (Fig. 6A, B), and list of genes was presented in Table 4 and 5, respectively. The functions of SAMD13 and interacted genes were classified into each functional group: molecular function, biological process, reactome pathway, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. Among functional and pathway enrichment analysis, SAMD13 was significantly associated with DNA binding on molecular function, including “chromatin binding”, “histone binding”, “methyl-CpG binding”, “general transcription initiation factor binding”; regulation of nucleic acid metabolic process on biological process, including “metaphase/anaphase transition of mitotic cell cycle”, “histone ubiquitination and/or methylation”, “regulation of transcription involved in G1/S transition of mitotic cell cycle”, “hepatocyte apoptotic process”; cell cycle regulation on reactome pathway, including “inactivation of APC/C via direct inhibition of the APC/C complex”, “aberrant regulation of mitotic exit in cancer due to RB1 defects”, “G1 phase”, “oncogene induced senescence”; proteolysis on KEGG pathway, including “ubiquitin mediated proteolysis” and “cell cycle” (Additional file 2: Table S1 and Additional file 3: Table S2). To further explore the potential network co-expressed with SAMD13 gene in HCC, the GEPIA2 database was employed to evaluate the relationship between SAMD13 and co-expressed genes. In the TCGA-HCC dataset, 13 genes were identified with significant correlation with SAMD13 from the tumor group only and a total of six genes including FOXM1, JUN, JARID2, BRE, BUB1B, and PHC2 were shown to significantly predict poor overall survival in log-rank tests in HCC (Fig. 6C). According to the high/low expression groups among these genes combined with SAMD13, we divided the TCGA-HCC samples into four combinations based on the median of gene expression for all six pairs of genes. In survival analysis of the six pairs of gene, the group with high expression of SAMD13 and FOXM1, JUN, JARID2, BRE, BUB1B, or PHC2 had the poorest prognosis, while the groups of low expression of SAMD13 with low expression of all six genes had the best prognosis, respectively (Fig. 6D).
Association of SAMD13 expression with the acquired drug resistance and immune checkpoint blockade (ICB) genes
Since the biological role of SAMD13 enriches the “cell cycle” and “nucleic acid metabolic process”, we further examined the chemosensitivity to cisplatin, doxorubicin, sorafenib, and JNJ-28841072 which have been shown to potently inhibit the cell cycle used for the treatment of HCC from three independent GEO data set, including GSE54175 (identification of chemo-resistant genes in human metastatic HCC), GSE121153 (sorafenib-resistant HCC), GSE125180 (doxorubicin-resistant HCC), and GSE93595 (anti-angiogenic drug resistant HCC). As shown Fig. 7A, SAMD13 mRNA level was elevated in almost all cell lines with conventional chemotherapy, based on analyses from the GEO databases. Typically, SAMD13 was positively correlated with the elevation of gene expression for cisplatin-resistant MHCC97L cells (GSE54175), sorafenib-resistant Huh7 cells (GSE121153), doxorubicin-resistant 834 cells (GSE125180), and JNJ-28841072-resistant Huh7 cells which had acquired resistance under long-term anti-angiogenic therapy from serial transplantation in immunocompromised mice (GSE93595). To further verify the clinical significance of SAMD13 expression for assessing response of chemotherapy, drug-cancer response data set was investigated. One cohort (GSE109211) treated with placebo or sorafenib was retrieved, which contained 140 patients randomized to receive 98 non-responders and 42 responders. Compared to placebo- or sorafenib-treatment, there was no significant differences between responders and non-responders (Fig. 7B), as well as recurrence of HCC (Additional file 1: Figure S4). Since the role of SAMD13 in ICB treatment of HCC has not been reported, we further assessed the TCGA–LIHC data to explore the correlation between SAMD13 and some ICB related genes. The results demonstrated that SAMD13 was significantly positive correlation with ICB genes, including CD276, CD86, ICOS, HAVCR2, LAIR1, CD44, TIGIT, CD80, CTLA4 (Fig. 7C). Collectively, these results suggest that SAMD13 expression might not be involved in chemotherapy response to cancer treatment, but rather a target for ICB treatment that can counter acquired chemical resistance.
HCC accounts for more than 90% of primary liver tumors and the 5-year survival rate of HCC is only 18% [26, 27]. The main reason for rapid progression and high mortality of HCC leads to treatment failure caused by resistance to conventional chemotherapy, recurrence, and metastasis in patients with HCC, and eventually shows a worse prognosis. Hence, the discovery of promising biomarkers for clinical diagnosis and treatment of patients with HCC is urgently needed.
In this study, we evaluated the expression, methylation, drug-resistance, and prognostic significance of SAMD13 as a new prognostic biomarker in HCC based on bioinformatics approaches. Our result revealed that SAMD13 was significantly overexpressed in TCGA including HCC and we validated using three independent GEO data sets (GSE22058, GSE25097, and GSE45436). Although SAMD13 is highly expressed in various types of cancer, it was significantly associated with shorter OS, DSS, DFI, and PFI in patients with HCC, even poorly differentiated HCC cell lines which tend to be more aggressive and higher grade. SAMD13 could be attributed to various immune infiltration events, such as B cells, T cells, macrophages, neutrophil, and DCs. In analysis of methylation profiling, CpGs of SAMD13 gene containing cg15103960, cg02041547, cg23694882, cg23086720, cg15089272, and cg23925111 was associated with prognosis, and SAMD13 gene showed hypomethylation. In networks and functional enrichment analyses, six hub genes (FOXM1, JUN, JARID2, BRE, BUB1B, and PHC2) which were up-regulated in HCC tissues in comparison with normal tissues. Meanwhile, high-level of six-hub gene expressions and high-level SAMD13 expression were significantly associated with the worst OS in correlation analysis. More importantly, SAMD13 expression could be a potential for immuno-checkpoint inhibitor (ICI) treatment, although it exhibits acquired diverse anti-cancer drug resistant in HCC. Taken together, these results suggested that SAMD13 play an important role in poor clinical outcome and additional genetic and/or epigenetic alterations of SAMD13 gene promotes the development of HCC. To our knowledge, this is report to demonstrate that SAMD13 is the novel biomarker to further clarify the prognostic factors for HCC.
SAMD13 is known a protein coding gene which enables chromatin and histone binding activity involved in negative regulation of transcription. Only one study showed that SAMD13 is negatively associated with invasive micropapillary carcinoma of breast cancer, especially micropapillary area , but its clinical impact on cancers remains totally unclear, even contrary. In this study, our findings revealed that the level of SAMD13 expression was significantly up-regulated in the tumor group compared to the non-tumor group, while high level of SAMD13 expression was significantly associated with the worst prognosis in HCC. As high-grade (poorly differentiated) cancer cells tended to be more motile and aggressive than well-differentiated cancer cells , we confirmed that mRNA expression of SAMD13 in well- and poorly differentiated HCC cell lines. We also indicated that SAMD13 was positively associated with immune cell infiltration. Since immune cells and stromal cells represent typical TME, their interactions are related to the tumor invasion, recurrence, metastasis, and the effect of immunotherapy, and prognosis [29,30,31]. Besides, Pappas et al. revealed that SAMD9L plays a role in the migration of INF-beta affected T cells . In the present study, the expression levels of SAMD13 exhibited a positive correlation with diverse immune cells, including B cell, T cell, TAM, macrophage, neutrophil, and DCs. Interestingly, low SAMD13 expression was associated with low immune cell infiltration of neutrophil, macrophage, macrophage M0/M2, and MDSC, it showed better prognosis (Fig. 3D). In patients with HCC, high levels of neutrophil, macrophage (tumor-derived macrophage or macrophage M2), and MDSC are reported to be associated with poor survival [33,34,35], it is hypothesized that SAMD13 may play an important role on the TME in HCC and further study of the immune microenvironment and HCC feature should be addressed. In epigenetic analysis, methylation in SAMD13 affected prognosis and copy number alterations, especially for some specific CpG sites. Methylation of cytosine residues in CpG dinucleotide pairs could serve as not only prognostic biomarkers but also therapeutic targets [36,37,38,39,40]. Because aberrant DNA methylation contribute directly to dysregulation of gene expression, such as oncogenes or tumor suppressor genes, respectively [41,42,43]. Since newly discovered six CpGs located in genes were able to predict HCC OS, CpG methylation of SAMD13 gene could provide a new prognostic clue in HCC.
In addition, using Pathway Commons and STING, we found that thirteen genes had significantly interacted with SAMD13 and finally identified six genes positively correlated with SAMD13 expression of HCC–FOXM1, JUN, JARID2, BRE, BUB1B, and PHC2. Among these six genes, FOXM1, which is a transcription factor fork-head box M1, has been proved to play a role in tumor progression and poor prognosis in HCC [44, 45]. JUN has been reported to be overexpressed in diverse cancers including HCC and can play an oncogenic transcription factor [46, 47]. JARID2 is a DNA-binding protein which plays an essential role in transcriptional regulation during embryonic development [48, 49]. Recently some researchers reported that JARID2 can promote invasion and metastasis and was correlated with worse OS with HCC [50, 51]. In HCC, the over expressed BRE which is core component of the deubiquitin complex. BRCA1-A and BUB1B belong to core elements of the spindle assembly checkpoint (SAC) correlated with tumor progression, worse OS and DFS in HCC [52,53,54,55]. The role of PHC2, which is a transcriptional repressor of target genes by mainly modulating histone methylation, has not yet been fully understood in HCC. Of note, SAMD13 and its associated six-hub gene identified as having transcriptional and epigenetic regulation in networks analyses and functional enrichment analyses resulted in poor progression in patients with HCC. In summary, our study provides clues regarding the biological roles and molecular characterization of these genes with SAMD13, and it could act as new prognostic biomarkers for HCC.
Furthermore, we demonstrated that upregulation of SAMD13 expression could act as acquired drug resistance, as well as a potential for immune checkpoint targeted therapy. Although conventional chemotherapies including targeted therapies based on tyrosine kinase inhibitors (TKIs) and conventional chemotherapy are in clinical use, the reasons for the high resistance in HCC are not fully understood. However, the proposed mechanisms of drug resistance are associated with cell cycle regulation, DNA damage repair, drug detoxification, and apoptotic pathway, and epigenetic functions such as DNA methylation which includes certain genes that play a role [56, 57]. Since we observed high SAMD13 expression in cisplatin-, sorafenib-, doxorubin- and JNK-28841072 resistant HCC sublines in the GEO database, this suggests that SAMD13 is likely to mediate acquired drug resistance. All these indicate that the biological role of SAMD13 enriched in “regulation of nucleic acid metabolic process” and “cell cycle regulation” participates in acquired chemotherapy resistance and translational process, further supporting a role for CpG methylation of SAMD13 linked with worse prognosis in patients with HCC. We also evaluated the relationship between SAMD13 and ICB related genes in HCC. As a result, SAMD13 showed significantly positive correlations with multiple immunotherapy-predicted pathways and checkpoints, which could affect the efficacy of ICB. Therefore, the SAMD13-based prediction model might provide individualized precision treatment to patients and further study should be addressed to confirm these results.
Our bioinformatics analysis identified that SAMD13 may be critical in the development and prognosis of HCC. However, our study was performed based on bioinformatic analysis, further experimental validations should be addressed to confirm the results of this prediction in HCC. We hope present study can provide some evidence for the potential use as cancer therapeutic and prognostic biomarker of HCC from new insights.
Availability of data and materials
Bray F, Ferlay J, Soerjomataram I, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Villanueva A. Hepatocellular carcinoma. N Engl J Med. 2019;380(15):1450–62.
Liu PH, Hsu CY, Hsia CY, et al. Surgical resection versus radiofrequency ablation for single hepatocellular carcinoma </= 2 cm in a propensity score model. Ann Surg. 2016;263(3):538–45.
Ponziani FR, Bhoori S, Castelli C, et al. Hepatocellular carcinoma is associated with gut microbiota profile and inflammation in nonalcoholic fatty liver disease. Hepatology. 2019;69(1):107–20.
Zhao, K., Xu, L. & Li, F. et al. Identification of hepatocellular carcinoma prognostic markers based on 10-immune gene signature, Biosci Rep. 2020;40(8). https://doi.org/10.1042/BSR20200894.
Peterson AJ, Kyba M, Bornemann D, et al. A domain shared by the Polycomb group proteins Scm and ph mediates heterotypic and homotypic interactions. Mol Cell Biol. 1997;17(11):6683–92.
Seidel JJ, Graves BJ. An ERK2 docking site in the Pointed domain distinguishes a subset of ETS transcription factors. Genes Dev. 2002;16(1):127–37.
Oberstrass FC, Lee A, Stefl R, et al. Shape-specific recognition in the structure of the Vts1p SAM domain with RNA. Nat Struct Mol Biol. 2006;13(2):160–7.
Aviv T, Lin Z, Lau S, et al. The RNA-binding SAM domain of Smaug defines a new family of post-transcriptional regulators. Nat Struct Biol. 2003;10(8):614–21.
Li H, Fung KL, Jin DY, et al. Solution structures, dynamics, and lipid-binding of the sterile alpha-motif domain of the deleted in liver cancer 2. Proteins. 2007;67(4):1154–66.
Knight MJ, Leettola C, Gingery M, Li H, Bowie JU. A human sterile alpha motif domain polymerizome. Protein Sci. 2011;20(10):1697–706.
Schultz J, Ponting CP, Hofmann K, Bork P. SAM as a protein interaction domain involved in developmental regulation. Protein Sci. 1997;6(1):249–53.
Kanomata N, Kurebayashi J, Koike Y, Yamaguchi R, Moriya T. CD1d- and PJA2-related immune microenvironment differs between invasive breast carcinomas with and without a micropapillary feature. BMC Cancer. 2019;19(1):76.
Chandrashekar DS, Bashel B, Balasubramanya SAH, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.
Tang Z, Li C, Kang B, et al. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102.
Cerami E, Gao J, Dogrusoz U, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.
Gyorffy B, Lanczky A, Eklund AC, et al. An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients. Breast Cancer Res Treat. 2010;123(3):725–31.
Ma L, Wang L, Khatib SA, et al. Single-cell atlas of tumor cell evolution in response to therapy in hepatocellular carcinoma and intrahepatic cholangiocarcinoma. J Hepatol. 2021;75(6):1397–408.
Massalha H, Bahar Halpern K, Abu-Gazala S, et al. A single cell atlas of the human liver tumor microenvironment. Mol Syst Biol. 2020;16(12):e9682.
Li T, Fu J, Zeng Z, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509–14.
Modhukur V, Iljasenko T, Metsalu T, et al. MethSurv: a web tool to perform multivariable survival analysis using DNA methylation data. Epigenomics. 2018;10(3):277–88.
Li Y, Ge D, Lu C. The SMART App: an interactive web application for comprehensive DNA methylation analysis and visualization. Epigenetics Chromatin. 2019;12(1):71.
Rodchenkov I, Babur O, Luna A, et al. Pathway commons 2019 update: integration, analysis and exploration of pathway data. Nucleic Acids Res. 2020;48(D1):D489–97.
The Gene Ontology C. The gene ontology resource: 20 years and still going strong. Nucleic Acids Res. 2019;47(D1):D330–8.
Ponten F, Jirstrom K, Uhlen M. The human protein atlas–a tool for pathology. J Pathol. 2008;216(4):387–93.
Llovet JM, Zucman-Rossi J, Pikarsky E, et al. Hepatocellular carcinoma. Nat Rev Dis Primers. 2016;2:16018.
Kulik LM, Chokechanachaisakul A. Evaluation and management of hepatocellular carcinoma. Clin Liver Dis. 2015;19(1):23–43.
Yuzugullu H, Benhaj K, Ozturk N, et al. Canonical Wnt signaling is antagonized by noncanonical Wnt5a in hepatocellular carcinoma cells. Mol Cancer. 2009;8:90.
Ge Z, Ding S. The crosstalk between Tumor-Associated Macrophages (TAMs) and tumor cells and the corresponding targeted therapy. Front Oncol. 2020;10:590941.
Obradovic A, Chowdhury N, Haake SM, et al. Single-cell protein activity analysis identifies recurrence-associated renal tumor macrophages. Cell. 2021;184(11):2988–3005.
Guo Y, Yang J, Ren K, et al. The heterogeneity of immune cell infiltration landscape and its immunotherapeutic implications in hepatocellular carcinoma. Front Immunol. 2022;13:861525.
Pappas DJ, Coppola G, Gabatto PA, et al. Longitudinal system-based analysis of transcriptional responses to type I interferons. Physiol Genomics. 2009;38(3):362–71.
Margetts J, Ogle LF, Chan SL, et al. Neutrophils: driving progression and poor prognosis in hepatocellular carcinoma? Br J Cancer. 2018;118(2):248–57.
Arvanitakis K, Koletsa T, Mitroulis I, Germanidis G. Tumor-associated macrophages in hepatocellular carcinoma pathogenesis, prognosis and therapy. Cancers (Basel). 2022;14(1):226.
Zhang X, Fu X, Li T, Yan H. The prognostic value of myeloid derived suppressor cell level in hepatocellular carcinoma: a systematic review and meta-analysis. PLoS ONE. 2019;14(12):e0225327.
Hao X, Luo H, Krawczyk M, et al. DNA methylation markers for diagnosis and prognosis of common cancers. Proc Natl Acad Sci U S A. 2017;114(28):7414–9.
Issa JP. DNA methylation as a therapeutic target in cancer. Clin Cancer Res. 2007;13(6):1634–7.
Cheng Y, He C, Wang M, et al. Targeting epigenetic regulators for cancer therapy: mechanisms and advances in clinical trials. Signal Transduct Target Ther. 2019;4:62.
Hernandez-Vargas H, Lambert MP, Le Calvez-Kelm F, et al. Hepatocellular carcinoma displays distinct DNA methylation signatures with potential as clinical predictors. PLoS ONE. 2010;5(3):e9749.
Chen MZ, Moily NS, Bridgford JL, et al. A thiol probe for measuring unfolded protein load and proteostasis in cells. Nat Commun. 2017;8(1):474.
Hlady RA, Zhou D, Puszyk W, et al. Initiation of aberrant DNA methylation patterns and heterogeneity in precancerous lesions of human hepatocellular cancer. Epigenetics. 2017;12(3):215–25.
Casalino L, Verde P. Multifaceted roles of DNA methylation in neoplastic transformation, from tumor suppressors to EMT and metastasis. Genes (Basel). 2020;11(8):922.
Jones PA, Baylin SB. The fundamental role of epigenetic events in cancer. Nat Rev Genet. 2002;3(6):415–28.
Wang L, Shi C, Yu J, Xu Y. FOXM1-induced TYMS upregulation promotes the progression of hepatocellular carcinoma. Cancer Cell Int. 2022;22(1):47.
Fan W, Ma H, Jin B. Expression of FOXM1 and PLK1 predicts prognosis of patients with hepatocellular carcinoma. Oncol Lett. 2022;23(5):146.
Shaulian E. AP-1–the Jun proteins: oncogenes or tumor suppressors in disguise? Cell Signal. 2010;22(6):894–9.
Brennan A, Leech JT, Kad NM, Mason JM. Selective antagonism of cJun for cancer therapy. J Exp Clin Cancer Res. 2020;39(1):184.
Kinkel SA, Galeev R, Flensburg C, et al. Jarid2 regulates hematopoietic stem cell function by acting with polycomb repressive complex 2. Blood. 2015;125(12):1890–900.
Pasini D, Cloos PA, Walfridsson J, et al. JARID2 regulates binding of the Polycomb repressive complex 2 to target genes in ES cells. Nature. 2010;464(7286):306–10.
Lei X, Xu JF, Chang RM, et al. JARID2 promotes invasion and metastasis of hepatocellular carcinoma by facilitating epithelial-mesenchymal transition through PTEN/AKT signaling. Oncotarget. 2016;7(26):40266–84.
Li S, Wu Z, Li Q, et al. The prognostic value of AT-Rich Interaction Domain (ARID) family members in patients with hepatocellular carcinoma. Evid Based Complement Alternat Med. 2022;2022:1150390.
Chan BC, Ching AK, To KF, et al. BRE is an antiapoptotic protein in vivo and overexpressed in human hepatocellular carcinoma. Oncogene. 2008;27(9):1208–17.
Chui YL, Ching AK, Chen S, et al. BRE over-expression promotes growth of hepatocellular carcinoma. Biochem Biophys Res Commun. 2010;391(3):1522–5.
Qiu J, Zhang S, Wang P, et al. BUB1B promotes hepatocellular carcinoma progression via activation of the mTORC1 signaling pathway. Cancer Med. 2020;9(21):8159–72.
Zhuang L, Yang Z, Meng Z. Upregulation of BUB1B, CCNB1, CDC7, CDC20, and MCM3 in tumor tissues predicted worse overall survival and disease-free survival in hepatocellular carcinoma patients. Biomed Res Int. 2018;2018:7897346.
Marin JJG, Macias RIR, Monte MJ, et al. Molecular bases of drug resistance in hepatocellular carcinoma. Cancers (Basel). 2020;12(6):1663.
Romero-Garcia S, Prado-Garcia H, Carlos-Reyes A. Role of DNA methylation in the resistance to therapy in solid tumors. Front Oncol. 2020;10:1152.
This work was supported by grants from National Research Foundation (NRF-2022R1F1A107455812) in the Republic of Korea.
Ethics approval and consent to participate
This manuscript is a not involving human or animal subjects were performed. Therefore, an institutional or ethical review was not considered necessary by the authors.
The authors have no competing interest to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Pan-cancer view of SAMD13 gene expression in GENT2 database. Figure S2. Scatterplot of DNA methylation levels of SAMD13 in samples of normal and HCC. The probes located in CpG N_Shore (A), S_Shore (B), S_Shelf (C), Open_Sea (D), and Island (E), respectively. Figure S3. Spearman’s correlation analysis between each methylation level and mRNA levels of SAMD13 gene in TCGA–LICH. The probes located in CpG N_Shore (A), S_Shore (B), S_Shelf (C), Open_Sea (D), and Island (E), respectively. Figure S4. Relevance of SAMD13 gene to patient HCC recurrence in GSE76427 (A) and TCGA–LIHC (B). The OSdream database (https://bioinfo.henu.edu.cn/OSdream/OSdream.html) was used to predict the SAMD13 gene and HCC recurrence.
Functional and pathway enrichment analysis of SAMD13 on Pathway Commons including molecular function, biological process, and reactome pathway.
Functional and pathway enrichment analysis of SAMD13 on STRING including molecular function, biological process, reactome pathway, and KEGG pathways.
About this article
Cite this article
Yoo, W., Kim, S. & Noh, K. SAMD13 serves as a useful prognostic biomarker for hepatocellular carcinoma. Eur J Med Res 28, 514 (2023). https://doi.org/10.1186/s40001-023-01347-5