Patient-derived organoid models to decode liver pathophysiology

Liver diseases represent a growing global health challenge, and the increasing prevalence of obesity and metabolic disorders is set to exacerbate this crisis. To meet evolving regulatory demands, patient-specific in vitro liver models are essential for understanding disease mechanisms and developing new therapeutic approaches. Organoid models, which faithfully recapitulate liver biology, can be established from both non-malignant and malignant liver tissues, offering insight into various liver conditions, from acute injuries to chronic diseases and cancer. Improved understanding of liver microenvironments, innovative biomaterials, and advanced imaging techniques now facilitate comprehensive and unbiased data analysis, paving the way for personalised medicine. In this review, we discuss state-of-the-art patient-derived liver organoid models, recent technological advancements, and strategies to enhance their clinical impact.

Diagnostic lncRNA high expression for liver patients prognosis and medication guidance: a systematic review and meta-analysis

BACKGROUND

The study of long non-coding RNAs (lncRNAs) has gained significant attention due to their roles in regulating gene expression and their potential as diagnostic biomarkers. This systematic review and meta-analysis aimed to evaluate the diagnostic value of high-expression lncRNAs in liver disease patients, including those with hepatitis, cirrhosis, and hepatocellular carcinoma (HCC).

METHODS

A comprehensive literature search was conducted across multiple electronic databases, including PubMed, Embase, Web of Science, and Cochrane Library, up to July 2024. Studies were included if they investigated the expression of lncRNAs in liver disease patients and evaluated their diagnostic performance. The Quality Assessment of Diagnostic Accuracy Studies-2 (QUADAS-2) tool was used to assess the quality of included studies. Pooled sensitivity, specificity, diagnostic odds ratios (DOR), and summary receiver operating characteristic (SROC) curves were calculated using a bivariate random-effects model.

RESULTS

Nine studies involving 888 samples were included in the meta-analysis. The pooled hazard ratio (HR) for overall survival (OS) was 2.01 (95% CI: 1.71-2.36), indicating a significant association between high lncRNA expression and poor liver disease outcomes. Subgroup analyses revealed a pooled odds ratio (OR) of 1.99 (95% CI: 1.53-2.60) for tissue samples and 8.62 (95% CI: 1.16-63.71) for blood samples, suggesting a stronger diagnostic value for blood-based lncRNAs. The funnel plots indicated minimal publication bias, and sensitivity analyses confirmed the robustness of the findings.

CONCLUSION

High-expression lncRNAs show significant potential as diagnostic biomarkers for liver diseases, offering non-invasive, accurate, and timely diagnostic information. Despite the promising results, further research is needed to standardize detection methods, elucidate the biological functions of lncRNAs, and validate their clinical utility in diverse patient populations. Integrating lncRNA biomarkers with traditional diagnostic approaches could enhance diagnostic accuracy and improve patient management and outcomes in liver disease.

Intratumoral CXCR4hi neutrophils display ferroptotic and immunosuppressive signatures in hepatoblastoma

Pediatric hepatoblastoma (HB) is the most common primary liver malignancy in infants and children. With great diversity and plasticity, tumor-infiltrating neutrophils were one of the most determining factors for poor prognosis in many malignant tumors. In this study, through bulk RNA sequencing for sorted blood and tumor-infiltrated neutrophils and comparison of neutrophils in tumor and para-tumor tissue by single-cell sequencing, we found that intratumoral neutrophils were composed of heterogenous functional populations at different development stages. Our study showed that terminally differentiated neutrophils with active ferroptosis prevailed in tumor tissue, whereas, in para-tumor, pre-fate naïve neutrophils were dominant and ferroptotic neutrophils dispersed in a broad spectrum of cell maturation. Gene profiling and in vitro T-cell coculture experiment confirmed that one of main functional intratumoral neutrophils was mainly immunosuppressive, which relied on the activation of ferroptosis. Combining the bulk RNA-seq, scRNA-seq data, and immunochemistry staining of tumor samples, CXCL12/CXCR4 chemotaxis pathway was suggested to mediate the migration of neutrophils in tumors as CXCR4 highly expressed by intratumoral neutrophils and its ligand CXCL12 expressed much higher level in tumor than that in para-tumor. Moreover, our study pinpointed that infiltrated CXCR4hi neutrophils, regardless of their differential distribution of cell maturation status in HB tumor and para-tumor regions, were the genuine perpetrators for immune suppression. Our data characterized the ferroptosis-dependent immunosuppression energized by intratumoral CXCR4 expression neutrophils and suggest a potential cell target for cancer immunotherapies.

