Genes & Cancer

Mechanistically based blood proteomic markers in the TGF-β pathway stratify risk of hepatocellular cancer in patients with cirrhosis

Xiyan Xiang1,*, Krishanu Bhowmick1,*, Kirti Shetty2, Kazufumi Ohshiro1, Xiaochun Yang1, Linda L. Wong3, Herbert Yu4, Patricia S. Latham5, Sanjaya K. Satapathy6, Christina Brennan7,8, Richard J. Dima7, Nyasha Chambwe9, Gulru Sharifova10, Fellanza Cacaj1, Sahara John1, James M. Crawford11, Hai Huang8, Srinivasan Dasarathy12, Adrian R. Krainer13, Aiwu R. He14, Richard L. Amdur1,15 and Lopa Mishra1,13,16

1The Institute for Bioelectronic Medicine, The Feinstein Institutes for Medical Research and Cold Spring Harbor Laboratory, Division of Gastroenterology and Hepatology, Northwell Health, Manhasset, NY 11030, USA
2Division of Gastroenterology and Hepatology, University of Maryland, Baltimore, MD 21201, USA
3Department of Surgery, University of Hawaii, Honolulu, HI 96813, USA
4Department of Epidemiology, University of Hawaii Cancer Center, Honolulu, HI 96813, USA
5Department of Pathology, The George Washington University, Washington, DC 20037, USA
6Department of Medicine, Sandra Atlas Bass Center for Liver Diseases and Transplantation, North Shore University Hospital/Northwell Health, Manhasset, NY 11030, USA
7Office of Clinical Research, Northwell Health, Lake Success, NY 11042, USA
8The Feinstein Institutes for Medical Research, Manhasset, NY 11030, USA
9Institute of Molecular Medicine, The Feinstein Institutes for Medical Research, Manhasset, NY 11030, USA
10Department of Medicine, Donald and Barbara Zucker School of Medicine at Hofstra/Northwell, Hempstead, NY 11549, USA
11Hofstra Northwell School of Medicine, Hempstead, NY 11549, USA
12Division of Gastroenterology and Hepatology, Cleveland Clinic, Cleveland, OH 44106, USA
13Cold Spring Harbor Laboratory, Cold Spring Harbor, NY 11724, USA
14Georgetown Lombardi Comprehensive Cancer Center, Washington, DC 20007, USA
15Quantitative Intelligence, The Institutes for Health Systems Science, The Feinstein Institutes for Medical Research, Manhasset, NY 11030, USA
16Department of Surgery, The George Washington University, Washington, DC 20037, USA
*These authors contributed equally to this work

Correspondence to: Lopa Mishra, email: [email protected], [email protected]
Richard L. Amdur, email: [email protected]

Keywords: TGF-β; hepatocellular carcinoma; myostatin; pyruvate kinase M2; biomarker

Received: November 04, 2023

Accepted: December 05, 2023

Published: February 05, 2024

Copyright: © 2024 Xiang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


Hepatocellular carcinoma (HCC) is the third leading cause of death from cancer worldwide but is often diagnosed at an advanced incurable stage. Yet, despite the urgent need for blood-based biomarkers for early detection, few studies capture ongoing biology to identify risk-stratifying biomarkers. We address this gap using the TGF-β pathway because of its biological role in liver disease and cancer, established through rigorous animal models and human studies. Using machine learning methods with blood levels of 108 proteomic markers in the TGF-β family, we found a pattern that differentiates HCC from non-HCC in a cohort of 216 patients with cirrhosis, which we refer to as TGF-β based Protein Markers for Early Detection of HCC (TPEARLE) comprising 31 markers. Notably, 20 of the patients with cirrhosis alone presented an HCC-like pattern, suggesting that they may be a group with as yet undetected HCC or at high risk for developing HCC. In addition, we found two other biologically relevant markers, Myostatin and Pyruvate Kinase M2 (PKM2), which were significantly associated with HCC. We tested these for risk stratification of HCC in multivariable models adjusted for demographic and clinical variables, as well as batch and site. These markers reflect ongoing biology in the liver. They potentially indicate the presence of HCC early in its evolution and before it is manifest as a detectable lesion, thereby providing a set of markers that may be able to stratify risk for HCC.


Most HCC is diagnosed late, resulting in poor clinical outcomes [14]. Current screening for HCC involves serial ultrasound (US) and alpha-fetoprotein (AFP) that do not have a high sensitivity/specificity, which may be related to the limited non-mechanistic basis of these strategies [5]. Outcomes remain poor using these screening methods [6] and most patients are diagnosed at an advanced stage [1, 3], indicating a failure of early detection. Thus, there is an urgent need for reliable circulating biomarkers that can stratify the risk of developing early HCC. The recently described PLSec [7] and GALAD [8] have not been validated in extensive and well-characterized studies and do not capture ongoing biology in the liver. Collectively, these results support the need for improved screening and surveillance strategies that reflect pathology in the continuum of liver disease from cirrhosis to HCC.

Molecular heterogeneity, multiple etiologies, and divergent pathophysiology associated with HCC have delineated molecular pathways that control the initiation and progression of HCC [911]. Etiologies associated with the initiation and progression of HCC include cirrhosis resulting from hepatitis viral infection, exposure to toxins (for example, alcohol or aflatoxin), metabolic syndromes, and mutations in one or more genes (for example, TERT, SMAD4, TP53, and HGF) [1215]. Additionally, signaling pathways linked to HCC include those involving the mitogen-activated protein kinase (MAPK) cascade, the phosphoinositide 3-kinase (PI3K) cascade, the WNT-β-catenin pathway, the insulin-like growth factor-1 (IGF1) pathway, and the transforming growth factor-β (TGF-β) pathway [16, 17].

