- Open Access
TM4SF1 upregulates MYH9 to activate the NOTCH pathway to promote cancer stemness and lenvatinib resistance in HCC
Biology Direct volume 18, Article number: 18 (2023)
TM4SF1, a member of the transmembrane 4 superfamily, is crucial for both healthy and malignant human tissues. The significant function of TM4SF1 in the incidence and progression of cancer has been widely recognized in recent years. Although some achievements have been made in the study of TM4SF1, the effect of TM4SF1 on cancer stemness in hepatocellular carcinoma (HCC) and its molecular basis are yet to be reported. We found through abundant in vitro and in vivo experiments which the expression of TM4SF1 was positively correlated with the progression and cancer stemness of HCC. We identified the downstream protein MYH9 of TM4SF1 and its final regulatory target NOTCH pathway using bioinformatics analysis and protein mass spectrometry. We cultivated a Lenvatinib-resistant strain from HCC cells to examine the relationship between cancer stemness and tumor drug resistance. The study confirmed that TM4SF1 could regulate the NOTCH pathway by upregulating MYH9, thus promoting cancer stemness and Lenvatinib resistance in HCC. This study not only provided a new idea for the pathogenesis of HCC but also confirmed that TM4SF1 might become a new intervention point to improve the clinical efficacy of Lenvatinib in treating HCC.
With 905,677 individuals, liver cancer rose to the sixth-most prevalent cancer worldwide in 2020. The number of people dying from liver cancer also reached 830,180, ranking third among all cancers . The majority (80%) of initial liver cancers are hepatocellular carcinomas (HCC) . The mechanism of HCC's occurrence and progression still needs to be investigated, despite recent advancement in research and treatment. Previous studies showed that the occurrence and development of HCC in many cases were caused by the imbalance of signal pathways instigated by gene mutations , including upregulations of Wnt/β-catenin, Akt/mTOR, P53 signal pathway  and so on. An in-depth exploration of these processes can help us to treat HCC more specifically. Moreover, some scholars have found that cancer stem cells (CSCs) are a major factor affecting tumor recurrence and metastasis [5, 6], which has also gained our attention. CSCs are a kind of cells with multi-directional differentiation potential and unlimited proliferation. Their division and differentiation are not regulated by the body . They are also considered one of the principal sources of cancer , as tumors continue to grow due to the continuous self-renewal and differentiation of CSCs. Therefore, the study of cancer stemness has clinical value and profound clinical significance . The conclusion that CSCs can effectively promote drug resistance in tumors has been widely accepted [10,11,12]. Among the first-line drugs used for treating HCC [13, 14], Lenvatinib is the first targeted drug approved for advanced-stage HCC and plays a considerable role in clinical treatment. However, it still cannot prevent patients from developing drug resistance, resulting in a poor prognosis . At present, different opinions are available on the molecular mechanism of these phenomena, leading to increasing attention from scientists for further research. TM4SF1 is located in 3q25.1 and has four highly conserved transmembrane domains, two extracellular rings, and one intracellular small ring [16, 17]. Tumor-associated antigen was how TM4SF1 was initially classified. It stabilizes the cell signal complex and plays a role in cell proliferation, adhesion, and movement [18, 19]. Studies have confirmed that TM4SF1 is upregulated in a variety of cancers [20,21,22,23] including HCC. A study by Tang and colleagues in 2020 showed that TM4SF1 might also have a potential effect on cancer stemness . Using the gene set enrichment analysis (GSEA) technology, we enriched the signaling pathways significantly related to TM4SF1 and found that the NOTCH pathway was related to cancer stemness. After protein mass spectrometry analysis of TM4SF1, we screened out the downstream protein MYH9, which, as a widely expressed cytoskeletal protein, was closely related to a variety of cancers. Hence, we can assume that TM4SF1 may activate the NOTCH pathway via upregulation of MYH9, ultimately promoting tumor stemness in HCC. Moreover, TM4SF1 may also affect the resistance of HCC to Lenvatinib and become a potential target for the clinical treatment of HCC.
Materials and methods
Five HCC cell lines and human normal liver cell LO2 were employed in this study. All HCC cell lines (Hep-G2, Hep-3B, LM3, Huh7, and MHCC97H) were acquired from the Chinese Academy of Sciences' cell bank, while LO2 was acquired from Southern Medical University's Cancer Research Institute. All cells were grown in Dulbecco's modified Eagle's medium (DMEM) with 10% fetal bovine serum (FBS) from GIBCO, NY, USA, and antibiotics from BIOMYC-3 Antibiotic Solution and Penicillin–Streptomycin Amphotericin B Solution from Biological Industries, Israel, at 37 °C and 5% CO2.
Clinical specimens and immunohistochemistry (IHC)
We purchased a HCC tissue microarrays from Shanghai Outdo Biotech Company, containing 90 samples of HCC and 90 samples of paracancerous tissues. The clinical samples were approved by patients' informed consent and authoritative histopathologists, and passed the review of the ethics committee on the use of the samples (Ctl No.:YB M-05-02). After a series of steps such as dewaxing, hydration, block the activity of endogenous peroxidase and nonspecific antigen, antigen retrieval, incubation with first antibody, staining and sealing. The tissue microarray was observed under Inverted microscope and double-blind scored in accordance with staining intensity (0-Colorless, 1-light yellow, 2-light brown, 3-dark brown) and positive area ratio (0: < 5%, 1:5–25%, 2:25–50%, 3:50–75%, 4: > 75%) by The Department of Pathology of the Hospital of Integrated traditional Chinese and Western Medicine of Southern Medical University (Score calculation: staining intensity grade*positive area proportion grade).
SiRNA, plasmids and transfection conditions
The cells to be transfected were cultured in cell petri dishes for 48–72 h in advance. SiRNA (Ribo-Bio, Guangzhou, China) and overexpressing plasmid (GeneChem, Shanghai, China) was introduced into cells by Lipofectamine 3000 (Thermo Fisher Scientific Company, USA) to obtain cell lines that temporarily silenced TM4SF1 or MYH9 (Additional file 1: Table S1).
Lentivirus production and infection
Lentivirus human TM4SF1 gene was introduced into lentiviral vector GV367 to construct lentiviral TM4SF1 (LV-TM4SF1) (GeneChem, Shanghai, China) to infect hepatoma cells, and stable hepatoma cell lines with loss of TM4SF1 expression was constructed. On this basis, the empty vector lentivirus (LVNC) (GeneChem, Shanghai, China) was used as the control, and the cells were screened with puromycin, and the polyclonal cells with green fluorescence signal were selected for follow-up experiment (Additional file 1: Table S1).
Extraction of RNA and QPCR
The total RNA of cells was obtained by using RNA extraction kit (Foregene, Chengdu, China), and cDNA was amplified by reverse transcription kit produced by Takara company. Finally, quantitative polymerase chain reaction of cDNA was performed with the help of Bio-Rad T100 system using primers ordered from Guangzhou IGE Biology Company (Additional file 1: Table S2), and the changes of genes were detected by calculating the data.
Western blot analysis and antibody
The protein of the cell was obtained from the solution system of lytic buffer, protease inhibitor in the ratio of 100:1:1 and phosphatase inhibitor. The pure protein was collected by ultrasound and centrifugation and quantified by BCA protein analysis kit. SDS–polyacrylamide gel electrophoresis was used to separate the protein. It was then moved to a polyvinylidene fluoride membrane and treated at 4 °C with the desired primary antibody. The antibodies used in the experiment include TM4SF1, CD44, CD133, OCT4, SOX2, MYH9, NOTCH1, JAGGED1, HES1, β-Tubulin (Additional file 1: Table S3). After the protein was fully combined with the primary antibody, the chemiluminescence kit (Thermo Fisher Scientific Company, USA) was used to image in the ChemiDoc XRS + molecular imager (Bio-Rad, Hercules, CA, USA)to analyze the expression of the protein.
Tumour sphere formation
The cells (5 × 103) were inoculated into six-well ultra-low attachment plates per well (Corning, NY, USA) and cultured with special sphere formation medium (It is composed of 1 mL B27, 20 uL human EGF of 20 ng/mL and 20 uL human FGF of 20 ng/mL added to 50 mL DMEM/F12 culture medium) at 37 °C, 5% CO2 for 1–2 weeks. Tumour spheres were observed and photographed with inverted microscope. The proportions of tumour spheres with diameters of 50–100 μm, 100–150 μm and > 150 μm in the total visual field were calculated, and statistical charts were drawn.
Detection of side population cells
After the experimental cells were taken out for digestion and centrifugation, the supernatant was abandoned and the cells were re-suspended with 1 mL high-glucose medium containing 2% serum. After cell counting, the cell concentration was adjusted to 1 × 106 cells/mL in high-glucose medium containing 2% serum heated by 37 °C water bath, and then transferred into the flow tube. A blank control group was set for each experimental group, plus 50 μM verapamil. After covering the lid, it was first incubated in a 37 °C water bath for 30 min. The flow tube was taken out every 10 min, mixed up and down gently once, centrifuged at 12000 r/min and 4 °C for 5 min, then the supernatant was discarded, and the medium was added for re-suspension. One of the samples was added with DCV dye solution and the concentration was adjusted to 20 uM. After covering the sample, it was placed in a 37 °C water bath and incubated for 90 min to make it fully stained. The flow tube was taken out every 10 min, gently mixed up and down once, centrifuged at 12000 r/min and 4 °C for 5 min, and the supernatant was discarded. The PBS cells were precooled at 4 °C and detected by upflow cytometry after treatment.
Fluorescence-activated cell sorting
Consistent with the above procedure for “3.8 Detection of Side population cells”, after the cells were suspended by PBS, 5μLPI was added into the cell suspension under the condition of avoiding light. Note that the experimental process was carried out on the ice box. After the operation, upflow cytometry was used to detect the expression of side population cells, and then FACS staining was performed. Co-incubated with FITCCD133, PE-CD44 antibody, and AF647 sheep anti-mouse immunoglobulin antibody (AF647), staining cells were taken for LSRFortessa or AriaII (BD) analysis or classification.
Protein mass spectrometry
The protein extracted from the tested cells was entrusted to Fitgene Biotechnology company (Guangzhou, China)for protein spectrum analysis, and the genes whose expression degree changed after silencing TM4SF1 were screened. These genes were further screened according to the multiple of gene expression difference and bioinformatics database search, the appropriate genes were determined for further research.
The experimental cells were removed and lysed with mild lysate. The cell suspension was transferred to a 1.5mLEP tube, centrifuged at 4 °C and 12000r for 30 min, supernatant was taken, and a small amount of protein samples were taken to a new 1.5mLEP tube and stored at − 20 °C for protein western blotting analysis. The remaining protein samples were added to the protein primary antibody of the decoy protein and incubated overnight with slow shaking at 4 °C. Pretreated 10 μL proteinA agarose beads were added to the protein samples incubated with the antibody overnight. The samples were slowly shaken at 4 °C for 2–4 h to conjugate the antibody to the proteinA agarose beads. After the immunoprecipitation reaction, the agar–agar beads were centrifuged at 4 °C and 3000r for 3 min to the bottom of the tube, the supernatant was absorbed, and the agar–agar beads were washed with 1 mL cracking buffer for 3–4 times. Adding 15μL loadingbuffer sample buffer, boiling water for 5 min, finally by western blot analysis.
Establishment of Lenvatinib-resistant cell lines
The cells were cultured for a long time by the method of increasing concentration, and the model of Lenvatinib-resistant strain was established. At 37 °C and 5% CO2, the cells were grown in DMEM containing 10% fetal bovine serum (FBS). Lenvatinib of 1 mg was administered to the culture media once the cells had stabilized. After 1 week of culture, if the cells could grow stably, they would be changed to a higher concentration of Lenvatinib (increasing 1 mg each time), until the cells could grow stably in the culture medium containing 10 mg Lenvatinib, indicating that Lenvatinib-resistant cells were induced successfully.
Etermination of CCK-8 cell viability
The cell survival rate and tolerance to Lenvatinib were detected by CCK-8 method. According to previous studies , a 96-well plate was prepared, the experimental cells were taken out for digestion, centrifugation and suspension, 5000 cells were counted by cell counting method, diluted with 100 μL medium and inoculated into each pore, then placed back into the incubator for 24 h. After observing the cell growth status, 10 μL of different concentrations of Lenvatinib were added into each pore and placed back into the incubator for 48 h. The 96-well plates were taken out and added with 10 μL CCK-8 reagent in each well under the condition of hiding from light. The 96-well plates were wrapped with aluminum foil and put back into the incubator for further culture for 1–4 h. The 96-hole plate was removed, the lid was opened, and OD values were measured at 450 nm wavelength with enzyme marker. Finally, standard curves were drawn with the help of Graphpadprism software.
All animal experiments meet the requirements of the International guidelines for Animal Care and Conservation and the Animal Research Committee of the academic Medical Center of Southern Medical University. The tumor growth was evaluated by subcutaneous transplantation of xenogeneic oysters in mice.
LM3 with stable expression loss of TM4SF1(LV-TM4SF1) in different concentration gradients was subcutaneously injected into 3 week-old male nude mice weighing 16–25 g, and the corresponding concentration of LM3 infected with the empty vector lentivirus (LV-NC) was used as negative control. The growth of the tumor was observed and the size of the tumor was recorded with Vernier caliper every three days. 45 days later, the tumor was taken out, photographed and weighed, and the data of the tumor were collected for statistics and analysis.
The effect of TM4SF1 on drug resistance of human hepatocellular carcinoma cell line Lenvatinib was studied by subcutaneous transplantation in nude mice. Saline containing Lenvatinib was injected into nude mice by intraperitoneal injection in advance. The experimental group of nude mice received subcutaneous injections of various quantities of Lenvatinib-resistant LM3; untreated LM3 served as the negative control. Then, according to the previous method, nude mice were cultured and tumor data were collected for follow-up experiments.
Primary cell culture
Fresh tumor tissue was taken, cut into small pieces with a diameter of about 4 mm, and repeatedly rinsed with PBS until the blood clots and adipocytes on the tissue fragments were basically removed, then transferred to a centrifuge tube, digested in a 37 °C water bath with 0.25% trypsin for 15–20 min, rinsed with PBS for 2–3 times, and then gently blown off the digested cells with a small amount of cell medium. Cells are counted and cultured at appropriate concentrations. Note that a sterile environment should be carefully maintained during the experiment to avoid cell contamination. If there is still undigested tissue, the above steps can be repeated for multiple digestion.
The study used bioinformatics databases to obtain materials that can support the conclusions of the study. The databases used include: The Cancer Genome Atlas (TCGA), TIMER2.0, Gene Expression Profiling Interactive Analysis (GEPIA). Gene Set Enrichment Analysis (GSEA) technique was used to enrich the KEGG and GO pathways of TM4SF1, and the pathways related to its expression were obtained.
Each group of in vitro experiments were repeated three times to ensure the credibility and reproducibility of the study, and the average-standard deviation (x ± s) was used to express the measured results. SPSS13.0 software was used to analyze the statistical process and data in the study, and t-test was used to compare between groups, the difference was statistically significant (P < 0.05).
High expression of TM4SF1 was closely related to the progression and poor prognosis of HCC
We found via the bioinformatics analysis using the TIMER2.0 and GEPIA databases in the early stage of the experiment that TM4SF1 was much more expressed in human HCC tissues than in healthy liver tissues (Additional file 2: Fig. S1A and B). TM4SF1 and the cumulative survival rate, overall survival rate, and disease-free survival rate in patients with HCC were further examined. Patients with strong TM4SF1 expression had a generally unfavorable prognosis, we found. (Additional file 2: Fig. S1C and D). We also measured the expression of TM4SF1 in all liver cancer cell lines using the TCGA database, and the results revealed that the majority of these cell lines had high levels of TM4SF1 expression (Additional file 2: Fig. S1E). Then, we purchased the tissue microarray with 90 hepatocellular carcinoma and paracancerous tissues for immunohistochemical experiment (Fig. 1A) and scored (Fig. 1B). The findings demonstrated that HCC tissues had much greater levels of TM4SF1 expression than paracancerous tissues. Additionally, utilizing the prognostic data from 90 patients with hepatocellular carcinoma, we examined the overall survival rate. The findings demonstrated that patients with low expression of TM4SF1 had a higher overall survival rate than those with high expression of TM4SF1 (Fig. 1C). And than, using SPSS 13.0 statistical software to analyze our data, we found that the expression of TM4SF1 was positively correlated with the tumor size, Edmondson–Steiner grade, and cirrhosis, but not with sex, age of the patients and so on (Table 1). The expression of TM4SF1, tumor size, and Edmondson-Steiner grade were found to be linked with the overall survival percentage of 90 HCC patients in both univariate and multivariate analyses (Table 2, Fig. 1D). We cultured five hepatoma cell lines (Hep3B, HepG2, LM3, Huh7, and MHCC97H) and normal liver cells LO2, extracted RNA and protein from the cells, and detected the expression of TM4SF1 at the RNA and protein levels using qPCR and Western blot (WB), respectively (Fig. 1E and F).
TM4SF1 upregulated the stemness of tumor cells in HCC
In order to prove the correlation between TM4SF1 and cancer stemness, we used OCLR algorithm to evaluate the stemness index of TM4SF1 based on the tissues of 371 patients with HCC in TCGA database. In HCC tissues compared to normal tissues, it was discovered that TM4SF1 had a greater stemness index (Fig. 2A). Using the GEPIA database, we examined the relationship between TM4SF1 and molecules associated with cancer stem cells. We discovered that TM4SF1 expression was positively correlated with CD44, CD133, and OCT4 (Fig. 2B). We verified the high expression of TM4SF1 in tumor spheres via the PCR experiment using the RNA extracted from the tumor spheres of liver cancer provided by the Cancer Research Institute of Southern Medical University, Guangzhou, China (Fig. 2C). Then, we selected two HCC cell lines Hep3B and LM3 as the research objects and inserted siRNA into the cells to silence TM4SF1 with the help of a transfection reagent Lipofectamine 3000, and verify the efficiency of silencing. (Fig. 2D and E). The relationship between TM4SF1 and cancer stem–related molecules CD44, CD133, OCT4, and SOX2 was detected via qPCR and WB, which confirmed the decrease in the expression of CD44, CD133, OCT4, and SOX2 with the silencing of TM4SF1 (Fig. 2F and G). To investigate the impact of TM4SF1 on the stem cell-related phenotype of HCC, we infected Hep3B and LM3 cells with lentiviruses carrying silenced portions of TM4SF1. After verifying the efficiency (Fig. 2D and E), The side-population cell detection test revealed that the proportion of CSCs was also lower in HCC cells silenced by TM4SF1 than in the control group (Fig. 2H). We confirmed via tumor sphere formation and side-population cell detection tests that TM4SF1 could affect the cancer stemness in HCC. The results showed that the size and number of tumor spheres formed by TM4SF1-silenced HCC cells were inferior to those in the control group (Fig. 2I and J). In order to further confirm whether TM4SF1 is related to the cancer stem of HCC, we isolated CD133+, CD44+ and CD133+/CD44+ cells from LM3 cells with silenced TM4SF1 and LM3 cells from untreated group by flow cytometry. It showed that the proportion of CD133+, CD44+ and CD133+/CD44+ cells in the silent TM4SF1 group was inferior to which in the control group (Fig. 2K).
When carrying out the in vitro experiment, we also considered the subcutaneous tumorigenesis experiment of nude mice as the in vivo experiment to verify the influence of TM4SF1 on the occurrence and progression of HCC. The nude mice in the experimental group were administered by injection with different numbers of LM3 cells infected with LV-TM4SF1; the LM3 cells infected with no-load Lentivirus were used as the control. The changes in size during tumor growth were recorded, and the tumor was removed 45 days later (Fig. 3A). However, we found that the tumor size and weight in the experimental group were smaller than those in the control group regardless of the number of days or the number of cells (Fig. 3B and C). We also made tissue sections of subcutaneous tumors of nude mice in previous experiments and performed immunohistochemical staining. We chose TM4SF1, CD44, CD133, and OCT4 as the primary antibodies to combine with tumor tissues. It demonstrated that, compared to the control group, the staining score in the experimental group was much lower (Fig. 3D). We grew primary cells from the collected tumor tissues of nude mice. After the cells grew stably, RNA and protein were extracted and qPCR and WB were performed. The results showed that the expression of TM4SF1 and stem-related molecules in primary tumor cells of nude mice increased (Fig. 3E and F).
TM4SF1 activated the NOTCH pathway and promoted cancer stemness in HCC cells by upregulating MYH9
When considering the molecular mechanism by which TM4SF1 affected the cancer stemness in HCC cells, we first considered enriching the pathways of KEGG and GO for TM4SF1 using the GSEA technology. Considering 374 patients with liver cancer in the TCGA database as the samples, we enriched many TM4SF1-related pathways. Among the pathways related to cancer stemness, the NOTCH pathway attracted our attention because of its high correlation score (Fig. 4A). We detected the correlation between TM4SF1 and the common ligand JAGGED1, common receptor NOTCH1, and target gene HES1 of the NOTCH pathway via qPCR and WB, and found that the expression of JAGGED1, NOTCH1, and HES1 decreased after silencing TM4SF1 (Fig. 4B and C). RNA and protein extracted from primary tumor cells in nude mice were used for quantitative polymerase chain reaction (QPCR) and Western blotting (WB). The results showed that the expression of Notch pathway related molecules decreased in primary tumor cells of nude mice silenced with TM4SF1 (Fig. 4D and E). In order to further verify the regulatory effect of TM4SF1 on the NOTCH pathway, we purchased a targeted inhibitor RO4929097 of the NOTCH pathway. With the help of WB experiment, we found that after overexpression of TM4SF1, the inhibition effect of RO4929097 on the NOTCH pathway was reversed, and the NOTCH pathway was activated again (Fig. 4F). After confirming upregulation between TM4SF1 and the NOTCH pathway, we wondered how TM4SF1, as a single gene, activated the NOTCH pathway. Therefore, we detected using protein mass spectrometry the genes whose expression changed with the silencing of TM4SF1, and screened out some genes with significant correlation with the cancer stem–related molecules and NOTCH pathway molecules via the correlation analysis of the GEPIA database. When sequencing these differentially expressed genes from multiple species, we were surprised to find that the cytoskeleton gene MYH9 met our screening requirements and had a higher ranking (Fig. 5A and B). Hence, we hypothesized that TM4SF1 might affect the cancer stemness of HCC cells by regulating MYH9 and activating the NOTCH pathway. Through the analysis of TCGA database, we find that there is a high correlation between TM4SF1 and MYH9 (Fig. 5C). Meanwhile, through CO-IP and WB experiments, we found that TM4SF1, as a decoy protein, could fully and specifically bind to MYH9 in HCC cells. (Fig. 5D). In order to explore the specific regulatory relationship between TM4SF1 and MYH9 in liver cancer cells, we silenced TM4SF1 and MYH9 respectively, and observed the expression of another molecule by WB experiment. As a result, we found that when we silenced TM4SF1, the expression of MYH9 was significantly down-regulated, but when we silenced MYH9, the expression of TM4SF1 did not change significantly. Therefore, we believe that MYH9 can act as the downstream molecule of TM4SF1 in hepatocellular carcinoma cells and be subject to unilateral positive regulation from TM4SF1 (Fig. 5E). Then, we silenced MYH9 in LM3 cells, verify the silencing efficiency and detect the relationship between TM4SF1 and NOTCH pathway by qPCR and WB. The consequences indicated that the expression of JAGGED1, NOTCH1, and HES1 was downregulated after silencing MYH9 (Fig. 5F–I). Finally, we overexpressed MYH9 while knocking down TM4SF1. We confirmed via WB, detection of side-population cells, and tumor sphere formation that the overexpression of MYH9 could save the inhibition of cancer stemness in HCC cells caused by the silencing of TM4SF1 (Fig. 6A–D).
TM4SF1 enhanced the lenvatinib resistance in HCC by promoting the cancer stemness in HCC cells
Currently, Lenvatinib is the treatment of choice for treating advanced-stage liver cancer; nevertheless, many patients have a poor therapeutic effect or a poor prognosis because of the resistance to Lenvatinib. We speculated that TM4SF1 affected not only the cancer stemness of HCC but also the Lenvatinib resistance in HCC because of the close relationship between cancer stemness and tumor drug resistance. We cultured hepatoma cell lines Hep3B, HepG2, LM3, Huh7, and MHCC97H in a cell medium with different concentrations of Lenvatinib, and used CCK-8 assay to detect the cell survival rate. The results showed that untreated HepG2 and MHCC97H cells had strong resistance to Lenvatinib, while untreated Hep3B, LM3, and Huh7 cells had poor resistance to Lenvatinib (Fig. 7A). We selected LM3, which was poorly resistant to Lenvatinib, and cultivated the drug-resistant strain LM3-LR having strong resistance to Lenvatinib with the increase in its concentration. We verified via the CCK-8 experiment that the resistance of LM3-LR cells to Lenvatinib was stronger than that of the untreated LM3 cells (Fig. 7C). The qPCR and WB detection showed that the expression of TM4SF1, cancer stem–related molecules, the NOTCH pathway–related molecules, and MYH9 in LM3-LR cells was higher than that in the control group (Fig. 7D and E). The tumor sphere formation experiment is also verified that the cancer stem function of LM3-LR cells was stronger than that of the control group (Fig. 7F and G). We extracted RNA from the tumor spheres, and through qPCR test, we found that the expression of TM4SF1 in the tumor spheres formed by LM3-LR was abnormally increased compared with the untreated group (Fig. 7H). In particular, the expression of TM4SF1 in LM3 cells increased with the increase in the Lenvatinib concentration (Fig. 7B). In order to further verify the accuracy of the previous experiment, we silenced TM4SF1 in LM3-LR cells. The CCK-8 experiment's findings demonstrated that the resistance of LM3-LR cells to lenvatinib decreased (Fig. 7M). Similarly, we also silenced TM4SF1 in HepG2 cells with strong Lenvatinib resistance (Fig. 7N and O), and found that the resistance of HepG2 with silenced TM4SF1 to Lenvatinib also decreased compared with that of the untreated group (Fig. 7P).
We also performed the in vivo drug resistance test on subcutaneous tumor formation in nude mice based on the aforementioned findings (Fig. 7I), and found that the tumor size and weight of nude mice decreased gradually with the increase in the lenvatinib concentration injected through the abdominal cavity. However, when injected with the same Lenvatinib concentration, the volume and weight of subcutaneous tumors of nude mice injected with LM3-LR cells were larger than those of the control group injected with untreated LM3 cells (Fig. 6J and K). Finally, in order to prove the safety of Lenvatinib in nude mice and the stability of treatment, we weighed the nude mice regularly. The findings demonstrated that during Lenvatinib treatment, there was no discernible variation in body weight between the two groups (Fig. 7L).
Although various treatments are available for liver cancer, the main treatment to date is hepatectomy or liver transplantation via clinical surgery [26, 27]. Only when liver cancer is in its early stages may surgery be utilized to treat it effectively. Due to the absence of distinctive signs in the early stages of liver cancer, many individuals have advanced stages of the disease by the time they are diagnosed . Hence, treating patients with advanced-stage liver cancer are particularly difficult and quite challenging. Therefore, the occurrence and development of liver cancer and its molecular mechanism need to be urgently investigated to find the clinical treatment target of liver cancer. Numerous studies have shown that TM4SF1 is strongly associated with the occurrence and progression of cancer [20,21,22]. Zeng and colleagues showed in 2021 that TM4SF1 played a role in advancing the proliferation, invasion, and metastasis of HCC . CSCs are crucial in cancer, according to a vast number of studies. Inducing the creation and growth of primary tumors is a crucial biological activity, but it also plays other significant roles [5, 6]. but also a driving factor to promote cancer drug resistance [10,11,12]. However, few reports exist on the relationship between TM4SF1 and cancer stemness. A recent study showed that TM4SF1 promoted proliferation, invasion, metastasis, and stemness in colon cancer . These studies suggested that TM4SF1 played a promoting role in various cancers and might have a potential link with cancer stemness. However, the effect of TM4SF1 on the stemness of HCC and its related mechanism have not been reported. Since TM4SF1 is strongly expressed in liver cancer cells and tissues compared to normal liver cells and tissues, we investigated its impact on the stemness of HCC in this work. We also confirmed that the expression of TM4SF1 negatively correlates with patient prognosis. The experiment of subcutaneous tumorigenesis in nude mice and the detection of RNA and protein levels in primary cells prepared from tumors in nude mice also confirmed this finding. We also discovered that the expression of stemness-related molecules CD44, CD133, OCT4, and SOX2 also decreased with the decrease in the expression of TM4SF1. The tumor sphere formation test and side-population cell detection also showed that the expression level of TM4SF1 was positively correlated with the cancer stemness of HCC. These results suggested that TM4SF1 might be a potential factor affecting the stemness of HCC.
The NOTCH signaling pathway was first found in Drosophila melanogaster and plays an major role in embryonic development . In this pathway, the extracellular ligand JAGGED1 has the main role: when activated, it binds to the receptor NOTCH1–4 on the cell membrane to form the intracellular signal transduction molecule NICD and finally acts on the nuclear target gene HES1 . The NOTCH pathway has been associated to the occurrence and progression of cancer in recent years [31, 32]. In 2013, Jia Lu et al. found that the activation of the NOTCH pathway promoted cancer stemness in patients with rectal cancer . However, in a 2021 study, Guenter and his collaborators found that the NOTCH pathway positively mediates cancer stemness in patients with thyroid cancer . The NOTCH pathway is crucial for the cancer's stemness. In this research, we took TM4SF1 as the research object and enriched a series of related pathways using the GSEA technology. The NOTCH pathway aroused our interest because of its close and high correlation with cancer stemness. We verified via qPCR and WB experiments that the expression of TM4SF1 was positively correlated with the expression of the NOTCH pathway ligand JAGGED1, receptor NOTCH1, and target gene HES1 in HCC cells, meanwhile, with the help of the NOTCH pathway targeted inhibitor RO4929097, we found that overexpression of TM4SF1 could reverse the inhibitory effect of RO4929097 on the NOTCH pathway. We hypothesized that TM4SF1 might promote cancer stemness in HCC via activating the NOTCH signaling pathway.
The heavy chain of non-muscle myosin IIA, or MYH9, is a protein of 1960 amino acids, which is an important constituent of the non-muscle cytoskeleton . In recent years, many researchers have provided strong evidence that MYH9 plays a key role in the initiation and progression of a number of malignancies [36,37,38,39]. Lin X and his collaborators showed in a 2020 study that silencing MYH9 could inhibit cancer stemness in HCC . These studies demonstrated that the cancer stemness in patients with liver cancer may be significantly impacted by MYH9. We screened out some genes, using protein mass spectrometry, whose expression changed with the silencing of TM4SF1 in HCC cells. Also, we tested the correlation between these genes and cancer stem–related molecules and the NOTCH pathway–related molecules using the GEPIA database. After screening and arranging the data, we found that MYH9 had good performance in differentially expressed genes from multiple species and their correlation rankings. We verified via Co-IP and WB analysis that TM4SF1 could pull down MYH9, Meanwhile, WB experiment also verified that as a downstream molecule, M MYH9 is positively regulated by TM4SF1 and MYH9 expression in HCC cells was favorably connected with the NOTCH pathway's associated components. Therefore, we hypothesized that TM4SF1 activated the NOTCH pathway by positively regulating MYH9, which in turn promoted cancer stemness in HCC.
Lenvatinib is an oral small-molecule inhibitor of multi-receptor tyrosine kinase that has received approval for use in many nations as the first-line therapy for advanced-stage or irreversible HCC . Lenvatinib, as a traditional first-line treatment drug, has made an outstanding contribution to the treatment of patients with advanced-stage liver cancer . However, the tolerance and poor prognosis of some patients are also the disadvantages of Lenvatinib [25, 42, 43]. In previous studies, some scholars made many speculations about the principle of lenvatinib resistance, but still many mysteries remain unsolved. After verifying the promoting effect of TM4SF1 on the cancer stemness of HCC, we speculated based on the close relationship between cancer stemness and tumor drug resistance, that TM4SF1 could enhance the Lenvatinib resistance in HCC cells by promoting the cancer stemness in HCC. The strain LM3 with poor resistance to Lenvatinib was screened using the CCK-8 test, and the resistant strain LM3-LR with strong tolerance to Lenvatinib was cultivated. A series of tests were carried out to verify the effect of drug resistance. Subsequently, we found that the expression of TM4SF1, MYH9, and the NOTCH pathway–related molecules increased in LM3-LR cells. We silenced TM4SF1 in HepG2 and LM3-LR cells, which were more tolerant to Lenvatinib, and found that their resistance to Lenvatinib decreased again. We carried out a subcutaneous tumor formation experiment in nude mice to confirm the accuracy of this result. The results also supported our previous hypothesis.
We verified that TM4SF1 could activate the NOTCH signaling pathway by positively regulating MYH9, which ultimately promoted the cancer stemness in HCC and enhanced the drug resistance of HCC cells to Lenvatinib (Additional file 3: Fig. S2). As a result, TM4SF1 may represent a promising target for the clinical treatment of HCC and contribute significantly to the enhancement of Lenvatinib's therapeutic effects in patients with advanced-stage HCC. First, in order to offer a theoretical foundation for a more precise clinical targeted therapy for HCC, we should conduct a more in-depth investigation into the molecular mechanism by which TM4SF1 activates the NOTCH pathway. In order to improve the curative effect of Lenvatinib and develop new medications to improve the prognosis in patients with HCC, additional research should be done to explore the precise mechanism of TM4SF1 affecting lenvatinib resistance in HCC. This will increase the effectiveness of medications in clinical trials while maintaining safety.
Availability of data and materials
1. The “expression of TM4SF1 in HCC” and “the cumulative survival rate” data that support the findings of this study are available in “TIMER2.0”, “http://timer.cistrome.org/”. 2. The “expression of TM4SF1 in HCC” “intermolecular correlation analysis” and “the overall survival rate and disease-free survival rate” data that support the findings of this study are available in “GEPIA”, “http://gepia2.cancer-pku.cn/#index”. 3. The “expression of TM4SF1 in HCC cell lines” and “stemness index of TM4SF1” data that support the findings of this study are available from “TCGA”, “https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga”. 4. The “GSEA Pathway enrichment” data that support the findings of this study are available from “KEGG”and “GO”, “https://www.kegg.jp/” and “http://geneontology.org/”.
International Agency for Research on Cancer Globocan. 2020. Available from: https://gco.iarc.fr/today/data/factsheets/cancers/11-Liver-fact-sheet.pdf.
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. https://doi.org/10.3322/caac.21492.
Zhu C, Luo X, Wu J, Liu Y, Liu L, Ma S, Xie R, Wang S, Ji W. TM4SF1, a binding protein of DVL2 in hepatocellular carcinoma, positively regulates beta-catenin/TCF signalling. J Cell Mol Med. 2021;25(5):2356–64. https://doi.org/10.1111/jcmm.14787.
Rebouissou S, Nault JC. Advances in molecular classification and precision oncology in hepatocellular carcinoma. J Hepatol. 2020;72(2):215–29. https://doi.org/10.1016/j.jhep.2019.08.017.
Najafi M, Farhood B, Mortezaee K. Cancer stem cells (CSCs) in cancer progression and therapy. J Cell Physiol. 2019;234(6):8381–95. https://doi.org/10.1002/jcp.27740.
Lee TK, Guan XY, Ma S. Cancer stem cells in hepatocellular carcinoma—from origin to clinical implications. Nat Rev Gastroenterol Hepatol. 2022;19(1):26–44. https://doi.org/10.1038/s41575-021-00508-3.
Kuşoğlu A, Biray AÇ. Cancer stem cells: a brief review of the current status. Gene. 2019;10(681):80–5. https://doi.org/10.1016/j.gene.2018.09.052.
Pan Y, Ma S, Cao K, Zhou S, Zhao A, Li M, Qian F, Zhu C. Therapeutic approaches targeting cancer stem cells. J Cancer Res Ther. 2018;14(7):1469–75. https://doi.org/10.4103/jcrt.JCRT_976_17.
Walcher L, Kistenmacher AK, Suo H, Kitte R, Dluczek S, Strauß A, Blaudszun AR, Yevsa T, Fricke S, Kossatz-Boehlert U. Cancer stem cells-origins and biomarkers: perspectives for targeted personalized therapies. Front Immunol. 2020;7(11):1280. https://doi.org/10.3389/fimmu.2020.01280.
Kreso A, Dick JE. Evolution of the cancer stem cell model. Cell Stem Cell. 2014;14(3):275–91. https://doi.org/10.1016/j.stem.2014.02.006.
Shibue T, Weinberg RA. EMT, CSCs, and drug resistance: the mechanistic link and clinical implications. Nat Rev Clin Oncol. 2017;14(10):611–29. https://doi.org/10.1038/nrclinonc.2017.44.
Dean M, Fojo T, Bates S. Tumour stem cells and drug resistance. Nat Rev Cancer. 2005;5(4):275–84. https://doi.org/10.1038/nrc1590.
Kudo M, Finn RS, Qin S, Han KH, Ikeda K, Piscaglia F, Baron A, Park JW, Han G, Jassem J, Blanc JF, Vogel A, Komov D, Evans TRJ, Lopez C, Dutcus C, Guo M, Saito K, Kraljevic S, Tamai T, Ren M, Cheng AL. Lenvatinib versus sorafenib in first-line treatment of patients with unresectable hepatocellular carcinoma: a randomised phase 3 non-inferiority trial. Lancet. 2018;391(10126):1163–73. https://doi.org/10.1016/S0140-6736(18)30207-1.
Al-Salama ZT, Syed YY, Scott LJ. Lenvatinib: a review in hepatocellular carcinoma. Drugs. 2019;79(6):665–74. https://doi.org/10.1007/s40265-019-01116-x.
Lu Y, Shen H, Huang W, He S, Chen J, Zhang D, Shen Y, Sun Y. Genome-scale CRISPR-Cas9 knockout screening in hepatocellular carcinoma with lenvatinib resistance. Cell Death Discov. 202118;7(1):359. https://doi.org/10.1038/s41420-021-00747-y. Erratum in: Cell Death Discov. 2022;8(1):74. PMID: 34795217; PMCID: PMC8602346.
Xing P, Dong H, Liu Q, Zhao T, Yao F, Xu Y, Chen B, Zheng X, Wu Y, Jin F, Li J. Upregulation of transmembrane 4 L6 family member 1 predicts poor prognosis in invasive breast cancer: a STROBE-compliant article. Medicine (Baltimore). 2017;96(52):e9476. https://doi.org/10.1097/MD.0000000000009476.
Lin C-I, Merley A, Sciuto TE, Li D, Dvorak AM, Melero-Martin JM, et al. TM4SF1: a new vascular therapeutic target in cancer. Angiogenesis. 2014;17(4):897–907. https://doi.org/10.1007/s10456-014-9437-2.
Roffler SR, Kao Y-R, Shih J-Y, Chen B-M, Chu Y-W, Yang P-C, et al. P-330 tumor-associated antigen L6 (TAL6) is involved in the invasion of human lung cancer cells. Lung Cancer. 2003;41:S175–S175.
Mazzocca A, Carloni V, Sciammetta SC, Cordella C, Pinzani M. Expression of transmembrane 4 superfamily (TM4SF) proteins and their role in hepatic stellate cell motility and wound healing migration. J Hepatol. 2002;37(3):322–30. https://doi.org/10.1016/S0168-8278(02)00175-7.
Fu XY, Zhou WB, Xu J. TM4SF1 facilitates non-small cell lung cancer progression through regulating YAP-TEAD pathway. Eur Rev Med Pharmacol Sci. 2020;24(4):1829–40. https://doi.org/10.26355/eurrev_202002_20361.
Gao C, Yao H, Liu H, Feng Y, Yang Z. TM4SF1 is a potential target for anti-invasion and metastasis in ovarian cancer. BMC Cancer. 2019;19(1):237. https://doi.org/10.1186/s12885-019-5417-7.
Yang JC, Zhang Y, He SJ, Li MM, Cai XL, Wang H, Xu LM, Cao J. TM4SF1 promotes metastasis of pancreatic cancer via regulating the expression of DDR1. Sci Rep. 2017;3(7):45895. https://doi.org/10.1038/srep45895.
Zeng Z, Shi Z, Liu Y, Zhao J, Lu Q, Guo J, Liu X, Huang D, Xu Q. HIF-1α-activated TM4SF1-AS1 promotes the proliferation, migration, and invasion of hepatocellular carcinoma cells by enhancing TM4SF1 expression. Biochem Biophys Res Commun. 2021;20(566):80–6. https://doi.org/10.1016/j.bbrc.2021.06.011.
Tang Q, Chen J, Di Z, Yuan W, Zhou Z, Liu Z, Han S, Liu Y, Ying G, Shu X, Di M. TM4SF1 promotes EMT and cancer stemness via the Wnt/β-catenin/SOX2 pathway in colorectal cancer. J Exp Clin Cancer Res. 2020;39(1):232. https://doi.org/10.1186/s13046-020-01690-z.
Jin H, Shi Y, Lv Y, Yuan S, Ramirez CFA, Lieftink C, Wang L, Wang S, Wang C, Dias MH, Jochems F, Yang Y, Bosma A, Hijmans EM, de Groot MHP, Vegna S, Cui D, Zhou Y, Ling J, Wang H, Guo Y, Zheng X, Isima N, Wu H, Sun C, Beijersbergen RL, Akkari L, Zhou W, Zhai B, Qin W, Bernards R. EGFR activation limits the response of liver cancer to lenvatinib. Nature. 2021;595(7869):730–4. https://doi.org/10.1038/s41586-021-03741-7.
Lian H, Han YP, Zhang YC, Zhao Y, Yan S, Li QF, Wang BC, Wang JJ, Meng W, Yang J, Wang QH, Mao WW, Ma J. Integrative analysis of gene expression and DNA methylation through one-class logistic regression machine learning identifies stemness features in medulloblastoma. Mol Oncol. 2019;13(10):2227–45.
Malta TM, Sokolov A, Gentles AJ, Burzykowski T, Poisson L, Weinstein JN, Kamińska B, Huelsken J, Omberg L, Gevaert O, Colaprico A, Czerwińska P, Mazurek S, Mishra L, Heyn H, Krasnitz A, Godwin AK, Lazar AJ, Cancer Genome Atlas Research Network, Stuart JM, Hoadley KA, Laird PW, Noushmehr H, Wiznerowicz M. Machine learning identifies stemness features associated with oncogenic dedifferentiation. Cell. 2018;173(2):338–54.
Cauchy F, Fuks D, Belghiti J. HCC: current surgical treatment concepts. Langenbecks Arch Surg. 2012;397(5):681–95. https://doi.org/10.1007/s00423-012-0911-2.
Morise Z, Kawabe N, Tomishige H, Nagata H, Kawase J, Arakawa S, Yoshida R, Isetani M. Recent advances in the surgical treatment of hepatocellular carcinoma. World J Gastroenterol. 2014;20(39):14381–92. https://doi.org/10.3748/wjg.v20.i39.14381.
Anwanwan D, Singh SK, Singh S, Saikam V, Singh R. Challenges in liver cancer and possible treatment approaches. Biochim Biophys Acta Rev Cancer. 2020;1873(1):188314. https://doi.org/10.1016/j.bbcan.2019.188314.
Mohr OL. Character changes caused by mutation of an entire region of a chromosome in drosophila. Genetics. 1919;4(3):275–82. https://doi.org/10.1093/genetics/4.3.275.
Andersson ER, Sandberg R, Lendahl U. Notch signaling: simplicity in design, versatility in function. Development. 2011;138(17):3593–612. https://doi.org/10.1242/dev.063610.
Capaccione KM, Pine SR. The Notch signaling pathway as a mediator of tumor survival. Carcinogenesis. 2013;34(7):1420–30. https://doi.org/10.1093/carcin/bgt127.
Aster JC, Pear WS, Blacklow SC. The varied roles of notch in cancer. Annu Rev Pathol. 2017;12:245–75. https://doi.org/10.1146/annurev-pathol-052016-100127.
Lu J, Ye X, Fan F, Xia L, Bhattacharya R, Bellister S, Tozzi F, Sceusi E, Zhou Y, Tachibana I, Maru DM, Hawke DH, Rak J, Mani SA, Zweidler-McKay P, Ellis LM. Endothelial cells promote the colorectal cancer stem cell phenotype through a soluble form of Jagged-1. Cancer Cell. 2013;23(2):171–85. https://doi.org/10.1016/j.ccr.2012.12.021.
Pecci A, Ma X, Savoia A, Adelstein RS. MYH9: structure, functions and role of non-muscle myosin IIA in human disease. Gene. 2018;664:152–67. https://doi.org/10.1016/j.gene.2018.04.048.
Zhong Y, Long T, Gu CS, Tang JY, Gao LF, Zhu JX, Hu ZY, Wang X, Ma YD, Ding YQ, Li ZG, Wang XY. MYH9-dependent polarization of ATG9B promotes colorectal cancer metastasis by accelerating focal adhesion assembly. Cell Death Differ. 2021;28(12):3251–69. https://doi.org/10.1038/s41418-021-00813-z.
Liu L, Yi J, Deng X, Yuan J, Zhou B, Lin Z, Zeng Z. MYH9 overexpression correlates with clinicopathological parameters and poor prognosis of epithelial ovarian cancer. Oncol Lett. 2019;18(2):1049–56. https://doi.org/10.3892/ol.2019.10406.
Savoia A, Pecci A. MYH9-Related Disease. 2008 Nov 20 [updated 2021 Feb 18]. In: Adam MP, Everman DB, Mirzaa GM, Pagon RA, Wallace SE, Bean LJH, Gripp KW, Amemiya A, editors. GeneReviews® [Internet]. Seattle (WA): University of Washington, Seattle; 1993–2022. PMID: 20301740.
Lin X, Li AM, Li YH, Luo RC, Zou YJ, Liu YY, Liu C, Xie YY, Zuo S, Liu Z, Liu Z, Fang WY. Silencing MYH9 blocks HBx-induced GSK3β ubiquitination and degradation to inhibit tumor stemness in hepatocellular carcinoma. Signal Transduct Target Ther. 2020;5(1):13. https://doi.org/10.1038/s41392-020-0111-4.
Liu Y, Sun L, Guo H, Zhou S, Wang C, Ji C, Meng F, Liang S, Zhang B, Yuan Y, Ma K, Li X, Guo X, Cui T, Zhang N, Wang J, Liu Y, Liu L. Targeting SLP2-mediated lipid metabolism reprograming restricts proliferation and metastasis of hepatocellular carcinoma and promotes sensitivity to Lenvatinib. Oncogene. 2023;42(5):374–88.
Xia S, Pan Y, Liang Y, Xu J, Cai X. The microenvironmental and metabolic aspects of sorafenib resistance in hepatocellular carcinoma. EBioMedicine. 2020;51: 102610.
Mendez-Blanco C, Fondevila F, Garcia-Palomo A, Gonzalez-Gallego J, Mauriz JL. Sorafenib resistance in hepatocarcinoma: role of hypoxia-inducible factors. Exp Mol Med. 2018;50(10):1–9.
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Ethics approval and consent to participate
The ethics committee approval from the authors’ institution and patients’ informed consent have been obtained for this study.
All authors have no competing interests. On behalf of all authors, the corresponding author states that there is no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
The sequences used in this study. Table S2. The primers used in this study.
Additional file 2: Fig. S1.
The bioinformatics databases showed that TM4SF1 is highly expressed in HCC and correlated with poor prognosis. A Research on the TIMER2.0 database showed that hepatocellular cancer had significant levels of TM4SF1 expression. B Information from the GEPIA database revealed that the expression of TM4SF1 was higher in tissues with hepatocellular carcinoma than in healthy tissues. C Analysis of the GEPIA database and TIMER2.0 database showed that the cumulative survival rate, overall survival rate and disease-free survival rate of patients with low TM4SF1 expression were better than those with high TM4SF1 expression. D Expression of TM4SF1 in hepatocellular carcinoma cell lines was analyzed using the TCGA database.
Additional file 3: Fig. S2.
Schematic diagram of the molecular mechanism of TM4SF1 upregulates MYH9 to activate the NOTCH pathway to promote cancer stemness and Lenvatinib resistance in HCC.
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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Yang, Sb., Zhou, Zh., Lei, J. et al. TM4SF1 upregulates MYH9 to activate the NOTCH pathway to promote cancer stemness and lenvatinib resistance in HCC. Biol Direct 18, 18 (2023). https://doi.org/10.1186/s13062-023-00376-8