Development and validation of prognostic and diagnostic model for pancreatic ductal adenocarcinoma based on scRNA-seq and bulk-seq datasets.

The 5-year overall survival (OS) of pancreatic ductal adenocarcinoma (PDAC) is only 10%, partly owing to the lack of reliable diagnostic and prognostic biomarkers. The raw gene-cell matrix for single-cell RNA-seq (scRNA-seq) analysis was downloaded from the GSA database. We drew cell atlas for PDAC and normal pancreatic tissues. The inferCNV analysis was used to distinguish tumor cells from normal ductal cells. We identified differential expression genes (DEGs) by comparing tumor cells and normal ductal cells. The common DEGs were used to conduct prognostic and diagnostic model using univariate and multivariate Cox or logistic regression analysis. Four genes, MET, KLK10, PSMB9 and ITGB6, were utilized to create risk score formula to predict OS and to establish diagnostic model for PDAC. Finally, we drew an easy-to-use nomogram to predict 2-year and 3-year OSs. In conclusion, we developed and validated the prognostic and diagnostic model for PDAC based on scRNA-seq and bulk-seq datasets.

Cross-species scRNA-seq reveals the cellular landscape of retina and early alterations in type 2 diabetes mice.

Single-cell RNA sequencing (scRNA-seq) analysis have provided an unprecedented resolution for the studies on diabetic retinopathy (DR). However, the early changes in the retina in diabetes remain unclear. A total of 8 human and mouse scRNA-seq datasets, containing 276,402 cells were analyzed individually to comprehensively delineate the retinal cell atlas. The neural retinas were isolated from the type 2 diabetes (T2D) and control mice, and scRNA-seq analysis was conducted to evaluate the early effects of diabetes on the retina. Bipolar cell (BC) heterogeneity were identified. We found some stable BCs across multiple datasets, and explored their biological functions. A new RBC subtype (Car8_RBC) in the mouse retina was validated using the multi-color immunohistochemistry. AC149090.1 was significantly upregulated in the rod cells, ON cone BCs (CBCs), OFF CBCs, and RBCs in T2D mice. Additionally, the interneurons, especially BCs, were the most vulnerable cells to diabetes by integrating scRNA-seq and genome-wide association studies (GWAS) analyses. In conclusion, this study delineated a cross-species retinal cell atlas and uncovered the early pathological alterations in the retina of T2D mice.

Immune profiling and prognostic model of pancreatic cancer using quantitative pathology and single-cell RNA sequencing.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) has a complex tumor immune microenvironment (TIME), the clinical value of which remains elusive. This study aimed to delineate the immune landscape of PDAC and determine the clinical value of immune features in TIME. METHODS: Univariable and multivariable Cox regression analyses were performed to evaluate the clinical value of immune features and establish a new prognostic model. We also conducted single-cell RNA sequencing (scRNA-seq) to further characterize the immune profiles of PDAC and explore cell-to-cell interactions. RESULTS: There was a significant difference in the immune profiles between PDAC and adjacent noncancerous tissues. Several novel immune features were captured by quantitative pathological analysis on multiplex immunohistochemistry (mIHC), some of which were significantly correlated with the prognosis of patients with PDAC. A risk score-based prognostic model was established based on these immune features. We also constructed a user-friendly nomogram plot to predict the overall survival (OS) of patients by combining the risk score and clinicopathological features. Both mIHC and scRNA-seq analysis revealed PD-L1 expression was low in PDAC. We found that PD1 + cells were distributed in different T cell subpopulations, and were not enriched in a specific subpopulation. In addition, there were other conserved receptor-ligand pairs (CCL5-SDC1/4) besides the PD1-PD-L1 interaction between PD1 + T cells and PD-L1 + tumor cells. CONCLUSION: Our findings reveal the immune landscape of PDAC and highlight the significant value of the combined application of mIHC and scRNA-seq for uncovering TIME, which might provide new clues for developing immunotherapy combination strategies.

Clinical characteristics of liver injury induced by immune checkpoint inhibitors in patients with advanced biliary tract carcinoma.

