Mining the interpretable prognostic features from pathological image of intrahepatic cholangiocarcinoma using multi-modal deep learning.

BACKGROUND: The advances in deep learning-based pathological image analysis have invoked tremendous insights into cancer prognostication. Still, lack of interpretability remains a significant barrier to clinical application. METHODS: We established an integrative prognostic neural network for intrahepatic cholangiocarcinoma (iCCA), towards a comprehensive evaluation of both architectural and fine-grained information from whole-slide images. Then, leveraging on multi-modal data, we conducted extensive interrogative approaches to the models, to extract and visualize the morphological features that most correlated with clinical outcome and underlying molecular alterations. RESULTS: The models were developed and optimized on 373 iCCA patients from our center and demonstrated consistent accuracy and robustness on both internal (n = 213) and external (n = 168) cohorts. The occlusion sensitivity map revealed that the distribution of tertiary lymphoid structures, the geometric traits of the invasive margin, the relative composition of tumor parenchyma and stroma, the extent of necrosis, the presence of the disseminated foci, and the tumor-adjacent micro-vessels were the determining architectural features that impacted on prognosis. Quantifiable morphological vector extracted by CellProfiler demonstrated that tumor nuclei from high-risk patients exhibited significant larger size, more distorted shape, with less prominent nuclear envelope and textural contrast. The multi-omics data (n = 187) further revealed key molecular alterations left morphological imprints that could be attended by the network, including glycolysis, hypoxia, apical junction, mTORC1 signaling, and immune infiltration. CONCLUSIONS: We proposed an interpretable deep-learning framework to gain insights into the biological behavior of iCCA. Most of the significant morphological prognosticators perceived by the network are comprehensible to human minds.

Identification and validation of a lactate metabolism-related six-gene prognostic signature in intrahepatic cholangiocarcinoma.

PURPOSE: Intrahepatic cholangiocarcinoma (iCCA) is a highly malignant and fatal liver tumor with increasing incidence worldwide. Lactate metabolism has been recently reported as a crucial contributor to tumor progression and immune regulation in the tumor microenvironment. However, it remains poorly identified about the biological functions of lactate metabolism in iCCA, which hinders the development of prognostic tools and therapeutic interventions. METHODS: The univariate Cox regression analysis and Boruta algorithm were utilized to identify key lactate metabolism-related genes (LMRGs), and a prognostic signature was constructed based on LMRG scores. Genomic variations and immune cell infiltration were evaluated in the high and low LMRG score groups. Finally, the biological functions of key LMRGs were verified with in vitro and in vivo experiments. RESULTS: Patients in the high LMRG score group exhibit a poor prognosis compared to those in the low LMRG score group, with a high frequency of TP53 and KRAS mutations. Moreover, the infiltration and function of NK cells were compromised in the high LMRG score group, consistent with the results from two independent single-cell RNA sequencing datasets and immunohistochemistry of tissue microarrays. Experimental data revealed that lactate dehydrogenase A (LDHA) knockdown inhibited proliferation and migration in iCCA cell lines and tumor growth in immunocompetent mice. CONCLUSION: Our study revealed the biological roles of LDHA in iCCA and developed a reliable lactate metabolism-related prognostic signature for iCCA, offering promising therapeutic targets for iCCA in the clinic.

Whole-Genome DNA Methylation Profiling of Intrahepatic Cholangiocarcinoma Reveals Prognostic Subtypes with Distinct Biological Drivers.