Chromosome 8p engineering reveals increased metastatic potential targetable by patient-specific synthetic lethality in liver cancer

Large-scale chromosomal aberrations are prevalent in human cancer, but their function remains poorly understood. We established chromosome-engineered hepatocellular carcinoma cell lines using CRISPR-Cas9 genome editing. A 33-mega-base pair region on chromosome 8p (chr8p) was heterozygously deleted, mimicking a frequently observed chromosomal deletion. Using this isogenic model system, we delineated the functional consequences of chr8p loss and its impact on metastatic behavior and patient survival. We found that metastasis-associated genes on chr8p act in concert to induce an aggressive and invasive phenotype characteristic for chr8p-deleted tumors. Genome-wide CRISPR-Cas9 viability screening in isogenic chr8p-deleted cells served as a powerful tool to find previously unidentified synthetic lethal targets and vulnerabilities accompanying patient-specific chromosomal alterations. Using this target identification strategy, we showed that chr8p deletion sensitizes tumor cells to targeting of the reactive oxygen sanitizing enzyme Nudix hydrolase 17. Thus, chromosomal engineering allowed for the identification of novel synthetic lethalities specific to chr8p loss of heterozygosity.

MicroRNA-27a-3p targets FoxO signalling to induce tumour-like phenotypes in bile duct cells

BACKGROUND & AIMS

Cholangiocarcinoma (CCA) is a heterogeneous and lethal malignancy, the molecular origins of which remain poorly understood. MicroRNAs (miRs) target diverse signalling pathways, functioning as potent epigenetic regulators of transcriptional output. We aimed to characterise miRNome dysregulation in CCA, including its impact on transcriptome homeostasis and cell behaviour.

METHODS

Small RNA sequencing was performed on 119 resected CCAs, 63 surrounding liver tissues, and 22 normal livers. High-throughput miR mimic screens were performed in three primary human cholangiocyte cultures. Integration of patient transcriptomes and miRseq together with miR screening data identified an oncogenic miR for characterization. MiR-mRNA interactions were investigated by a luciferase assay. MiR-CRISPR knockout cells were generated and phenotypically characterized in vitro (proliferation, migration, colony, mitochondrial function, glycolysis) and in vivo using subcutaneous xenografts.

RESULTS

In total, 13% (140/1,049) of detected miRs were differentially expressed between CCA and surrounding liver tissues, including 135 that were upregulated in tumours. CCA tissues were characterised by higher miRNome heterogeneity and miR biogenesis pathway expression. Unsupervised hierarchical clustering of tumour miRNomes identified three subgroups, including distal CCA-enriched and IDH1 mutant-enriched subgroups. High-throughput screening of miR mimics uncovered 71 miRs that consistently increased proliferation of three primary cholangiocyte models and were upregulated in CCA tissues regardless of anatomical location, among which only miR-27a-3p had consistently increased expression and activity in several cohorts. FoxO signalling was predominantly downregulated by miR-27a-3p in CCA, partially through targeting of FOXO1. MiR-27a knockout increased FOXO1 levels in vitro and in vivo, impeding tumour behaviour and growth.

CONCLUSIONS

The miRNomes of CCA tissues are highly remodelled, impacting transcriptome homeostasis in part through regulation of transcription factors like FOXO1. MiR-27a-3p arises as an oncogenic vulnerability in CCA.

IMPACT AND IMPLICATIONS

Cholangiocarcinogenesis entails extensive cellular reprogramming driven by genetic and non-genetic alterations, but the functional roles of these non-genetic events remain poorly understood. By unveiling global miRNA upregulation in patient tumours and their functional ability to increase proliferation of cholangiocytes, these small non-coding RNAs are implicated as critical non-genetic alterations promoting biliary tumour initiation. These findings identify possible mechanisms for transcriptome rewiring during transformation, with potential implications for patient stratification.

