Skip to main content

Proteomics-based identification of TMED9 is linked to vascular invasion and poor prognoses in patients with hepatocellular carcinoma



Due to the difficulties in early diagnosing and treating hepatocellular carcinoma (HCC), prognoses for patients remained poor in the past decade. In this study, we established a screening model to discover novel prognostic biomarkers in HCC patients.


Candidate biomarkers were screened by liquid chromatography with tandem mass spectrometry (LC-MS/MS) analyses of five HCC normal (N)/tumor (T) paired tissues and preliminarily verified them through several in silico database analyses. Expression levels and functional roles of candidate biomarkers were respectively evaluated by immunohistochemical staining in N/T paired tissue (n = 120) and MTS, colony formation, and transwell migration/invasion assays in HCC cell lines. Associations of clinicopathological features and prognoses with candidate biomarkers in HCC patients were analyzed from GEO and TCGA datasets and our recruited cohort.


We found that the transmembrane P24 trafficking protein 9 (TMED9) protein was elevated in HCC tissues according to a global proteomic analysis. Higher messenger (m)RNA and protein levels of TMED9 were observed in HCC tissues compared to normal liver tissues or pre-neoplastic lesions. The TMED9 mRNA expression level was significantly associated with an advanced stage and a poor prognosis of overall survival (OS, p = 0.00084) in HCC patients. Moreover, the TMED9 protein expression level was positively correlated with vascular invasion (p = 0.026), OS (p = 0.044), and disease-free survival (p = 0.015) in our recruited Taiwanese cohort. In vitro, manipulation of TMED9 expression in HCC cells significantly affected cell migratory, invasive, proliferative, and colony-forming abilities.


Ours is the first work to identify an oncogenic role of TMED9 in HCC cells and may provide insights into the application of TMED9 as a novel predictor of clinical outcomes and a potential therapeutic target in patients with HCC.


Hepatocellular carcinoma (HCC) is the most fatal malignancy disease among primary liver cancers and is one of the leading causes of cancer-related deaths worldwide [1]. HCC develops from underlying chronic liver disease and cirrhosis, which are usually caused by hepatitis B and C viral (HBV and HCV) infections, alcohol consumption, or diabetes and nonalcoholic fatty liver disease [2,3,4]. The incidence of HCC is highest in Asia, including Taiwan, and Africa, and it occurs in males more often than in females [5]. The median survival time of HCC in Taiwan is significantly higher than that in Africa, where the median survival is only 2.5 months [6]. To the present, several therapeutic strategies have been applied as HCC treatment, including chemoradiotherapy, surgery, and systemic therapy with sorafenib [7, 8], but low therapeutic and diagnostic efficiencies still lead to severe mortality of patients with this disease. Thus, we wanted to illuminate novel biomarkers for HCC diagnoses or therapeutic targets.

Images from computed tomography (CT) or magnetic resonance imaging (MRI) are key approaches for evaluating the progression of HCC and formulating future treatment plans [9]. In addition, a high serum alpha-fetoprotein (AFP) level has also been widely used for HCC diagnoses and recurrence predictions. However, due to its low sensitivity and specificity, serum AFP is not a highly accurate biomarker for HCC diagnoses [10]. Recently, systemic proteomic analyses have shown promise in exploring proper biomarkers for predicting cancer prognoses [11]. For example, matrix-assisted laser desorption/ionization (MALDI)-imaging mass spectrometry (IMS), multiple reaction monitoring (MRM) MS and liquid chromatography-tandem MS (LC-MS/MS) have been used to discover diagnostic, prognostic and survival biomarkers in liver cancer [12,13,14,15]. Although MS protein expression data in HCC are available from recently published works, the prognostic values and functional roles of these proteins in clinical HCC samples and HCC cells have not yet been investigated [12, 14, 16].

Herein, we mined candidate biomarkers of HCC from five paired normal and HCC tissues using LC-MS/MS analyses. Among the candidate proteins, overexpression of transmembrane P24 trafficking protein 9 (TMED9) was observed in HCC samples and pre-neoplastic lesions compared to normal tissues. HCC patients with TMED9high tumors had a higher frequency of developing an advanced stage and vascular invasion, and had shorter overall survival (OS) and disease-free survival (DFS) times compared to patients with TMED9low tumors. Moreover, TMED9-knockdown significantly inhibited the growth and motility of HCC cell lines. Taken together, this study first identified that TMED9 might be a novel prognostic biomarker and therapeutic target for HCC.

Material and methods

Cell culture

The human Mahlavu and HCC36 cell lines were respectively maintained in the Dulbecco's modified Eagle and minimum essential medium (DMEM and MEM; Gibco, Waltham, MA, USA) with 10% fetal bovine serum (FBS), 2 mM glutamine, 100 units/mL penicillin, and 100 mg/mL streptomycin, and was incubated in a 5% CO2 humidified atmosphere at 37 °C. The TMED9 short hairpin (sh)RNA constructs were purchased from the National RNAi core Facility at Academic Sinica (Taipei, Taiwan). The target sequences were GCC AGT CTT CTG TCT TCC TTT and CGG CAC CTC AAG AGC TTC TTT. The HA-TMED9 construct was purchased form Applied Biological Materials Inc.

In-gel digestion

Coomassie dye staining gels were cut into small pieces and washed with distilled water (d2H2O). Gel pieces were destained with 25 mM ammonium bicarbonate (ABC) overnight at 4 °C, which was then replaced with reduction buffer (10 mM dithiothreitol (DTT)/25 mM ABC) for 1 h at 56 °C. Next, the reduction buffer was discarded, and alkylation buffer (55 mM iodoacetamide/25 mM ABC) was added for 1 h at room temperature in the dark and washed twice with 40% acetonitrile (ACN) for 10 min and then dehydrated by treatment with 100% ACN. Dry pieces were re-swollen with 0.12 μg of modified trypsin in 25 mM ABC and digested overnight at 37 °C. Peptides were harvested using 60% ACN/0.1% trifluoroacetic acid (TFA) and sonication. After the solution was dried in a vacuum, the peptides were re-dissolved in 0.1% TFA for the LC-MS/MS analysis.

