- Open Access
TGF beta1 and related-Smads contribute to pulmonary metastasis of hepatocellular carcinoma in mice model
Journal of Experimental & Clinical Cancer Research volume 31, Article number: 93 (2012)
Recent studies indicate that Transforming Growth Factor beta (TGF β) correlated with pulmonary metastasis of cancers. However, the correlation between TGF β and pulmonary metastasis of hepatocellular carcinoma (HCC) is till unknown.
We detected the in vitro and in vivo expression levels of TGF β1/Smads by Real-time PCR and Western blot in MHCC97-H and MHCC97–L cell lines, which are HCC cell lines and have higher and lower pulmonary metastatic potential respectively.
TGF β1 mRNA level in MHCC97-L tumors were higher than that in MHCC97-H tumors, (2.81±1.61 vs. 1.24±0.96, P=0.002), TGF β1 protein level in MHCC97-L tumors were also higher than that in MHCC97-H tumors (1.37±0.95 vs. 0.32±0.22, P<0.001). In addition, the TGF β1 mRNA level positively correlated with pulmonary metastasis, and the relations between TGF β1 and Smads were also found (R2=0.12 and 0.40, respectively).
Our results suggest that TGF β/ Smads promote pulmonary metastasis of HCC.
Hepatocellular carcinoma (HCC) is one of the most common cancers in the world. The overall five-year survival rate following resection has remained as poor as 35–50%[1–3]. The extremely poor prognosis of HCC is largely the result of a high rate of recurrence after surgery and of metastasis[4, 5]. Lung is the most common site for extrahepatic recurrence of HCC. The incidence of pulmonary metastasis after hepatic resection for HCC ranges from 37% to 58%. Therefore, to reduce the pulmonary metastasis could ameliorate the prognosis of HCC.
Transforming growth factor beta (TGF β) is a known regulator of epithelial cell, autonomous tumor initiation, progression and metastasis[7–9]. There are three kinds of molecular in TGF β family, while, TGFβ1 is predominantly and importantly expressed in liver cells, whereas two other members, TGFβ2 and TGFβ3, are present in a little amounts and its roles are even ignored in many studies. Signaling of TGF β1 play a role mainly through Smad proteins. Recently, a report indicates that transient exposure of breast cancer cells to TGF β which produced in the primary tumor microenvironment promotes cancer cells to extravagate from blood vessels and entry into the lung by upregulation of the adipokine angiopoietin-like 4.
In HCC, TGF β is a useful serologic marker for diagnosis because it shows higher sensitivity than AFP in earlier stage of cancer. In addition, the role of TGF β1 in HCC metastasis is emphasized. In a study by Giannelli et al. Laminin-5 (Ln-5) and TGF β1 cooperatively induce epithelial mesenchymal transition (EMT) and cancer invasion in HCC. However, although a multitude of studies have presented evidence for TGF β changes in HCC tumors, the direction of the changes is not always consistent. In several studies, TGF β1 levels are demonstrated to be lower[16, 17], while, in other studies, the levels are demonstrated to be higher versus healthy individuals[18, 19].
In this study, by comparing the different expression of TGF β/Smads in HCC cell lines, we tried to investigate the correlation between TGF β/Smads levels and potential of pulmonary metastasis in HCC.
Materials and methods
MHCC97-L and MHCC97-H, were human HCC cell lines, and which have a lower and higher metastatic potential respectively. These cell lines were clonally selected from the same parent cell lines, MHCC97, they have an identical genetic background[20, 21]. Both cell lines were cultured in high glucose Dulbecco’s modified Eagle's medium (H-DMEM, Gibco) and supplemented with 10% fetal calf serum (Gibco) at 37°C in a humidified incubator that contained 5% CO2.
31 samples and observed data were selected randomly from our previous experiment, which were tissues of MHCC97-H models (n=20) and MHCC97-L models (n=11). The models were established as follow: 6×106 MHCC97-H and 6×106 MHCC97-L cells were inoculated subcutaneously into the right side backs of the nude mice (average weight 25g). After tumor formed, the tumor size was estimated according to the formula: volume (mm3) = 0.5 a2×b, in which “a” is the major diameter of tumor and “b” is the minor diameter perpendicular to the major one . According to our experience, to guarantee enough tumor size and pulmonary metastasis, the MHCC97-L models were feed longer (40days) than MHCC97-H models (35days). In the end of feeding, animals were sacrificed. The tumor and lung tissues were removed and partly cryopreserved in -70°C for real-time PCR analysis, and partly paraffin embedded for immunohistochemstry or H&E (hematoxylin and eosin) staining.
These experiments were approved by the Shanghai Medical Experimental Animal Care Commission, and were in accordance with the Helsinki Declaration of 1975.
Analysis of pulmonary metastasis
Each lung tissues were sliced for 20 sections with 5μm in thickness, and 50μm interval between two successive sections. After stained with HE, sections were independently observed under microscopic to evaluate pulmonary metastasis by two pathologists.
RNA extraction and Real-time PCR
Total RNA of MHCC97-H, MHCC97-L cell lines and tumor tissues were extracted by TRIZOL Reagent (Invitrogen corp, USA) according instruction of the product. Real-time RT-PCR analysis was performed to identify the expression level of TGF β1, smad2 and smad7 by using SYBR Green mix(ToYoBo Co, Japan). The primers were designed by software (premier premier 5.0) as follow: TGF β1 (sense 5′ GGCGATACCTCAGCAACCG 3′; antisense, 5′ CTAAGGCGAAAGCCCTCAAT 3′), Smad2 (sense, 5′ TACTACTCTTTCCCTGT 3′; antisense, 5′ TTCTTGTCATTTCTACCG 3′), Smad7 (sense, 5′ CAACCGCAGCAGTTACCC 3′; antisense, 5′ CGAAAGCCTTGATGGAGA 3′), β-actins (sense, 5′ -TCGTGCGTGACATTAAGGAG-3′; antisense, 5′ - ATGCCAGGGTACATGGTAAT-3′). Amplification conditions were: 95°C for 9 min, followed by 45 cycles of 95°C for 30s, 57°C for 30s and 72°C for 15s, and followed by an extension at 72°C for 5 min. β-actins was used as a control for the presence of amplifiable cDNA. The mRNA expression level was assessed by 2-△△Ct in brief, the Ct value for target gene was subtracted from the Ct value of β-actins to yield a △Ct value. The average △Ct was calculated for the control group and this value was subtracted from the △Ct of all other samples (including the control group). This resulted in a △△Ct value for all samples which was then used to calculate the fold-induction of mRNA expression of target gene using the formula 2-△△Ct, as recommended by the manufacturer (Bio-Rad, Hercules, CA, USA). In this study, we used MHCC97-H model samples as control group. The detection about mRNA expression in MHCC97-H and MHCC97-L cell lines was repeated for 10 times.
Protein extraction and western blot analysis
1×106 MHCC97-H, MHCC97-L cells and parts of freeze tumor sample (n=12) were lysed in RIPA buffer (50 mM Tris–HCl pH7.5; 150 mM NaCl; 0.5% NaDOC; 1% NP40; 0.1% SDS) plus protease inhibitors. Protein was extracted by micro centrifugation for 30 minutes, Protein concentration was determined using Bradford Reagent. 20ul equal amount of samples and 10ul markers were run onto 10% SDS-PAGE gel and electro-transferred onto PVDF membrane using Mini-Genie blotting system (Bio-Rad). The membranes were incubated with primary antibody, Mouse anti-human TGF β1 antibody (Chemicon, 1:1000 diluted) and Mouse anti-human β-actins antibody (Chemicon, 1:2000 diluted), and HRP-conjugated goat anti-mouse IgG secondary antibody (SIGMA, 1:2000 diluted), The membranes were washed and incubated with 10ml LumiGLO and exposed to film. The blot bands intensity was quantitatively analyzed using FURI Smart View 2000 software (Shanghai). The ratio of TGF β1 to β-actin bands intensity was assessed.
Cytoimmunochemistry and Immunohistochemistry
2×105 MHCC97-H and MHCC97-L cell were plated and cultured in six-well plate respectively, when reached to 60% confluent, the cells were fixed with 100% methanol, permeabilized with 0.5% Triton X-100, and sequentially incubated with the primary anti- TGF β1 monoclonal antibodies and anti-mouse immunoglobulin (Ig) coupled to Horseradish peroxidase (HRP), then, the cells were stained with DAB (3, 3′-diaminobenzidine) and counterstained with hematoxylin. Paraffin-embedded tumor tissues were sliced as 5μm sections in thickness and mounted on glass. Slides were deparaffinated and rehydrated over 10 min through a graded alcohol series to deionized water; 1% Antigen Unmasking Solution (Vector Laboratories) and microwaved were used to enhance antigen retrieval; the slide were incubated with anti-TGF β1 monoclonal antibodies and HRP-conjugated secondary antibody, and then, stained with DAB.
Total protein of all tumor tissues were extracted as described above. TGF β1 protein levels in tumors were determined using the Quantikine TGF β1 Immunoassay (R&D, Minneapolis, MN,USA). The operational approach was performed according to manufacture specification.
Statistical analysis was performed using SPSS 11.5 software (SPSS Inc, USA). The data were analyzed by Students’ t test, one-way analysis of variance and covariance analysis. All statistical tests were two-sided; a P value of less than 0.05 was considered statistically significant.
The tumor weight and pulmonary metastatic rate
The tumors of MHCC9-H model grew fast than that of MHCC97-L, and especially in early stage of tumor formation, MHCC9-H spent shorter time (days) than MHCC97-L getting to the size of 500mm3 (21.93±3.67 vs. 30.83±1.94, P<0.001) (Figure1A), however, the growth speed became similar from the size of 500mm3 to 1500 mm3 (9.00±2.69 vs.10.83±1.47, P=0.14 ) (Figure1B). MHCC9-H model had bigger pulmonary metastatic loci than MHCC97-L model (Figure1C,D). The mean tumor weight (g) in MHCC9-H and MHCC97-L were 1.75±0.75 and 1.26±0.51, and the pulmonary metastatic rate were 55% and 36.36%; and the average number of metastatic cell in lung were 119.25±177.39 and 43.36±47.80 respectively (Table1).
The MHCC97-H cells have lower mRNA expression level of TGF β1 and Smads than MHCC97-L in vitro and in vivo
As shown in Table2, mRNA levels of TGF β1 and Smad2 in MHCC97-H cell line were lower than that of MHCC97-L cell line (0.18±0.15 vs. 0.40±0.19, P=0.011; 0.99±0.17 vs. 2.56±0.66, P=0.047), and TGF β1 in MHCC97-H model was also lower than that of MHCC97-L models (1.24±0.96 vs. 2.81±1.61, P=0.002). Compared with MHCC97-L cells, the expression of TGF β1 protein in MHCC97-H was also lower by western blot analysis (Figure2A), and in mice models, According to quantitative band-intensity analysis of Western blots, the average ratio of TGF β1 to β-actin bands intensity in MHCC97-L models, MHCC97-H models were 0.75±0.45 and 0.57±0.37 (Figure2B).
By cytoimmunochemistry (Figure1Ca, b) and immunohistochemistry method (Figure2Da, b), we found MHCC97-L cell lines and MHCC97-L models have higher expression level of TGF β1 than MHCC97-H cell lines and MHCC97-H models.
The TGF β1 protein levels correlated with metastasis
Compared with MHCC97-H models, MHCC97-L models have a higher TGF β1 protein level by ELASA (0.32±0.22 vs. 1.37±0.95, P<0.001) (Figure3A). And in MHCC97-H and MHCC97-L models, we divided all samples (31cases) into two groups according to metastasis, and we found the TGF β1 protein level in metastasis group was higher than in none metastasis group by covariance analysis (0.16±0.15 vs. 0.12±0.10, P<0.001) (Figure3B). In addition, in mRNA levels, the relations between TGF β1 and Smad2, Smad7 were also found (R2=0.12, P=0.059 and R2=0.40, P<0.001, respectively) (Figure3C,D), but none of them correlated to tumor size.
Although MHCC97-L cell line and MHCC97-H cell line have an identical genetic background, in this study, we observed the expression of TGF β1, Smad2 and Smad7 in MHCC97-L cell lines was higher than that in MHCC97-H cell lines both in vitro and in vivo, in addition, MHCC97-L have a slower growth speed in early stage of tumor formation. Our results were in agreement with other documents, which demonstrate TGF β can induce apoptosis of human hepatoma cell line in vitro, and enhance tumor formation by transfection of an antisense TGF-β1 expression vector into cancer cells[24, 25]. Our results suggest that the basic level of TGF β in cell line could affect on its growth, and TGF β1/Smads play an inhibitory role in the course of tumorigenensis.
We also found the TGF β1 protein were positively correlated with pulmonary metastasis in the models, and in mRNA levels, TGF β1 correlated with that of Smad2 and Smad7. Our results were consistent with other studies regarding the association between TGF β1/Smads and HCC metastasis[7, 15, 26], and these results support the veiw that TGF β1/Smads promote pulmonary metastasis of HCC.
The contradict results in this study, inhibitory role in tumorgenesis and promoting role in tumor metastasis, may arise from the dual role of TGF β1 in different stage of cancer development. It has reported during the early stages of tumor formation, TGF β1 acts as a tumor suppressor, inhibiting proliferation and inducing apoptosis of tumor cells. However, during later stages of tumorigenesis, many tumor cells become unresponsive to the growth inhibitory functions of TGF β1, and get more motile, more invasive, and more resistant to apoptosis. In addition, TGF β can stimulate non-invasive HCC cells to acquire invasive phenotypes. Our results support the view that TGF β1/Smads play a dual role in the development of HCC. We also observed MHCC97-L cell lines have a higher TGF β1/Smads levels but a lower metastasis than MHCC97-H cell lines, and both cell lines have an upregulated levels of TGF β1 during the course of metastasis. These results reflected the basic levels of TGF β1 were not the only factor for metastasis, and highlight that the role of TGF β1/Smads should be decided in an active course.
The result that TGF β correlate with pulmonary metastasis in our study will give a new insight to investigate the metastatic mechanism of HCC. The cells in the tumor tissue communicate through the secretion of growth factors, chemokines, and cytokines during tumor progression, and TGF β is unique in its ability to both promote and inhibit tumorigenesis, depending on the cell type it is acting on. Moreover, TGFβ1 could affect various molecular expression, such as P160ROCK, Integrin and Matrix Metalloproteinases,and all of these molecules relate to HCC invasion.
Collectively, our results suggest that TGF β1 play an important role in the process of tumor growth and pulmonary metastasis of HCC, and the role were time-dependent and based on cell type itself. Strategies to modulate expression levels of TGF β1 could provide a better approach for the treatment of pulmonary metastasis in HCC.
This work was supported in part by China National Natural Science Foundation for distinguished Young Scholars (30325041), the China National '863' R & D High-tech Key Project.
Polymerase chain reaction
Dulbecoo's modified Eagle's medium
Phosphate buffered saline
Alpha fetal protein.
Ono T, Yamanoi A, Nazmy E, Assal O, Kohno H, Nagasue N: Adjuvant chemotherapy after resection of hepatocellular carcinoma causes deterioration of long-term prognosis in cirrhotic patients: meta analysis of three randomized controlled trials. Cancer. 2001, 91: 2378-2385. 10.1002/1097-0142(20010615)91:12<2378::AID-CNCR1271>3.0.CO;2-2.
Kurokawa Y, Matoba R, Takemasa I, Nagano H, Dono K, Nakamori S, Umeshita K, Sakon M, Ueno N, Oba S, et al: Monden MMolecular-based prediction of early recurrence in hepatocellular carcinoma. J Hepatol. 2004, 41: 284-291. 10.1016/j.jhep.2004.04.031.
Lai EC, Fan ST, Lo CM, Chu KM, Liu CL, Wong J: Hepatic resection for hepatocellular carcinoma. An audit of 343 patients. Ann Surg. 1995, 221: 291-298. 10.1097/00000658-199503000-00012.
Ye QH, Qin LX, Forgues M, He P, Kim JW, Peng AC, Simon R, Li Y, Robles AI, Chen Y, et al: Predicting hepatitis B virus–positive metastatic hepatocellular carcinomas using gene expression profiling and supervised machine learning. Nat Med. 2003, 9: 416-423. 10.1038/nm843.
Genda T, Sakamoto M, Ichida T, Asakura H, Hirohashi S: Cell motility mediated by rho and rho-associated protein kinase plays a critical role in intrahepatic metastasis of human hepatocellular carcinoma. Hepatology. 1999, 30: 1027-1036. 10.1002/hep.510300420.
Nakamura T, Kimura T, Umehara Y, Suzuki K, Okamoto K, Okumura T, Morizumi S, Kawabata T, Komiyama A: Long-term survival after report resection of pulmonary metastases from hepatocellular carcinoma: report of two cases. Surg Today. 2005, 35: 890-892. 10.1007/s00595-005-3022-7.
Giannelli G, Fransvea E, Marinosci F, Bergamini C, Colucci S, Schiraldi O, Antonaci S: Transforming growth factor-beta1 triggers hepatocellular carcinoma invasiveness via alpha3beta1 integrin. Am J Pathol. 2002, 161: 183-193. 10.1016/S0002-9440(10)64170-3.
Lyer S, Wang ZG, Akhtari M, Zhao W, Seth P: TargetingTGFbeta signaling for cancer therapy. Cancer Biol Ther. 2005, 4: 261-266. 10.4161/cbt.4.3.1566.
Oft M, Peli J, Rudaz C, Schwarz H, Beug H, Reichmann E: TGFbeta1 and Ha-Ras collaborate in modulating the phenotypic plasticity and invasiveness of epithelial tumor cells. Genes Dev. 1996, 10: 2462-2477. 10.1101/gad.10.19.2462.
Kinnman N, Andersson U, Hultcrantz C: In situ expression of transforming growth factor-beta 1–3, latent transforming growth factor-beta binding protein and tumor necrosis factor-alpha in liver tissue from patients with chronic hepatitis C. Scand J Gastroenterol. 2000, 35: 1294-1300. 10.1080/003655200453656.
Rubtsov YP, Rudensky AY: TGFβ signalling in control of T-cell-mediated self-reactivity. Nature Immunol. 2007, 7: 443-453. 10.1038/nri2095.
Itoh S, Itoh F, Goumans MJ,PTD: Signaling of transforming growth factor-b family members through Smad proteins. Eur J Biochem. 2000, 267: 6954-6967. 10.1046/j.1432-1327.2000.01828.x.
Welm AL: TGFβ primes breast tumor cells for metastasis. Cell. 2008, 133: 27-28. 10.1016/j.cell.2008.03.012.
Song BC, Chung YH, Kim JA, Choi WB, Suh DD, Pyo SI, Shin JW, Lee HC, Lee YS, Suh DJ: Transforming growth factor-beta1 as a useful serologic marker of small hepatocellular carcinoma. Cancer. 2002, 94: 175-180. 10.1002/cncr.10170.
Giannelli G, Bergamini C, Fransvea E, Sgarra C, Antonaci S: Laminin-5 With Transforming Growth Factor-β1 Induces Epithelial to Mesenchymal Transition in Hepatocellular Carcinoma. Gastroenterology. 2005, 129: 1375-1383. 10.1053/j.gastro.2005.09.055.
Grasl-Kraupp B, Rossmanith W, Ruttkay-Nedecky B, Mullauer L, Kammerer B, Bursch W, Schulte-Hermann R: Levels of transforming growth factor β and transforming growth factor β receptors in rat liver during growth, regression by apoptosis and neoplasia. Hepatology. 1998, 28: 717-726.
Jaskiewicz K, Chasen MR: Differential expression of transforming growth factor β, adhesions molecules and integrins in primary, metastatic liver tumors and in liver cirrhosis. Anticancer Res. 1995, 15: 559-562.
Yuen MF, Norris S, Evans LW, Langley PG, Hughes RD: Transforming growth factor-β1, activin and follistatin in patients with hepatocellular carcinoma and patients with alcoholic cirrhosis. Scand J Gastroenterol. 2002, 37: 233-238. 10.1080/003655202753416939.
Kim YJ, Lee HS, Im JP, Min BH, Kim HD, Jeong JB, Yoon JH, Kim CY, Kim MS, Kim JY, et al: Association of transforming growth factor-β1 gene polymorphisms with a hepatocellular carcinoma risk in patients with chronic hepatitis B virus infection. Exp Mol Med. 2003, 35: 196-202.
Li Y, Tang Y, Ye L, Liu B, Liu K, Chen J, Xue Q: Establishment of a hepatocellular carcinoma cell line with unique metastatic characteristics through in vivo selection and screening for metastasis-related genes through cDNA microarray. J Cancer Res Clin Oncol. 2003, 129: 43-51.
Li Y, Tang ZY, Ye SL, Liu YK, Chen J, Xue Q, Chen J, Gao DM, Bao WH: Establishment of cell clones with different metastatic potential from the metastatic hepatocellular carcinoma cell line MHCC97. World J Gastroenterol. 2001, 7: 630-636.
Carey KD, Garton AJ, Romero MS, Kahler J, Thomson S, Ross S, Park F, Haley JD, Gibson N, Sliwkowski MX: Kinetic analysis of epidermal growth factor receptor somatic mutant proteins shows increased sensitivity to the epidermal growth factor receptor tyrosine kinase inhibitor, erlotinib. Cancer Res. 2006, 66: 8163-8171. 10.1158/0008-5472.CAN-06-0453.
Lin JK, Chou CK: In Vitro apoptosis in the human hepatoma cell line induced by Transforming Growth Factor beta1. Cancer Res. 1992, 52: 385-388.
Wu SP, Sun LZ, Willson JK, Humphrey L, Kerbel R, Brattain MG: Repression of autocrine transforming growth factor beta 1 and beta 2 in quiescent CBS colon carcinoma cells leads to progression of tumorigenic properties. Cell Growth Diff. 1993, 4: 115-123.
Wu SP, Theodorescu D, Kerbel RS, Willson JK, Mulder KM, Humphrey LE, Brattain MG: TGF-beta 1 is an autocrine-negative growth regulator of human colon carcinoma FET cells in vivo as revealed by transfection of an antisense expression vector. J Cell Biol. 1992, 116: 187-196. 10.1083/jcb.116.1.187.
Fransvea E, Angelotti U, Antonaci S, Giannelli G: Blocking transforming growth factor-beta up-regulates E-cadherin and reduces migration and invasion of hepatocellular carcinoma cells. Hepatology. 2008, 47: 1557-1566. 10.1002/hep.22201.
Derynck R, Akhurst RJ, Balmain A: TGF-β signaling in tumor suppression and cancer progression. Nat Genet. 2001, 29: 117-129. 10.1038/ng1001-117.
Katabami K, Mizuno H, Sano R, Saito Y, Ogura M, Itoh S, Tsuji T: Transforming growth factor-β1 upregulates transcription of a3 integrin gene in hepatocellular carcinoma cells via Ets-transcription factor-binding motif in the promoter region. Clin Exp Metastas. 2005, 22: 539-548. 10.1007/s10585-005-5260-x.
Littlepage LE, Egeblad M, Werb Z: Coevolution of cancer and stromal cellular responses. Cancer Cell. 2005, 7: 499-500. 10.1016/j.ccr.2005.05.019.
Bhowmick NA, Ghiassi M, Aakre M, Brown K, Singh V, Moses HL: TGF-beta-induced RhoA and p160ROCK activation is involved in the inhibition of Cdc25A with resultant cell-cycle arrest. PNAS. 2003, 100: 15548-15553. 10.1073/pnas.2536483100.
Wahl SM, Allen JB, Weekst BS HLW, Klotmant PE: Transforming growth factor 1–3 enhances integrin expression and type IV collagenase secretion in human monocytes. PNAS. 1993, 90: 15548-15553.
Li GC, Ye QH, Xue YH, Sun HJ, Zhou HJ, Ren N, Jia HL, Shi J, Wu JC, Dai C, et al: Human mesenchymal stem cells inhibit metastasis of a hepatocellular carcinoma model using the MHCC97-H cell line. Cancer Sci. 2010, 101: 2546-2553. 10.1111/j.1349-7006.2010.01738.x.
We would like to thank Mrs. Qiong Xue, Dong-Mei Gao, Rui-Xia Sun and Jie Chen, Drs. Hai-Ying Zeng, Teng-Fang Zhu and Jun Chen for their help in the animal experiments and cell culture.
The authors declare that they have no competing interests.
GCL wrote the first draft of the manuscript, performed cell culture and contributed to the experimental design; QHY analyzed data and contributed to writing and editing of the manuscript; DQZ and NR analyzed data; LXQ designed experiments and wrote and revised the manuscript; HLJ performed western blot and wrote and edited the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Li, G., Ye, Q., Dong, Q. et al. TGF beta1 and related-Smads contribute to pulmonary metastasis of hepatocellular carcinoma in mice model. J Exp Clin Cancer Res 31, 93 (2012). https://doi.org/10.1186/1756-9966-31-93
- Hepatocellular carcinoma
- Transforming growth factor beta