[The role of molecular genetic factors in the development of cholangiocellular carcinoma]

The article lists the main inducers of cholangiocarcinogenesis. The main inflammatory mediators (IL-6, nitric oxide, COX2) have been considered. Data on the study of gene mutations in cholangiocarcinomas are presented. The spectrum of genetic mutations depends on the biliary cancer origin (FGFR2 with intrahepatic cholangiocarcinoma, PRKACA, PRKACB with extrahepatic cholangiocarcinoma). Mutations in the KRAS, TP53, ARIAD1A genes are common in extrahepatic bile duct cancer. The role of epigenetic changes such as DNA hypermethylation, histone modifications, chromatin remodeling, as well as disturbances in miRNA expression is presented. A number of epigenetic features, such as the presence of a TP53 mutations with hypermethylation of p14ARF, DAPK, and/or ASC, correlate with a more aggressive course of the disease. The role of the SOX17 gene in the development of drug resistance is highlighted. The study of the molecular genetic features of extrahepatic bile duct cancer can help to better understand the pathogenesis of this type of tumor, to establish new prognostic and diagnostic markers of the disease.

Mutational signatures and processes in hepatobiliary cancers

The evolutionary history of hepatobiliary cancers is embedded in their genomes. By analysing their catalogue of somatic mutations and the DNA sequence context in which they occur, it is possible to infer the mechanisms underpinning tumorigenesis. These mutational signatures reflect the exogenous and endogenous origins of genetic damage as well as the capacity of hepatobiliary cells to repair and replicate DNA. Genomic analysis of thousands of patients with hepatobiliary cancers has highlighted the diversity of mutagenic processes active in these malignancies, highlighting a prominent source of the inter-cancer-type, inter-patient, intertumour and intratumoural heterogeneity that is observed clinically. However, a substantial proportion of mutational signatures detected in hepatocellular carcinoma and biliary tract cancer remain of unknown cause, emphasizing the important contribution of processes yet to be identified. Exploiting mutational signatures to retrospectively understand hepatobiliary carcinogenesis could advance preventative management of these aggressive tumours as well as potentially predict treatment response and guide the development of therapies targeting tumour evolution.

The Consumption of Cholesterol-Enriched Diets Conditions the Development of a Subtype of HCC with High Aggressiveness and Poor Prognosis

Non-alcoholic fatty liver disease (NAFLD) and progression to non-alcoholic steatohepatitis (NASH) result as a consequence of diverse conditions, mainly unbalanced diets. Particularly, high-fat and cholesterol content, as well as carbohydrates, such as those commonly ingested in Western countries, frequently drive adverse metabolic alterations in the liver and promote NAFLD development. Lipid liver overload is also one of the main risk factors for initiation and progression of hepatocellular carcinoma (HCC), but detailed knowledge on the relevance of high nutritional cholesterol remains elusive. We were aimed to characterize HCC development in mice fed with a Western diet (high in lipids and cholesterol) and to identify molecular alterations that define a subtype of liver cancer induced by lipid overload. Mice under western or high cholesterol diets more frequently developed tumors with a more aggressive phenotype than animals fed with a chow diet. Associated changes involved macrophage infiltration, angiogenesis, and stemness features. RNA-seq revealed a specific gene expression signature (Slc41a; Fabp5; Igdcc4 and Mthfd1l) resembling the adverse phenotypic features and poor clinical outcomes seen in patients with HCC. In conclusion; consumption of lipid enriched diets; particularly cholesterol; could accelerate HCC development with an aggressive phenotype and poor prognosis.

Patient-Derived Cancer Organoids as Predictors of Treatment Response

Patient-derived cancer organoids have taken a prominent role in pre-clinical and translational research and have been generated for most common solid tumors. Cancer organoids have been shown to retain key genetic and phenotypic characteristics of their tissue of origin, tumor subtype and maintain intratumoral heterogeneity and therefore have the potential to be used as predictors for individualized treatment response. In this review, we highlight studies that have used cancer organoids to compare the efficacy of standard-of-care and targeted combination treatments with clinical patient response. Furthermore, we review studies using cancer organoids to identify new anti-cancer treatments using drug screening. Finally, we discuss the current limitations and improvements needed to understand the full potential of cancer organoids as avatars for clinical management of cancer therapy.

