PLK1 inhibition leads to mitotic arrest and triggers apoptosis in cholangiocarcinoma cells.

Cholangiocarcinoma (CCA) is a lethal cancer originating from the epithelial cells within the bile duct and ranks as the second most prevalent form of liver cancer in Thailand. Polo-like kinase 1 (PLK1), a protein serine/threonine kinase, regulates a number of steps in cell mitosis and is upregulated in several types of cancer, including CCA. Our previous study identified PLK1 as a biomarker of the C1 subtype, correlating with poor prognosis in intrahepatic CCA. The present study aimed to examine the effect of PLK1 inhibition on CCA cells. Different CCA cell lines developed from Thai patients, HuCCA1, KKU055, KKU100 and KKU213A, were treated with two PLK1 inhibitors, BI2536 and BI6727, and were transfected with small interfering RNA, followed by analysis of cell proliferation, cell cycle distribution and cell apoptosis. It was discovered that BI2536 and BI6727 inhibited cell proliferation and caused G2/M-phase arrest in CCA cells. Furthermore, the number of total apoptotic cells was increased in PLK1 inhibitor-treated CCA cells. The expression levels of mitotic proteins, aurora kinase A, phosphorylated PLK1 (T210) and cyclin B1, were augmented in PLK1-inhibited CCA cells. Additionally, inhibition of PLK1 led to increased DNA damage, as determined by the upregulated levels of γH2AX and increased cleavage of poly (ADP-ribose) polymerase, an apoptotic marker. These results suggested that inhibiting PLK1 prolonged mitotic arrest and subsequently triggered cell apoptosis. Validation of the antiproliferative effects of PLK1 inhibition was accomplished through silencing of the PLK1 gene. In conclusion, targeting PLK1 provided promising results for further study as a potential candidate for targeted therapy in CCA.

Characterization of tumor evolution by functional clonality and phylogenetics in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is a molecularly heterogeneous solid malignancy, and its fitness may be shaped by how its tumor cells evolve. However, ability to monitor tumor cell evolution is hampered by the presence of numerous passenger mutations that do not provide any biological consequences. Here we develop a strategy to determine the tumor clonality of three independent HCC cohorts of 524 patients with diverse etiologies and race/ethnicity by utilizing somatic mutations in cancer driver genes. We identify two main types of tumor evolution, i.e., linear, and non-linear models where non-linear type could be further divided into classes, which we call shallow branching and deep branching. We find that linear evolving HCC is less aggressive than other types. GTF2IRD2B mutations are enriched in HCC with linear evolution, while TP53 mutations are the most frequent genetic alterations in HCC with non-linear models. Furthermore, we observe significant B cell enrichment in linear trees compared to non-linear trees suggesting the need for further research to uncover potential variations in immune cell types within genomically determined phylogeny types. These results hint at the possibility that tumor cells and their microenvironment may collectively influence the tumor evolution process.

Spatial proximity of tumor-immune interactions predicts patient outcome in hepatocellular carcinoma.

BACKGROUND AND AIMS: The fitness and viability of a tumor ecosystem are influenced by the spatial organization of its cells. We aimed to study the structure, architecture, and cell-cell dynamics of the heterogeneous liver cancer tumor microenvironment using spatially resolved multiplexed imaging. APPROACH AND RESULTS: We performed co-detection by indexing multiplexed immunofluorescence imaging on 68 HCC biopsies from Thai patients [(Thailand Initiative in Genomics and Expression Research for Liver Cancer (TIGER-LC)] as a discovery cohort, and then validated the results in an additional 190 HCC biopsies from Chinese patients [Liver Cancer Institute (LCI)]. We segmented and annotated 117,270 and 465,632 cells from the TIGER-LC and LCI cohorts, respectively. We observed 4 patient groups of TIGER-LC (IC1, IC2, IC3, and IC4) with distinct tumor-immune cellular interaction patterns. In addition, patients from IC2 and IC4 had much better overall survival than those from IC1 and IC3. Noticeably, tumor and CD8 + T-cell interactions were strongly enriched in IC2, the group with the best patient outcomes. The close proximity between the tumor and CD8 + T cells was a strong predictor of patient outcome in both the TIGER-LC and the LCI cohorts. Bulk transcriptomic data from 51 of the 68 HCC cases were used to determine tumor-specific gene expression features of our classified subtypes. Moreover, we observed that the presence of immune spatial neighborhoods in HCC as a measure of overall immune infiltration is linked to better patient prognosis. CONCLUSIONS: Highly multiplexed imaging analysis of liver cancer reveals tumor-immune cellular heterogeneity within spatial contexts, such as tumor and CD8 + T-cell interactions, which may predict patient survival.

Genome-wide profiling of transcription factor activity in primary liver cancer using single-cell ATAC sequencing.

Primary liver cancer (PLC) consists of two main histological subtypes; hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA). The role of transcription factors (TFs) in malignant hepatobiliary lineage commitment between HCC and iCCA remains underexplored. Here, we present genome-wide profiling of transcription regulatory elements of 16 PLC patients using single-cell assay for transposase accessible chromatin sequencing. Single-cell open chromatin profiles reflect the compositional diversity of liver cancer, identifying both malignant and microenvironment component cells. TF motif enrichment levels of 31 TFs strongly discriminate HCC from iCCA tumors. These TFs are members of the nuclear/retinoid receptor, POU, or ETS motif families. POU factors are associated with prognostic features in iCCA. Overall, nuclear receptors, ETS and POU TF motif families delineate transcription regulation between HCC and iCCA tumors, which may be relevant to development and selection of PLC subtype-specific therapeutics.