Immune-related liver injuries are closely associated with the liver's fundamental state. Patients with advanced biliary tract carcinoma (BTC) have poor liver function. We evaluated the clinical data of immune-related liver injury in patients with advanced BTC and gastric cancer (GC) during immune checkpoint inhibitor (ICI) treatment between February 2019 and July 2022 at Peking University First Hospital. Twenty-five patients with advanced BTC were identified. Fifteen patients (60%) experienced immune-related liver injury during ICI treatment. We also evaluated the clinical status of patients with GC in another group receiving immunotherapy. The results demonstrated that the incidence of immune-related liver injury was higher in patients with BTC than in GC cancer (p=0.040). Multivariate analysis suggested that the type of malignant tumor and baseline liver function status were high-risk factors for grade 2 and higher immune-related liver injuries. Two patients were diagnosed with immune-related cholangitis. Both biliary enzymes can be decreased to a certain degree by corticosteroid and ursodeoxycholic acid (UDCA) therapy but are difficult to reduce to normal levels. Liver function normalized, and symptoms improved after local treatment for cholestasis (stent implantation or PTBD). We observed a higher incidence of immune-related liver injury after ICI treatment in patients with advanced BTC. Effect of baseline liver function on the incidence of liver injury associated with immunotherapy. Interventional therapy provides rapid relief from cholestasis and is an indispensable and effective approach to the treatment of immune-related cholangitis.

Efficacy and safety of early drain removal following pancreatic resections: a meta-analysis.

BACKGROUND: No consensus was reached with regard to the effect of EDR on postoperative outcomes after pancreatic surgery. The meta-analysis was designed to explore the efficacy and safety of early drain removal (EDR). METHODS: Systematic literature search was performed. Data extraction and correction were performed by three researchers. For dichotomous and continuous outcomes, we calculated the pooled risk difference and mean difference with 95% confidence intervals, respectively. The heterogeneity of included studies was evaluated using Cochran's Q and I2 test. The stratified analyses of pancreaticoduodenectomy (PD) and distal pancreatectomy (DP) were performed. RESULTS: A total of 10 studies including 3 RCTs and 7 non RCTs were included for meta-analysis, among which 1780 patients with EDR and 5613 patients with late drain removal (LDR) were enrolled. The meta-analysis of both all the available studies and studies only with selected low risk patients indicated that EDR group had significantly lower incidences of Grade B/C postoperative pancreatic fistula (POPF) and total complications for both PD and DP. However, no advantages of EDR were observed in the meta-analysis of the 3 RCTs. In addition, EDR was associated with a lower incidence of intra-abdominal infection after PD. While for DP, EDR group had decreased risk of delayed gastric emptying and re-operation, and shorter postoperative in-hospital stay. CONCLUSIONS: The meta-analysis demonstrates that EDR is effective and safe for both PD and DP considering POPF and total complications, especially for patients with low concentration of postoperative drain fluid amylase.

Early Drain Removal is Safe in Patients With Low or Intermediate Risk of Pancreatic Fistula After Pancreaticoduodenectomy: A Multicenter, Randomized Controlled Trial.

OBJECTIVE: This multicenter randomized controlled trial was designed to test the hypothesis that early drain removal (EDR) could decrease the incidence of grade 2 to 4 complications for patients undoing pancreaticoduodenectomy (PD) with low or intermediate risk of postoperative pancreatic fistula (POPF). BACKGROUND: The safety and effects of EDR on postoperative complications after PD are still controversial. METHODS: A multicenter randomized controlled trial at 6 tertiary referral hospitals was carried out (NCT03055676). Patients who met the inclusion criteria, including drain amylase level less than 5000 U/L on postoperative day (POD) 1 and POD 3, and drain output less than 300 mL per day within 3 days after surgery, were enrolled. Patients were then randomized to the EDR group or the routine drain removal (RDR) group. In the EDR group, all drainage tubes were removed on POD3. In the RDR group, drainage tubes were removed on POD 5 or beyond. Primary outcome was the incidence of Clavien-Dindo grade 2 to 4 complications. Secondary outcomes were comprehensive complication index, grade B/C POPF, total medical expenses and postoperative in-hospital stay etc, within 90 days after surgery. RESULTS: A total of 692 patients were screened, and 312 patients were eligible for randomization. Baseline characteristics were well balanced between the 2 groups and 96.8% of these 312 patients had low or intermediate risk of POPF, according to the 10-point fistula risk score. A total of 20.5% of the patients in the EDR group suffered at least 1 grade 2 to 4 complication, versus 26.3% in the RDR group (P = 0.229). Multi-variate analysis showed older age (>65 years old) and blood transfusion were independent risk factors for grade 2 to 4 complications. The rate of grade B/C POPF was low in either group (3.8% vs 6.4%, P = 0.305). The comprehensive complication index of the 2 groups was also comparable (20.9 vs 20.9, P = 0.253). Total medical expenses were not significantly different. Postoperative in-hospital stay was clinically similar (15 days vs 16 days, P = 0.010). CONCLUSIONS: Nearly half of the patients undergoing PD met the inclusion criteria, predicting low incidence of grade B/C POPF and major complications. EDR was safe in these patients but did not significantly decrease major complications.