Molecular Targets in Cholangiocarcinoma

Cholangiocarcinoma (CCA) encompasses a heterogeneous collection of malignancies for which diagnostic biomarkers are lacking and population screening is infeasible because of its status as a rare disease. Coupled with high postsurgical recurrence rates among the minority of patients diagnosed at resectable stages, systemic clinical management will inevitably be required for the majority of patients with CCA with recurrent and advanced disease. In this review, we discuss the therapeutic potential of different classes of molecular targets at various stages of development in CCA, including those targeted to the tumor epithelia (oncogenic, developmental, metabolic, epigenomic) and tumor microenvironment (angiogenesis, checkpoint regulation). Furthermore, we discuss the successes and failures of CCA-targeted therapies, emphasizing key lessons learned that should pave the way for future molecular target evaluation in this uncommon yet bona fide target-rich disease.

GDF11 restricts aberrant lipogenesis and changes in mitochondrial structure and function in human hepatocellular carcinoma cells

Growth differentiation factor 11 (GDF11) has been characterized as a key regulator of differentiation in cells that retain stemness features. Recently, it has been reported that GDF11 exerts tumor-suppressive properties in hepatocellular carcinoma cells, decreasing clonogenicity, proliferation, spheroid formation, and cellular function, all associated with a decrement in stemness features, resulting in mesenchymal to epithelial transition and loss of aggressiveness. The aim of the present work was to investigate the mechanism associated with the tumor-suppressive properties displayed by GDF11 in liver cancer cells. Hepatocellular carcinoma-derived cell lines were exposed to GDF11 (50 ng/ml), RNA-seq analysis in Huh7 cell line revealed that GDF11 exerted profound transcriptomic impact, which involved regulation of cholesterol metabolic process, steroid metabolic process as well as key signaling pathways, resembling endoplasmic reticulum-related functions. Cholesterol and triglycerides determination in Huh7 and Hep3B cells treated with GDF11 exhibited a significant decrement in the content of these lipids. The mTOR signaling pathway was downregulated, and this was associated with a reduction in key proteins involved in the mevalonate pathway. In addition, real-time metabolism assessed by Seahorse technology showed abridged glycolysis as well as glycolytic capacity, closely related to an impaired oxygen consumption rate and decrement in adenosine triphosphate production. Finally, transmission electron microscopy revealed mitochondrial abnormalities, such as cristae disarrangement, consistent with metabolic changes. Results provide evidence that GDF11 impairs cancer cell metabolism targeting lipid homeostasis, glycolysis, and mitochondria function and morphology.

From genome sequencing to the discovery of potential biomarkers in liver disease

Chronic liver disease progresses through several stages, fatty liver, steatohepatitis, cirrhosis, and eventually, it leads to hepatocellular carcinoma (HCC) over a long period of time. Since a large proportion of patients with HCC are accompanied by cirrhosis, it is considered to be an important factor in the diagnosis of liver cancer. This is because cirrhosis leads to an irreversible harmful effect, but the early stages of chronic liver disease could be reversed to a healthy state. Therefore, the discovery of biomarkers that could identify the early stages of chronic liver disease is important to prevent serious liver damage. Biomarker discovery at liver cancer and cirrhosis has enhanced the development of sequencing technology. Next generation sequencing (NGS) is one of the representative technical innovations in the biological field in the recent decades and it is the most important thing to design for research on what type of sequencing methods are suitable and how to handle the analysis steps for data integration. In this review, we comprehensively summarized NGS techniques for identifying genome, transcriptome, DNA methylome and 3D/4D chromatin structure, and introduced framework of processing data set and integrating multi-omics data for uncovering biomarkers.

BAX Redistribution Induces Apoptosis Resistance and Selective Stress Sensitivity in Human HCC