Intrahepatic cholangiocarcinoma (iCCA) is the second most prevalent primary liver cancer. Although the genetic characterization of iCCA has led to targeted therapies for treating tumors with FGFR2 alterations and IDH1/2 mutations, only a limited number of patients can benefit from these strategies. Epigenomic profiles have emerged as potential diagnostic and prognostic biomarkers for improving the treatment of cancers. In this study, we conducted whole-genome bisulfite sequencing on 331 iCCAs integrated with genetic, transcriptomic, and proteomic analyses, demonstrating the existence of four DNA methylation subtypes of iCCAs (S1-S4) that exhibited unique postoperative clinical outcomes. The S1 group was an IDH1/2 mutation-specific subtype with moderate survival. The S2 subtype was characterized by the lowest methylation level and the highest mutational burden among the four subtypes and displayed upregulation of a gene-expression pattern associated with cell cycle/DNA replication. The S3 group was distinguished by high interpatient heterogeneity of tumor immunity, a gene-expression pattern associated with carbohydrate metabolism, and an enrichment of KRAS alterations. Patients with the S2 and S3 subtypes had the shortest survival among the four subtypes. Tumors in the S4 subtype, which had the best prognosis, showed global methylation levels comparable to normal controls, increased FGFR2 fusions/BAP1 mutations, and the highest copy-number variant burdens. Further integrative and functional analyses identified GBP4 demethylation, which is highly prevalent in the S2 and S3 groups, as an epigenetic oncogenic factor that regulates iCCA proliferation, migration, and invasion. Together, this study identifies prognostic methylome alterations and epigenetic drivers in iCCA. SIGNIFICANCE: Characterization of the DNA methylome of intrahepatic cholangiocarcinoma integrated with genomic, transcriptomic, and proteomic analyses uncovers molecular mechanisms affected by genome-wide DNA methylation alterations, providing a resource for identifying potential therapeutic targets.

Proteomics of adjacent-to-tumor samples uncovers clinically relevant biological events in hepatocellular carcinoma.

Normal adjacent tissues (NATs) of hepatocellular carcinoma (HCC) differ from healthy liver tissues and their heterogeneity may contain biological information associated with disease occurrence and clinical outcome that has yet to be fully evaluated at the proteomic level. This study provides a detailed description of the heterogeneity of NATs and the differences between NATs and healthy livers and revealed that molecular features of tumor subgroups in HCC were partially reflected in their respective NATs. Proteomic data classified HCC NATs into two subtypes (Subtypes 1 and 2), and Subtype 2 was associated with poor prognosis and high-risk recurrence. The pathway and immune features of these two subtypes were characterized. Proteomic differences between the two NAT subtypes and healthy liver tissues were further investigated using data-independent acquisition mass spectrometry, revealing the early molecular alterations associated with the progression from healthy livers to NATs. This study provides a high-quality resource for HCC researchers and clinicians and may significantly expand the knowledge of tumor NATs to eventually benefit clinical practice.

An Inflammatory Checkpoint Generated by IL1RN Splicing Offers Therapeutic Opportunity for KRAS-Mutant Intrahepatic Cholangiocarcinoma.

KRAS mutations are causally linked to protumor inflammation and are identified as driving factors in tumorigenesis. Here, using multiomics data gathered from a large set of patients, we showed that KRAS mutation was associated with a specific landscape of alternative mRNA splicing that connected to myeloid inflammation in intrahepatic cholangiocarcinoma (iCCA). Then, we identified a negative feedback mechanism in which the upregulation of interleukin 1 receptor antagonist (IL1RN)-201/203 due to alternative splicing confers vital anti-inflammatory effects in KRAS-mutant iCCA. In KRAS-mutant iCCA mice, both IL1RN-201/203 upregulation and anakinra treatment ignited a significant antitumor immune response by altering neutrophil recruitment and phenotypes. Furthermore, anakinra treatment synergistically enhanced anti-PD-1 therapy to activate intratumoral GZMB+ CD8+ T cells in KRAS-mutant iCCA mice. Clinically, we found that high IL1RN-201/203 levels in patients with KRAS-mutant iCCA were significantly associated with superior response to anti-PD-1 immunotherapy. SIGNIFICANCE: This work describes a novel inflammatory checkpoint mediated by IL1RN alternative splicing variants that may serve as a promising basis to develop therapeutic options for KRAS-mutant iCCA and other cancers. This article is featured in Selected Articles from This Issue, p. 2109.

Pharmaco-proteogenomic characterization of liver cancer organoids for precision oncology.