Predictive Model of Early Death of Resectable Pancreatic Ductal Adenocarcinoma After Curative Resection: A SEER-Based Study.

OBJECTIVE: This study aims to determine the factors that predict early death and establish a predictive model for early death by analyzing clinical characteristics of patients with resectable pancreatic ductal adenocarcinoma (R-PDAC) who die early after radical surgery. MATERIALS AND METHODS: This was a retrospective study of patients who underwent radical surgical resection for R-PDAC in the Surveillance, Epidemiology, and End Results (SEER) database. Patients with overall survival ≤ 12 months were assigned as early death group and above 1 year as the late death group. Univariate and multivariate logistic regression was conducted to identify factors significantly associated with early death. An early death predictive model was constructed based on the identified independent risk factors. RESULTS: A total of 9695 patients were analyzed, and the total incidence of early death was 30.72%. Multivariable analysis showed that factors significantly associated with early death included age at diagnosis, race, marital status, tumor location, tumor size, tumor grade, number of positive lymph nodes, number of examined lymph nodes, positive lymph node ratio, chemotherapy, and radiotherapy. The predictive model showed good discrimination with a C-index of 0.722 (95% confidence interval: 0.711-0.733) and convincing calibration. CONCLUSIONS: We developed a predictive model that may be easily applied to patients with R-PDAC after radical resection to predict the chance of death within 1 year. For patients with high risk of early death, neoadjuvant therapy should be considered. Even after radical resection, more aggressive adjuvant chemotherapy (with or without combined radiotherapy) must be used to minimize the chance of early death.

Crosstalk between Pancreatic Cancer Cells and Cancer-Associated Fibroblasts in the Tumor Microenvironment Mediated by Exosomal MicroRNAs.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most malignant digestive tumors, characterized by a low rate of early diagnosis, strong invasiveness, and early metastasis. The abundant stromal cells, dense extracellular matrix, and lack of blood supply in PDAC limit the penetration of chemotherapeutic drugs, resulting in poor efficacy of the current treatment regimens. Cancer-associated fibroblasts (CAFs) are the major stromal cells in the tumor microenvironment. Tumor cells can secrete exosomes to promote the generation of activated CAFs, meanwhile exosomes secreted by CAFs help promote tumor progression. The aberrant expression of miRNAs in exosomes is involved in the interaction between tumor cells and CAFs, which provides the possibility for the application of exosomal miRNAs in the diagnosis and treatment of PDAC. The current article reviews the mechanism of exosomal miRNAs in the crosstalk between PDAC cells and CAFs in the tumor microenvironment, in order to improve the understanding of TME regulation and provide evidence for designing diagnostic and therapeutic targets against exosome miRNA in human PDAC.

High mobility group AT-hook 2 promotes tumorigenicity of pancreatic cancer cells via upregulating ANLN.

HMGA2 is associated with the regulation of cellular biological processes in various human disorders and cancer progression, yet little is known about how HMGA2 controls tumorigenesis. This study uncovered the mechanism of HMGA2-mediated regulation of tumorigenicity in pancreatic cancer. We showed that HMGA2 was highly expressed in pancreatic cancer cells and correlated with poor prognosis. HMGA2 expression knockdown inhibited the tumorigenicity of pancreatic cancer cells. Conversely, overexpression of HMGA2 promoted tumorigenicity. Combination of ChIP-Seq, RNA-Seq and dual-luciferase reporter assays revealed HMGA2 could directly regulate ANLN expression. Furthermore, we found ANLN could mediate the HMGA2-induced effects on pancreatic cancer cells. The identification of the regulatory mechanism of HMGA2 and ANLN will provide insights into the progression for human pancreatic cancer.