LC-MS/MS analysis

LC-MS/MS analyses were performed on an LTQ-FT ion trap mass spectrometer (LTQ FT Ultra, Thermo Scientific, Waltham, MA, USA). The Mascot software package (Matrix Science, Boston, MA, USA) was used for protein identification according to SWISS-PROT, an annotated protein sequence database [17]. The peptide mass tolerance and fragment mass tolerance were set at 100 ppm and 0.25 Da, respectively. Matches with scores higher than the 95% confidence level were regarded as significant.

Cell-proliferation and colony-formation assays

Mahlavu and HCC36 cells were stably infected with virus carrying either shTMED9, HA-TMED9, or their respective controls. For the proliferation assay, Mahlavu (5 × 103) and HCC36 (6 × 103) cells were seeded into 96-well dishes. After 48 (for HCC36) or 72 h (for Mahlavu), Alamar blue dye was added to each well and incubated at 37 °C for 3 h. The absorbance was measured at 570 nm using a microplate reader (MQX200; Bio-Tek Instruments, Winooski, VT, USA). Values are the mean ± standard deviation (SD) of triplicate wells and were normalized to that of the control group to determine the multiples of proliferative ability.

For the colony-formation assay, 1.5 × 103 Mahlavu cells were seeded into 6-well dishes. After 7–10 days, cells were fixed with methanol and then stained with 1% crystal violet. Numbers of colonies were further quantified by Image J software (National Institutes of Health, Bethesda, MD, USA).

Transwell migration and Matrigel-invasion assays

For the cell-migration analysis, 105 HCC36 cells infected with virus carrying the Luc-control or HA-TMED9 were plated in an uncoated top chamber (24-well insert; pore size, 8 μm; Corning Costar, Corning, NY). For the cell-invasion analysis, 1.5 × 104 Mahlavu cells infected with virus carrying the sh-control or shTMED9 were seeded in a 20 μg Matrigel (BD Biosciences, Bedford, MA)-coated top chamber. In both assays, top chamber all contained serum-free medium and medium supplemented with serum was used as a chemoattractant in the lower chamber. After 48 and 24 h of incubation, migrated and invaded cells in the lower chamber were fixed with methanol and stained with 0.1% crystal violet. The number of cells migrating through or invading through the membrane was counted under a light microscope (× 100, five random fields per well).

Bioinformatics analysis