Among these pathways, as a regulator of hepatocyte proliferation and tissue regeneration, the TGF-β pathway is an essential link between fibrosis/cirrhosis and the development of HCC. Alterations in the TGF-β pathway members have been observed in nearly 40% of liver and gastrointestinal cancers [1821]. TGF-β signaling plays a complex context-driven role in liver disease and cancer. It regulates not only hepatocyte proliferation but also fibrosis, functional activity of the stroma including cells of the immune system and stellate cells, genomic stability in the context of genotoxic injury, as well as hepatocarcinogenesis and the tumor microenvironment [12, 2225]. Furthermore, multiple genetically engineered mouse models with impaired TGF-β signaling are susceptible to HCC [16, 2628]. Mice heterozygous for Sptbn1, encoding the multifunctional Smad3 adaptor protein βII-spectrin (β2SP), alone or when combined with the heterozygosity of Smad3 (Sptbn1+/−, Smad3+/−Sptbn1+/−) spontaneously develop GI cancers and hepatocellular cancer (HCC) in the C57BL/6J background [13, 21]. Furthermore, haploinsufficiency phenotypes include severe fatty liver disease, fibrosis, inflammation, and early adenomas of the colon, consistent with a premalignant state. The remaining Smad3 or Smad3 wild-type allele is also lost in the late stages of tumor progression. These data suggest that gene dosage of Smad3 and Sptbn1 drives tumorigenesis and that progression involves loss of the wildtype Smad allele. Interestingly, βII-spectrin has a liver-specific role in response to the diet-induced fatty liver through interacting with the lipogenic transcription factor SREBP1 [29]. A caspase-mediated cleavage product of βII-spectrin contributes to acetaminophen toxicity in the liver [30], and mice treated with shRNA targeting β2SP show significantly reduced acetaminophen-induced hepatotoxicity [30]. Furthermore, liver-specific knockout of Sptbn1 (βS-LSKO) or siRNA targeting βII-spectrin reduces diet-induced Metabolic Dysfunction Associated Steatohepatitis (MASH) and protects against HCC [29], suggesting that targeting βII-spectrin in experimental studies might provide new insights into MASH and steatosis-associated HCC. Moreover, these mouse models and mechanistic studies provide an ongoing approach towards identifying biologically important biomarkers that could stratify risk for HCC.

Alterations in the TGF-β signaling pathway could reflect a continuum of fibrosis to cirrhosis to cancer in the liver. Thus, we hypothesize that the TGF-β pathway-enriched biomarkers may serve as biomarkers in the evolution of HCC and stratify patients at risk for HCC. In addition, we hypothesized that the integrated animal model-to-human studies program would yield new TGF-β driven mechanistic biomarkers that could be valuable in yielding additional biomarkers that could stratify the risk of HCC. Here, we identify a 31-biomarker panel, which we refer to as TGF-β based Protein Markers for Early Detection of HCC (TPEARLE), that could potentially stratify cirrhotic patients at high risk for developing HCC. In addition, we found two other biologically relevant markers, Myostatin and Pyruvate Kinase M2 (PKM2), which were significantly associated with HCC.


Selection of potential HCC protein markers

Our cohort for the main analysis included 216 patients with cirrhosis, of whom 59 had HCC (18 early stage defined as one tumor <5 cm or </= 3 tumors, each <3 cm, with no extrahepatic spread). Demographics and lab values differed between groups (Table 1). Patients with HCC were older, more often male, more often Black or Asian, more likely to be Child-Pugh class A, and with a cirrhosis etiology more likely to be HCV or HBV. The HCC patients had higher levels of AFP and albumin and lower levels of bilirubin and Alkaline phosphatase (ALP) than those without HCC.

Our examination of 108 proteins from the TGF-β family members [4, 19] and some known biomarkers for HCC [7, 3134], after unsupervised clustering of markers and patients using our batch-1 data, revealed a promising proteomic pattern that differentiates cirrhosis from HCC, using an early cohort of N = 170 (Figure 1A). Markers positively associated with HCC included INHBC, HN1, DKK1, RAN, SMAD3, SMAD2, and YWHAZ. Markers negatively associated with HCC included COL3a1, BMP7, BMP4, BMP6, LUM, MMP19, TGFBR3, and MMP7. Notably, 20 patients with cirrhosis presented an HCC-like pattern, suggesting that they may be a group at high risk for developing HCC. We refer to the proteomic signature differentiating HCC from non-HCC in patients with cirrhosis as TGF-β based Protein Markers for Early Detection of HCC (TPEARLE). PKM2 and myostatin (Figure 1B), markers known to play a role in HCC, were also significantly associated with HCC in a later cohort of N = 467. Including the batch 1 and 2 cohort, the pattern was similar (Figure 2). SMAD2, SMAD3, INHBC, MSTN, HN1, PKM2, and DKK were positive HCC markers (dark red indicating the highest marker quintile, blue indicating the lowest quintile). BMP7, BMP4, COL5a1, MMP19, TGFBR3, MMP7, LUM were negatively associated with HCC.

We and others have focused on refining untargeted circulating proteins unique to HCC as potential biomarkers. In preliminary studies, we identified circulating proteins that are components of the TGF-β signaling pathway in HCC.

Additional markers PKM2 and myostatin

Lastly, we explored potential new biomarkers to continue examining our mouse models and query our mechanistic studies. We observed that core metabolic pathways, such as glycolysis and gluconeogenesis, were significantly altered, and critical metabolites, such as serine, pyruvate, aconitate, and malate, were substantially increased in obese mouse tissues -which we refer to as Obese MT [29, 35] -that were corrected with blocking SPTBN1 (βS-LSKO) [29]. Transcriptomic analysis shows increased expression of glucose-related metabolic genes involved in the TCA cycle, pyruvate metabolism, pentose phosphate pathway, and serine metabolism in the liver tissues of obese mice, compared to the normal WT group (Supplementary Figure 1). We observed that a key enzyme in glycolysis- PKM2 expression is substantially increased in DEN-induced HCC liver tissues. In contrast, PKM2 expression is significantly reduced in liver tissues of βS-LSKO mice fed on normal chow diet (NC) and Western diet (WD) (Figure 3A). Furthermore, we observed increased PKM2 expression in hepatic Kupffer cells of the MASH-associated HCC mouse model, and the βS-LSKO mice showed a significant reduction in PKM2 expression in these cells (Figure 3A). Indeed, we observed that PKM2 levels are significantly increased in HCC patients compared to cirrhotic individuals (Figure 3B).