A modified alternative fistula risk score (a-FRS) obtained from the computed tomography enhancement pattern of the pancreatic parenchyma predicts pancreatic fistula after pancreatoduodenectomy.

BACKGROUND: Alternative fistula risk score (a-FRS) is useful to predict clinically relevant postoperative pancreatic fistula (CR-POPF) after pancreatoduodenectomy (PD). METHODS: Clinical data from 239 patients undergoing PD were collected. The CT value of the pancreatic parenchyma was measured in the nonenhanced (N), arterial (A), portal venous (P), and late (L) phases. The A/N, A/P, P/L and A/L ratios were calculated and their correlation with CR-POPF were analyzed. By replacing pancreatic texture with the best CT attenuation ratio, a modified a-FRS was developed. RESULTS: Forty-seven patients developed CR-POPF. The A/P ratio (P < 0.001), P/L ratio (P = 0.002) and A/L ratio (P < 0.001) were significantly higher in the CR-POPF group. The A/L ratio performed best in predicting CR-POPF (AUC: 0.803) and the cut-off value is 1.36. A/L ratio >1.36 (P < 0.001), body mass index (P = 0.005) and duct diameter (P = 0.037) were independently associated with CR-POPF. By replacing soft texture with an A/L ratio >1.36, a modified a-FRS was developed and performed better than the a-FRS (AUC: 0.823 vs 0.748, P = 0.006) in predicting CR-POPF. CONCLUSIONS: The modified a-FRS is an objective and preoperative model for predicting the occurrence of CR-POPF after PD.

Validation and modification of the AJCC 8th TNM staging system for pancreatic ductal adenocarcinoma in a Chinese cohort: A nationwide pancreas data center analysis.

OBJECTIVE: To validate the 8th edition of the American Joint Committee on Cancer (AJCC) staging system for pancreatic ductal adenocarcinoma (PDAC) in a Chinese cohort of radically resected patients and to develop a refined staging system for PDAC. METHODS: Data were collected from the China Pancreas Data Center (CPDC) for patients with resected PDAC in 2016 and 2017, and cancer-specific survival (CSS) was evaluated using the Kaplan-Meier method and log-rank test. Univariate and multivariate analyses based on Cox regression were performed to identify prognostic factors. The recursive partitioning analysis (RPA), Kaplan-Meier method, and log-rank test were performed on the training dataset to generate a proposed modification for the 8th TNM staging system utilizing the preoperative carbohydrate antigen (CA)19-9 level. Validation was performed for both staging systems in the validation cohort. RESULTS: A total of 1,676 PDAC patients were retrieved, and the median CSS was significantly different between the 8th TNM groupings, with no significant difference in survival between stage IB and IIA. The analysis of T and N stages demonstrated a better prognostic value in the N category. Multivariate analysis showed that the preoperative serum CA19-9 level was the strongest prognostic indicator among all the independent risk factors. All patients with CA19-9 >500 U/mL had similar survival, and we proposed a new staging system by combining IB and IIA and stratifying all patients with high CA19-9 into stage III. The modified staging system had a better performance for predicting CSS than the 8th AJCC staging scheme. CONCLUSIONS: The 8th AJCC staging system for PDAC is suitable for a Chinese cohort of resected patients, and the N category has a better prognostic value than the T category. Our modified staging system has superior accuracy in predicting survival than the 8th AJCC TNM staging system.

Features and Risk Factors for Paraneoplastic Autoimmune Multiorgan Syndrome in 145 Chinese Patients.

Paraneoplastic autoimmune multiorgan syndrome is a complex and deadly disease. We retrospectively reviewed the clinical features and risk factors for paraneoplastic autoimmune multiorgan syndrome in 145 Chinese patients. The most common neoplasm was Castleman disease (56%), and patients with Castleman disease tended to be younger (≤ 42 years old: 83% vs. 29%) and to have a greater proportions of lichen planus-like lesions (47% vs. 27%) and bronchiolitis obliterans (49% vs. 29%), compared to other neoplasm-associated patients. Among all 145 patients in the study, the survival rates were 84% at 1 year, 65% at 3 years, and 54% at 5 years. Kaplan-Meier curve analysis revealed that mortality was associated with older age (> 42 years), neoplasm type, labial lesions, and larger skin lesion area (> 17.5% of the body surface area). However, only older age and larger skin lesion area were independent factors associated with mortality in multivariate analysis. We suggest that patients with Castleman disease and paraneoplastic autoimmune multiorgan syndrome have many unique characteristics and the underlying risk factors for death require further exploration.

Chaperonin-GroEL as a Smart Hydrophobic Drug Delivery and Tumor Targeting Molecular Machine for Tumor Therapy.

The targeted delivery of hydrophobic therapeutic drugs to tumors is one of the major challenges in drug development. The use of natural proteins as drug delivery vehicles holds great promise due to various functionalities of proteins. In the current study, we exploited a natural protein, GroEL, which possesses a double layer cage structure, as a hydrophobic drug container, which is switchable by ATP binding to a hydrophilic status, to design a novel and intelligent hydrophobic drug delivery molecular machine with a controlled drug release profile. When loaded with the hydrophobic antitumor drug, Doxorubicin (Dox), GroEL was able to shield the drug from the aqueous phase of blood, releasing the drug once in the presence of a critical concentration of ATP at the tumor site. Unexpectedly, we found that GroEL has a specific affinity for the cell structural protein, plectin, which is expressed at abnormally elevated levels on the membranes of tumor cells but not in normal cells. This finding, in combination with the ATP sensitivity, makes GroEL a superior natural tumor targeting nanocarrier. Our data show that GroEL-Dox is able to effectively, and highly selectively, deliver the hydrophobic drug to fast growing tumors without overt adverse effects on the major organs. GroEL is therefore a promising drug delivery platform that can overcome the obstacles to hydrophobic drug targeting and delivery.

Demethylzeylasteral (ZST93) inhibits cell growth and enhances cell chemosensitivity to gemcitabine in human pancreatic cancer cells via apoptotic and autophagic pathways.

The overall 5-year survival rate of patients with human pancreatic cancer remains less than 8% because of its aggressive growth, early metastasis and resistance to conventional chemoradiotherapy. It is essential to develop innovative and effective therapeutic agents to improve its prognosis. Demethylzeylasteral (ZST93) is a novel triterpenoid monomer extracted from the xylem of Tripterygium roots. Our study aimed to assess the effects of ZST93 on cell proliferation and its role in the chemosensitivity to gemcitabine in human pancreatic cancer cells. The effects of ZST93 on cancer cell proliferation, cell cycle distribution, apoptosis and autophagy were evaluated in various human pancreatic cancer cell lines, and the antitumor effects of ZST93 alone and in combination with gemcitabine were identified in a xenograft mouse model. The results showed that ZST93 could inhibit the proliferation of pancreatic cancer cells and arrest cell cycle at G0/G1 phase by regulating the expression of Cyclin D1 and Cyclin A2. Moreover, ZST93 killed pancreatic cancer cells through two different mechanisms: inducing autophagic cell death at low concentrations and apoptotic cell death at high concentrations. Furthermore, ZST93 could enhance the chemosensitivity of pancreatic cancer cells to gemcitabine both in vitro and in vivo through modulation of the cross talk between autophagy and apoptosis. ZST93 is a potential therapeutic agent for developing novel therapeutic strategies in human pancreatic cancer.

[Risk factors of postoperative pancreatic fistula after pancreaticoduodenectomy and its predictive score].

OBJECTIVE: To develop and test a scoring system to predict the risks of postoperative pancreatic fistula (POPF) after pancreaticoduodenectomy(PD). METHODS: Clinic data and postoperative complications of the 445 consecutive patients who underwent a PD procedure between January 2008 and April 2015 in Peking University First Hospital were retrospectively collected and analyzed.The patients were randomly selected to modelling and validation sets at a ratio of 3∶1, respectively.The patient data were subjected to univariate and multivariate analysis in the modelling set of patients.A score predictive of POPF was designed and tested in the validation set. RESULTS: POPF occurred in 88 of 334 patients(26.4%) in the modelling set.The multivariate analysis showed that body mass index (BMI, P<0.01) and pancreatic duct width(P=0.001) are associated with POPF independently.A risk score to predict POPF was constructed based on these factors and successfully tested.The area under the receiver operating characteristic curve were 0.829(95% CI: 0.777-0.881) on the modelling set and 0.885(95% CI: 0.825-0.945) on the validation set, respectively. CONCLUSIONS: BMI and pancreatic duct width were associated with POPF after PD. The preoperative assessment of a patient's risk for POPF is feasible.The present risk score is a valid tool to predict POPF in patients undergoing PD, to make the selection on anastomosis types, and to take precautions against POPF.

Feasibility of novel intraoperative navigation for anatomical liver resection using real-time virtual sonography combined with indocyanine green fluorescent imaging technology.

To analyze the feasibility and clinical effect of novel intraoperative navigation of real-time virtual sonography (RVS) combined with indocyanine green (ICG) fluorescent imaging technology in anatomical liver resection (ALR) for hepatocellular carcinoma. The clinical data of 41 patients who underwent ALR using RVS intraoperative navigation combined with ICG fluorescent imaging technology in the Department of Hepatobiliary Surgery of Peking University International Hospital from January 2020 to May 2022 were retrospectively analyzed. RVS was applied to guide the surgical plane through fusing real-time intraoperative ultrasound images with corresponding preoperative CT or MRI images. Operation methods, operation time, intraoperative blood loss, operative margin, hospital stay and postoperative complications were analyzed. The 1-year overall survival rate and tumor-free survival rate of patients were followed up by outpatient review or telephone calls. ALR surgery was performed on each of 41 patients. There were no deaths during perioperative period and postoperative complications occurred in 7 cases (17.1%). The postoperative pathological examinations demonstrated all cases of hepatocellular carcinoma and negative operative margins. The 41 patients were followed up for 12 to 20 months, with a median follow-up time of 14 months. The overall survival rate 1 year after surgery was 100.0% (41/41), 3 patients (7.3%) experienced tumor recurrence, and the tumor-free survival rate of 1 year after surgery was 92.7% (38/41). In conclusion, novel intraoperative navigation of RVS combined with ICG fluorescent imaging technology is safe and feasible in anatomical segmental hepatectomy of hepatocellular carcinoma.

Single-cell transcriptome profiling of primary tumors and paired organoids of pancreatobiliary cancer.

Single-cell RNA-seq (scRNA-seq) and cancer organoid model have shown promise in investigating tumor microenvironment heterogeneity and facilitating chemotherapeutic drug testing to inform treatment selection. It is still unknown whether the scRNA-seq results based on organoid can faithfully reflect the heterogeneity of primary pancreatobiliary cancer. To reveal the similarities and differences between primary tumors and their matched organoids at transcriptome level, we conducted scRNA-seq for paired primary tumors and organoids from one cholangiocarcinoma (CCA) and two pancreatic ductal adenocarcinoma (PDAC) patients. We identified inter-patient and intra-tumor heterogeneity and found that the organoids retained copy number variation (CNV) patterns of primary tumors. There was no significant difference in cancer stem cell (CSC) properties between the primary tumors and the organoids, whereas organoid from one PDAC case had increased mesenchymal-score and decreased epithelial-score compared with the primary tumors. All organoids showed a transition tendency from the classical subtype to the basal-like subtype in the transcriptional level. Organoids and primary tumors differed in metabolic and unfolded protein response (UPR) signatures. In addition, we revealed the heterogeneity of cancer associated fibroblasts (CAFs) and T cells, and explored the developmental trajectory of T cells. Our findings facilitate further understanding of organoid model and confirm its application prospects in pancreatobiliary cancer.

HDAC Inhibition Sensitize Hepatocellular Carcinoma to Lenvatinib via Suppressing AKT Activation.

Hepatocellular carcinoma (HCC) is a deadly malignancy with limited treatment options. As a first-line treatment for advanced HCC, Lenvatinib has been applicated in clinic since 2018. Resistance to Lenvatinib, however, has severely restricted the clinical benefits of this drug. Therefore, it is urgent to explore the potential resistance mechanisms of Lenvatinib and identify appropriate methods to reduce resistance for the treatment of HCC. We identified SAHA, a HDAC inhibitor, to have effective anti-tumor activity against Lenvatinib-resistant HCC organoids by screening a customized drug library. Mechanism analysis revealed that SAHA upregulates PTEN expression and suppresses AKT signaling, which contributes to reversing Lenvatinib resistance in liver cancer cells. Furthermore, combinational application of Lenvatinib and HDAC inhibitor or AKT inhibitor synergistically inhibits HCC cell proliferation and induces cell apoptosis. Finally, we confirmed the synergistic effects of Lenvatinib and SAHA, or AZD5363 in primary liver cancer patient derived organoids. Collectively, these findings may enable the development of Lenvatinib combination therapies for HCC.