We collected the GSE6764, GSE76311, and GSE102079 microarray datasets, which contain transcriptome profiles of normal tissues, dysplastic tissues, and HCC tissues, from the Gene Expression Omnibus (GEO) database [18,19,20]. Representative figures for the RNA sequencing analysis were obtained from the visualization platform Gene Expression Profiling Interactive Analysis (GEPIA) ( [21]. Further, the website STRING ( was utilized to explore protein–protein interaction networks of TMED9-regulated genes, and the cBioPortal ( was used to investigate correlations of TMED9 with other genes. Correlations of TMED3 with TMED9 and with patient prognosis data were calculated and downloaded from UCSC Xena ( [22].

Patients and tissue microarray (TMA) construction

This study is a hospital-based case–control study. A total of 182 patients who had undergone surgical resection of HCC at Taipei Veterans General Hospital in Taipei, Taiwan were recruited as a case group between 1990 and 2006, with an average follow-up time of 60.02 months. The patients were diagnosed with HCC according to the characteristic criteria of the national guidelines for HCC. Collection of patient data was approved by the Ethics Committee of Taipei Veterans General Hospital (201010021IC). We constructed a formalin-fixed, paraffin-embedded TMA composed of 182 HCC tissue cores. The clinicopathologic data of 182 patients with HCC were obtained by a retrospective review of the medical records. Before 2002, patients with tumor recurrence could receive transarterial chemoembolization, surgery, and a percutaneous ethanol injection. After 2002, other therapeutic choices, such as radiofrequency tumor ablation and liver transplantation, were provided for recurrent patients. Because approval of targeted therapy or immunotherapy was not available in 1990 to 2006, no patients in this study received that drug. Collection of our Taiwanese HCC cohort was described in detail in our previous study [23].

Immunohistochemical (IHC) staining

IHC staining was automatically performed with an immunostainer (Ventana Discovery XT autostainer, Tucson, AZ, USA). TMA slides were stained with a polyclonal rabbit anti-human TMED9 antibody (GeneTex, Hsinchu, Taiwan), and cytoplasmic expression in HCC cells was evaluated. The staining intensity and percentage of TMED9 were recorded for further analysis. According to the intensity of staining, samples were scored from 0 (no staining) to 3 (heavy staining). The percentage of positive cells was calculated to range 0–100%. The intensity score was multiplied by the percentage of positive cells as the final IHC score (0–300). The high or low expression of TMED9 in HCC patients was divided by a cutoff value of 60.

Western blotting

Protein extraction and Western blot preparation were described in a previous study [24]. The primary antibodies used in this study were TMED9 (1:1000; GeneTex, Hsinchu, Taiwan); TGF-α (1:250; Santa Cruz, CA, USA), p-ERK (1:1000; Santa Cruz, CA, USA), GLI1 (1:1000; Santa Cruz, CA, USA), β-catenin (1:1000; Cell Signaling, MA, USA), p-STAT3 (1:1000; Cell Signaling, MA, USA), t-STAT3 (1:1000; Cell Signaling, MA, USA), t-ERK (1:1000; Cell Signaling, MA, USA), and GAPDH (1:5000, Sigma, St. Louis, MO, USA).

Dot blot analysis

Briefly, 5 × 105 HCC cells were seeded in a 6-cm Petridish and infected with virus carrying either shTMED9, HA-TMED9, or their respective controls. After infection, the cell culture medium was replaced by serum-free medium for another 24 h and harvested for dot blotting. Each dot was presented the proteins in 300 µL soup. The hybridization and detection were followed by Western blot protocol.

Statistical analysis

All spectra were processed by baseline subtraction, peak detection, and peak-area calculation according to default settings. Discrimination between control and tumor samples was performed using an algorithm. Values from in vitro studies are presented as the mean ± standard deviation (SD). Data were analyzed using Student’s t-test when two groups were compared. Correlations of TMED9 with clinicopathologic features of HCC were examined by Pearson’s Chi-squared test. Cumulative survival was analyzed by the Kaplan–Meier method. Risk factors affecting survival were assessed by a Cox proportional hazards regression model. A p value of < 0.05 was considered a statistically significant difference.


Potential biomarkers of HCC are discovered by the LC–MS/MS analysis

In this study, we first applied sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) coupled with LC-MS/MS to identify candidate biomarkers from five sets of paired HCC N/T tissues. Fresh protein lysates from paired N/T tissues were separated by SDS-PAGE and analyzed by LC/MS/MS. Then, protein IDs were identified by a Mascot engine through searching the SWISS-PROT database, and protein expression levels were calculated according to Mascot scores. Hematoxylin and eosin (H&E) staining of samples was used by a pathologist to define the cancer and non-cancerous parts (Fig. 1a). By interpretation of the protein ID analysis and a literature search, we selected several candidate proteins which had a > 1 and those with a < 1 score ratio of HCC tumor versus normal samples (Table 1, Additional file 1: Table S1).

Fig. 1
figure 1

Clinical values of candidate genes derived from an LC–MS/MS analysis of paired normal (N)/tumorous (T) hepatocellular carcinoma (HCC) samples. a Representative H&E staining pictures of the histomorphology of normal and HCC samples. b Forest plot showing hazard ratios (HRs) and 95% confidence intervals for the association of candidate genes and overall survival in patients with HCC

Table 1 List of protein score ratios (> 1) of hepatocellular carcinoma (HCC) tumor (T) versus normal (N) samples

We next validated the prognostic values of these 20 HCC tissue-upregulated candidates by online databases, including SurvExpress ( and TCGA ( [22, 25]. Most of these 20 candidate proteins exhibited higher mRNA levels in HCC tissues compared to normal liver tissues according to the GEPIA analysis ( (Table 2). According to the OS analysis by the SurvExpress website, we further found that elevated mRNA expression of some candidates correlated with a poor prognosis in HCC patients, especially Ezrin, Rab18, ILEU (SERPINB1), PSD7 (PSMD7), TMED9, and Rab7 (Fig. 1b). Among these, Rab1a and Rab18 had previously been reported to promote HCC cell proliferation or act as a poor prognostic marker in HCC patients [26, 27]. Ezrin overexpression was associated with vascular invasion and a poor prognosis in patients with HCC [28, 29]. These results indicated the high credibility and applicability of our screening system to search for candidate biomarkers in HCC patients.

Table 2 mRNA expression levels of candidate proteins in hepatocellular carcinoma (HCC) tissues

mRNA expression level of TMED9 is significantly higher in HCC tumor tissues and is correlated with an advanced stage and poor patient prognoses

Among these 20 candidate proteins, the role of TMED9, a protein secretion modulator, in HCC is largely unexplored. According to LC–MS/MS data, several peptides belonging to TMED9 were identified in patient samples (Additional file 1: Table S2). For instance, the LC–MS/MS fragment of R.QLVEQVEQIQK.E was identified as TMED9_HUMAN (Fig. 2). Thus, we hypothesized that TMED9 could be a putative prognostic biomarker in HCC. To further identify the potential prognostic value of TMED9 in HCC, we first analyzed its expression levels in 30 pre-neoplastic lesions (including 13 cirrhotic and 17 dysplastic samples) and 35 HCC samples (including 18 early and 17 advanced HCC) from the GSE6764 dataset of the GEO database. Significantly higher TMED9 transcripts were observed in tumors samples, especially in advanced HCC, than in pre-neoplastic lesions (Fig. 3a, b). In addition, we also examined TMED9 levels in paired N/T tissues in the GSE76311 dataset (Additional file 2: Figure S1a) and normal tissue versus tumor tissue in GSE102079 dataset (Additional file 2: Figure S1b). These results indicated significantly elevated TMED9 transcripts in malignant tissues. In addition to the GEO database, we then explored TMED9 in HCC via the GEPIA website ( based on the TCGA and GTEx databases. We also observed that TMED9 transcripts were more highly expressed in HCC tumor samples compared to normal liver tissues (Fig. 3c). Moreover, the prognostic value of TMED9 was also verified by GEPIA, and results showed that the OS of HCC patients was adversely affected by a higher TMED9 expression level (p = 0.00084) (Fig. 3d).

Fig. 2
figure 2

Results of an LC–MS/MS analysis for TMED9. a Representative data showing the sequence of the trypsinized peptide, R.QLVEQVEQIQK.E, which was matched to TMED9 in hepatocellular carcinoma (HCC) patient samples. b Matching B-ions and Y-ions (red) fragments were recorded to identify the sequence of TMED9. The retained charge on the amino-terminal part of the peptide is shown as B-ions, and the retained charge on the carboxyl-terminal part of the peptide is Y-ions

Fig. 3
figure 3

In silico analysis of TMED9 expression in hepatocellular carcinoma (HCC). a, b Relative TMED9 mRNA levels (208757_at) in 30 pre-neoplastic lesions (including 13 cirrhotic samples and 17 dysplastic samples) and 35 HCC (including 18 early-stage and 17 advanced-stage HCC) samples from GSE6764 microarray datasets. Expression data were normalized to the median intensity of all probes and scaled to log2 transformation. c Expression levels of TMED9 transcripts in normal and HCC samples. d Kaplan–Meier curves for overall survival of patients with HCC, as categorized according to high or low expression of TMED9. The p value indicates a comparison between patients with TMED9high and TMED9low. Data in both c and d are from GEPIA online available databases (

Expression of TMED9 protein is significantly elevated in HCC tissues, and high TMED9 expression is correlated with vascular invasion and poor prognoses in HCC patients

In addition to analyzing mRNA levels of TMED9 from datasets available online, we further verified protein levels of TMED9 by IHC staining from a Taiwanese HCC tissue array cohort. Representative examples of tumors showing overall negative (score 0), low (score 1), moderate (score 2), and high TMED9 (score 3) intensities are illustrated in Fig. 4a, and all positive cases revealed a diffuse cytoplasmic TMED9 distribution in cancer cells. According to our analysis, the score of TMED9 in HCC was significantly higher than that in corresponding non-cancerous tissues (p < 0.0001) (Fig. 4b). We thus stratified patients into high- and low-TMED9 groups, and high expression levels of TMED9 were significantly correlated with worse OS (p = 0.044; Fig. 4c) and DFS (p = 0.015; Fig. 4d) in our recruited cohort.

Fig. 4
figure 4

TMED9 expression correlates with overall and disease-free survival in patients with hepatocellular carcinoma (HCC) after surgical resection. a Tissue microarray of HCC was immunohistochemically analyzed for TMED9 and showed different intensity scores for TMED9 expression. b TMED9 expression scores were significantly higher in tumor than non-tumor tissues. Statistical significance was analyzed by a paired t test (p < 0.0001). c, d Kaplan–Meier plot of overall c and disease-free d survival of 182 patients with HCC stratified by TMED9 expression levels. A log-rank test was used to show differences between groups

Clinicopathological parameters of HCC patients are shown in Table 3. There were 182 patients, including 155 males and 27 females. Survival follow-up was available in all cases and ranged 1 ~ 172 months (median, 52 months; mean, 60.02 months). Among these patients, 43 cases (23.6%) experienced no recurrence, and 139 cases (76.4%) had disease recurrence during the follow-up period. A Chi-squared analysis was conducted to explore relationships between selected clinicopathologic features and TMED9 expression. Table 3 shows that the TMED9 expression level was significantly correlated with hepatic vein invasion (p = 0.026), but not other clinicopathologic features such as the AFP level, HBV, liver cirrhosis, tumor size, and so on.

Table 3 Relationship between TMED9 expression and clinicopathological features in 182 patients with hepatocellular carcinoma (HCC)

Moreover, we utilized a univariate analysis to examine the prognostic significance of clinicopathologic variables in our HCC cohort. In addition to TMED9 expression (hazard ratio (HR), 1.463; p = 0.047), the tumor size (HR, 1.543; p = 0.023), portal vein invasion (HR, 3.463; p < 0.001), and hepatic vein invasion (HR, 3.971; p < 0.001) were all shown to have adverse impacts on OS (Table 4). Furthermore, DFS was also adversely affected by high TMED9 expression (HR, 1.506; p = 0.017), portal vein invasion (HR, 3.662; p < 0.001), and hepatic vein invasion (HR, 11.333; p < 0.001) (Table 5).

Table 4 Cox univariate and multivariate regression analysis of prognostic factors and TMED9 expression for overall survival (OS) in 182 hepatocellular carcinoma (HCC) patients
Table 5 Cox univariate and multivariate regression analysis of prognostic factors and TMED9 expression for disease-free survival (DFS) in 182 hepatocellular carcinoma (HCC) patients

TMED9 modulates the cell-proliferative, colony-forming, and invasive/migratory abilities of HCC cells

To further investigate the functional role of TMED9 in HCC cells, we knocked-down TMED9 by two lentiviral-based shRNAs in highly invasive Mahlavu HCC cells (Fig. 5a) and examined the cell-invasive, proliferative, and colony-forming abilities of HCC cells. As shown in Fig. 5b, a significantly lower proliferative ability of HCC cells was observed in TMED9-knockdown groups than in control groups. Moreover, we examined the effect of TMED9 on the long-term growth (10 days) of Mahlavu cells using a colony-formation assay and found that the number of colonies was reduced following TMED9 knockdown (Fig. 5c). Furthermore, results from the transwell invasion assay showed that significant attenuation of the invasive ability was observed in TMED9-depleted HCC cells (Fig. 5d). In comparison to Mahlavu cells, we ectopically expressed TMED9 in low-invasive HCC36 cells (Additional file 2: Figure S2a) and observed that overexpression of TMED9 significantly increased cell proliferation (Additional file 2: Figure S2b) and motility (Additional file 2: Figure S2c). These results indicated that TMED9 may be deeply involved in HCC progression, and suggested that targeting TMED9 might be a potential therapeutic strategy for arresting HCC progression.

Fig. 5
figure 5

Depletion of TMED9 decreases cell growth and invasion abilities of hepatocellular carcinoma (HCC) cells. a Mahlavu cells were infected with a lentivirus carrying control shRNA (sh-Ctrl) or two TMED9 shRNAs (shTMED9). After 72 h, knockdown efficiencies of the two TMED9 shRNAs were determined by Western blotting. bd Cell-proliferative, colony-forming, and invasive abilities of TMED9-knockdown Mahlavu cells were respectively determined by Alamar blue cell-viability (b), colony-formation (c), and transwell-invasion d assays. Upper panel of C and left panel of D: representative photomicrographs. Values are presented as the mean ± SD from three independent experiments. ***p < 0.001 compared to the sh-Ctrl group


The main finding of the present study was to identify TMED9, one of the transport proteins located in the endoplasmic reticulum (ER) and Golgi network, as a prognostic biomarker of HCC progression. Reliable mass spectrometry-based analyses have been used for drug discovery, metabolomic analysis, and cancer biomarker investigations [30, 31]. Herein, we interpreted LC–MS/MS data to identify elevated candidate proteins in HCC tissues (Table 1) for diagnostic and prognostic predictions in HCC. Similar to TMED9, one of the candidate markers, aldehyde dehydrogenase 7 family member A1 (ALDH7A1), was previously reported to be more highly expressed in HCC tissues than in non-tumorous liver tissues, but its expression in tumor tissues was not correlated with OS, recurrence, or microvascular invasion [32]. Although some candidate markers, such as Rab1, Rab18, and Ezrin which are elevated in HCC, were reported to promote HCC progression [26,27,28,29], the oncogenic or tumor-suppressive roles of other candidate proteins we found in HCC tissues remain unclear. For example, Rab7 expression regulated the migratory ability of lung cancer H1299 cells through Ras-related C3 botulinum toxin substrate 1 (RAC1) and vimentin [33]. In MCF-7 breast cancer cells, Rab7 is responsible for Akt survival signal maintenance to promote cells anoikis resistance [34]. In contrast, Rab7 was reported to be a tumor suppressor in prostate cancer based on the fact that it inhibited ligand-induced c-Met signaling [35]. In addition to Rab7, it was revealed that overexpression of serine protease inhibitor B1 (SERPINB1; also known as ILEU) promotes the motility of oral squamous cell carcinoma [36], but the epigenetic suppression of SERPINB1 predicts a poorer prognosis in prostate cancer patients [37]. Thus, further investigations of the roles of candidate proteins which we found in HCC tissues are still needed. After a literature review and our preliminary in silico analysis, we found that TMED9 showed a poor prognostic impact on HCC, but has rarely been investigated in HCC.

TMED9 belongs to the transmembrane emp24 domain-containing protein (TMED)/p24 family, including TMED1–7 and TMED9–11 which are responsible for selecting cargo in the processing ER-Golgi network of proteins (secretory pathway) and innate immune signaling [38, 39]. Given that a huge diversity exists among the proteins of this family, TMED proteins can influence the stability of each other while forming monomers or dynamic complexes [40, 41]. Abnormalities in the controlled transport of proteins in the secretory pathway contribute to diseases such as cancer. For example, a reduction in TMED2 was reported to increase the possibility of liver tumorigenesis [39]. In contrast, elevated levels of TMED2 were observed in patients with ovarian and breast cancers. TMED2 was identified as a promoter of cell proliferation and migration in ovarian cancer cells through activating Akt and as a poor prognostic factor in patients with breast cancer [42, 43]. TMED3 was reported to be an oncoprotein in HCC, prostate cancers, breast cancers, and kidney cancers. [44,45,46,47]. For example, TMED3 enhanced cell migration by activating interleukin (IL)-11/signal transducer and activator of transcription 3 (STAT3) signaling and promoted HCC progression [47]. TMED3 was also reported to promote proliferation and motility of breast cancer cells by activating WNT/β-catenin pathway [48]. In contrast, TMED3 was identified as a tumor suppressor in colon cancer and proposed to inhibit metastasis by repressing TMED9. TMED3 depletion can trigger a TMED9-induced metastasis via inducing transforming growth factor (TGF)-α secretion and upregulating CNIH4/TGF-α/GLI1 signaling [49]. According to those results, the malignant properties of TMED2 and TMED3 in cancer are cell type-specific.

TMED2 and TMED3 were respectively reported to be a tumor suppressor and an oncogene in HCC [39, 47]. However, the oncogenic role of TMED9 was only identified in colon cancer [49], but not in other cancer types. In this study, we first reported the prognostic and functional roles of TMED9 in HCC. We further dissected interacting neighbors of TMED9 using the STRING database ( and found ten potential TMED9-associated genes, including TMED10, TMED2, TMED3, TMED7, STX5, KDELR2, SEC22B, SURF4, GORASP1, and NAGLU (Fig. 6a), which are necessary for Golgi vesicle transport and related vesicle-trafficking functions (Fig. 6b). We then verified correlations of TMED9 with these genes in 373 HCC human samples using the cBioportal platform and observed that TMED9 expression was only significantly correlated with STX5, TMED3, and SURF4 (Fig. 6c). Furthermore, TCGA database showed that a poorer prognosis was observed in HCC patients with high levels of TMED9 and TMED3 in tumors compared to patients with low levels of TMED9 and TMED3 in tumors (Fig. 6d). Opposite roles of TMED9 and TMED3 in regulating the progression of colon cancer were reported [49]. Although the roles of TMED9 and TMED3 in HCC differed from those in colon cancer, how the interaction between TMED9 and TMED3 promotes HCC progression remains for further evaluation. Herein, we first evaluated the effect of TMED9 knockdown on TMED9- and TMED3-modulated signaling molecules which have been reported in other cancer types [48, 49]. We observed that TMED9 knockdown in Mahlavu cells decreased expression of β-catenin and GLI1, phosphorylation of ERK and STAT3, and secretion of TGF-α compared to control cells (Fig. 6e). In comparison, β-catenin, GLI1, p-ERK, p-STAT3, and secreted TGF-α were all upregulated in TMED9-overexpressed HCC36 cells (Fig. 6f). These results suggested that CNIH4/TGF-α/GLI1 and WNT/β-catenin pathways might be involved in TMED9-induced HCC progression and these issues should be further investigated in future research.

Fig. 6
figure 6

Prediction of potential interactions between cargo-transport proteins with TMED9 and their prognostic values in hepatocellular carcinoma (HCC). a TMED9 protein–protein interaction network of 10 differentially expressed genes from the STRING database. b Gene ontology analysis revealed that TMED9 participates in endoplasmic reticular (ER)-Golgi cargo trafficking. c Correlation analysis of The Cancer Genome Atlas (TCGA) Liver Hepatocellular Carcinoma database using the cBioPortal revealed positive correlations between TMED9 and mRNA levels of STX5, TMED3, and SURF4. d Effects of TMED9 and TMED3 mRNA levels on patient overall survival from the TCGA cohort. e, f Mahlavu (e) and HCC36 (f) cells respectively expressed shTMED9 and HA-TMED9 as indicated. Cell lysates and conditioned media were harvested to determine the expression of β-catenin and GLI1, phosphorylation of ERK and STAT3, and secretion of TGF-α in both cells by Western blotting and dot blot analyses

Our current study is not without limitations. First, since the HCC samples were collected from 1990 to 2006 in our study, the examination of HBV/HCV viral levels was not been determined and no anti-viral therapies were available at that period. Thus, we can't analyze whether the baseline viral level affecting HCC patient outcome. Moreover, the sample size of the current study was still not large enough and might lead to a limited statistical impact on the accuracy and precision of the results. For example, TMED9 expression level has a negative impact on OS and DFS from univariate, but not multivariate analysis. Larger HCC cohort should be further recruited to confirm this result in future. Furthermore, the clinical data and follow-up conditions were prospectively collected. Even we did our best to collect detail clinical information, there were still few patients lacking the status of cirrhosis or hepatic/portal vein invasion.


According to our study, TMED9 is present in fresh-frozen HCC tissues and can be successfully detected by an LC-MS/MS analysis. Using this technology, we first identified that TMED9 could be a valuable prognostic biomarker for HCC and revealed the functional role of TMED9 in HCC cells. Whether TMED9 can interact with TMED3 or other cargo-transport proteins to promote HCC cell proliferation and motility needs be further addressed; yet our present findings strongly support the targeting of TMED9 as a novel strategy for HCC treatment.

Availability of data and materials

All data used during the current study are available from the corresponding author on reasonable request.



Ammonium bicarbonate




Disease-free survival




Gene Expression Omnibus


Gene Expression Profiling Interactive Analysis


Genotype-Tissue Expression


Hepatocellular carcinoma


Hazard ratio




Liquid chromatography-tandem mass spectrometry


Optimum cutting temperature


Overall survival


The Cancer Genome Atlas


Transmembrane P24 trafficking protein 9


Transcripts per million


  1. Forner A, Reig M, Bruix J. Hepatocellular carcinoma. Lancet. 2018;391:1301–14.

    Article  PubMed  Google Scholar 

  2. Balogh J, Victor D 3rd, Asham EH, Burroughs SG, Boktour M, Saharia A, Li X, Ghobrial RM, Monsour HP Jr. Hepatocellular carcinoma: a review. J Hepatocell Carcinoma. 2016;3:41–53.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Batey RG, Burns T, Benson RJ, Byth K. Alcohol consumption and the risk of cirrhosis. Med J Aust. 1992;156:413–6.

    Article  CAS  PubMed  Google Scholar 

  4. Rinella ME. Nonalcoholic fatty liver disease: a systematic review. JAMA. 2015;313:2263–73.

    Article  CAS  PubMed  Google Scholar 

  5. Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM. Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer. 2010;127:2893–917.

    Article  CAS  PubMed  Google Scholar 

  6. Yang JD, Hainaut P, Gores GJ, Amadou A, Plymoth A, Roberts LR. A global view of hepatocellular carcinoma: trends, risk, prevention and management. Nat Rev Gastroenterol Hepatol. 2019;16:589–604.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Heimbach JK, Kulik LM, Finn RS, Sirlin CB, Abecassis MM, Roberts LR, Zhu AX, Murad MH, Marrero JA. AASLD guidelines for the treatment of hepatocellular carcinoma. Hepatology. 2018;67:358–80.

    Article  PubMed  Google Scholar 

  8. Gosalia AJ, Martin P, Jones PD. Advances and future directions in the treatment of hepatocellular carcinoma. Gastroenterol Hepatol (N Y). 2017;13:398–410.

    Google Scholar 

  9. Lee YJ, Lee JM, Lee JS, Lee HY, Park BH, Kim YH, Han JK, Choi BI. Hepatocellular carcinoma: diagnostic performance of multidetector CT and MR imaging-a systematic review and meta-analysis. Radiology. 2015;275:97–109.

    Article  PubMed  Google Scholar 

  10. Oka H, Tamori A, Kuroki T, Kobayashi K, Yamamoto S. Prospective study of alpha-fetoprotein in cirrhotic patients monitored for development of hepatocellular carcinoma. Hepatology. 1994;19:61–6.

    Article  CAS  PubMed  Google Scholar 

  11. Sallam RM. Proteomics in cancer biomarkers discovery: challenges and applications. Dis Markers. 2015;2015:321370.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Marquardt C, Tolstik T, Bielecki C, Kaufmann R, Crecelius AC, Schubert US, Settmacher U, Stallmach A, Dirsch O. MALDI imaging-based classification of hepatocellular carcinoma and non-malignant lesions in fibrotic liver tissue. Z Gastroenterol. 2015;53:33–9.

    Article  CAS  PubMed  Google Scholar 

  13. Zhang JY, Liang R, Wei JZ, Ye JX, He Q, Chunling Y, Ye JZ, Li YQ, Liu ZH, Lin Y. Identification of candidate biomarkers in malignant ascites from patients with hepatocellular carcinoma by iTRAQ-based quantitative proteomic analysis. Biomed Res Int. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Baniasadi H, Gowda GAN, Gu HW, Zeng A, Zhuang S, Skill N, Maluccio M, Raftery D. Targeted metabolic profiling of hepatocellular carcinoma and hepatitis C using LC-MS/MS. Electrophoresis. 2013;34:2910–7.

    CAS  PubMed  Google Scholar 

  15. Ahn YH, Shin PM, Oh NR, Park GW, Kim H, Yoo JS. A lectin-coupled, targeted proteomic mass spectrometry (MRM MS) platform for identification of multiple liver cancer biomarkers in human plasma. J Proteomics. 2012;75:5507–15.

    Article  CAS  PubMed  Google Scholar 

  16. Chen XL, Zhou L, Yang J, Shen FK, Zhao SP, Wang YL. Hepatocellular carcinoma-associated protein markers investigated by MALDI-TOF MS. Mol Med Rep. 2010;3:589–96.

    CAS  PubMed  Google Scholar 

  17. Bairoch A, Apweiler R. The SWISS-PROT protein sequence data bank and its supplement TrEMBL in 1999. Nucleic Acids Res. 1999;27:49–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Wurmbach E, Chen YB, Khitrov G, Zhang W, Roayaie S, Schwartz M, Fiel I, Thung S, Mazzaferro V, Bruix J, et al. Genome-wide molecular profiles of HCV-induced dysplasia and hepatocellular carcinoma. Hepatology. 2007;45:938–47.

    Article  CAS  PubMed  Google Scholar 

  19. Chaisaingmongkol J, Budhu A, Dang H, Rabibhadana S, Pupacdi B, Kwon SM, Forgues M, Pomyen Y, Bhudhisawasdi V, Lertprasertsuke N, et al. Common molecular subtypes among asian hepatocellular carcinoma and cholangiocarcinoma. Cancer Cell. 2017;32(57–70):e53.

    Google Scholar 

  20. Chiyonobu N, Shimada S, Akiyama Y, Mogushi K, Itoh M, Akahoshi K, Matsumura S, Ogawa K, Ono H, Mitsunori Y, et al. Fatty acid binding protein 4 (FABP4) overexpression in intratumoral hepatic stellate cells within hepatocellular carcinoma with metabolic risk factors. Am J Pathol. 2018;188:1213–24.

    Article  CAS  PubMed  Google Scholar 

  21. 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:W98–102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Goldman MJ, Craft B, Hastie M, Repecka K, McDade F, Kamath A, Banerjee A, Luo YH, Rogers D, Brooks AN, et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. 2020;38:675–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Chen MH, Jan YH, Chang PM, Chuang YJ, Yeh YC, Lei HJ, Hsiao M, Huang SF, Huang CY, Chau GY. Expression of GOLM1 correlates with prognosis in human hepatocellular carcinoma. Ann Surg Oncol. 2013;20:S616–24.

    Article  PubMed  Google Scholar 

  24. Yang YC, Chien MH, Liu HY, Chang YC, Chen CK, Lee WJ, Kuo TC, Hsiao M, Hua KT, Cheng TY. Nuclear translocation of PKM2/AMPK complex sustains cancer stem cell populations under glucose restriction stress. Cancer Lett. 2018;421:28–40.

    Article  CAS  PubMed  Google Scholar 

  25. Aguirre-Gamboa R, Gomez-Rueda H, Martinez-Ledesma E, Martinez-Torteya A, Chacolla-Huaringa R, Rodriguez-Barrientos A, Tamez-Pena JG, Trevino V. SurvExpress: an online biomarker validation tool and database for cancer gene expression data using survival analysis. PLoS ONE. 2013;8:e74250.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Dong Z, Qi R, Guo X, Zhao X, Li Y, Zeng Z, Bai W, Chang X, Hao L, Chen Y, et al. MiR-223 modulates hepatocellular carcinoma cell proliferation through promoting apoptosis via the Rab1-mediated mTOR activation. Biochem Biophys Res Commun. 2017;483:630–7.

    Article  CAS  PubMed  Google Scholar 

  27. Gong T, Zhou B, Liu M, Chen X, Huang S, Xu Y, Luo R, Chen Z. RAB18 promotes proliferation and metastasis in hepatocellular carcinoma. Am J Transl Res. 2019;11:1009–19.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Yeh CN, Pang ST, Chen TW, Wu RC, Weng WH, Chen MF. Expression of ezrin is associated with invasion and dedifferentiation of hepatitis B related hepatocellular carcinoma. BMC Cancer. 2009;9:233.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Okamura D, Ohtsuka M, Kimura F, Shimizu H, Yoshidome H, Kato A, Miyazaki M. Ezrin expression is associated with hepatocellular carcinoma possibly derived from progenitor cells and early recurrence after surgical resection. Mod Pathol. 2008;21:847–55.

    Article  CAS  PubMed  Google Scholar 

  30. Cazares LH, Troyer D, Mendrinos S, Lance RA, Nyalwidhe JO, Beydoun HA, Clements MA, Drake RR, Semmes OJ. Imaging mass spectrometry of a specific fragment of mitogen-activated protein kinase/extracellular signal-regulated kinase kinase kinase 2 discriminates cancer from uninvolved prostate tissue. Clin Cancer Res. 2009;15:5541–51.

    Article  CAS  PubMed  Google Scholar 

  31. Djidja MC, Francese S, Loadman PM, Sutton CW, Scriven P, Claude E, Snel MF, Franck J, Salzet M, Clench MR. Detergent addition to tryptic digests and ion mobility separation prior to MS/MS improves peptide yield and protein identification for in situ proteomic investigation of frozen and formalin-fixed paraffin-embedded adenocarcinoma tissue sections. Proteomics. 2009;9:2750–63.

    Article  CAS  PubMed  Google Scholar 

  32. Pommergaard HC, Hasselby JP, Willemoe GL, Ralbovska A, Rostved AA, Rasmussen A, Schultz NA, Hillingsø J, Nørgaard Larsen P, Kugler JM. Peroxisome proliferator-activated receptor activity correlates with poor survival in patients resected for hepatocellular carcinoma. J Hepatobiliary Pancreat Sci. 2020.

    Article  PubMed  Google Scholar 

  33. Margiotta A, Progida C, Bakke O, Bucci C. Rab7a regulates cell migration through Rac1 and vimentin. Biochim Biophys Acta Mol Cell Res. 2017;1864:367–81.

    Article  CAS  PubMed  Google Scholar 

  34. Wang T, Zhang M, Ma Z, Guo K, Tergaonkar V, Zeng Q, Hong W. A role of Rab7 in stabilizing EGFR-Her2 and in sustaining Akt survival signal. J Cell Physiol. 2012;227:2788–97.

    Article  CAS  PubMed  Google Scholar 

  35. Steffan JJ, Dykes SS, Coleman DT, Adams LK, Rogers D, Carroll JL, Williams BJ, Cardelli JA. Supporting a role for the GTPase Rab7 in prostate cancer progression. PLoS ONE. 2014;9:e87882.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Tseng MY, Liu SY, Chen HR, Wu YJ, Chiu CC, Chan PT, Chiang WF, Liu YC, Lu CY, Jou YS, et al. Serine protease inhibitor (SERPIN) B1 promotes oral cancer cell motility and is over-expressed in invasive oral squamous cell carcinoma. Oral Oncol. 2009;45:771–6.

    Article  CAS  PubMed  Google Scholar 

  37. Lerman I, Ma X, Seger C, Maolake A, Garcia-Hernandez ML, Rangel-Moreno J, Ackerman J, Nastiuk KL, Susiarjo M, Hammes SR. Epigenetic suppression of SERPINB1 promotes inflammation-mediated prostate cancer progression. Mol Cancer Res. 2019;17:845–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Strating JR, Martens GJ. The p24 family and selective transport processes at the ER-Golgi interface. Biol Cell. 2009;101:495–509.

    Article  CAS  PubMed  Google Scholar 

  39. Aber R, Chan W, Mugisha S, Jerome-Majewska LA. Transmembrane emp24 domain proteins in development and disease. Genet Res (Camb). 2019;101:e14.

    Article  Google Scholar 

  40. Jerome-Majewska LA, Achkar T, Luo L, Lupu F, Lacy E. The trafficking protein Tmed2/p24beta(1) is required for morphogenesis of the mouse embryo and placenta. Dev Biol. 2010;341:154–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Jenne N, Frey K, Brugger B, Wieland FT. Oligomeric state and stoichiometry of p24 proteins in the early secretory pathway. J Biol Chem. 2002;277:46504–11.

    Article  CAS  PubMed  Google Scholar 

  42. Shi-Peng G, Chun-Lin C, Huan W, Fan-Liang M, Yong-Ning C, Ya-Di Z, Guang-Ping Z, Ye-Ping C. TMED2 promotes epithelial ovarian cancer growth. Oncotarget. 2017;8:94151–65.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Lin X, Liu J, Hu SF, Hu X. Increased expression of TMED2 is an unfavorable prognostic factor in patients with breast cancer. Cancer Manag Res. 2019;11:2203–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Ha M, Moon H, Choi D, Kang W, Kim JH, Lee KJ, Park D, Kang CD, Oh SO, Han ME, et al. Prognostic role of TMED3 in clear cell renal cell carcinoma: a retrospective multi-cohort analysis. Front Genet. 2019;10:355.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Vainio P, Mpindi JP, Kohonen P, Fey V, Mirtti T, Alanen KA, Perälä M, Kallioniemi O, Iljin K. High-throughput transcriptomic and RNAi analysis identifies AIM1, ERGIC1, TMED3 and TPX2 as potential drug targets in prostate cancer. PLoS ONE. 2012;7:e39801.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Pei J, Zhang J, Yang X, Wu Z, Sun C, Wang Z, Wang B. TMED3 promotes cell proliferation and motility in breast cancer and is negatively modulated by miR-188-3p. Cancer Cell Int. 2019;19:75.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Zheng H, Yang Y, Han J, Jiang WH, Chen C, Wang MC, Gao R, Li S, Tian T, Wang J, et al. TMED3 promotes hepatocellular carcinoma progression via IL-11/STAT3 signaling. Sci Rep. 2016;6:37070.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Zhang X, Luo Y, Li Q. TMED3 promotes proliferation and migration in breast cancer cells by activating Wnt/β-catenin signaling. Onco Targets Ther. 2020;13:5819–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Mishra S, Bernal C, Silvano M, Anand S, Ruiz IAA. The protein secretion modulator TMED9 drives CNIH4/TGFalpha/GLI signaling opposing TMED3-WNT-TCF to promote colon cancer metastases. Oncogene. 2019;38:5817–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank Ms. Tracy Tsai for her assistance with the immunohistochemical work. We also thank Dr. Hao-Jan Lei for providing the HCC patient samples.


This study was supported by grant AS-SUMMIT-108 from Academia Sinica (to M. Hsiao).

Author information

Authors and Affiliations



M-H, M-HC, C-NY, Y-CY, and M-HC designed and conceived the study; Y-CY, T-CL, Y-HJ, and W-MC performed the research; M-CT and S-MJ participated in the results discussion and technical support; Y-CY, M-HC, T-CL, and M-H wrote the paper. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Ming-Hsien Chien, Ming-Huang Chen, Chun-Nan Yeh or Michael Hsiao.

Ethics declarations

Ethical approval and consent to participate

The procedures for collecting patients were approved by the Ethics Committee of Taipei Veterans General Hospital (201010021IC). This study was supported by the Taipei Medical University Research Center of Cancer Translational Medicine from The Featured Areas Research Center Program within the framework of the Higher Education Sprout Project by the Ministry of Education in Taiwan (to M.-H. Chien).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflicts of interest.

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

: Table S1. List of protein score ratio (< 1) of hepatocellular carcinoma (HCC) tumor versus normal samples. Table S2. Matched TMED9 peptides in hepatocellular carcinoma (HCC) tissues.

Additional file 2: Figure S1.

Expression of TMED9 transcripts in paired adjacent (GSE76311, probe: TC05001018.hg.1; a and unpaired normal and tumor tissues (GSE102079, probe: 208757_at; b derived from patients with hepatocellular carcinoma (HCC). Figure S2. TMED9 overexpression promotes cell growth and migration of HCC36 cells. a HCC36 cells were infected with a lentivirus carrying control vector or HA-TMED9. After 72 h, the expression of TMED9 was determined by Western blotting. b, c Cell-proliferative b and migratory c abilities of HCC36 cells expressing control vector or HA-TMED9. Values are presented as the mean ± SD from three independent experiments. *p < 0.05; **p < 0.01 versus the vehicle control group.

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

Yang, YC., Chien, MH., Lai, TC. et al. Proteomics-based identification of TMED9 is linked to vascular invasion and poor prognoses in patients with hepatocellular carcinoma. J Biomed Sci 28, 29 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: