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.

Lineage and ecology define liver tumor evolution in response to treatment.

A tumor ecosystem constantly evolves over time in the face of immune predation or therapeutic intervention, resulting in treatment failure and tumor progression. Here, we present a single-cell transcriptome-based strategy to determine the evolution of longitudinal tumor biopsies from liver cancer patients by measuring cellular lineage and ecology. We construct a lineage and ecological score as joint dynamics of tumor cells and their microenvironments. Tumors may be classified into four main states in the lineage-ecological space, which are associated with clinical outcomes. Analysis of longitudinal samples reveals the evolutionary trajectory of tumors in response to treatment. We validate the lineage-ecology-based scoring system in predicting clinical outcomes using bulk transcriptomic data of additional cohorts of 716 liver cancer patients. Our study provides a framework for monitoring tumor evolution in response to therapeutic intervention.

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.

Pan-viral serology uncovers distinct virome patterns as risk predictors of hepatocellular carcinoma and intrahepatic cholangiocarcinoma.

This study evaluates the pan-serological profiles of hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA) compared to several diseased and non-diseased control populations to identify risk factors and biomarkers of liver cancer. We used phage immunoprecipitation sequencing, an anti-viral antibody screening method using a synthetic-phage-displayed human virome epitope library, to screen patient serum samples for exposure to over 1,280 strains of pathogenic and non-pathogenic viruses. Using machine learning methods to develop an HCC or iCCA viral score, we discovered that both viral scores were positively associated with several liver function markers in two separate at-risk populations independent of viral hepatitis status. The HCC score predicted all-cause mortality over 8 years in patients with chronic liver disease at risk of HCC, while the viral hepatitis status was not predictive of survival. These results suggest that non-hepatitis viral infections may contribute to HCC and iCCA development and could be biomarkers in at-risk populations.

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.

Tumor biology and immune infiltration define primary liver cancer subsets linked to overall survival after immunotherapy.

Primary liver cancer is a rising cause of cancer deaths in the US. Although immunotherapy with immune checkpoint inhibitors induces a potent response in a subset of patients, response rates vary among individuals. Predicting which patients will respond to immune checkpoint inhibitors is of great interest in the field. In a retrospective arm of the National Cancer Institute Cancers of the Liver: Accelerating Research of Immunotherapy by a Transdisciplinary Network (NCI-CLARITY) study, we use archived formalin-fixed, paraffin-embedded samples to profile the transcriptome and genomic alterations among 86 hepatocellular carcinoma and cholangiocarcinoma patients prior to and following immune checkpoint inhibitor treatment. Using supervised and unsupervised approaches, we identify stable molecular subtypes linked to overall survival and distinguished by two axes of aggressive tumor biology and microenvironmental features. Moreover, molecular responses to immune checkpoint inhibitor treatment differ between subtypes. Thus, patients with heterogeneous liver cancer may be stratified by molecular status indicative of treatment response to immune checkpoint inhibitors.

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.

A compendium of co-regulated mitoribosomal proteins in pan-cancer uncovers collateral defective events in tumor malignancy.

Mitochondria are major organelles responsible for cellular energy and metabolism, and their dysfunction is tightly linked to cancer. The mitochondrial ribosome (mitoribosome) is a protein complex consisting of 82 mitoribosomal proteins (MRPs) encoded by nuclear genes and is essential for mitochondrial protein synthesis. However, their roles in tumorigenesis remain poorly understood. We performed pan-cancer analyses of 18,177 tumors representing 28 cancer types to determine somatic alterations of MRP genes as a genetic basis for tumorigenesis. We identified a set of 20 altered MRPs known to be involved in early assembly of the mitoribosome complex. We found that tumors with affected MRPs were associated with impaired mitochondrial functions and TP53 mutations accompanied by increased genomic instability and intra-tumor heterogeneity. MRP deletions were associated with poor survival. Our results reveal a key role for mitochondrial ribosome biogenesis in tumor malignancy across cancer types.

Single-cell biology uncovers apoptotic cell death and its spatial organization as a potential modifier of tumor diversity in HCC.

BACKGROUND AND AIMS: HCC is a highly aggressive and heterogeneous cancer type with limited treatment options. Identifying drivers of tumor heterogeneity may lead to better therapeutic options and favorable patient outcomes. We investigated whether apoptotic cell death and its spatial architecture is linked to tumor molecular heterogeneity using single-cell in situ hybridization analysis. APPROACH AND RESULTS: We analyzed 254 tumor samples from two HCC cohorts using tissue microarrays. We developed a mathematical model to quantify cellular diversity among HCC samples using two tumor markers, cyclin-dependent kinase inhibitor 3 and protein regulator of cytokinesis 1 as surrogates for heterogeneity and caspase 3 (CASP3) as an apoptotic cell death marker. We further explored the impact of potential dying-cell hubs on tumor cell diversity and patient outcome by density contour mapping and spatial proximity analysis. We also developed a selectively controlled in vitro model of cell death using CRISPR/CRISPR-associated 9 to determine therapy response and growth under hypoxic conditions. We found that increasing levels of CASP3+ tumor cells are associated with higher tumor diversity. Interestingly, we discovered regions of densely populated CASP3+ , which we refer to as CASP3+ cell islands, in which the nearby cellular heterogeneity was found to be the greatest compared to cells farther away from these islands and that this phenomenon was associated with survival. Additionally, cell culture experiments revealed that higher levels of cell death, accompanied by increased CASP3 expression, led to greater therapy resistance and growth under hypoxia. CONCLUSIONS: These results are consistent with the hypothesis that increased apoptotic cell death may lead to greater tumor heterogeneity and thus worse patient outcomes.

Genetic susceptibility to hepatocellular carcinoma in chromosome 22q13.31, findings of a genome-wide association study.

BACKGROUND AND AIM: Chronic hepatitis C virus (HCV) infection, long-term alcohol use, cigarette smoking, and obesity are the major risk factors for hepatocellular carcinoma (HCC) in the United States, but the disease risk varies substantially among individuals with these factors, suggesting host susceptibility to and gene-environment interactions in HCC. To address genetic susceptibility to HCC, we conducted a genome-wide association study (GWAS). METHODS: Two case-control studies on HCC were conducted in the United States. DNA samples were genotyped using the Illumian microarray chip with over 710 000 single nucleotide polymorphisms (SNPs). We compared these SNPs between 705 HCC cases and 1455 population controls for their associations with HCC and verified our findings in additional studies. RESULTS: In this GWAS, we found that two SNPs were associated with HCC at P < 5E-8 and six SNPs at P < 5E-6 after adjusting for age, sex, and the top three principal components (PCs). Five of the SNPs in chromosome 22q13.31, three in PNPLA3 (rs2281135, rs2896019, and rs4823173) and two in SAMM50 (rs3761472, rs3827385), were replicated in a small US case-control study and a cohort study in Singapore. The associations remained significant after adjusting for body mass index and HCV infection. Meta-analysis of multiple datasets indicated that these SNPs were significantly associated with HCC. CONCLUSIONS: SNPs in PNPLA3 and SAMM50 are known risk loci for nonalcoholic fatty liver disease (NAFLD) and are suspected to be associated with HCC. Our GWAS demonstrated the associations of these SNPs with HCC in a US population. Biological mechanisms underlying the relationship remain to be elucidated.

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.

Single-cell atlas of tumor cell evolution in response to therapy in hepatocellular carcinoma and intrahepatic cholangiocarcinoma.

BACKGROUND & AIMS: Intratumor molecular heterogeneity is a key feature of tumorigenesis and is linked to treatment failure and patient prognosis. Herein, we aimed to determine what drives tumor cell evolution by performing single-cell transcriptomic analysis. METHODS: We analyzed 46 hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA) biopsies from 37 patients enrolled in interventional studies at the NIH Clinical Center, with 16 biopsies collected before and after treatment from 7 patients. We developed a novel machine learning-based consensus clustering approach to track cellular states of 57,000 malignant and non-malignant cells including tumor cell transcriptome-based functional clonality analysis. We determined tumor cell relationships using RNA velocity and reverse graph embedding. We also studied longitudinal samples from 4 patients to determine tumor cellular state and its evolution. We validated our findings in bulk transcriptomic data from 488 patients with HCC and 277 patients with iCCA. RESULTS: Using transcriptomic clusters as a surrogate for functional clonality, we observed an increase in tumor cell state heterogeneity which was tightly linked to patient prognosis. Furthermore, increased functional clonality was accompanied by a polarized immune cell landscape which included an increase in pre-exhausted T cells. We found that SPP1 expression was tightly associated with tumor cell evolution and microenvironmental reprogramming. Finally, we developed a user-friendly online interface as a knowledge base for a single-cell atlas of liver cancer. CONCLUSIONS: Our study offers insight into the collective behavior of tumor cell communities in liver cancer as well as potential drivers of tumor evolution in response to therapy. LAY SUMMARY: Intratumor molecular heterogeneity is a key feature of tumorigenesis that is linked to treatment failure and patient prognosis. In this study, we present a single-cell atlas of liver tumors from patients treated with immunotherapy and describe intratumoral cell states and their hierarchical relationship. We suggest osteopontin, encoded by the gene SPP1, as a candidate regulator of tumor evolution in response to treatment.

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.

Intratumoral γδ T-Cell Infiltrates, Chemokine (C-C Motif) Ligand 4/Chemokine (C-C Motif) Ligand 5 Protein Expression and Survival in Patients With Hepatocellular Carcinoma.

BACKGROUND AND AIMS: Hepatocellular carcinoma (HCC) is an aggressive malignancy which is often associated with a complex tumor microenvironment attributable to etiology-induced cellular inflammation. γδ T cells are known to detect and react to chronic inflammation, which is linked to cancer development, progression, and metastasis. Our recent genomic study revealed an increased infiltration of several immune cell types, including γδ T cells, in tumor microenvironments of a Thai HCC subtype associated with a good prognosis. APPROACH AND RESULTS: Here, we quantified the amount of γδ T cells using a γδ T-cell-specific gene signature in 247 Chinese HCC patients. We also validated the γδ T-cell signature in American HCC patients. Additionally, such an association was only found in tumor transcriptomic data, but not in adjacent nontumor transcriptomic data, suggesting a selective enrichment of γδ T cells in the tumor microenvironment. Moreover, the γδ T-cell signature was positively correlated with the expression of natural killer cell receptor genes, such as NKG2D and cytolytic T-cell genes granzymes and perforin, suggesting a stronger T-cell-mediated cytotoxic activity. Furthermore, we found that the γδ T-cell-specific gene expression is positively correlated with the expression of chemokine (C-C motif) ligand 4 (CCL4)/chemokine (C-C motif) ligand 5 (CCL5) and C-C chemokine receptor type 1 (CCR1)/C-C chemokine receptor type 5 (CCR5), the receptors for γδ T cells. We validated these results using immunohistochemical analysis of formalin-fixed, paraffin-embedded tumor biopsies from 182 HCC patients. Moreover, we found evidence of CCL4/CCL5-mediated recruitment of γδ T cells both in vitro and in a murine orthotopic Hepa1-6 HCC model. CONCLUSIONS: We propose that CCL4/CCL5 may interact with their receptor, CCR1/CCR5, which may facilitate the recruitment of γδ T cells from peripheral blood or peritumor regions to the tumor regions. Consequently, an increasing infiltration of γδ T cells in tumors may enhance antitumor immunity and improve patients' prognosis.

The genomic landscape of Mongolian hepatocellular carcinoma.

Mongolia has the highest incidence of hepatocellular carcinoma (HCC) in the world, but its causative factors and underlying tumor biology remain unknown. Here, we describe molecular characteristics of HCC from 76 Mongolian patients by whole-exome and transcriptome sequencing. We present a comprehensive analysis of mutational signatures, driver genes, and molecular subtypes of Mongolian HCC compared to 373 HCC patients of different races and ethnicities and diverse etiologies. Mongolian HCC consists of prognostic molecular subtypes similar to those found in patients from other areas of Asia, Europe, and North America, as well as other unique subtypes, suggesting the presence of distinct etiologies linked to Mongolian patients. In addition to common driver mutations (TP53, CTNNB1) frequently found in pan-cancer analysis, Mongolian HCC exhibits unique drivers (most notably GTF2IRD2B, PNRC2, and SPTA1), the latter of which is associated with hepatitis D viral infection. These results suggest the existence of new molecular mechanisms at play in Mongolian hepatocarcinogenesis.