Cancer therapies induce differential cell responses, ranging from efficient cell death to complete stress resistance. The BCL-2 proteins BAX and BAK govern the cellular decision between survival and mitochondrial apoptosis. Therefore, the status of BAX/BAK regulation can predict the cellular apoptosis predisposition. Relative BAX/BAK localization was analyzed in tumor and corresponding non-tumor samples from 34 hepatocellular carcinoma (HCC) patients. Key transcriptome changes and gene expression profiles related to the status of BAX regulation were applied to two independent cohorts including over 500 HCC patients. The prediction of apoptotic response was tested using cell lines and polyclonal tumor isolates. Cellular protection from BAX was confirmed by challenging cells with mitochondrial BAX. We discovered a subgroup of HCC with selective protection from BAX-dependent apoptosis. BAX-protected tumors showed enrichment of signaling pathways associated with oxidative stress response and DNA repair as well as increased genetic heterogeneity. Gene expression profiles characteristic to BAX-specific protection are enriched in poorly differentiated HCCs and show significant association to the overall survival of HCC patients. Consistently, addiction to DNA repair of BAX-protected cancer cells caused selective sensitivity to PARP inhibition. Molecular characteristics of BAX-protected HCC were enriched in cells challenged with mitochondrial BAX. Our results demonstrate that predisposition to BAX activation impairs tumor biology in HCC. Selective BAX inhibition or lack thereof delineates distinct subgroups of HCC patients with molecular features and differential response pattern to apoptotic stimuli and inhibition of DNA repair mechanisms.

Severe metabolic alterations in liver cancer lead to ERK pathway activation and drug resistance

BACKGROUND

The extracellular signal-regulated kinase (ERK) pathway regulates cell growth, and is hyper-activated and associated with drug resistance in hepatocellular carcinoma (HCC). Metabolic pathways are profoundly dysregulated in HCC. Whether an altered metabolic state is linked to activated ERK pathway and drug response in HCC is unaddressed.

METHODS

We deprived HCC cells of glutamine to induce metabolic alterations and performed various assays, including metabolomics (with 13C-glucose isotope tracing), microarray analysis, and cell proliferation assays. Glutamine-deprived cells were also treated with kinase inhibitors (e.g. Sorafenib, Erlotinib, U0126 amongst other MEK inhibitors). We performed bioinformatics analysis and stratification of HCC tumour microarrays to determine upregulated ERK gene signatures in patients.

FINDINGS

In a subset of HCC cells, the withdrawal of glutamine triggers a severe metabolic alteration and ERK phosphorylation (pERK). This is accompanied by resistance to the anti-proliferative effect of kinase inhibitors, despite pERK inhibition. High intracellular serine is a consistent feature of an altered metabolic state and contributes to pERK induction and the kinase inhibitor resistance. Blocking the ERK pathway facilitates cell proliferation by reprogramming metabolism, notably enhancing aerobic glycolysis. We have identified 24 highly expressed ERK gene signatures that their combined expression strongly indicates a dysregulated metabolic gene network in human HCC tissues.

INTERPRETATION

A severely compromised metabolism lead to ERK pathway induction, and primes some HCC cells to pro-survival phenotypes upon ERK pathway blockade. Our findings offer novel insights for understanding, predicting and overcoming drug resistance in liver cancer patients. FUND: DFG, BMBF and Sino-German Cooperation Project.

Recreating Tumour Complexity in a Dish: Organoid Models to Study Liver Cancer Cells and their Extracellular Environment

Primary liver cancer, consisting predominantly of hepatocellular carcinoma (HCC) and cholangiocarcinoma (CCA), remains one of the most lethal malignancies worldwide. This high malignancy is related to the complex and dynamic interactions between tumour cells, stromal cells and the extracellular environment. Novel in vitro models that can recapitulate the tumour are essential in increasing our understanding of liver cancer. Herein, primary liver cancer-derived organoids have opened up new avenues due to their patient-specificity, self-organizing ability and potential recapitulation of many of the tumour properties. Organoids are solely of epithelial origin, but incorporation into co-culture models can enable the investigation of the cellular component of the tumour microenvironment. However, the extracellular component also plays a vital role in cancer progression and representation is lacking within current in vitro models. In this review, organoid technology is discussed in the context of liver cancer models through comparisons to other cell culture systems. In addition, the role of the tumour extracellular environment in primary liver cancer will be highlighted with an emphasis on its importance in in vitro modelling. Converging novel organoid-based models with models incorporating the native tumour microenvironment could lead to experimental models that can better recapitulate liver tumours in vivo.

Models for Understanding Resistance to Chemotherapy in Liver Cancer