Organoid models have the potential to recapitulate the biological and pharmacotypic features of parental tumors. Nevertheless, integrative pharmaco-proteogenomics analysis for drug response features and biomarker investigation for precision therapy of patients with liver cancer are still lacking. We established a patient-derived liver cancer organoid biobank (LICOB) that comprehensively represents the histological and molecular characteristics of various liver cancer types as determined by multiomics profiling, including genomic, epigenomic, transcriptomic, and proteomic analysis. Proteogenomic profiling of LICOB identified proliferative and metabolic organoid subtypes linked to patient prognosis. High-throughput drug screening revealed distinct response patterns of each subtype that were associated with specific multiomics signatures. Through integrative analyses of LICOB pharmaco-proteogenomics data, we identified the molecular features associated with drug responses and predicted potential drug combinations for personalized patient treatment. The synergistic inhibition effect of mTOR inhibitor temsirolimus and the multitargeted tyrosine kinase inhibitor lenvatinib was validated in organoids and patient-derived xenografts models. We also provide a user-friendly web portal to help serve the biomedical research community. Our study is a rich resource for investigation of liver cancer biology and pharmacological dependencies and may help enable functional precision medicine.

SSH1 promotes progression of intrahepatic cholangiocarcinoma via p38 MAPK-CXCL8 axis.

Protein tyrosine phosphatases (PTPs) are involved in malignant transformation and metastasis. According to one of our previous studies, Slingshot homolog 1 (SSH1), a member of PTPs, is significantly associated with the survival of intrahepatic cholangiocarcinoma (iCCA) patients. However, the underlying mechanisms of SSH1 in iCCA remain largely elusive. Here, the expression and clinical significance of SSH1 were assessed using the iCCA patient samples. The results showed that SSH1 was dramatically up-regulated in iCCA tissues and elevated SSH1 expression was associated with worse overall survival of iCCA patients. Overexpression of SSH1 accelerated the proliferation, migration, and invasion of iCCA cells, and also inhibited cell apoptosis. Furthermore, the downstream signaling pathway of SSH1 in iCCA was explored and it was revealed that the increased expression of SSH1 could activate the p38 mitogen-activated protein kinase (MAPK) pathway and enhance the expression of C-X-C motif chemokine ligand 8 (CXCL8). Notably, the high correlation of SSH1 with CXCL8 jointly indicated the poor prognosis in iCCA patients. Thus, our study suggests SSH1 as a potentially promising target for iCCA, which promoted iCCA progression through a potential p38 MAPK-CXCL8 axis.

Gene-based cancer-testis antigens as prognostic indicators in hepatocellular carcinoma.

Cancer/testis antigens (CTAs) are reproductive tissue-restricted genes, frequently ectopic expressed in tumors. CTA genes associate with a poor prognosis in some solid tumors, due to their potential roles in the tumorigenesis and progression. However, whether CTAs relate with hepatocellular carcinoma (HCC) remains unclear. In this study, the prognostic signatures based on CTA genes were investigated and validated in three cohorts including Chinese HCC patients with hepatitis B virus infection (CHCC-HBV), International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) cohorts. Univariate, LASSO, and multivariate Cox regression analyses were used to screen prognostic genes and develop the prognostic gene signature. A prognosis model was established with six CTA genes (SSX1, CTCFL, OIP5, CEP55, NOL4, and TPPP2) in CHCC-HBV cohort, and further validated in the ICGC and TCGA cohorts. The CTA signature was an essential prognostic predictor independent of other clinical pathological factors. High-risk group exhibited up-regulated cell cycle-related and tumor-related pathways and more M0 macrophage, activated mast cell, activated memory CD4+ T cell, and memory B cell infiltration. Furthermore, CTA signature correlated with the sensitivity to multiple chemotherapy drugs. Our results highlighted that the CTA gene profiling was a prognostic assessment tool for HCC patients.

Bispecific antibodies targeting immunomodulatory checkpoints for cancer therapy.