Gut dysbiosis in Thai intrahepatic cholangiocarcinoma and hepatocellular carcinoma.

Primary liver cancer (PLC), which includes intrahepatic cholangiocarcinoma (iCCA) and hepatocellular carcinoma (HCC), has the highest incidence of all cancer types in Thailand. Known etiological factors, such as viral hepatitis and chronic liver disease do not fully account for the country's unusually high incidence. However, the gut-liver axis, which contributes to carcinogenesis and disease progression, is influenced by the gut microbiome. To investigate this relationship, fecal matter from 44 Thai PLC patients and 76 healthy controls were subjected to whole-genome metagenomic shotgun sequencing and then analyzed by marker gene-based and assembly based methods. Results revealed greater gut microbiome heterogeneity in iCCA compared to HCC and healthy controls. Two Veillonella species were found to be more abundant in iCCA samples and could distinguish iCCA from HCC and healthy controls. Conversely, Ruminococcus gnavus was depleted in iCCA patients and could distinguish HCC from iCCA samples. High Veillonella genus counts in the iCCA group were associated with enriched amino acid biosynthesis and glycolysis pathways, while enriched phospholipid and thiamine metabolism pathways characterized the HCC group with high Blautia genus counts. These findings reveal distinct landscapes of gut dysbiosis among Thai iCCA and HCC patients and warrant further investigation as potential biomarkers.

Purine anabolism creates therapeutic vulnerability in hepatocellular carcinoma through m 6 A-mediated epitranscriptomic regulation.

BACKGROUND AND AIMS: Purines are building blocks for the cellular genome, and excessive purine nucleotides are seen in tumors. However, how purine metabolism is dysregulated in tumors, and impacting tumorigenesis remains elusive. APPROACH AND RESULTS: Transcriptomic and metabolomic analyses of purine biosynthesis and purine degradation pathways were performed in the tumor and associated nontumor liver tissues obtained from 62 patients with HCC, one of the most lethal cancers worldwide. We found that most genes in purine synthesis are upregulated, while genes in purine degradation are inhibited in HCC tumors. High purine anabolism is associated with unique somatic mutational signatures linked to patient prognosis. Mechanistically, we discover that increasing purine anabolism promotes epitranscriptomic dysregulation of DNA damage repairing (DDR) machinery through upregulating RNA N6-methyladenosine (m 6 A) modification. High purine anabolic HCC is sensitive to DDR-targeting agents but not to standard HCC treatments, correlating with the clinical outcomes in 5 independent HCC cohorts containing 724 patients. We further showed that high purine anabolism determines the sensitivity to DDR-targeting agents in 5 HCC cell lines in vitro and in vivo . CONCLUSIONS: Our results reveal a central role of purine anabolism in regulating DDR, which could be therapeutically exploited in HCC.

Multiregional single-cell dissection of tumor and immune cells reveals stable lock-and-key features in liver cancer.