The lack of response to pharmacological treatment constitutes a substantial limitation in the handling of patients with primary liver cancers (PLCs). The existence of active mechanisms of chemoresistance (MOCs) in hepatocellular carcinoma, cholangiocarcinoma, and hepatoblastoma hampers the usefulness of chemotherapy. A better understanding of MOCs is needed to develop strategies able to overcome drug refractoriness in PLCs. With this aim, several experimental models are commonly used. These include in vitro cell-free assays using subcellular systems; studies with primary cell cultures; cancer cell lines or heterologous expression systems; multicellular models, such as spheroids and organoids; and a variety of in vivo models in rodents, such as subcutaneous and orthotopic tumor xenografts or chemically or genetically induced liver carcinogenesis. Novel methods to perform programmed genomic edition and more efficient techniques to isolate circulating microvesicles offer new opportunities for establishing useful experimental tools for understanding the resistance to chemotherapy in PLCs. In the present review, using three criteria for information organization: (1) level of research; (2) type of MOC; and (3) type of PLC, we have summarized the advantages and limitations of the armamentarium available in the field of pharmacological investigation of PLC chemoresistance.

Dynamics and predicted drug response of a gene network linking dedifferentiation with beta-catenin dysfunction in hepatocellular carcinoma

BACKGROUND & AIMS

Alterations of individual genes variably affect the development of hepatocellular carcinoma (HCC). Thus, we aimed to characterize the function of tumor-promoting genes in the context of gene regulatory networks (GRNs).

METHODS

Using data from The Cancer Genome Atlas, from the LIRI-JP (Liver Cancer - RIKEN, JP project), and from our transcriptomic, transfection and mouse transgenic experiments, we identify a GRN which functionally links LIN28B-dependent dedifferentiation with dysfunction of β-catenin (CTNNB1). We further generated and validated a quantitative mathematical model of the GRN using human cell lines and in vivo expression data.

RESULTS

We found that LIN28B and CTNNB1 form a GRN with SMARCA4, Let-7b (MIRLET7B), SOX9, TP53 and MYC. GRN functionality is detected in HCC and gastrointestinal cancers, but not in other cancer types. GRN status negatively correlates with HCC prognosis, and positively correlates with hyperproliferation, dedifferentiation and HGF/MET pathway activation, suggesting that it contributes to a transcriptomic profile typical of the proliferative class of HCC. The mathematical model predicts how the expression of GRN components changes when the expression of another GRN member varies or is inhibited by a pharmacological drug. The dynamics of GRN component expression reveal distinct cell states that can switch reversibly in normal conditions, and irreversibly in HCC. The mathematical model is available via a web-based tool which can evaluate the GRN status of HCC samples and predict the impact of therapeutic agents on the GRN.

CONCLUSIONS

We conclude that identification and modelling of the GRN provide insights into the prognosis of HCC and the mechanisms by which tumor-promoting genes impact on HCC development.

LAY SUMMARY

Hepatocellular carcinoma (HCC) is a heterogeneous disease driven by the concomitant deregulation of several genes functionally organized as networks. Here, we identified a gene regulatory network involved in a subset of HCCs. This subset is characterized by increased proliferation and poor prognosis. We developed a mathematical model which uncovers the dynamics of the network and allows us to predict the impact of a therapeutic agent, not only on its specific target but on all the genes belonging to the network.

Epigenome Remodeling in Cholangiocarcinoma

Cholangiocarcinoma (CCA) comprises a heterogeneous collection of malignancies arising within the biliary tract, characterized by late diagnosis, innate chemoresistance, and abysmal prognosis. Sequencing data have uncovered recurrent mutations in diverse epigenetic regulators, implicating epigenetic destabilization at the root of these tumors. However, few studies have characterized biliary tumor epigenomes. In this Opinion article, we argue that an epigenome-oriented approach to CCA could establish diverse interconnections between many key aspects of research on this disease, including molecular heterogeneity, diverse cells of origin, and prominent tumor microenvironments. Moreover, we discuss plausible causes of epigenome dysregulation in biliary tumors, including genetic, epigenetic, metabolic, microenvironmental, and physiological factors. Lastly, we assess the translational potential of epigenomics in CCA to uncover robust biomarkers and therapeutic opportunities for this growing group of patients with limited treatment options.