A Viral Exposure Signature Defines Early Onset of Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is an aggressive malignancy with its global incidence and mortality rate continuing to rise, although early detection and surveillance are suboptimal. We performed serological profiling of the viral infection history in 899 individuals from an NCI-UMD case-control study using a synthetic human virome, VirScan. We developed a viral exposure signature and validated the results in a longitudinal cohort with 173 at-risk patients who had long-term follow-up for HCC development. Our viral exposure signature significantly associated with HCC status among at-risk individuals in the validation cohort (area under the curve: 0.91 [95% CI 0.87-0.96] at baseline and 0.98 [95% CI 0.97-1] at diagnosis). The signature identified cancer patients prior to a clinical diagnosis and was superior to alpha-fetoprotein. In summary, we established a viral exposure signature that can predict HCC among at-risk patients prior to a clinical diagnosis, which may be useful in HCC surveillance.

AGO2 Mediates MYC mRNA Stability in Hepatocellular Carcinoma.

Deregulated RNA-binding proteins (RBP), such as Argonaute 2 (AGO2), mediate tumor-promoting transcriptomic changes during carcinogenesis, including hepatocellular carcinoma (HCC). While AGO2 is well characterized as a member of the RNA-induced silencing complex (RISC), which represses gene expression through miRNAs, its role as a bona fide RBP remains unclear. In this study, we investigated AGO2's role as an RBP that regulates the MYC transcript to promote HCC. Using mRNA and miRNA arrays from patients with HCC, we demonstrate that HCCs with elevated AGO2 levels are more likely to have the mRNA transcriptome deregulated and are associated with poor survival. Moreover, AGO2 overexpression stabilizes the MYC transcript independent of miRNAs. These observations provide a novel mechanism of gene regulation by AGO2 and provide further insights into the potential functions of AGO2 as an RBP in addition to RISC. IMPLICATIONS: Authors demonstrate that the RBP Argonaute 2 stabilizes the MYC transcript to promote HCC.

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.

Urinary Metabolites Diagnostic and Prognostic of Intrahepatic Cholangiocarcinoma.

BACKGROUND: Liver cancer is the second leading cause of cancer-related deaths worldwide. With a predicted 2.4-fold rise in liver cancer incidence by 2020, there is an urgent need for early, inexpensive diagnostic biomarkers to deploy in the clinic. METHODS: We employed ultraperformance liquid chromatography tandem mass-spectrometry (UPLC/MS-MS) for the quantitation of four metabolites, creatine riboside (CR), N-acetylneuraminic acid (NANA), cortisol sulfate, and a lipid molecule designated as 561+, in urine samples from the NCI-MD cohort comprising 98 hepatocellular carcinoma (HCC) cases, 101 high-risk subjects, and 95 controls. Validation was carried out in the TIGER-LC cohort [n = 370 HCC and intrahepatic cholangiocarcinoma (ICC) cases, 471 high-risk subjects, 251 controls], where ICC, the second most common primary hepatic malignancy, is highly prevalent. Metabolite quantitation was also conducted in TIGER-LC tissue samples (n = 48 ICC; n = 51 HCC). RESULTS: All profiled metabolites were significantly increased in liver cancer when compared with high-risk subjects and controls in the NCI-MD study. In the TIGER-LC cohort, the four-metabolite profile was superior at classifying ICC than a clinically utilized marker, CA19-9, and their combination led to a significantly improved model (AUC = 0.88, P = 4E-8). Metabolites CR and NANA were significantly elevated in ICC when compared with HCC cases in both urine and tissue samples. High levels of CR were associated with poorer prognosis in ICC. CONCLUSIONS: Four metabolites are significantly increased in HCC and ICC and are robust at classifying ICC in combination with the clinically utilized marker CA19-9. IMPACT: Noninvasive urinary metabolite biomarkers hold promise for diagnostic and prognostic evaluation of ICC.