Intratumor heterogeneity may result from the evolution of tumor cells and their continuous interactions with the tumor microenvironment which collectively drives tumorigenesis. However, an appearance of cellular and molecular heterogeneity creates a challenge to define molecular features linked to tumor malignancy. Here we perform multiregional single-cell RNA sequencing (scRNA-seq) analysis of seven liver cancer patients (four hepatocellular carcinoma, HCC and three intrahepatic cholangiocarcinoma, iCCA). We identify cellular dynamics of malignant cells and their communication networks with tumor-associated immune cells, which are validated using additional scRNA-seq data of 25 HCC and 12 iCCA patients as a stable fingerprint embedded in a malignant ecosystem representing features of tumor aggressiveness. We further validate the top ligand-receptor interaction pairs (i.e., LGALS9-SLC1A5 and SPP1-PTGER4 between tumor cells and macrophages) associated with unique transcriptome in additional 542 HCC patients. Our study unveils stable molecular networks of malignant ecosystems, which may open a path for therapeutic exploration.

Integration of adeno-associated virus (AAV) into the genomes of most Thai and Mongolian liver cancer patients does not induce oncogenesis.

BACKGROUND: Engineered versions of adeno-associated virus (AAV) are commonly used in gene therapy but evidence revealing a potential oncogenic role of natural AAV in hepatocellular carcinoma (HCC) has raised concerns. The frequency of potentially oncogenic integrations has been reported in only a few populations. AAV infection and host genome integration in another type of liver cancer, cholangiocarcinoma (CCA), has been studied only in one cohort. All reported oncogenic AAV integrations in HCC come from strains resembling the fully sequenced AAV2 and partly sequenced AAV13. When AAV integration occurs, only a fragment of the AAV genome is detectable in later DNA or RNA sequencing. The integrated fragment is typically from the 3' end of the AAV genome, and this positional bias has been only partly explained. Three research groups searched for evidence of AAV integration in HCC RNAseq samples in the Cancer Genome Atlas (TCGA) but reported conflicting results. RESULTS: We collected and analyzed whole transcriptome and viral capture DNA sequencing in paired tumor and non-tumor samples from two liver cancer Asian cohorts from Thailand (N = 147, 47 HCC and 100 intrahepatic cholangiocarcinoma (iCCA)) and Mongolia (N = 70, all HCC). We found only one HCC patient with a potentially oncogenic integration of AAV, in contrast to higher frequency reported in European patients. There were no oncogenic AAV integrations in iCCA patients. AAV genomic segments are present preferentially in the non-tumor samples of Thai patients. By analyzing the AAV genome positions of oncogenic and non-oncogenic integrated fragments, we found that almost all the putative oncogenic integrations overlap the X gene, which is present and functional only in the strain AAV2 among all fully sequenced strains. This gene content difference could explain why putative oncogenic integrations from other AAV strains have not been reported. We resolved the discrepancies in previous analyses of AAV presence in TCGA HCC samples and extended it to CCA. There are 12 TCGA samples with an AAV segment and none are in Asian patients. AAV segments are present in preferentially in TCGA non-tumor samples, like what we observed in the Thai patients. CONCLUSIONS: Our findings suggest a minimal AAV risk of hepatocarcinogenesis in Asian liver cancer patients. The partial genome presence and positional bias of AAV integrations into the human genome has complicated analysis of possible roles of AAV in liver cancer.

Tumor metabolism and associated serum metabolites define prognostic subtypes of Asian hepatocellular carcinoma.

Treatment effectiveness in hepatocellular carcinoma (HCC) depends on early detection and precision-medicine-based patient stratification for targeted therapies. However, the lack of robust biomarkers, particularly a non-invasive diagnostic tool, precludes significant improvement of clinical outcomes for HCC patients. Serum metabolites are one of the best non-invasive means for determining patient prognosis, as they are stable end-products of biochemical processes in human body. In this study, we aimed to identify prognostic serum metabolites in HCC. To determine serum metabolites that were relevant and representative of the tissue status, we performed a two-step correlation analysis to first determine associations between metabolic genes and tissue metabolites, and second, between tissue metabolites and serum metabolites among 49 HCC patients, which were then validated in 408 additional Asian HCC patients with mixed etiologies. We found that certain metabolic genes, tissue metabolites and serum metabolites can independently stratify HCC patients into prognostic subgroups, which are consistent across these different data types and our previous findings. The metabolic subtypes are associated with ß-oxidation process in fatty acid metabolism, where patients with worse survival outcome have dysregulated fatty acid metabolism. These serum metabolites may be used as non-invasive biomarkers to define prognostic tumor molecular subtypes for HCC.