Advances in antibody engineering have led to the generation of more innovative antibody drugs, such as bispecific antibodies (bsAbs). Following the success associated with blinatumomab, bsAbs have attracted enormous interest in the field of cancer immunotherapy. By specifically targeting two different antigens, bsAbs reduce the distance between tumor and immune cells, thereby enhancing tumor killing directly. There are several mechanisms of action upon which bsAbs have been exploited. Accumulating experience on checkpoint-based therapy has promoted the clinical transformation of bsAbs targeting immunomodulatory checkpoints. Cadonilimab (PD-1 × CTLA-4) is the first approved bsAb targeting dual inhibitory checkpoints, which confirms the feasibility of bsAbs in immunotherapy. In this review we analyzed the mechanisms by which bsAbs targeting immunomodulatory checkpoints and their emerging applications in cancer immunotherapy.

Genome-Scale CRISPR screen identifies LAPTM5 driving lenvatinib resistance in hepatocellular carcinoma.

ABBREVIATIONS: cld-CASP3: cleaved caspase 3; cld-PARP: cleaved PARP; DTP: drug tolerant persister; GO: Gene Ontology; GTEx: The Genotype-Tissue Expression; HCC: hepatocellular carcinoma; HCQ: hydroxychloroquine; IC50: half maximal inhibitory concentration value; KEGG: Kyoto Encyclopedia of Genes and Genomes; LAPTM5: lysosomal protein transmembrane 5; NT: non-targeting; PDC: patient-derived primary cell lines; PDO: patient-derived primary organoid; TCGA: The Cancer Genome Atlas.

The influence of 18F-fluorodeoxyglucose positron emission tomography/computed tomography on the N- and M-staging and subsequent clinical management of intrahepatic cholangiocarcinoma.

Background: Intrahepatic cholangiocarcinoma (ICC) is a highly metastatic cancer. 18F-fluorodeoxyglucose positron emission tomography/computed tomography (18F-FDG PET/CT) enables sensitive tumor and metastasis detection. Our aim is to evaluate the influence of pre-treatment PET/CT on the N- and M-staging and subsequent clinical management in ICC patients. Methods: Between August 2010 and August 2018, 660 consecutive ICC patients, without prior anti-tumor treatments nor other malignancies, were enrolled. The diagnostic performance of PET/CT on the N- and M-staging was compared with conventional imaging, and the preoperative staging accuracy and treatment re-allocation by PET/CT were retrospectively calculated. Survival difference was compared between patients receiving PET/CT or not after propensity score matching. Results: Patients were divided into group A (n=291) and group B (n=369) according to whether PET/CT was performed. Among 291 patients with both PET/CT and conventional imaging for staging in group A, PET/CT showed significantly higher sensitivity (83.0% vs. 70.5%, P=0.001), specificity (88.3% vs. 74.9%, P<0.001) and accuracy (86.3% vs. 73.2%, P<0.001) than conventional imaging in diagnosing regional lymph node metastasis, as well as higher sensitivity (87.8% vs. 67.6%, P<0.001) and accuracy (93.5% vs. 89.3%, P=0.023) in diagnosing distant metastasis. Overall, PET/CT improved the accuracy of preoperative staging from 60.1% to 71.8% (P<0.001), and modified clinical treatment strategy in 5.8% (17/291) of ICC patients, with unique roles in different tumor-node-metastasis (TNM) stages. High tumor-to-non-tumor ratio (TNR) predicted poor overall survival [hazard ratio (HR) = 2.17; 95% confidence interval (CI): 1.49-3.15; P<0.001]. Furthermore, patients performing PET/CT had longer overall survival compared with those without PET/CT (HR =0.74; 95% CI: 0.58-0.93; P=0.011) after propensity score matching. Conclusions: PET/CT was valuable for diagnosing regional lymph node metastasis and distant metastasis in ICC patients, and facilitated accurate tumor staging and optimal treatment allocation.

Multi-omics characterization reveals the pathogenesis of liver focal nodular hyperplasia.

The molecular landscape and pathogenesis of focal nodular hyperplasia (FNH) have yet to be elucidated. We performed multi-omics approaches on FNH and paired normal liver tissues from 22 patients, followed by multi-level bioinformatic analyses and experimental validations. Generally, FNH had low mutation burden with low variant allele frequencies, and the mutation frequency significantly correlated with proliferation rate. Although no recurrently deleterious genomic events were found, some putative tumor suppressors or oncogenes were involved. Mutational signatures indicated potential impaired mismatch function and possible poison contact. Integrated analyses unveiled a group of FNH specific endothelial cells that uniquely expressed SOST and probably had strong interaction with fibroblasts through PDGFB/PDGFRB pathway to promote fibrosis. Notably, in one atypical FNH (patient No.11) with pronounced copy number variations, we observed a unique immune module. Most FNH are benign, but molecularly atypical FNH still exist; endothelial cell derived PDGFB probably promotes the fibrogenic process in FNH.

Geospatial Immune Heterogeneity Reflects the Diverse Tumor-Immune Interactions in Intrahepatic Cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (iCCA) exhibits extensive intratumoral heterogeneity and an extremely high mortality rate. Here, we performed whole-exome sequencing, RNA sequencing, T-cell receptor (TCR) sequencing, and multiplexed immunofluorescence on 207 tumor regions from 45 patients with iCCA. Over half of iCCA displayed intratumoral heterogeneity of immune infiltration, and iCCA were classified into sparsely, heterogeneously, and highly infiltrated subgroups with distinct immunogenomic characteristics. Sparsely infiltrated tumors displayed active copy-number loss of clonal neoantigens, and heterogeneous immune infiltration played an important role in the subclonal evolution across tumor subregions. Highly infiltrated tumors were characterized by extensive immune activation and a similar TCR repertoire across tumor subregions, but counteracted with T-cell exhaustion and pervasive antigen presentation defects. Notably, FGFR2 mutations and fusions correlated with low mutation burden and reduced immune infiltration. Our work delineated the dynamic tumor-immune interactions and developed a robust classification system to divide patients with iCCA into high and low immune evasion groups with different prognoses. SIGNIFICANCE: This study elucidates the impact of spatial immune heterogeneity upon tumor evolution of iCCA and reveals distinct immune evasion mechanisms developed in different immune microenvironments, which can be exploited for the development of personalized immunotherapy strategies. This article is highlighted in the In This Issue feature, p. 2221.

Single-cell transcriptomic analysis suggests two molecularly subtypes of intrahepatic cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (iCCA) is a highly heterogeneous cancer with limited understanding of its classification and tumor microenvironment. Here, by performing single-cell RNA sequencing on 144,878 cells from 14 pairs of iCCA tumors and non-tumor liver tissues, we find that S100P and SPP1 are two markers for iCCA perihilar large duct type (iCCAphl) and peripheral small duct type (iCCApps). S100P + SPP1- iCCAphl has significantly reduced levels of infiltrating CD4+ T cells, CD56+ NK cells, and increased CCL18+ macrophages and PD1+CD8+ T cells compared to S100P-SPP1 + iCCApps. The transcription factor CREB3L1 is identified to regulate the S100P expression and promote tumor cell invasion. S100P-SPP1 + iCCApps has significantly more SPP1+ macrophage infiltration, less aggressiveness and better survival than S100P + SPP1- iCCAphl. Moreover, S100P-SPP1 + iCCApps harbors tumor cells at different status of differentiation, such as ALB + hepatocyte differentiation and ID3+ stemness. Our study extends the understanding of the diversity of tumor cells in iCCA.

Proteogenomic characterization identifies clinically relevant subgroups of intrahepatic cholangiocarcinoma.

We performed proteogenomic characterization of intrahepatic cholangiocarcinoma (iCCA) using paired tumor and adjacent liver tissues from 262 patients. Integrated proteogenomic analyses prioritized genetic aberrations and revealed hallmarks of iCCA pathogenesis. Aflatoxin signature was associated with tumor initiation, proliferation, and immune suppression. Mutation-associated signaling profiles revealed that TP53 and KRAS co-mutations may contribute to iCCA metastasis via the integrin-FAK-SRC pathway. FGFR2 fusions activated the Rho GTPase pathway and could be a potential source of neoantigens. Proteomic profiling identified four patient subgroups (S1-S4) with subgroup-specific biomarkers. These proteomic subgroups had distinct features in prognosis, genetic alterations, microenvironment dysregulation, tumor microbiota composition, and potential therapeutics. SLC16A3 and HKDC1 were further identified as potential prognostic biomarkers associated with metabolic reprogramming of iCCA cells. This study provides a valuable resource for researchers and clinicians to further identify molecular pathogenesis and therapeutic opportunities in iCCA.

A pre-screening strategy to assess resected tumor margins by imaging cytoplasmic viscosity and hypoxia.

To assure complete tumor removal, frozen section analysis is the most common procedure for intraoperative pathological assessment of resected tumor margins. However, during one operation, multiple biopsies may be sent for examination, but only few of them are made into cryosections because of the complex preparation protocols and time-consuming pathological analysis, which potentially increases the risk of overlooking tumor involvement. Here, we propose a fluorescence-based pre-screening strategy that allows high-throughput, convenient, and fast gross assessment of resected tumor margins. A dual-activatable cationic fluorescent molecular rotor was developed to specifically illuminate live tumor cells' cytoplasm by emitting two different fluorescence signals in response to elevations in hypoxia-induced nitroreductase (a biochemical marker) and cytoplasmic viscosity (a biophysical marker), two characteristics of cancer cells. The ability of the fluorescent molecular rotor in detecting tumor cells was evaluated in mouse and human specimens of multiple tissues by comparing with hematoxylin and eosin staining. Importantly, the fluorescent molecular rotor achieved 100 % specificity in discriminating lung and liver cancers from normal tissue, allowing pre-screening of the tumor-free surgical margins and promoting clinical decision. Altogether, this type of fluorescent molecular rotor and the proposed strategy may serve as a new option to facilitate intraoperative assessment of resected tumor margins.

Global immune characterization of HBV/HCV-related hepatocellular carcinoma identifies macrophage and T-cell subsets associated with disease progression.

Diverse immune cells in the tumor microenvironment form a complex ecosystem, but our knowledge of their heterogeneity and dynamics within hepatocellular carcinoma (HCC) still remains limited. To assess the plasticity and phenotypes of immune cells within HBV/HCV-related HCC microenvironment at single-cell level, we performed single-cell RNA sequencing on 41,698 immune cells from seven pairs of HBV/HCV-related HCC tumors and non-tumor liver tissues. We combined bio-informatic analyses, flow cytometry, and multiplex immunohistochemistry to assess the heterogeneity of different immune cell subsets in functional characteristics, transcriptional regulation, phenotypic switching, and interactions. We identified 29 immune cell subsets of myeloid cells, NK cells, and lymphocytes with unique transcriptomic profiles in HCC. A highly complex immunological network was shaped by diverse immune cell subsets that can transit among different states and mutually interact. Notably, we identified a subset of M2 macrophage with high expression of CCL18 and transcription factor CREM that was enriched in advanced HCC patients, and potentially participated in tumor progression. We also detected a new subset of activated CD8+ T cells highly expressing XCL1 that correlated with better patient survival rates. Meanwhile, distinct transcriptomic signatures, cytotoxic phenotypes, and evolution trajectory of effector CD8+ T cells from early-stage to advanced HCC were also identified. Our study provides insight into the immune microenvironment in HBV/HCV-related HCC and highlights novel macrophage and T-cell subsets that could be further exploited in future immunotherapy.

Systemic Therapy for Hepatocellular Carcinoma: Advances and Hopes.

The majority of patients with hepatocellular carcinoma (HCC) are diagnosed at an advanced stage that can only benefit from systemic treatments. Although HCC is highly treatmentresistant, significant achievements have been made in the molecular targeted therapy and immunotherapy of HCC. In addition to regorafenib, cabozantinib and ramucirumab were approved for the second- line targeted treatment by the FDA after disease progression on sorafenib. Nivolumab failed to demonstrate remarkable benefit in overall survival (OS) as first-line therapy, while pembrolizumab did not achieve pre-specified statistical significance in both OS and progression-free survival (PFS) as second-line treatment. Combinations of targeted agents, immune checkpoint inhibitors and other interventions showed favorable results. In this review, we summarized the progress of systemic therapy in HCC and discussed the future directions of the treatment of HCC.