Another critical characteristic of cirrhosis and HCC, as well as in our mouse models of HCC, is severe muscle loss [3638]. We observed that similar to the human fatty liver disease, our obese mice show marked sarcopenia that is corrected in the liver-specific knockout of SPTBN1. Because of its biological functional role as a TGF-β superfamily member and its role in muscle loss, we also examined Myostatin levels in our samples. Compared to the WT mice, the obese mice (obese MT) showed significant muscle loss (Figure 3C) with increased Myostatin expression (Figure 3D, 3E), and decreased levels of p-RPS6/total-RPS6 ratio (indicating mTOR signaling inhibition, Figure 3F) in striated muscle tissue. These changes are associated with the development of sarcopenia during liver cirrhosis and cancer [39, 40]. In contrast, we observed that Myostatin expression is decreased in βS-LSKO tissues (Figure 3D, 3E).Using the 108 markers as well as PKM2 and myostatin (MSTN), in batch 1 and 2 cohorts, we found 31 markers with false discovery rate (FDR) p values < 0.10 from univariable Kruskal-Wallis tests for the associations of markers with HCC vs. cirrhosis, which we retained for further analysis (Table 2). Markers positively associated with HCC included PKM2, MSTN, GPC3, INHBC, HN1, DKK1, and SMAD3 (Figure 4A). Markers negatively associated with HCC included COL28A1, TGFBR3, BMP7, LUM, and MMP7 (Figure 4B).

Potential confounding variables that were significantly associated with HCC in the clinical multivariable model included age (p = .002), HCV (p = .008), HBV (p = .0005), log10AFP (p < .0001), log10ALT (p = .04), site (p < .0001), and batch (p = .003). After adjusting for these, we found 16 markers with significant independent associations with HCC (Table 2). Based on adjusted OR (aOR), the strongest positive HCC markers adjusted for confounding factors were: DKK1 (aOR 1.73), HN1 (aOR 1.68), MSTN (aOR 1.64), SMAD3 (aOR 1.60), and PKM2 (aOR 1.58). The strongest negative markers for HCC were COL28a1 (aOR 0.48), BMP7 (aOR 0.52), MMP7 (aOR 0.52), TGFBR3 (aOR 0.53), and LUM (aOR 0.53) (Table 2). Other significant markers included MMP19, MMP2, BMP6, BMP4, RSPO3, and COL15a1.


A fundamental hypothesis we sought to test was whether biomarkers from the TGF-β signaling pathway might be of novel value in risk stratification of HCC in the clinical cirrhotic setting. We observe that in this cohort of patients with cirrhosis, alterations in the serum levels of members of TGF-β protein family may be promising markers for early detection of HCC [41]. We found 16 such markers significantly associated with HCC after adjusting for potentially confounding variables. These proteins are promising as potential additional markers that could be used in the surveillance of patients at risk for HCC along with AFP and ultrasound. These markers can be compiled with other published markers to develop a panel of markers predictive of HCC using small volumes of serum. TGF-β signaling initiates context-dependent responses and is critical for protein homeostasis (proteostasis) [22, 23].

Metabolic alterations are a major contributing factor in the progression of diseases in fatty liver disease to cancer. Earlier, we had observed a liver-specific role of the Smad3 adaptor βII-spectrin in response to the diet-induced fatty liver through its interactions with the lipogenic transcription factor SREBP1 [29]. Moreover, liver-specific knockout of Sptbn1 (βS-LSKO) or siRNA targeting βII-spectrin reduces diet-induced MASH and protects against HCC [29], suggesting that targeting βII-spectrin may be a viable strategy to better understand mechanisms, as well as identify new biomarkers of MASH and steatosis-associated HCC and effects on restoration of metabolic function [29]. Therefore, we examined multiple enzymes involved in glycolysis. Interestingly, we observed that PKM2 expression is substantially increased in DEN-induced HCC liver tissues. In contrast, PKM2 expression is significantly reduced in liver tissues of βS-LSKO mice. Enzymes such as PKM2 that catalyze the last step of glycolysis to generate pyruvate and ATP serve as “gatekeepers” of metabolic flux from glycolysis to the tricarboxylic acid cycle (TCA cycle). Importantly, mutually exclusive alternative splicing of exons 9 and 10 in the PKM gene leads to the expression of either the PKM1 or PKM2 isoform. PKM2 is overexpressed in human HCC and also in chemically induced mouse HCC [42, 43]. These studies suggest that PKM2 could represent an additional TGF-β driven mechanistic marker for risk stratification of HCC.

Cirrhosis and HCC are characterized by severe muscle loss- sarcopenia [3638]. Myostatin, a TGF-β superfamily member expressed primarily in skeletal muscle, inhibits muscle growth and causes sarcopenia in cirrhosis [39, 4449]. Muscle-specific deletion of the myostatin (myostatinMSD) gene prevents hepatic steatosis with high-fat diet (HFD) in mice [50]. Because of its biological functional role as a TGF-β superfamily member, we also examined Myostatin levels in our samples. We also observed that, similar to human fatty liver disease, obese mice show marked sarcopenia with increased Myostatin levels that are corrected in the liver-specific knockout of SPTBN1. These studies reflect that these are biologically functional circuits that could provide new serum markers for risk stratification of HCC.

A limitation of this study is the size of the patient cohort and the number of cases with HCC. The major strengths of this study, however, include the use of a set of markers with known mechanistic roles in HCC tumor formation and an approach that examines their association with HCC by first using unsupervised machine learning methods, and then by use of using methods that account for skewed marker distributions, multiple testing, and the potential for confounding by other variables that are associated with HCC. These results reflect the discovery of new and potentially important blood-based markers that show promise in their ability to provide early detection and risk stratification of HCC.


Cohort analysis