Tumor methionine metabolism drives T-cell exhaustion in hepatocellular carcinoma.

T-cell exhaustion denotes a hypofunctional state of T lymphocytes commonly found in cancer, but how tumor cells drive T-cell exhaustion remains elusive. Here, we find T-cell exhaustion linked to overall survival in 675 hepatocellular carcinoma (HCC) patients with diverse ethnicities and etiologies. Integrative omics analyses uncover oncogenic reprograming of HCC methionine recycling with elevated 5-methylthioadenosine (MTA) and S-adenosylmethionine (SAM) to be tightly linked to T-cell exhaustion. SAM and MTA induce T-cell dysfunction in vitro. Moreover, CRISPR-Cas9-mediated deletion of MAT2A, a key SAM producing enzyme, results in an inhibition of T-cell dysfunction and HCC growth in mice. Thus, reprogramming of tumor methionine metabolism may be a viable therapeutic strategy to improve HCC immunity.

Functional Genomic Complexity Defines Intratumor Heterogeneity and Tumor Aggressiveness in Liver Cancer.

Chronic inflammation and chromosome aneuploidy are major traits of primary liver cancer (PLC), which represent the second most common cause of cancer-related death worldwide. Increased cancer fitness and aggressiveness of PLC may be achieved by enhancing tumoral genomic complexity that alters tumor biology. Here, we developed a scoring method, namely functional genomic complexity (FGC), to determine the degree of molecular heterogeneity among 580 liver tumors with diverse ethnicities and etiologies by assessing integrated genomic and transcriptomic data. We found that tumors with higher FGC scores are associated with chromosome instability and TP53 mutations, and a worse prognosis, while tumors with lower FGC scores have elevated infiltrating lymphocytes and a better prognosis. These results indicate that FGC scores may serve as a surrogate to define genomic heterogeneity of PLC linked to chromosomal instability and evasion of immune surveillance. Our findings demonstrate an ability to define genomic heterogeneity and corresponding tumor biology of liver cancer based only on bulk genomic and transcriptomic data. Our data also provide a rationale for applying this approach to survey liver tumor immunity and to stratify patients for immune-based therapy.

Gene signature predictive of hepatocellular carcinoma patient response to transarterial chemoembolization.

Transarterial chemoembolization (TACE) is a commonly used treatment modality in hepatocellular carcinoma (HCC). The ability to identify patients who will respond to TACE represents an important clinical need, and tumor gene expression patterns may be associated with TACE response. We investigated whether tumor transcriptome is associated with TACE response in patients with HCC. We analyzed transcriptome data of treatment-naïve tumor tissues from a Chinese cohort of 191 HCC patients, including 105 patients who underwent TACE following resection with curative intent. We then developed a gene signature, TACE Navigator, which was associated with improved survival in patients that received either adjuvant or post-relapse TACE. To validate our findings, we applied our signature in a blinded manner to three independent cohorts comprising an additional 130 patients with diverse ethnic backgrounds enrolled in three different hospitals who received either adjuvant TACE or palliative TACE. TACE Navigator stratified patients into Responders and Non-Responders which was associated with improved survival following TACE in our test cohort (Responders: 67 months vs Non-Responders: 39.5 months, p<0.0001). In addition, multivariable Cox model demonstrates that TACE Navigator was independently associated with survival (HR: 9.31, 95% CI: 3.46-25.0, p<0.001). In our validation cohorts, the association between TACE Navigator and survival remained robust in both Asian patients who received adjuvant TACE (Hong Kong: 60 months vs 25.6 months p=0.007; Shandong: 61.3 months vs 32.1 months, p=0.027) and European patients who received TACE as primary therapy (Mainz: 60 months vs 41.5 months, p=0.041). These results indicate that a TACE-specific molecular classifier is robust in predicting TACE response. This gene signature can be used to identify patients who will have the greatest survival benefit after TACE treatment and enable personalized treatment modalities for patients with HCC.

DNA methylation at an enhancer of the three prime repair exonuclease 2 gene (TREX2) is linked to gene expression and survival in laryngeal cancer.

BACKGROUND: Genetic aberrations in DNA repair genes are linked to cancer, but less is reported about epigenetic regulation of DNA repair and functional consequences. We investigated the intragenic methylation loss at the three prime repair exonuclease 2 (TREX2) locus in laryngeal (n = 256) and colorectal cancer cases (n = 95) and in pan-cancer data from The Cancer Genome Atlas (TCGA). RESULTS: Significant methylation loss at an intragenic site of TREX2 was a frequent trait in both patient cohorts (p = 0.016 and < 0.001, respectively) and in 15 out of 22 TCGA studies. Methylation loss correlated with immunohistochemically staining for TREX2 (p < 0.0001) in laryngeal tumors and improved overall survival of laryngeal cancer patients (p = 0.045). Chromatin immunoprecipitation, demethylation experiments, and reporter gene assays revealed that the region of methylation loss can function as a CCAAT/enhancer binding protein alpha (CEBPA)-responsive enhancer element regulating TREX2 expression. CONCLUSIONS: The data highlight a regulatory role of TREX2 DNA methylation for gene expression which might affect incidence and survival of laryngeal cancer. Altered TREX2 protein levels in tumors may affect drug-induced DNA damage repair and provide new tailored therapies.

Common Molecular Subtypes Among Asian Hepatocellular Carcinoma and Cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (ICC) and hepatocellular carcinoma (HCC) are clinically disparate primary liver cancers with etiological and biological heterogeneity. We identified common molecular subtypes linked to similar prognosis among 199 Thai ICC and HCC patients through systems integration of genomics, transcriptomics, and metabolomics. While ICC and HCC share recurrently mutated genes, including TP53, ARID1A, and ARID2, mitotic checkpoint anomalies distinguish the C1 subtype with key drivers PLK1 and ECT2, whereas the C2 subtype is linked to obesity, T cell infiltration, and bile acid metabolism. These molecular subtypes are found in 582 Asian, but less so in 265 Caucasian patients. Thus, Asian ICC and HCC, while clinically treated as separate entities, share common molecular subtypes with similar actionable drivers to improve precision therapy.

Reduced promoter methylation and increased expression of CSPG4 negatively influences survival of HNSCC patients.

Proteoglycans are often overexpressed in tumors and can be found on several normal and neoplastic stem cells. In this study, we analyzed in-depth the role of CSPG4 in head and neck squamous cell carcinomas (HNSCC). Analysis of CSPG4 in a homogeneous study sample of HPV-negative stage IVa HNSCCs revealed overexpression of protein and mRNA levels in a subgroup of HNSCC tumors and a significant association of high CSPG4 protein levels with poor survival. This could be validated in three publicly available microarray datasets. As a potential cause for upregulated CSPG4 expression, we identified DNA hypomethylation in a CpG-island of the promoter region. Accordingly, we found an inverse correlation of methylation and patient outcome. Finally, CSPG4 re-expression was achieved by demethylating treatment of highly methylated HNSCC cell lines establishing a direct link between methylation and CSPG4 expression. In conclusion, we identified CSPG4 as a novel biomarker in HNSCC on several biological levels and established a causative link between DNA methylation and CSPG4 protein and mRNA expression.

Epigenetically mediated downregulation of the differentiation-promoting chaperon protein CRABP2 in astrocytic gliomas.

Impairment of endogenous differentiation pathways like retinoic acid (RA) signaling seems to be a central pathogenetic event in astrocytic gliomas. Among others, expression of the differentiation-promoting RA chaperon protein cellular retinoic acid binding protein 2 (CRABP2) is extenuated in high-grade gliomas. Against this background, we aimed at identifying potential pathomechanisms underlying reduced CRABP2 expression in these tumors. Using MassARRAY methylation analysis, we detected extensive CpG methylation upstream of the CRABP2 gene locus in a study sample comprising 100 astrocytic gliomas of WHO Grade II to IV. Compared to nontumorous control samples, tumors revealed increased CpG methylation and methylation levels were inversely correlated to CRABP2 mRNA expression. Substantiating our in situ findings, CRABP2 mRNA levels increased in glioma cell lines after exposure to the demethylating agent 5-aza-2'-deoxycytidine. Finally, a distinct CpG methylation signature distinguished between primary glioblastoma on the one hand and the group of astrocytoma WHO II-III and secondary glioblastoma on the other hand. Altogether, our observations suggest that epigenetic silencing of CRABP2 might contribute to an immature phenotype in glioma cells.