The cohort of cirrhosis patients (216 patients) was recruited from five medical centers in the United States (George Washington University-GW, University of Maryland -UM, University of Hawaii -UH, Northwell Health -NH and University of California at Los Angeles -UCLA), using two batches of samples, designated batch 1 and 2. Cirrhosis was diagnosed using either FibroSure/FibroTest >0.74, AST to Platelet Ratio Index (APRI) >2 [51], or Fibrosis-4 (FIB-4) >3.25 [52], or with histological/ imaging evidence of cirrhosis, or vibration controlled transient elastography (VCTE) >12.5 kPA or clinical evidence of portal hypertension. HCC diagnosis was confirmed by histology or using contrast-enhanced imaging (CT or MRI) showing Liver Reporting and Data System (LI-RADS®) 5/Organ Procurement Transplant Network 5 lesions ( All HCC patients were treatment naïve. Controls were designated as those with cirrhosis, without HCC at enrollment, and on an imaging study 6 months after enrollment. Because we focused on detecting models that might be clinically useful in a patient population with cirrhosis under surveillance for HCC, we chose not to include a normal control group. In support of this approach, several studies have shown that models using normal controls are over-optimistic, can be misleading, and require different controls to define a more realistic clinical setting [5355].

We compared patients with vs. without HCC demographic and clinical variables using chi-square for categorical variables and either a t-test or Kruskal-Wallis test for continuous variables, depending on whether they were normally distributed or skewed.

SomaScan analysis

We collaborated with SomaLogic, a Biotech company focused on precision health, for serum proteome analyses [5659]. SomaScan assay utilizes a small volume of 130 µl serum or plasma with a readout of ~5000 proteins, potentially including multiple liver cancer biomarkers. Samples were analyzed in two batches. In this study, we focused on 108 markers that included members of the TGF-β pathway. Raw proteomic data were used in units of relative fluorescence units. Only one aptamer per protein was used for proteins with multiple aptamers, selecting the one with the strongest association with HCC.

We used unsupervised k-means clustering of subjects and proteins to produce a heat map showing clusters of patients with similar patterns of protein levels, using the R pheatmap and ComplexHeatmap packages. When the patient clustering corresponded to those patients with vs. without HCC, it suggested that the protein pattern differentiating this patient group was predictive of HCC. This group of proteins was then explored further, as described below.

To examine the distribution of marker levels in patients with vs. without HCC, we examined boxplots. Since most markers had skewed distributions, we used the non-parametric Kruskal-Wallis test to associate marker level with HCC status. Then we used false-discovery-rate (FDR) adjusted p values to account for multiple testing. Markers with FDR p < 0.10 were considered for further study. We coded these markers into quintiles to eliminate skewness and to examine the functional relationship of marker level with the incidence of HCC using chi-square test.

We next examined each marker’s association with HCC after adjusting for the demographic and clinical variables that differed between HCC and non-HCC patients. These baseline variables could have acted as confounding variables, leading to differences in marker levels due to the baseline variables, rather than HCC. Coding the proteomic markers as quintiles also had the effect of equating the variances of the markers so that their effect sizes (i.e., their adjusted odds ratio with HCC after adjusting for covariates in the multivariate logistic regression model) could be compared with each other to determine which had the strongest effect. To examine the independent association of each marker with HCC after adjusting for potential confounding variables, we used one multiple logistic regression model per marker with demographic and clinical variables that had univariable p < 0.10 used as covariates. The markers that had significant independent effects in this model were considered to be potentially useful markers that deserved further study. We used SAS (version 9.4, Cary, NC) and R (ComplexHeatmap and pheatmap packages) for data analysis, with p < .05 considered significant, unless otherwise specified.

Animal models

All mice utilized in the study originated from a C57BL/6 background and are considered suitable for studies in obesity and cancer [29, 6062]. To induce liver steatosis through a high-fat diet (HFD), male and female mice aged 10 to 12 weeks were subjected to either the normal chow diet (NC) or HFD (ENVIGO, catalog no. TD.06414) for durations spanning 12 to 28 weeks. In the case of MASH-associated HCC induced by a combination of Diethylnitrosamine (DEN) and a Western diet (WD), mice were administered weekly injections of DEN (50 mg/kg) for two consecutive weeks, followed by a one-week pause, amounting to 6 cycles over an 18-week period. The WD diet was introduced one week after the initial DEN injection, continuing for 21 to 22 weeks in total.

Immunoblotting analysis

For Western blot analysis of Myostatin levels in muscle, skeletal muscle from WT, obese mouse (Obese MT), and βS-LSKO mouse tissues were homogenized with 250 µL of Tissue cell lysis buffer (GOLDBIO, Cat. GB-181) containing protease inhibitor and phosphatase inhibitor cocktails (Millipore-Sigma) on ice using a pellet pestle motor homogenizer. After homogenization, samples were gently rotated on a rotator at 4°C for 40 min. Samples were vortexed and then centrifuged at 16000 g in a microcentrifuge at 4°C for 20 min to pellet insoluble material. Protein concentration was determined using the Bio-Rad protein assay dye reagent, following the manufacturer’s instructions (DC protein assay reagent A and B; #5000113, #5000114). Fifty micrograms of protein were then separated by SDS-PAGE and transferred to nitrocellulose Membrane (Bio-Rad, #1620112). Immunoblotting was performed using the following antibodies: Myostatin (Proteintech, Cat. 19142-1-AP) and GAPDH (Santa Cruz, sc-32233). ImageJ software (version 1.53t, National Institutes of Health) evaluated and quantified the protein bands.

Immunohistochemical analysis

For immunohistochemical analyses, mouse tissues were fixed in 10% formalin and embedded in paraffin (FFPE) following standard procedures. FFPE sections were deparaffinized, hydrated, and pretreated for antigen retrieval in 10 mM citrate buffer, pH 6.0 (Vector Laboratories, Cat. H-330). Endogenous peroxidase activity was blocked with 3% H2O2 solution for 30 min. Nonspecific binding was blocked with 10% normal goat serum (Vector Laboratories) in 0.1% Triton X-100 in a 1× PBS solution (PBST) for 1 hr at room temperature in a humidifying chamber. Sections were labeled with antibodies specific to PKM2 (Cell signaling, #4053), Myostatin (Proteintech, Cat. 19142-1-AP), phospho-S6 Ribosomal (Cell signaling, #2215S), and S6 Ribosomal Protein (Cell signaling, #2217T) overnight at 4°C in the humidifying chamber. After being washed with 1× PBST three times for 3 min, sections were incubated with biotinylated secondary antibodies (Vectastain Elite ABC kit, PK6106) for 30 min at room temperature in the humidifying chamber. Finally, the sections were stained using DAB kit (Dako EnVision Dual link System HRT, K4065) for signal amplification and detection.


HCC: hepatocellular carcinoma; TGF-β: transforming growth factor β; PKM2: Pyruvate Kinase M2; βS-LSKO: liver-specific knockout of Sptbn1; TPEARLE: TGF-β based Protein Markers for Early Detection of HCC; NC: normal chow diet; HFD: high-fat diet; WD: Western diet; DEN: Diethylnitrosamine.

Author contributions

LM and RA conceived the study and hypotheses, designed the experiments, supervised the research, and wrote the manuscript. XX and RA performed the statistical analyses. KB, KO and SJ performed immunohistochemical labeling and provided the essential input to the manuscript. XY and KB did Western blotting. PL provided intellectual input and analyzed the immunohistochemical labeling slides. RA, HY, LW, and KS performed pre-processing of patient data. SS, CB, RD, NC, GS, FC, JC, HH, SD, AK, AH critically reviewed and revised the manuscript. LM and RLA have full access to all data in the study and take responsibility for data integrity and the accuracy of data analysis.


We acknowledge all the patients who participated in this study.


The authors declare that there are no conflicts of interest.


This study was approved by the local Ethics Committee; the principles of the Declaration of Helsinki were followed. The protocol/approval number for the human study is Northwell Health IRB #20-0150.

All animal experiments were performed according to the guidelines for the care and use of laboratory animals and were approved by the Institutional Research Ethics Committee of the Feinstein Institute of Medical Research (The protocol/approval number #2021-020).


All patients gave their informed consent for the use of clinical/omics data and samples for research purposes.


This research work was supported by National Institute of Health grants R01AA023146 (L. Mishra), R01CA236591 (L. Mishra), and U01 CA230690 (L. Mishra).

Table 1
Patient characteristics in full sample and stratified by presence of HCC
Table 1: Patient characteristics in full sample and stratified by presence of HCC
Table 2
Univariable and multivariable p values for 31 proteomic markers
Table 2: Univariable and multivariable p values for 31 proteomic markers
Figure 1
TGF-β pathway-related serum markers stratify HCC risk in cirrhotic patients (TPEARLE).
(A) Unsupervised clustering analysis of TGF-β pathway-associated proteins in cirrhosis (n = 104) and HCC patients (n = 66). (B) In a larger cohort, serum Myostatin is higher in HCC (n = 157), compared to cirrhotic patients (n = 380).
Figure 1: TGF-β pathway-related serum markers stratify HCC risk in cirrhotic patients (TPEARLE).
Figure 2
Heatmap showing TGF-β pathway-related serum proteomic markers associated with HCC in cirrhotic patients.
Unsupervised k-means clustering of markers and patients N = 216 patients, 59 with HCC, using batch 1 and 2 data. Markers are scored as quintiles (0–4). Grey = missing. There are distinct HCC-like and non-HCC-like patterns of marker levels and batch effects.
Figure 2: Heatmap showing TGF-β pathway-related serum proteomic markers associated with HCC in cirrhotic patients.
Figure 3
Expression of TGF-β driven mechanistic marker Myostatin and PKM2 in MASH/HCC mouse model.
(A) Increased PKM2 expression in hepatic Kupffer cells of liver tissue from the MASH/HCC mouse model (Western diet (WD) + Diethylnitrosamine (DEN)), and (B) in human HCC. PKM2 expression is significantly reduced in βS-LSKO mice, which block HCC. Arrows indicate positive staining. (CF) Muscle loss and increased expression of Myostatin in muscle tissue from obese mouse tissues (Obese MT): (C) muscle/body weight ratio decreased, *p < 0.05; (D, E) Increased Myostatin; (F) decreased p-RPS6/total-RPS6 ratio in muscle tissue from obese mouse.
Figure 3: Expression of TGF-β driven mechanistic marker Myostatin and PKM2 in MASH/HCC mouse model.
Figure 4
Association of selected marker levels by HCC status.
(A) Markers positively associated with HCC. (B) Markers are negatively associated with HCC. In box plots, the black square is the mean, the horizontal line is the median, the box shows the intra-quartile range, and ‘whiskers’ show the range. In the quintile chart, error bars show the 95% confidence interval. aFDR-adjusted Kruskal-Wallis p value, bchi-square p value.
Figure 4: Association of selected marker levels by HCC status.
  • 1. Siegel RL, Miller KD, Wagle NS, Jemal A. Cancer statistics, 2023. CA Cancer J Clin. 2023; 73:17–48. [Pubmed]
  • 2. Mattiuzzi C, Lippi G. Current Cancer Epidemiology. J Epidemiol Glob Health. 2019; 9:217–22. [Pubmed]
  • 3. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018; 68:394–424. [Pubmed]
  • 4. Zaidi S, Amdur R, Xiang X, Yu H, Wong LL, Rao S, He AR, Amin K, Zaheer D, Narayan RK, Satapathy SK, Latham PS, Shetty K, et al. Using quantitative immunohistochemistry in patients at high risk for hepatocellular cancer. Genes Cancer. 2022; 13:9–20. [Pubmed]
  • 5. Yarchoan M, Agarwal P, Villanueva A, Rao S, Dawson LA, Llovet JM, Finn RS, Groopman JD, El-Serag HB, Monga SP, Wang XW, Karin M, Schwartz RE, et al. Recent Developments and Therapeutic Strategies against Hepatocellular Carcinoma. Cancer Res. 2019; 79:4326–30. [Pubmed]
  • 6. Tzartzeva K, Obi J, Rich NE, Parikh ND, Marrero JA, Yopp A, Waljee AK, Singal AG. Surveillance Imaging and Alpha Fetoprotein for Early Detection of Hepatocellular Carcinoma in Patients With Cirrhosis: A Meta-analysis. Gastroenterology. 2018; 154:1706–18.e1. [Pubmed]
  • 7. Fujiwara N, Kobayashi M, Fobar AJ, Hoshida A, Marquez CA, Koneru B, Panda G, Taguri M, Qian T, Raman I, Li QZ, Hoshida H, Sezaki H, et al. A blood-based prognostic liver secretome signature and long-term hepatocellular carcinoma risk in advanced liver fibrosis. Med. 2021; 2:836–50.e10. [Pubmed]
  • 8. Yang JD, Addissie BD, Mara KC, Harmsen WS, Dai J, Zhang N, Wongjarupong N, Ali HM, Ali HA, Hassan FA, Lavu S, Cvinar JL, Giama NH, et al. GALAD Score for Hepatocellular Carcinoma Detection in Comparison with Liver Ultrasound and Proposal of GALADUS Score. Cancer Epidemiol Biomarkers Prev. 2019; 28:531–38. [Pubmed]
  • 9. Chan LK, Tsui YM, Ho DW, Ng IO. Cellular heterogeneity and plasticity in liver cancer. Semin Cancer Biol. 2022; 82:134–49. [Pubmed]
  • 10. Lee DH, Lee JM. Primary malignant tumours in the non-cirrhotic liver. Eur J Radiol. 2017; 95:349–61. [Pubmed]
  • 11. Sanyal AJ, Van Natta ML, Clark J, Neuschwander-Tetri BA, Diehl A, Dasarathy S, Loomba R, Chalasani N, Kowdley K, Hameed B, Wilson LA, Yates KP, Belt P, et al. Prospective Study of Outcomes in Adults with Nonalcoholic Fatty Liver Disease. N Engl J Med. 2021; 385:1559–69. [Pubmed]
  • 12. Hagiwara S, Nishida N, Ueshima K, Minami Y, Komeda Y, Aoki T, Takita M, Morita M, Chishina H, Yoshida A, Ida H, Kudo M. Accumulation of Genetic and Epigenetic Alterations in the Background Liver and Emergence of Hepatocellular Carcinoma in Patients with Non-Alcoholic Fatty Liver Disease. Cells. 2021; 10:3257. [Pubmed]
  • 13. Chen J, Yao ZX, Chen JS, Gi YJ, Muñoz NM, Kundra S, Herlong HF, Jeong YS, Goltsov A, Ohshiro K, Mistry NA, Zhang J, Su X, et al. TGF-β/β2-spectrin/CTCF-regulated tumor suppression in human stem cell disorder Beckwith-Wiedemann syndrome. J Clin Invest. 2016; 126:527–42. [Pubmed]
  • 14. Zhu R, Baker SS, Moylan CA, Abdelmalek MF, Guy CD, Zamboni F, Wu D, Lin W, Liu W, Baker RD, Govindarajan S, Cao Z, Farci P, et al. Systematic transcriptome analysis reveals elevated expression of alcohol-metabolizing genes in NAFLD livers. J Pathol. 2016; 238:531–42. [Pubmed]
  • 15. Li S, Mao M. Next generation sequencing reveals genetic landscape of hepatocellular carcinomas. Cancer Lett. 2013; 340:247–53. [Pubmed]
  • 16. Tang Y, Kitisin K, Jogunoori W, Li C, Deng CX, Mueller SC, Ressom HW, Rashid A, He AR, Mendelson JS, Jessup JM, Shetty K, Zasloff M, et al. Progenitor/stem cells give rise to liver cancer due to aberrant TGF-beta and IL-6 signaling. Proc Natl Acad Sci U S A. 2008; 105:2445–50. [Pubmed]
  • 17. Heldin CH, Moustakas A. Signaling Receptors for TGF-β Family Members. Cold Spring Harb Perspect Biol. 2016; 8:a022053. [Pubmed]
  • 18. Rao S, Zaidi S, Banerjee J, Jogunoori W, Sebastian R, Mishra B, Nguyen BN, Wu RC, White J, Deng C, Amdur R, Li S, Mishra L. Transforming growth factor-β in liver cancer stem cells and regeneration. Hepatol Commun. 2017; 1:477–93. [Pubmed]
  • 19. Korkut A, Zaidi S, Kanchi RS, Rao S, Gough NR, Schultz A, Li X, Lorenzi PL, Berger AC, Robertson G, Kwong LN, Datto M, Roszik J, et al. A Pan-Cancer Analysis Reveals High-Frequency Genetic Alterations in Mediators of Signaling by the TGF-β Superfamily. Cell Syst. 2018; 7:422–37.e7. [Pubmed]
  • 20. Chen J, Zaidi S, Rao S, Chen JS, Phan L, Farci P, Su X, Shetty K, White J, Zamboni F, Wu X, Rashid A, Pattabiraman N, et al. Analysis of Genomes and Transcriptomes of Hepatocellular Carcinomas Identifies Mutations and Gene Expression Changes in the Transforming Growth Factor-β Pathway. Gastroenterology. 2018; 154:195–210. [Pubmed]
  • 21. Yao ZX, Jogunoori W, Choufani S, Rashid A, Blake T, Yao W, Kreishman P, Amin R, Sidawy AA, Evans SR, Finegold M, Reddy EP, Mishra B, et al. Epigenetic silencing of beta-spectrin, a TGF-beta signaling/scaffolding protein in a human cancer stem cell disorder: Beckwith-Wiedemann syndrome. J Biol Chem. 2010; 285:36112–20. [Pubmed]
  • 22. Morikawa M, Derynck R, Miyazono K. TGF-β and the TGF-β Family: Context-Dependent Roles in Cell and Tissue Physiology. Cold Spring Harb Perspect Biol. 2016; 8:a021873. [Pubmed]
  • 23. David CJ, Massagué J. Contextual determinants of TGFβ action in development, immunity and cancer. Nat Rev Mol Cell Biol. 2018; 19:419–35. [Pubmed]
  • 24. Baek HJ, Pishvaian MJ, Tang Y, Kim TH, Yang S, Zouhairi ME, Mendelson J, Shetty K, Kallakury B, Berry DL, Shin KH, Mishra B, Reddy EP, et al. Transforming growth factor-β adaptor, β2-spectrin, modulates cyclin dependent kinase 4 to reduce development of hepatocellular cancer. Hepatology. 2011; 53:1676–84. [Pubmed]
  • 25. Matsuura I, Denissova NG, Wang G, He D, Long J, Liu F. Cyclin-dependent kinases regulate the antiproliferative function of Smads. Nature. 2004; 430:226–31. [Pubmed]
  • 26. Liu S, Huang F, Ru G, Wang Y, Zhang B, Chen X, Chu L. Mouse Models of Hepatocellular Carcinoma: Classification, Advancement, and Application. Front Oncol. 2022; 12:902820. [Pubmed]
  • 27. Gough NR, Xiang X, Mishra L. TGF-β Signaling in Liver, Pancreas, and Gastrointestinal Diseases and Cancer. Gastroenterology. 2021; 161:434–52.e15. [Pubmed]
  • 28. Baek HJ, Lim SC, Kitisin K, Jogunoori W, Tang Y, Marshall MB, Mishra B, Kim TH, Cho KH, Kim SS, Mishra L. Hepatocellular cancer arises from loss of transforming growth factor beta signaling adaptor protein embryonic liver fodrin through abnormal angiogenesis. Hepatology. 2008; 48:1128–37. [Pubmed]
  • 29. Rao S, Yang X, Ohshiro K, Zaidi S, Wang Z, Shetty K, Xiang X, Hassan MI, Mohammad T, Latham PS, Nguyen BN, Wong L, Yu H, et al. β2-spectrin (SPTBN1) as a therapeutic target for diet-induced liver disease and preventing cancer development. Sci Transl Med. 2021; 13:eabk2267. [Pubmed]
  • 30. Baek HJ, Lee YM, Kim TH, Kim JY, Park EJ, Iwabuchi K, Mishra L, Kim SS. Caspase-3/7-mediated Cleavage of β2-spectrin is Required for Acetaminophen-induced Liver Damage. Int J Biol Sci. 2016; 12:172–83. [Pubmed]
  • 31. Zhan Z, Chen Y, Duan Y, Li L, Mew K, Hu P, Ren H, Peng M. Identification of key genes, pathways and potential therapeutic agents for liver fibrosis using an integrated bioinformatics analysis. PeerJ. 2019; 7:e6645. [Pubmed]
  • 32. Singh G, Yoshida EM, Rathi S, Marquez V, Kim P, Erb SR, Salh BS. Biomarkers for hepatocellular cancer. World J Hepatol. 2020; 12:558–73. [Pubmed]
  • 33. Wang W, Wei C. Advances in the early diagnosis of hepatocellular carcinoma. Genes Dis. 2020; 7:308–19. [Pubmed]
  • 34. Beudeker BJB, Boonstra A. Circulating biomarkers for early detection of hepatocellular carcinoma. Therap Adv Gastroenterol. 2020; 13:1756284820931734. [Pubmed]
  • 35. Chang YC, Lee HL, Yang W, Hsieh ML, Liu CC, Lee TY, Huang JY, Nong JY, Li FA, Chuang HL, Ding ZZ, Su WL, Chueh LY, et al. A common East-Asian ALDH2 mutation causes metabolic disorders and the therapeutic effect of ALDH2 activators. Nat Commun. 2023; 14:5971. [Pubmed]
  • 36. Fujiwara N, Nakagawa H, Kudo Y, Tateishi R, Taguri M, Watadani T, Nakagomi R, Kondo M, Nakatsuka T, Minami T, Sato M, Uchino K, Enooku K, et al. Sarcopenia, intramuscular fat deposition, and visceral adiposity independently predict the outcomes of hepatocellular carcinoma. J Hepatol. 2015; 63:131–40. [Pubmed]
  • 37. Iritani S, Imai K, Takai K, Hanai T, Ideta T, Miyazaki T, Suetsugu A, Shiraki M, Shimizu M, Moriwaki H. Skeletal muscle depletion is an independent prognostic factor for hepatocellular carcinoma. J Gastroenterol. 2015; 50:323–32. [Pubmed]
  • 38. Tantai X, Liu Y, Yeo YH, Praktiknjo M, Mauro E, Hamaguchi Y, Engelmann C, Zhang P, Jeong JY, van Vugt JLA, Xiao H, Deng H, Gao X, et al. Effect of sarcopenia on survival in patients with cirrhosis: A meta-analysis. J Hepatol. 2022; 76:588–99. [Pubmed]
  • 39. Dasarathy S, Merli M. Sarcopenia from mechanism to diagnosis and treatment in liver disease. J Hepatol. 2016; 65:1232–44. [Pubmed]
  • 40. Dasarathy S. Myostatin and beyond in cirrhosis: all roads lead to sarcopenia. J Cachexia Sarcopenia Muscle. 2017; 8:864–69. [Pubmed]
  • 41. Lu WQ, Qiu JL, Huang ZL, Liu HY. Enhanced circulating transforming growth factor beta 1 is causally associated with an increased risk of hepatocellular carcinoma: a mendelian randomization meta-analysis. Oncotarget. 2016; 7:84695–704. [Pubmed]
  • 42. Lu DH, Lv WW, Li WX, Gao YD. High PKM2 expression is independently correlated with decreased overall survival in hepatocellular carcinoma. Oncol Lett. 2018; 16:3603–10. [Pubmed]
  • 43. Wong CC, Au SL, Tse AP, Xu IM, Lai RK, Chiu DK, Wei LL, Fan DN, Tsang FH, Lo RC, Wong CM, Ng IO. Switching of pyruvate kinase isoform L to M2 promotes metabolic reprogramming in hepatocarcinogenesis. PLoS One. 2014; 9:e115036. [Pubmed]
  • 44. McPherron AC, Lawler AM, Lee SJ. Regulation of skeletal muscle mass in mice by a new TGF-beta superfamily member. Nature. 1997; 387:83–90. [Pubmed]
  • 45. McPherron AC, Lee SJ. Double muscling in cattle due to mutations in the myostatin gene. Proc Natl Acad Sci U S A. 1997; 94:12457–61. [Pubmed]
  • 46. Tandon P, Montano-Loza AJ, Lai JC, Dasarathy S, Merli M. Sarcopenia and frailty in decompensated cirrhosis. J Hepatol. 2021 (Suppl 1); 75:S147–62. [Pubmed]
  • 47. Nishikawa H, Enomoto H, Ishii A, Iwata Y, Miyamoto Y, Ishii N, Yuri Y, Hasegawa K, Nakano C, Nishimura T, Yoh K, Aizawa N, Sakai Y, et al. Elevated serum myostatin level is associated with worse survival in patients with liver cirrhosis. J Cachexia Sarcopenia Muscle. 2017; 8:915–25. [Pubmed]
  • 48. Kim JH, Kang SH, Lee M, Youn GS, Kim TS, Jun BG, Kim MY, Kim YD, Cheon GJ, Kim DJ, Baik SK, Choi DH, Suk KT. Serum Myostatin Predicts the Risk of Hepatocellular Carcinoma in Patients with Alcoholic Cirrhosis: A Multicenter Study. Cancers (Basel). 2020; 12:3347. [Pubmed]
  • 49. García PS, Cabbabe A, Kambadur R, Nicholas G, Csete M. Brief-reports: elevated myostatin levels in patients with liver disease: a potential contributor to skeletal muscle wasting. Anesth Analg. 2010; 111:707–9. [Pubmed]
  • 50. Burgess K, Xu T, Brown R, Han B, Welle S. Effect of myostatin depletion on weight gain, hyperglycemia, and hepatic steatosis during five months of high-fat feeding in mice. PLoS One. 2011; 6:e17090. [Pubmed]
  • 51. Wai CT, Greenson JK, Fontana RJ, Kalbfleisch JD, Marrero JA, Conjeevaram HS, Lok AS. A simple noninvasive index can predict both significant fibrosis and cirrhosis in patients with chronic hepatitis C. Hepatology. 2003; 38:518–26. [Pubmed]
  • 52. Sterling RK, Lissen E, Clumeck N, Sola R, Correa MC, Montaner J, S Sulkowski M, Torriani FJ, Dieterich DT, Thomas DL, Messinger D, Nelson M, and APRICOT Clinical Investigators. Development of a simple noninvasive index to predict significant fibrosis in patients with HIV/HCV coinfection. Hepatology. 2006; 43:1317–25. [Pubmed]
  • 53. Johnson PJ, Pirrie SJ, Cox TF, Berhane S, Teng M, Palmer D, Morse J, Hull D, Patman G, Kagebayashi C, Hussain S, Graham J, Reeves H, Satomura S. The detection of hepatocellular carcinoma using a prospectively developed and validated model based on serological biomarkers. Cancer Epidemiol Biomarkers Prev. 2014; 23:144–53. [Pubmed]
  • 54. Wang M, Devarajan K, Singal AG, Marrero JA, Dai J, Feng Z, Rinaudo JA, Srivastava S, Evans A, Hann HW, Lai Y, Yang H, Block TM, Mehta A. The Doylestown Algorithm: A Test to Improve the Performance of AFP in the Detection of Hepatocellular Carcinoma. Cancer Prev Res (Phila). 2016; 9:172–79. [Pubmed]
  • 55. Mehta AS, Lau DT, Wang M, Aslam A, Nasir B, Javaid A, Poongkunran M, Block TM. Application of the Doylestown algorithm for the early detection of hepatocellular carcinoma. PLoS One. 2018; 13:e0203149. [Pubmed]
  • 56. Kraemer S, Vaught JD, Bock C, Gold L, Katilius E, Keeney TR, Kim N, Saccomano NA, Wilcox SK, Zichi D, Sanders GM. From SOMAmer-based biomarker discovery to diagnostic and clinical applications: a SOMAmer-based, streamlined multiplex proteomic assay. PLoS One. 2011; 6:e26332. [Pubmed]
  • 57. Gold L, Ayers D, Bertino J, Bock C, Bock A, Brody EN, Carter J, Dalby AB, Eaton BE, Fitzwater T, Flather D, Forbes A, Foreman T, et al. Aptamer-based multiplexed proteomic technology for biomarker discovery. PLoS One. 2010; 5:e15004. [Pubmed]
  • 58. Lehallier B, Gate D, Schaum N, Nanasi T, Lee SE, Yousef H, Moran Losada P, Berdnik D, Keller A, Verghese J, Sathyan S, Franceschi C, Milman S, et al. Undulating changes in human plasma proteome profiles across the lifespan. Nat Med. 2019; 25:1843–50. [Pubmed]
  • 59. Williams SA, Kivimaki M, Langenberg C, Hingorani AD, Casas JP, Bouchard C, Jonasson C, Sarzynski MA, Shipley MJ, Alexander L, Ash J, Bauer T, Chadwick J, et al. Plasma protein patterns as comprehensive indicators of health. Nat Med. 2019; 25:1851–57. [Pubmed]
  • 60. Lang P, Hasselwander S, Li H, Xia N. Effects of different diets used in diet-induced obesity models on insulin resistance and vascular dysfunction in C57BL/6 mice. Sci Rep. 2019; 9:19556. [Pubmed]
  • 61. Zhang BT, Xu JY, Wang W, Zeng Y, Jiang J. Obesity and cancer: Mouse models used in studies. Front Oncol. 2023; 13:1125178. [Pubmed]
  • 62. Grohmann M, Wiede F, Dodd GT, Gurzov EN, Ooi GJ, Butt T, Rasmiena AA, Kaur S, Gulati T, Goh PK, Treloar AE, Archer S, Brown WA, et al. Obesity Drives STAT-1-Dependent NASH and STAT-3-Dependent HCC. Cell. 2018; 175:1289–306.e20. [Pubmed]