Comprehensive and Integrative Genomic Characterization of Hepatocellular Carcinoma

Immune landscape of hepatocellular carcinoma microenvironment: Implications for prognosis and therapeutic applications

The development of immunotherapy for solid tumours has boosted interest in the contexture of tumour immune microenvironment (TIME). Several lines of evidence indicate that the interplay between tumour cells and TIME components is a key factor for the evolution of hepatocellular carcinoma (HCC) and for the likelihood of response to immunotherapeutics. The availability of high-resolution methods will be instrumental for a better definition of the complexity and diversity of TIME with the aim of predicting disease outcome, treatment response and possibly new therapeutic targets. Here, we review current knowledge about the immunological mechanisms involved in shaping the clinical course of HCC. Effector cells, regulatory cells and soluble mediators are discussed for their role defining TIME and as targets for immune modulation, together with possible immune signatures for optimization of immunotherapeutic strategies.

Development and validation of a CIMP-associated prognostic model for hepatocellular carcinoma

BACKGROUND

CpG island methylator phenotype (CIMP), a common biological phenomenon characterized by a subset of concurrently methylated genes, can have an influence on the progression of multiple cancers. However, the potential mechanism of CIMP in hepatocarcinogenesis and its clinical relevance remains only partially understood.

METHODS

We used a methylation array from the cancer genome atlas (TCGA) to stratify HCC patients into different CIMP subtypes, and evaluated their correlation with clinical characteristics. In addition, mutation, CNV, and transcriptome profiles were also utilized to evaluate the distinctive genomic patterns correlated with CIMP. Finally, a CIMP-associated prognostic model (CPM) was trained and validated using four independent datasets.

FINDINGS

A subgroup of patients was identified as having CIMP-H, which was associated with worse OS and DFS. Gene enrichment analysis indicated that the terms "liver cancer with EPCAM up", "tumor invasiveness up", "methyltransferase complex", and "translational initiation" were enriched in CIMP-H subgroup. Notably, somatic mutation analysis indicated that CIMP-H patients presented with a higher mutation burden of BRD4, DDIAS and NOX1. Moreover, four CPM associated genes could significantly categorize patients into low- and high-risk groups in the training dataset and another 3 independent validation datasets. Finally, a nomogram incorporating a classifier based on four mRNAs, pathological M stage and CIMP status was established, which showed a favorable discriminating ability and might contribute to clinical decision-making for HCC.

INTERPRETATION

Our work highlights the potential clinical application value of CPM in predicting the overall survival of HCC patients and the mechanisms underlying the role of CIMP in hepatocarcinogenesis. FUND: This work was supported by the State Key Project on Infectious Diseases of China (2018ZX10723204-003), the National Nature Science Foundation of China (Nos. 81874065, 81500565, 81874149, 81572427, and 81401997), the Hepato-Biliary-Pancreatic Malignant Tumor Investigation Fund of Chen Xiao-ping Foundation for the Development of Science and Technology of Hubei Province (CXPJJH11800001-2018356).

Analysis of Liver Cancer Cell Lines Identifies Agents With Likely Efficacy Against Hepatocellular Carcinoma and Markers of Response

BACKGROUND AND AIMS

Hepatocellular carcinomas (HCCs) are heterogeneous aggressive tumors with low rates of response to treatment at advanced stages. We screened a large panel of liver cancer cell lines (LCCLs) to identify agents that might be effective against HCC and markers of therapeutic response.

METHODS

We performed whole-exome RNA and microRNA sequencing and quantification of 126 proteins in 34 LCCLs. We screened 31 anticancer agents for their ability to decrease cell viability. We compared genetic, RNA, and protein profiles of LCCLs with those of primary HCC samples and searched for markers of response.

RESULTS

The protein, RNA and mutational signatures of the LCCLs were similar to those of the proliferation class of HCC, which is the most aggressive tumor type. Cell lines with alterations in genes encoding members of the Ras-MAPK signaling pathway and that required fibroblast growth factor (FGF)19 signaling via FGF receptor 4 for survival were more sensitive to trametinib than to FGF receptor 4 inhibitors. Amplification of FGF19 resulted in increased activity of FGF19 only in tumor cells that kept a gene expression pattern of hepatocyte differentiation. We identified single agents and combinations of agents that reduced viability of cells with features of the progenitor subclass of HCC. LCCLs with inactivating mutations in TSC1 and TSC2 were sensitive to the mammalian target of rapamycin inhibitor rapamycin, and cells with inactivating mutations in TP53 were sensitive to the Aurora kinase A inhibitor alisertib. Amplification of MET was associated with hypersensitivity to cabozantinib and the combination of sorafenib and inhibitors of MAP kinase 1 and MAP kinase2 had a synergistic antiproliferative effect.

CONCLUSION

LCCLs can be screened for drugs and agents that might be effective for treatment of HCC. We identified genetic alterations and gene expression patterns associated with response to these agents. This information might be used to select patients for clinical trials.

Growth arrest-specific gene 2 suppresses hepatocarcinogenesis by intervention of cell cycle and p53-dependent apoptosis

BACKGROUND

Growth arrest-specific gene 2 (GAS2) plays a role in modulating in reversible growth arrest cell cycle, apoptosis, and cell survival. GAS2 protein is universally expressed in most normal tissues, particularly in the liver, but is depleted in some tumor tissues. However, the functional mechanisms of GAS2 in hepatocellular carcinoma (HCC) are not fully defined.

AIM

To investigate the function and mechanism of GAS2 in HCC.

METHODS

GAS2 expression in clinic liver and HCC specimens was analyzed by real-time PCR and western blotting. Cell proliferation was analyzed by counting, MTS, and colony formation assays. Cell cycle analysis was performed by flow cytometry. Cell apoptosis was investigated by Annexin V apoptosis assay and western blotting.

RESULTS

GAS2 protein expression was lower in HCC than in normal tissues. Overexpression of GAS2 inhibited the proliferation of HCC cells with wide-type p53, while knockdown of GAS2 promoted the proliferation of hepatocytes (P < 0.05). Furthermore, GAS2 overexpression impeded the G1-to-S cell cycle transition and arrested more G1 cells, particularly the elevation of sub G1 (P < 0.01). Apoptosis induced by GAS2 was dependent on p53, which was increased by etoposide addition. The expression of p53 and apoptosis markers was further enhanced when GAS2 was upregulated, but became diminished upon downregulation of GAS2. In the clinic specimen, GAS2 was downregulated in more than 60% of HCCs. The average fold changes of GAS2 expression in tumor tissues were significantly lower than those in paired non-tumor tissues (P < 0.05).

CONCLUSION

GAS2 plays a vital role in HCC cell proliferation and apoptosis, possibly by regulating the cell cycle and p53-dependent apoptosis pathway.

Hippo signaling and epithelial cell plasticity in mammalian liver development, homeostasis, injury and disease

A traditional view of cellular differentiation is unidirectional: progenitor cells adopt specific fates in response to environmental cues resulting in deployment of cell-specific gene expression programs and acquisition of unique differentiated cellular properties such as production of structural and functional proteins that define individual cell types. In both development and in tissue repair stem and progenitor cells are thought to both self-renew to maintain the pool of precursors and to expand to give rise to transient amplifying and differentiated cell types. Recently, however, it has become appreciated that differentiated cell types can be reprogrammed to adopt progenitor and stem cell properties. In the case of epithelial cells in the mammalian liver, hepatocytes and biliary epithelial cells there is a significant degree of plasticity between these lineages that has been implicated in mechanisms of tissue repair and in liver pathologies such as cancer. Recent studies have highlighted the role of Hippo signaling, an emerging growth control and tumor suppressor pathway, in regulating epithelial cell plasticity in the mammalian liver and in this review, the role of cellular plasticity and Hippo signaling in regulating normal and abnormal tissue responses in the mammalian liver will be discussed.

Tumor-Infiltrating Leukocyte Composition and Prognostic Power in Hepatitis B- and Hepatitis C-Related Hepatocellular Carcinomas

Background: Tumor-infiltrating leukocytes (TILs) are immune cells surrounding tumor cells, and several studies have shown that TILs are potential survival predictors in different cancers. However, few studies have dissected the differences between hepatitis B- and hepatitis C-related hepatocellular carcinoma (HBV−HCC and HCV−HCC). Therefore, we aimed to determine whether the abundance and composition of TILs are potential predictors for survival outcomes in HCC and which TILs are the most significant predictors. Methods: Two bioinformatics algorithms, ESTIMATE and CIBERSORT, were utilized to analyze the gene expression profiles from 6 datasets, from which the abundance of corresponding TILs was inferred. The ESTIMATE algorithm examined the overall abundance of TILs, whereas the CIBERSORT algorithm reported the relative abundance of 22 different TILs. Both HBV−HCC and HCV−HCC were analyzed. Results: The results indicated that the total abundance of TILs was higher in non-tumor tissue regardless of the HCC type. Alternatively, the specific TILs associated with overall survival (OS) and recurrence-free survival (RFS) varied between subtypes. For example, in HBV−HCC, plasma cells (hazard ratio [HR] = 1.05; 95% CI 1.00–1.10; p = 0.034) and activated dendritic cells (HR = 1.08; 95% CI 1.01–1.17; p = 0.03) were significantly associated with OS, whereas in HCV−HCC, monocytes (HR = 1.21) were significantly associated with OS. Furthermore, for RFS, CD8+ T cells (HR = 0.98) and M0 macrophages (HR = 1.02) were potential biomarkers in HBV−HCC, whereas neutrophils (HR = 1.01) were an independent predictor in HCV−HCC. Lastly, in both HBV−HCC and HCV−HCC, CD8+ T cells (HR = 0.97) and activated dendritic cells (HR = 1.09) had a significant association with OS, while γ delta T cells (HR = 1.04), monocytes (HR = 1.05), M0 macrophages (HR = 1.04), M1 macrophages (HR = 1.02), and activated dendritic cells (HR = 1.15) were highly associated with RFS. Conclusions: These findings demonstrated that TILs are potential survival predictors in HCC and different kinds of TILs are observed according to the virus type. Therefore, further investigations are warranted to elucidate the role of TILs in HCC, which may improve immunotherapy outcomes.

Degenerate minigene library analysis enables identification of altered branch point utilization by mutant splicing factor 3B1 (SF3B1)

Cancer-associated mutations of the core splicing factor 3 B1 (SF3B1) result in selection of novel 3′ splice sites (3′SS), but precise molecular mechanisms of oncogenesis remain unclear. SF3B1 stabilizes the interaction between U2 snRNP and branch point (BP) on the pre-mRNA. It has hence been speculated that a change in BP selection is the basis for novel 3′SS selection. Direct quantitative determination of BP utilization is however technically challenging. To define BP utilization by SF3B1-mutant spliceosomes, we used an overexpression approach in human cells as well as a complementary strategy using isogenic murine embryonic stem cells with monoallelic K700E mutations constructed via CRISPR/Cas9-based genome editing and a dual vector homology-directed repair methodology. A synthetic minigene library with degenerate regions in 3′ intronic regions (3.4 million individual minigenes) was used to compare BP usage of SF3B1^K700E and SF3B1^WT. Using this model, we show that SF3B1^K700E spliceosomes utilize non-canonical sequence variants (at position −1 relative to BP adenosine) more frequently than wild-type spliceosomes. These predictions were confirmed using minigene splicing assays. Our results suggest a model of BP utilization by mutant SF3B1 wherein it is able to utilize non-consensus alternative BP sequences by stabilizing weaker U2-BP interactions.

Mechanistically detailed systems biology modeling of the HGF/Met pathway in hepatocellular carcinoma

Hepatocyte growth factor (HGF) signaling through its receptor Met has been implicated in hepatocellular carcinoma tumorigenesis and progression. Met interaction with integrins is shown to modulate the downstream signaling to Akt and ERK (extracellular-regulated kinase). In this study, we developed a mechanistically detailed systems biology model of HGF/Met signaling pathway that incorporated specific interactions with integrins to investigate the efficacy of integrin-binding peptide, AXT050, as monotherapy and in combination with other therapeutics targeting this pathway. Here we report that the modeled dynamics of the response to AXT050 revealed that receptor trafficking is sufficient to explain the effect of Met-integrin interactions on HGF signaling. Furthermore, the model predicted patient-specific synergy and antagonism of efficacy and potency for combination of AXT050 with sorafenib, cabozantinib, and rilotumumab. Overall, the model provides a valuable framework for studying the efficacy of drugs targeting receptor tyrosine kinase interaction with integrins, and identification of synergistic drug combinations for the patients.

Multiple “Omics” data-based biomarker screening for hepatocellular carcinoma diagnosis

The huge prognostic difference between early and late stage hepatocellular carcinoma (HCC) is a challenging diagnostic problem. Alpha-fetoprotein is the mostly widely used biomarker for HCC used in the clinic, however it’s sensitivity and specificity of is not optimal. The development and application of multiple biotechnologies, including next generation sequencing, multiple “omics” data, that include genomics, epigenomics, transcriptomics, proteomics, metabolomics, metagenomics has been used for HCC diagnostic biomarker screening. Effective biomarkers/panels/models have been identified and validated at different clinical levels. A large proportion of these have a good diagnostic performance for HCC, especially for early HCC. In this article, we reviewed the various HCC biomarkers derived from “omics” data and discussed the advantages and disadvantages for diagnosis HCC.

TumorFusions: an integrative resource for cancer-associated transcript fusions

Gene fusion represents a class of molecular aberrations in cancer and has been exploited for therapeutic purposes. In this paper we describe TumorFusions, a data portal that catalogues 20 731 gene fusions detected in 9966 well characterized cancer samples and 648 normal specimens from The Cancer Genome Atlas (TCGA). The portal spans 33 cancer types in TCGA. Fusion transcripts were identified via a uniform pipeline, including filtering against a list of 3838 transcript fusions detected in a panel of 648 non-neoplastic samples. Fusions were mapped to somatic DNA rearrangements identified using whole genome sequencing data from 561 cancer samples as a means of validation. We observed that 65% of transcript fusions were associated with a chromosomal alteration, which is annotated in the portal. Other features of the portal include links to SNP array-based copy number levels and mutational patterns, exon and transcript level expressions of the partner genes, and a network-based centrality score for prioritizing functional fusions. Our portal aims to be a broadly applicable and user friendly resource for cancer gene annotation and is publicly available at http://www.tumorfusions.org.

A Pharmacogenomic Landscape in Human Liver Cancers

Liver cancers are highly heterogeneous with poor prognosis and drug response. A better understanding between genetic alterations and drug responses would facilitate precision treatment for liver cancers. To characterize the landscape of pharmacogenomic interactions in liver cancers, we developed a protocol to establish human liver cancer cell models at a success rate of around 50% and generated the Liver Cancer Model Repository (LIMORE) with 81 cell models. LIMORE represented genomic and transcriptomic heterogeneity of primary cancers. Interrogation of the pharmacogenomic landscape of LIMORE discovered unexplored gene-drug associations, including synthetic lethalities to prevalent alterations in liver cancers. Moreover, predictive biomarker candidates were suggested for the selection of sorafenib-responding patients. LIMORE provides a rich resource facilitating drug discovery in liver cancers.

TCF12 promotes the tumorigenesis and metastasis of hepatocellular carcinoma via upregulation of CXCR4 expression

TCF12, which is known to be involved in the regulation of cell growth and differentiation, has been reported to function as an oncogene or a tumor suppressor gene in the progression of various malignant tumors. However, its function and molecular mechanism in hepatocellular carcinoma (HCC) remain unclear.

Methods: Stable ectopic TCF12 expression or knockdown in HCC cell lines was established by lentiviral infection. Then, MTT, colony formation, migration, invasion and HUVECs tube formation assays as well as an orthotopic xenograft model were used to investigate the biologic function of TCF12 in HCC cells in vitro and in vivo. Subsequently, RNA-Seq analysis was utilized to explore the target genes regulated by TCF12. RT-qPCR, western blotting, a dual-luciferase reporter assay, Ch-IP, CHIP-Seq and functional rescue experiments were used to confirm the target gene regulated by TCF12. Finally, RT-qPCR, western blot and immunohistochemical (IHC) staining were performed to detect the expression level of TCF12 and to analyze the correlation of TCF12 with downstream genes as well as the clinical significance of TCF12 in human primary HCC.

Results: Our functional studies revealed that stable overexpression of TCF12 in human HCC cells enhanced cell proliferation, migration and invasion in vitro and in vivo, whereas knockdown of TCF12 showed opposing effects. Mechanistically, CXCR4 was a downstream target of TCF12, and TCF12 directly bound to the CXCR4 promoter to regulate its expression. Moreover, CXCR4, with its ligand CXCL12, played a critical role in tumor progression induced by TCF12 via activation of the MAPK/ERK and PI3K/AKT signaling pathways. Clinically, IHC analysis revealed that TCF12 was significantly associated with poor survival of HCC patients and that TCF12 expression was closely correlated with CXCR4 expression in primary HCC tissues.

Conclusion: Our findings are the first to indicate that TCF12 could promote the tumorigenesis and progression of HCC mainly by upregulating CXCR4 expression and is a prognostic indicator for patients with HCC.

SLMP53-2 Restores Wild-Type-Like Function to Mutant p53 through Hsp70: Promising Activity in Hepatocellular Carcinoma

Half of human cancers harbor TP53 mutations that render p53 inactive as a tumor suppressor. In these cancers, reactivation of mutant p53 (mutp53) through restoration of wild-type-like function constitutes a valuable anticancer therapeutic strategy. In order to search for mutp53 reactivators, a small library of tryptophanol-derived oxazoloisoindolinones was synthesized and the potential of these compounds as mutp53 reactivators and anticancer agents was investigated in human tumor cells and xenograft mouse models. By analysis of their anti-proliferative effect on a panel of p53-null NCI-H1299 tumor cells ectopically expressing highly prevalent mutp53, the compound SLMP53-2 was selected based on its potential reactivation of multiple structural mutp53. In mutp53-Y220C-expressing hepatocellular carcinoma (HCC) cells, SLMP53-2-induced growth inhibition was mediated by cell cycle arrest, apoptosis, and endoplasmic reticulum stress response. In these cells, SLMP53-2 restored wild-type-like conformation and DNA-binding ability of mutp53-Y220C by enhancing its interaction with the heat shock protein 70 (Hsp70), leading to the reestablishment of p53 transcriptional activity. Additionally, SLMP53-2 displayed synergistic effect with sorafenib, the only approved therapy for advanced HCC. Notably, it exhibited potent antitumor activity in human HCC xenograft mouse models with a favorable toxicological profile. Collectively, SLMP53-2 is a new mutp53-targeting agent with promising antitumor activity, particularly against HCC.

Comprehensive transcriptomic analysis of cell lines as models of primary tumors across 22 tumor types

Cancer cell lines are a cornerstone of cancer research but previous studies have shown that not all cell lines are equal in their ability to model primary tumors. Here we present a comprehensive pan-cancer analysis utilizing transcriptomic profiles from The Cancer Genome Atlas and the Cancer Cell Line Encyclopedia to evaluate cell lines as models of primary tumors across 22 tumor types. We perform correlation analysis and gene set enrichment analysis to understand the differences between cell lines and primary tumors. Additionally, we classify cell lines into tumor subtypes in 9 tumor types. We present our pancreatic cancer results as a case study and find that the commonly used cell line MIA PaCa-2 is transcriptionally unrepresentative of primary pancreatic adenocarcinomas. Lastly, we propose a new cell line panel, the TCGA-110-CL, for pan-cancer studies. This study provides a resource to help researchers select more representative cell line models.

Cell lines are used ubiquitously in cancer research but how well they represent the tumor type they were derived from is variable. Here, the authors compare transcriptomic profiles of 22 tumor types and cell lines and propose a new comprehensive cell line panel for pan-cancer studies.

Identification of crucial genes based on expression profiles of hepatocellular carcinomas by bioinformatics analysis

Hepatocellular carcinoma (HCC) is one of the most heterogeneous malignant cancers with no effective targets and treatments. However, the molecular pathogenesis of HCC remains largely uncertain. The aims of our study were to find crucial genes involved in HCC through multidimensional methods and revealed potential molecular mechanisms. Here, we reported the gene expression profile GSE121248 findings from 70 HCC and 37 adjacent normal tissues, all of which had chronic hepatitis B virus (HBV) infection, we were seeking to identify the dysregulated pathways, crucial genes and therapeutic targets implicated in HBV-associated HCC. We found 164 differentially expressed genes (DEGs) (92 downregulated genes and 72 upregulated genes). Gene ontology (GO) analysis of DEGs revealed significant functional enrichment of mitotic nuclear division, cell division, and the epoxygenase P450 pathway. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the DEGs were mainly enriched in metabolism, cell cycle regulation and the p53 signaling pathway. The Mcode plugin was calculated to construct a module complex of DEGs, and the module was mainly enriched in cell cycle checkpoints, RHO GTPase effectors and cytochrome P450. Considering a weak contribution of each gene, gene set enrichment analysis (GSEA) was performed, revealing results consistent with those described above. Six crucial proteins were selected based on the degree of centrality, including NDC80, ESR1, ZWINT, NCAPG, ENO3 and CENPF. Real-time quantitative PCR analysis validated the six crucial genes had the same expression trend as predicted. Furthermore, the methylation data of The Cancer Genome Atlas (TCGA) with HCC showed that mRNA expression of crucial genes was negatively correlated with methylation levels of their promoter region. The overall survival reflected that high expression of NDC80, CENPF, ZWINT, and NCAPG significantly predicted poor prognosis, whereas ESR1 high expression exhibited a favorable prognosis. The identification of the crucial genes and pathways would contribute to the development of novel molecular targets and biomarker-driven treatments for HCC.

Regional Chemotherapy for Biliary Tract Tumors and Hepatocellular Carcinoma

Locally advanced hepatocellular carcinoma and intrahepatic cholangiocarcinoma are associated with a grim prognosis. The development of highly effective systemic therapies for these tumors has been challenging; however, numerous locoregional treatment alternatives have emerged, including transarterial hepatic embolization (TAE), transarterial chemoembolization (TACE), drug-eluting bead TACE (DEB-TACE), hepatic arterial infusion chemotherapy (HAI), radioembolization, and stereotactic body radiation therapy. Although there is potential for long-term disease control for these therapies, the evidence to guide adequate patient selection and choose among different treatment alternatives is still limited. This review focuses on the rationale and data supporting TAE, TACE, DEB-TACE, and HAI in hepatobiliary cancers.

CAGI5: Objective performance assessments of predictions based on the Evolutionary Action equation

Many computational approaches estimate the effect of coding variants, but their predictions often disagree with each other. These contradictions confound users and raise questions regarding reliability. Performance assessments can indicate the expected accuracy for each method and highlight advantages and limitations. The Critical Assessment of Genome Interpretation (CAGI) community aims to organize objective and systematic assessments: They challenge predictors on unpublished experimental and clinical data and assign independent assessors to evaluate the submissions. We participated in CAGI experiments as predictors, using the Evolutionary Action (EA) method to estimate the fitness effect of coding mutations. EA is untrained, uses homology information, and relies on a formal equation: The fitness effect equals the functional sensitivity to residue changes multiplied by the magnitude of the substitution. In previous CAGI experiments (between 2011 and 2016), our submissions aimed to predict the protein activity of single mutants. In 2018 (CAGI5), we also submitted predictions regarding clinical associations, folding stability, and matching genomic data with phenotype. For all these diverse challenges, we used EA to predict the fitness effect of variants, adjusted to specifically address each question. Our submissions had consistently good performance, suggesting that EA predicts reliably the effects of genetic variants.

Identification of Coding and Long Noncoding RNAs Differentially Expressed in Tumors and Preferentially Expressed in Healthy Tissues

The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) datasets allow unprecedented gene expression analyses. Here, using these datasets, we performed pan-cancer and pan-tissue identification of coding and long noncoding RNA (lncRNA) transcripts differentially expressed in tumors and preferentially expressed in healthy tissues and/or tumors. Pan-cancer comparison of mRNAs and lncRNAs showed that lncRNAs were deregulated in a more tumor-specific manner. Given that lncRNAs are more tissue-specific than mRNAs, we identified healthy tissues that preferentially express lncRNAs upregulated in tumors and found that testis, brain, the digestive tract, and blood/spleen were the most prevalent. In addition, specific tumors also upregulate lncRNAs preferentially expressed in other tissues, generating a unique signature for each tumor type. Most tumors studied downregulated lncRNAs preferentially expressed in their tissue of origin, probably as a result of dedifferentiation. However, the same lncRNAs could be upregulated in other tumors, resulting in "bimorphic" transcripts. In hepatocellular carcinoma (HCC), the upregulated genes identified were expressed at higher levels in patients with worse prognosis. Some lncRNAs upregulated in HCC and preferentially expressed in healthy testis or brain were predicted to function as oncogenes and were significantly associated with higher tumor burden, and poor prognosis, suggesting their relevance in hepatocarcinogenesis and/or tumor evolution. Taken together, therapies targeting oncogenic lncRNAs should take into consideration the healthy tissue, where the lncRNAs are preferentially expressed, to predict and decrease unwanted secondary effects and increase potency. SIGNIFICANCE: Comprehensive analysis of coding and noncoding genes expressed in different tumors and normal tissues, which should be taken into account to predict side effects from potential coding and noncoding gene-targeting therapies.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/79/20/5167/F1.large.jpg.

MERIT: Systematic Analysis and Characterization of Mutational Effect on RNA Interactome Topology

The interaction between RNA-binding proteins (RBPs) and RNA plays an important role in regulating cellular function. However, decoding genome-wide protein-RNA regulatory networks as well as how cancer-related mutations impair RNA regulatory activities in hepatocellular carcinoma (HCC) remains mostly undetermined. We explored the genetic alteration patterns of RBPs and found that deleterious mutations are likely to occur on the surface of RBPs. We then constructed protein-RNA interactome networks by integration of target binding screens and expression profiles. Network analysis highlights regulatory principles among interacting RBPs. In addition, somatic mutations selectively target functionally important genes (cancer genes, core fitness genes, or conserved genes) and perturb the RBP-gene regulatory networks in cancer. These regulatory patterns were further validated using independent data. A computational method (Mutational Effect on RNA Interactome Topology) and a web-based, user-friendly resource were further proposed to analyze the RBP-gene regulatory networks across cancer types. Pan-cancer analysis also suggests that cancer cells selectively target "vulnerability" genes to perturb protein-RNA interactome that is involved in cancer hallmark-related functions. Specifically, we experimentally validated four pairs of RBP-gene interactions perturbed by mutations in HCC, which play critical roles in cell proliferation. Based on the expression of perturbed RBP and target genes, we identified three subtypes of HCC with different survival rates. Conclusion: Our results provide a valuable resource for characterizing somatic mutation-perturbed protein-RNA regulatory networks in HCC, yielding valuable insights into the genotype-phenotype relationships underlying human cancer, and potential biomarkers for precision medicine.

Gene Expression and DNA Methylation Alterations in the Glycine N-Methyltransferase Gene in Diet-Induced Nonalcoholic Fatty Liver Disease-Associated Carcinogenesis

Nonalcoholic fatty liver disease (NAFLD) is becoming a major etiological risk factor for hepatocellular carcinoma (HCC) in the United States and other Western countries. In this study, we investigated the role of gene-specific promoter cytosine DNA methylation and gene expression alterations in the development of NAFLD-associated HCC in mice using (1) a diet-induced animal model of NAFLD, (2) a Stelic Animal Model of nonalcoholic steatohepatitis-derived HCC, and (3) a choline- and folate-deficient (CFD) diet (CFD model). We found that the development of NAFLD and its progression to HCC was characterized by down-regulation of glycine N-methyltransferase (Gnmt) and this was mediated by progressive Gnmt promoter cytosine DNA hypermethylation. Using a panel of genetically diverse inbred mice, we observed that Gnmt down-regulation was an early event in the pathogenesis of NAFLD and correlated with the extent of the NAFLD-like liver injury. Reduced GNMT expression was also found in human HCC tissue and liver cancer cell lines. In in vitro experiments, we demonstrated that one of the consequences of GNMT inhibition was an increase in genome methylation facilitated by an elevated level of S-adenosyl-L-methionine. Overall, our findings suggest that reduced Gnmt expression caused by promoter hypermethylation is one of the key molecular events in the development of NAFLD-derived HCC and that assessing Gnmt methylation level may be useful for disease stratification.

Molecular portrait of high alpha-fetoprotein in hepatocellular carcinoma: implications for biomarker-driven clinical trials

The clinical utility of serum alpha-fetoprotein (AFP) in patients with hepatocellular carcinoma (HCC) is widely recognised. However, a clear understanding of the mechanisms of AFP overexpression and the molecular traits of patients with AFP-high tumours are not known. We assessed transcriptome data, whole-exome sequencing data and DNA methylome profiling of 520 HCC patients from two independent cohorts to identify distinct molecular traits of patients with AFP-high tumours (serum concentration > 400 ng/ml), which represents an accepted prognostic cut-off and a predictor of response to ramucirumab. Those AFP-high tumours (18% of resected cases) were characterised by significantly lower AFP promoter methylation (p < 0.001), significant enrichment of progenitor-cell features (CK19, EPCAM), higher incidence of BAP1 oncogene mutations (8.5% vs 1.6%) and lower mutational rates of CTNNB1 (14% vs 30%). Specifically, AFP-high tumours displayed significant activation of VEGF signalling (p < 0.001), which might provide the rationale for the reported benefit of ramucirumab in this subgroup of patients.

Feasibility of circulating tumor DNA testing in hepatocellular carcinoma

Background

Advancement of technologies enabling clinical assessment of circulating tumor DNA (ctDNA) are allowing for assessment of tumor specific genetic alterations in patients. This holds incredible promise for early detection of disease, serial monitoring of tumor heterogeneity, elucidation of therapeutic targets, and evaluation of treatment response and mechanisms of resistance. Hepatocellular carcinoma (HCC) is the most common primary malignancy of the liver and is often diagnosed late, recurs commonly, and is often diagnosed based upon imaging features alone. A comprehensive evaluation of real-time evaluation of ctDNA in patients with HCC has thus far not been undertaken.

Methods

From January 2015 to February 2018, 35 patients with biliary tract cancer (BTC) at the Mayo Clinic Comprehensive Cancer Center underwent ctDNA testing using a clinically available assay. The majority of samples were tested utilizing the 73-gene panel which includes somatic genomic targets, including complete or critical exon coverage in 30 and 40 genes, respectively; and in some, amplifications, fusions, and indels.

Results

A total of 44 samples were collected on these 35 patients, with >70% having stage 3 or 4 disease. Among all samples the median number of alterations per sample, excluding variants of undetermined significance (VUS), was 3.5, with a median allele frequency of 0.65%. A total of 122 unique genetic alterations, excluding VUS or synonymous alterations, were seen. The overall landscape of alterations is described. The top 10 genes altered in this cohort of patients, excluding VUS or synonymous alterations, were TP53 (18%), TERT (14%), CTNNB1 (13%), ARID1A (9%), MYC (5%), BRAF (4%), CCND1 (4%), CDK6 (4%), and MET (4%), and EGFR (3%).

Conclusions

Herein, we describe feasibility of ctDNA testing and results from such testing in HCC patients undergoing ctDNA testing in a real-time clinical context. Patients with these cancers stand to benefit immensely from the use of ctDNA technologies, and concerted efforts at further investigation of such are critically needed.

Liver-derived cell-free nucleic acids in plasma: Biology and applications in liquid biopsies

There is much global research interest surrounding the use of cell-free DNA (cfDNA) for liquid biopsies. cfDNA-based non-invasive prenatal testing for foetal chromosomal aneuploidies was the first successful application of cfDNA technology that transformed clinical practice - it has since been rapidly adopted in dozens of countries and is used by millions of pregnant women every year. Prompted by such developments, efforts to use cfDNA in other fields, especially for cancer detection and monitoring have been actively pursued in recent years. Cancer-associated aberrations including single nucleotide mutations, copy number aberrations, aberrations in methylation and alterations in DNA fragmentation patterns have been detected in the cfDNA of patients suffering from a wide variety of cancers. In addition, the analysis of methylation and fragmentomic patterns has enabled the tissue origin of cfDNA to be determined. In this review, different approaches for detecting circulating liver-derived nucleic acids and cancer-associated aberrations, as well as their potential clinical applications for the detection, monitoring and management of hepatocellular carcinoma, will be discussed.

Current State of Liver-Directed Therapies and Combinatory Approaches with Systemic Therapy in Hepatocellular Carcinoma (HCC)

The increasing set of liver-directed therapies (LDT) have become an integral part of hepatocellular carcinoma (HCC) treatment. These range from percutaneous ablative techniques to arterial embolization, and varied radiotherapy strategies. They are now used for local disease control, symptom palliation, and bold curative strategies. The big challenge in the face of these innovative and sometimes overlapping technologies is to identify the best opportunity of use. In real practice, many patients may take benefit from LDT used as a bridge to curative treatment such as resection and liver transplantation. Varying trans-arterial embolization strategies are used, and comparison between established and developing technologies is scarce. Also, radioembolization utilizing yttrium-90 (Y-90) for locally advanced or intermediate-stage HCC needs further evidence of clinical efficacy. There is increasing interest on LDT-led changes in tumor biology that could have implications in systemic therapy efficacy. Foremost, additional to its apoptotic and necrotic properties, LDT could warrant changes in vascular endothelial growth factor (VEGF) expression and release. However, trans-arterial chemoembolization (TACE) used alongside tyrosine-kinase inhibitor (TKI) sorafenib has had its efficacy contested. Most recently, interest in associating Y-90 and TKI has emerged. Furthermore, LDT-led differences in tumor immune microenvironment and immune cell infiltration could be an opportunity to enhance immunotherapy efficacy for HCC patients. Early attempts to coordinate LDT and immunotherapy are being made. We here review LDT techniques exposing current evidence to understand its extant reach and future applications alongside systemic therapy development for HCC.

lncRNA XIST regulates proliferation and migration of hepatocellular carcinoma cells by acting as miR-497-5p molecular sponge and targeting PDCD4

BACKGROUND

MicroRNAs (miRNAs) play a pivotal role in hepatocellular carcinoma (HCC) progression and have been confirmed to participate in the carcinogenesis and development of HCC. However, the relationship between miR-497-5p and HCC remains unclear.

METHODS

Kaplan-Meier curve analysis and the log-rank test were used to investigate the efficacy of miR-497-5p on overall survival (OS) and disease-free survival (DFS) in patients with HCC. According to in vitro experiments, programmed cell death 4 (PDCD4) was a target of miR-497-5p by the dual-luciferase activity assay. The efficacy of PDCD4 on cell proliferation and metastasis in HCC was examined by transwell assays, CCK-8 assays and reverse transcription quantitative PCR (RT-qPCR). Additionally, we conducted a luciferase activity reporter assay to confirm the interaction between lncRNA XIST and miR-49-5p. Then, to evaluate the relationship between lncRNA XIST and miR-497-5p, several mechanistic experiments, including qRT-PCR, Western blotting, transwell assays and tumor xenograft assays, were performed.

RESULTS

miR-497-5p was upregulated in HCC tissues, and high expression of miR-497-5p resulted in increases in tumor size and tumor number and a higher tumor-node-metastasis (TNM) stage and Edmondson grade in patients with HCC. Silencing miR-497-5p inhibited the proliferation and migration of HCC cells. PDCD4, which was downregulated in HCC tissues, was shown to be a target of miR-497-5p and was negatively correlated with the expression of miR-497-5p. lncRNA XIST was found to act as a miR-497-5p sponge and to regulate the level of PDCD4, which is targeted by miR-497-5p. lncRNA XIST was observed to be downregulated in the HCC tissues and positively correlated with the expression of PDCD4.

CONCLUSIONS

Our findings reveal that the XIST/miR-497-5p/PDCD4 axis participates in HCC development and that XIST could be used as a biomarker of HCC.

Immunomodulatory TGF-β Signaling in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is an inflammation-induced and chemotherapy-resistant cancer. Dysregulated signaling in the transforming growth factor beta (TGF-β) pathway plays a central role in inflammation, fibrogenesis, and immunomodulation in the HCC microenvironment. This review dissects the genetic landscape of the TGF-β superfamily genes in HCC and discusses the essential effects of this pathway on the tumor immune microenvironment. We highlight the TGF-β signature as a potential biomarker for identifying individualized immunotherapeutic approaches in HCC. An improved understanding of the detailed mechanisms of liver cancer immunogenicity and the specific role of TGF-β in mediating immunotherapy resistance in HCC will provide important insights into HCC immune escape and promote the development of biomarker-derived combination immunotherapies for HCC.

Signatures of gene expression, DNA methylation and microRNAs of hepatocellular carcinoma with vascular invasion

Aim: Micro and macro vascular invasion (VI) are known as independent predictors of tumor recurrence and poor survival after surgical treatment of hepatocellular carcinoma (HCC). Here, we aimed to re-analyze The Cancer Genome Atlas of liver hepatocellular carcinoma datasets to identify the VI-expression signatures. Materials & methods: We filtered The Cancer Genome Atlas liver hepatocellular carcinoma (LIHC) datasets into three groups: no VI (NVI = 198); micro VI (MIVI = 89) and macro VI (MAVI = 16). We performed differential gene expression, methylation and microRNA analyses. Results & conclusion: We identified 12 differentially expressed genes and 55 differentially methylated genes in MAVI compared with no VI. The GPD1L gene appeared in all of the comparative analyses. Higher GPD1L expression was associated with VI and poor outcomes in the HCC patients.

Immunotherapy for hepatocellular carcinoma: recent advances and future perspectives

The introduction of immunotherapies has been a major development in the treatment of many advanced cancers, including hepatocellular carcinoma (HCC). We are entering a new era of systemic therapy for advanced HCC associated with an explosion of clinical trial activity. Data from phase I/II studies of checkpoint inhibitors in advanced HCC have been promising, with durable objective response rates of approximately 20% seen (in both first- and second-line settings) and acceptable safety profiles (including immune-mediated hepatitis). Phase III studies evaluating anti-programmed cell death protein 1 (anti-PD-1) and anti-programmed cell death ligand 1 (anti-PD-L1) antibodies compared with sorafenib are already underway. The potential synergistic effects of anti-PD-1/anti-PD-L1 when used in combination with agents against other checkpoint molecules, systemic therapies, as well as conventional surgical and locoregional therapies are also being explored in upcoming clinical trials. Aside from this, other strategies to harness the immune system, including chimeric antigen receptor-engineered T cells, natural killer cell therapies, and peptide vaccines directed against HCC antigens have entered phase I/II studies. Current limitations of immunotherapies and areas of future research include the accurate assessment and prediction of tumor response, overcoming the immunosuppressive effects of a hypoxic microenvironment, and the management of immune-related hepatitis in patients who already have limited liver reserve.

Genomic landscape of lymphoepithelioma-like hepatocellular carcinoma

Lymphoepithelioma-like hepatocellular carcinoma (LEL-HCC) is a distinct variant of HCC that is characterized by dense tumor-infiltrating lymphocytes (TILs). Patients with LEL-HCC also show better clinical outcomes compared to conventional HCC (c-HCC), which is commonly presented with low TIL. Emerging evidence has begun to highlight tumor-intrinsic genetic abnormalities in the tumor-host immune interfaces. However, genome-wide characterization of LEL-HCC remains largely unexplored. Here, we defined the genomic landscape of 12 LEL-HCC using whole-exome sequencing, and further underpinned those genetic alterations related to an immune active microenvironment by comparing findings to 15 c-HCC that were sequenced in parallel. Overall, the mutational load between LEL-HCC and c-HCC was similar. Interestingly, SNV incidences of specific genes (CTNNB1, AXIN1, NOTCH1, and NOTCH2) were significantly higher in c-HCC than LEL-HCC, suggesting a plausible link between activated Wnt/β-catenin and Notch signaling pathways and immune avoidance. Marked focal amplification of chromosome 11q13.3 was prevalent in LEL-HCC. Using The Cancer Genome Atlas dataset, we further established oncogenes expressed from chromosome 11q13.3 (CCND1, FGF19, and FGF4) to be strongly associated with the immune checkpoint signature (CD274, PDCD1, BTLA, CTLA4, HAVCR2, IDO1, and LAG3). Our results have illustrated for the first time the somatic landscape of LEL-HCC, and highlighted molecular alterations that could be exploited in combinatory therapy with checkpoint inhibitors in targeting HCC. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Identification of hedgehog signaling as a potential oncogenic driver in an aggressive subclass of human hepatocellular carcinoma: A reanalysis of the TCGA cohort

Hepatocellular carcinoma (HCC) is a heterogeneous disease and the second most common cause of cancer-related death worldwide. Marked developments in genomic technologies helped scientists to understand the heterogeneity of HCC and identified multiple HCC-related molecular subclasses. An integrative analysis of genomic datasets including 196 patients from The Cancer Genome Atlas (TCGA) group has recently reported a new HCC subclass, which contains three subgroups (iCluster1, iCluster2, and iCluster3). However, the transcriptional molecular characteristics underlying the iClusters have not been thoroughly investigated. Herein, we identified a more aggressive subset of HCC patients in the iCluster1, and re-clustered the TCGA samples into novel HCC subclasses referred to as aggressive (Ag), moderate-aggressive (M-Ag), and less-aggressive (L-Ag) subclasses. The Ag subclass had a greater predictive power than the TCGA iCluster1, and a higher level of alpha fetoprotein, microscopic vascular invasion, immune infiltration, isocitrate dehydrogenase 1/2 mutation status, and a worse survival than M-Ag and L-Ag subclasses. Global transcriptomic analysis showed that activation of hedgehog signaling in the Ag subclass may play key roles in tumor development of aggressive HCC. GLI1, a key transcriptional regulator of hedgehog signaling upregulated in the Ag subclass, was correlated with poor prognosis of HCC, and may be a potential prognostic biomarker and therapeutic target for Ag subclass HCC patients.

Metabolic Regulation of Redox Balance in Cancer

Reactive oxygen species (ROS) are chemically active free radicals produced by partial reduction of oxygen that can activate discrete signaling pathways or disrupt redox homeostasis depending on their concentration. ROS interacts with biomolecules, including DNA, and can cause mutations that can transform normal cells into cancer cells. Furthermore, certain cancer-causing mutations trigger alterations in cellular metabolism that can increase ROS production, resulting in genomic instability, additional DNA mutations, and tumor evolution. To prevent excess ROS-mediated toxicity, cancer-causing mutations concurrently activate pathways that manage this oxidative burden. Hence, an understanding of the metabolic pathways that regulate ROS levels is imperative for devising therapies that target tumor cells. In this review, we summarize the dual role of metabolism as a generator and inhibitor of ROS in cancer and discuss current strategies to target the ROS axis.

Loss of tumor suppressor IGFBP4 drives epigenetic reprogramming in hepatic carcinogenesis

Genomic sequencing of hepatocellular carcinoma (HCC) uncovers a paucity of actionable mutations, underscoring the necessity to exploit epigenetic vulnerabilities for therapeutics. In HCC, EZH2-mediated H3K27me3 represents a major oncogenic chromatin modification, but how it modulates the therapeutic vulnerability of signaling pathways remains unknown. Here, we show EZH2 acts antagonistically to AKT signaling in maintaining H3K27 methylome through epigenetic silencing of IGFBP4. ChIP-seq revealed enrichment of Ezh2/H3K27me3 at silenced loci in HBx-transgenic mouse-derived HCCs, including Igfbp4 whose down-regulation significantly correlated with EZH2 overexpression and poor survivals of HCC patients. Functional characterizations demonstrated potent growth- and invasion-suppressive functions of IGFBP4, which was associated with transcriptomic alterations leading to deregulation of multiple signaling pathways. Mechanistically, IGFBP4 stimulated AKT/EZH2 phosphorylation to abrogate H3K27me3-mediated silencing, forming a reciprocal feedback loop that suppressed core transcription factor networks (FOXA1/HNF1A/HNF4A/KLF9/NR1H4) for normal liver homeostasis. Consequently, the in vivo tumorigenicity of IGFBP4-silenced HCC cells was vulnerable to pharmacological inhibition of EZH2, but not AKT. Our study unveils chromatin regulation of a novel liver tumor suppressor IGFBP4, which constitutes an AKT-EZH2 reciprocal loop in driving H3K27me3-mediated epigenetic reprogramming. Defining the aberrant chromatin landscape of HCC sheds light into the mechanistic basis of effective EZH2-targeted inhibition.

Hepatocellular carcinoma: Therapeutic advances in signaling, epigenetic and immune targets

Hepatocellular carcinoma (HCC) remains a global medical burden with rising incidence due to chronic viral hepatitis and non-alcoholic fatty liver diseases. Treatment of advanced disease stages is still unsatisfying. Besides first and second generation tyrosine kinase inhibitors, immune checkpoint inhibitors have become central for the treatment of HCC. New modalities like epigenetic therapy using histone deacetylase inhibitors (HDACi) and cell therapy approaches with chimeric antigen receptor T cells (CAR-T cells) are currently under investigation in clinical trials. Development of such novel drugs is closely linked to the availability and improvement of novel preclinical and animal models and the identification of predictive biomarkers. The current status of treatment options for advanced HCC, emerging novel therapeutic approaches and different preclinical models for HCC drug discovery and development are reviewed here.

Genome-wide promoter DNA methylation profiling of hepatocellular carcinomas arising either spontaneously or due to chronic exposure to Ginkgo biloba extract (GBE) in B6C3F1/N mice

Epigenetic modifications, such as DNA methylation, play an important role in carcinogenesis. In a recent NTP study, chronic exposure of B6C3F1/N mice to Ginkgo biloba extract (GBE) resulted in a high incidence of hepatocellular carcinomas (HCC). Genome-wide promoter methylation profiling on GBE-exposed HCC (2000 mg/kg group), spontaneous HCC (vehicle-control group), and age-matched vehicle control liver was performed to identify differentially methylated genes in GBE-exposed HCC and spontaneous HCC. DNA methylation alterations were correlated to the corresponding global gene expression changes. Compared to control liver, 1296 gene promoters (719 hypermethylated, 577 hypomethylated) in GBE-exposed HCC and 738 (427 hypermethylated, 311 hypomethylated) gene promoters in spontaneous HCC were significantly differentially methylated, suggesting an impact of methylation on GBE-exposed HCC. Differential methylation of promoter regions in relevant cancer genes (cMyc, Spry2, Dusp5) and their corresponding differential gene expression was validated by quantitative pyrosequencing and qRT-PCR, respectively. In conclusion, we have identified differentially methylated promoter regions of relevant cancer genes altered in GBE-exposed HCC compared to spontaneous HCC. Further study of unique sets of differentially methylated genes in chemical-exposed mouse HCC could potentially be used to differentiate treatment-related tumors from spontaneous-tumors in cancer bioassays and provide additional understanding of the underlying epigenetic mechanisms of chemical carcinogenesis.

Cullin-3/KCTD10 complex is essential for K27-polyubiquitination of EIF3D in human hepatocellular carcinoma HepG2 cells

Eukaryotic translation initiation factor 3 subunit D (EIF3D) binds to the 5'-cap of specific mRNAs, initiating their translation into polypeptides. From a pathological standpoint, EIF3D has been observed to be essential for cell growth in various cancer types, and cancer patients with high EIF3D mRNA levels exhibit poor prognosis, indicating involvement of EIF3D in oncogenesis. In this study, we found, by mass spectrometry, that Cullin-3 (CUL3)/KCTD10 ubiquitin (Ub) ligase forms a complex with EIF3D. We also demonstrated that EIF3D is K27-polyubiquitinated at the lysine 153 and 275 residues in a KCTD10-dependent manner in human hepatocellular carcinoma HepG2 cells. Similar to other cancers, high expression of EIF3D significantly correlated with poor prognosis in hepatocellular carcinoma patients, and depletion of EIF3D drastically suppressed HepG2 cell proliferation. These results indicate that EIF3D is a novel substrate of CUL3/KCTD10 Ub ligase and suggest involvement of K27-polyubiquitinated EIF3D in the development of hepatocellular carcinoma.

Novel significant stage-specific differentially expressed genes in hepatocellular carcinoma

BACKGROUND

Liver cancer is among top deadly cancers worldwide with a very poor prognosis, and the liver is a vulnerable site for metastases of other cancers. Early diagnosis is crucial for treatment of the predominant liver cancers, namely hepatocellular carcinoma (HCC). Here we developed a novel computational framework for the stage-specific analysis of HCC.

METHODS

Using publicly available clinical and RNA-Seq data of cancer samples and controls and the AJCC staging system, we performed a linear modelling analysis of gene expression across all stages and found significant genome-wide changes in the log fold-change of gene expression in cancer samples relative to control. To identify genes that were stage-specific controlling for confounding differential expression in other stages, we developed a set of six pairwise contrasts between the stages and enforced a p-value threshold (< 0.05) for each such contrast. Genes were specific for a stage if they passed all the significance filters for that stage. The monotonicity of gene expression with cancer progression was analyzed with a linear model using the cancer stage as a numeric variable.

RESULTS

Our analysis yielded two stage-I specific genes (CA9, WNT7B), two stage-II specific genes (APOBEC3B, FAM186A), ten stage-III specific genes including DLG5, PARI, NCAPG2, GNMT and XRCC2, and 35 stage-IV specific genes including GABRD, PGAM2, PECAM1 and CXCR2P1. Overexpression of DLG5 was found to be tumor-promoting contrary to the cancer literature on this gene. Further, GABRD was found to be signifincantly monotonically upregulated across stages. Our work has revealed 1977 genes with significant monotonic patterns of expression across cancer stages. NDUFA4L2, CRHBP and PIGU were top genes with monotonic changes of expression across cancer stages that could represent promising targets for therapy. Comparison with gene signatures from the BCLC staging system identified two genes, HSP90AB1 and ARHGAP42. Gene set enrichment analysis indicated overrepresented pathways specific to each stage, notably viral infection pathways in HCC initiation.

CONCLUSIONS

Our study identified novel significant stage-specific differentially expressed genes which could enhance our understanding of the molecular determinants of hepatocellular carcinoma progression. Our findings could serve as biomarkers that potentially underpin diagnosis as well as pinpoint therapeutic targets.

Characterization of pediatric hepatocellular carcinoma reveals genomic heterogeneity and diverse signaling pathway activation

BACKGROUND

Pediatric hepatocellular carcinoma (HCC) is a rare liver tumor in children with a poor prognosis. Comprehensive molecular profiling to understand the underlying genomic drivers of this tumor has not been completed, and it is unclear whether nonfibrolamellar pediatric HCC is more genomically similar to hepatoblastoma or adult HCC.

PROCEDURE

To characterize the molecular landscape of these tumors, we analyzed a cohort of 15 pediatric non-FL-HCCs by sequencing a panel of cancer-associated genes and conducting copy-number and gene-expression analyses.

RESULTS

We detected multiple types of molecular alterations in Wnt signaling genes, including APC inversion, AMER1 somatic mutation, and most commonly CTNNB1 intragenic deletions. There were multiple alterations to the telomerase pathway via TERT activation or ATRX mutation. Therapeutically targetable activating mutations in MAPK/ERK signaling pathway genes, including MAPK1 and BRAF, were detected in 20% of tumors. TP53 mutations occurred far less frequently in our pediatric HCC cohort than reported in adult cohorts. Tumors arising in children with underlying liver disease were found to be molecularly distinct from the remainder and lacking detectable oncogenic drivers, as compared with those arising in patients without a history of underlying liver disease; the majority of both types were positive for glypican-3, another potential therapeutic target.

CONCLUSION

Our study revealed pediatric HCC to be a molecularly heterogeneous group of tumors. Those non-FL-HCC tumors arising in the absence of underlying liver disease harbor genetic alterations affecting multiple cancer pathways, most notably Wnt signaling, and share some characteristics with adult HCC.

Sorafenib for Treatment of Hepatocellular Carcinoma: A Survival Analysis From the South American Liver Research Network

GOALS

We aim to describe the efficacy, safety profile, and variables associated with survival in patients with hepatocellular carcinoma (HCC) treated with sorafenib in South America.

BACKGROUND

Sorafenib has been shown to improve survival in patients with advanced HCC. There are few data on sorafenib use for HCC in South America.

STUDY

We performed a retrospective analysis of HCC cases treated with sorafenib from 8 medical centers in 5 South American countries, between January 2010 and June 2017. The primary endpoint was overall survival (OS), which was defined as time from sorafenib initiation to death or last follow-up. Risk factors for decreased OS were assessed using Cox proportional hazard regression and log-rank tests.

RESULTS

Of 1336 evaluated patients, 127 were treated with sorafenib and were included in the study. The median age of individuals was 65 years (interquartile range, 55 to 71) and 70% were male individuals. Median OS in all patients was 8 months (interquartile range, 2 to 17). Variables associated with survival on multivariate analysis were platelets >/<250,000 mm (2 vs. 8 mo, P=0.01) and Barcelona Clinic Liver Cancer (BCLC) stage (A/B, 13 vs. C/D, 6 mo; P=0.04). In a subanalysis of patients with BCLC stage C, platelets >/<250,000 mm were also independently associated with survival (2 vs. 5.5 mo, P=0.03). Patients lived longer if they experienced any side effects from sorafenib use (11 vs. 2 mo, P=0.009). Patients who stopped sorafenib because of side effects had shorter survival compared with patients who were able to tolerate side effects and continue treatment (7.5 vs. 13 mo, P=0.01).

CONCLUSIONS

Pretreatment elevation of platelets and advanced BCLC stage were independently associated with poor survival on sorafenib in a South American cohort.

Integrative analysis of DNA methylation and gene expression identify a six epigenetic driver signature for predicting prognosis in hepatocellular carcinoma

DNA methylation is a crucial regulator of gene transcription in the etiology and pathogenesis of hepatocellular carcinoma (HCC). Thus, it is reasonable to identify DNA methylation-related prognostic markers. Currently, we aimed to make an integrative epigenetic analysis of HCC to identify the effectiveness of epigenetic drivers in predicting prognosis for HCC patients. By the software pipeline TCGA-Assembler 2, RNA-seq, and methylation data were downloaded and processed from The Cancer Genome Atlas. A bioconductor package MethylMix was utilized to incorporate gene expression and methylation data on all 363 samples and identify 589 epigenetic drivers with transcriptionally predictive. By univariate survival analysis, 72 epigenetic drivers correlated with overall survival (OS) were selected for further analysis in our training cohort. By the robust likelihood-based survival model, six epi-drivers (doublecortin domain containing 2, flavin containing monooxygenase 3, G protein-coupled receptor 171, Lck interacting transmembrane adaptor 1, S100 calcium binding protein P, small nucleolar RNA host gene 6) serving as prognostic markers was identified and then a DNA methylation signature for HCC (MSH) predicting OS was identified to stratify patients into low-risk and high-risk groups in the training cohort (p < 0.001). The capability of MSH was also assessed in the validation cohort (p = 0.002). Furthermore, a receiver operating characteristic curve confirmed MSH as an effective prognostic model for predicting OS in HCC patients in training area under curve (AUC = 0.802) and validation (AUC = 0.691) cohorts. Finally, a nomogram comprising MSH and pathologic stage was generated to predict OS in the training cohort, and it also operated effectively in the validation cohort (concordance index: 0.674). In conclusion, MSH, a six epi-drivers based signature, is a potential model to predict prognosis for HCC patients.

Comprehensive analysis of genetic aberrations linked to tumorigenesis in regenerative nodules of liver cirrhosis

BACKGROUND

Hepatocellular carcinoma (HCC) recurrently develops in cirrhotic liver containing a number of regenerative nodules (RNs). However, the biological tumorigenic potential of RNs is still unclear. To uncover the molecular bases of tumorigenesis in liver cirrhosis, we investigated the genetic aberrations in RNs of cirrhotic tissues using next-generation sequencing.

METHODS

We isolated 205 RNs and 7 HCC tissues from the whole explanted livers of 10 randomly selected patients who had undergone living-donor liver transplantation. Whole-exome sequencing and additional targeted deep sequencing on 30 selected HCC-related genes were conducted to reveal the mutational landscape of RNs and HCCs.

RESULTS

Whole-exome sequencing demonstrated that RNs frequently harbored relatively high-abundance genetic alterations, suggesting a clonal structure of each RN in cirrhotic liver. The mutation signature observed in RNs was similar to those determined in HCC, characterized by a predominance of C>T transitions, followed by T>C and C>A mutations. Targeted deep sequencing analyses of RNs identified nonsynonymous low-abundance mutations in various tumor-related genes, including TP53 and ARID1A. In contrast, TERT promoter mutations were not detected in any of the RNs examined. Consistently, TERT expression levels in RNs were comparable to those in normal livers, whereas every HCC tissue demonstrated an elevated level of TERT expression.

CONCLUSION

Analyses of RNs constructing cirrhotic liver indicated that a variety of genetic aberrations accumulate in the cirrhotic liver before the development of clinically and histologically overt HCC. These aberrations in RNs could provide the basis of tumorigenesis in patients with liver cirrhosis.

Ccne1 Overexpression Causes Chromosome Instability in Liver Cells and Liver Tumor Development in Mice

BACKGROUND & AIMS

The CCNE1 locus, which encodes cyclin E1, is amplified in many types of cancer cells and is activated in hepatocellular carcinomas (HCCs) from patients infected with hepatitis B virus or adeno-associated virus type 2, due to integration of the virus nearby. We investigated cell-cycle and oncogenic effects of cyclin E1 overexpression in tissues of mice.

METHODS

We generated mice with doxycycline-inducible expression of Ccne1 (Ccne1^T mice) and activated overexpression of cyclin E1 from age 3 weeks onward. At 14 months of age, livers were collected from mice that overexpress cyclin E1 and nontransgenic mice (controls) and analyzed for tumor burden and by histology. Mouse embryonic fibroblasts (MEFs) and hepatocytes from Ccne1^T and control mice were analyzed to determine the extent to which cyclin E1 overexpression perturbs S-phase entry, DNA replication, and numbers and structures of chromosomes. Tissues from 4-month-old Ccne1^T and control mice (at that age were free of tumors) were analyzed for chromosome alterations, to investigate the mechanisms by which cyclin E1 predisposes hepatocytes to transformation.

RESULTS

Ccne1^T mice developed more hepatocellular adenomas and HCCs than control mice. Tumors developed only in livers of Ccne1^T mice, despite high levels of cyclin E1 in other tissues. Ccne1^T MEFs had defects that promoted chromosome missegregation and aneuploidy, including incomplete replication of DNA, centrosome amplification, and formation of nonperpendicular mitotic spindles. Whereas Ccne1^T mice accumulated near-diploid aneuploid cells in multiple tissues and organs, polyploidization was observed only in hepatocytes, with losses and gains of whole chromosomes, DNA damage, and oxidative stress.

CONCLUSIONS

Livers, but not other tissues of mice with inducible overexpression of cyclin E1, develop tumors. More hepatocytes from the cyclin E1-overexpressing mice were polyploid than from control mice, and had losses or gains of whole chromosomes, DNA damage, and oxidative stress; all of these have been observed in human HCC cells. The increased risk of HCC in patients with hepatitis B virus or adeno-associated virus type 2 infection might involve activation of cyclin E1 and its effects on chromosomes and genomes of liver cells.

Macrotrabecular Hepatocellular Carcinoma: An Aggressive Subtype of Hepatocellular Carcinoma

The macrotrabecular (MT) pattern of hepatocellular carcinoma (HCC) has been suggested to represent a distinct HCC subtype. We retrospectively reviewed 231 HCC cases. Detailed pathologic evaluation for histologic patterns, including MT-pattern, was performed for each case and recorded as percentage involved at 10% intervals. MT-pattern was defined as having trabeculae >6 cells thick. After excluding all recognized HCC subtypes, remaining cases were deemed conventional HCC (CV-HCC) and served as controls. HCCs with a component of ≥10%, ≥30% and ≥50% MT-pattern were identified in 41 (17.7%), 24 (10.4%) and 4 (1.7%) cases, respectively. The clinicopathologic features of HCCs with 10% to 29% MT-pattern (n=17, 7.4%) were largely similar to CV-HCC. No significant difference was observed between the 30% and 49% (n=20) and ≥50% (n=4) MT groups, hence these were combined for further analysis as MT-HCC. MT-HCCs (≥30% MT-pattern) were larger tumors (5.5 vs. 3.1 cm), were less likely to be associated with cirrhosis (54% vs. 79%), were more likely to have hepatitis B (21% vs. 5%) and less likely hepatitis C infection (33% vs. 58%) compared with CV-HCC. MT-HCC was associated with the presence of anaplastic tumor cells (42% vs. 14%), higher alpha-fetoprotein level, higher AJCC stage, and higher histologic grade. Compared with patients with CV-HCC, patients with MT-HCC had poorer overall survival. Patients with MT-HCC who underwent primary resection or transplantation had a higher recurrence rate and worse recurrence-free survival. Our findings suggest that ≥30% MT-pattern could be used as the more appropriate cut-off for defining MT-HCC, which represents a unique and aggressive HCC histologic subtype.

Alpha-1 Antitrypsin Deficiency Liver Disease, Mutational Homogeneity Modulated by Epigenetic Heterogeneity With Links to Obesity

Alpha-1 antitrypsin deficiency (AATD) liver disease is characterized by marked heterogeneity in presentation and progression, despite a common underlying gene mutation, strongly suggesting the involvement of other genetic and/or epigenetic modifiers. Variation in clinical phenotype has added to the challenge of detection, diagnosis, and testing of new therapies in patients with AATD. We examined the contribution of DNA methylation (5-methylcytosine [5mC]) to AATD liver disease heterogeneity because 5mC responds to environmental and genetic cues and its deregulation is a major driver of liver disease. Using liver biopsies from adults with early-stage AATD and the ZZ genotype, genome-wide 5mC patterns were interrogated. We compared DNA methylation among patients with early AATD, and among patients with normal liver, cirrhosis, and hepatocellular carcinoma derived from multiple etiologic exposures, and linked patient clinical/demographic features. Global analysis revealed significant genomic hypomethylation in AATD liver-impacting genes related to liver cancer, cell cycle, and fibrosis, as well as key regulatory molecules influencing growth, migration, and immune function. Further analysis indicated that 5mC changes are localized, with hypermethylation occurring within a background of genome-wide 5mC loss and with patients with AATD manifesting distinct epigenetic landscapes despite their mutational homogeneity. By integrating clinical data with 5mC landscapes, we observed that CpGs differentially methylated among patients with AATD disease are linked to hallmark clinical features of AATD (e.g., hepatocyte degeneration and polymer accumulation) and further reveal links to well-known sex-specific effects of liver disease progression. Conclusion: Our data reveal molecular epigenetic signatures within this mutationally homogeneous group that point to ways to stratify patients for liver disease risk.

Immunogenomic landscape of hepatocellular carcinoma with immune cell stroma and EBV-positive tumor-infiltrating lymphocytes

BACKGROUND & AIMS

The immunogenomic characteristics of hepatocellular carcinomas (HCCs) with immune cell stroma (HCC-IS), defined histologically, have not been clarified. We investigated the clinical and molecular features of HCC-IS and the prognostic impact of Epstein-Barr virus (EBV) infection.

METHODS

We evaluated 219 patients with conventional HCC (C-HCC) and 47 with HCC-IS using in situ hybridization for EBV, immunohistochemistry, multiplex immunofluorescence staining, and whole exome and transcriptome sequencing. Human leukocyte antigen types were also extracted from the sequencing data. Genomic and prognostic parameters were compared between HCC-IS and C-HCC.

RESULTS

CD8 T cell infiltration was more frequent in HCC-IS than C-HCC (mean fraction/sample, 22.6% vs. 8.9%, false discovery rate q <0.001), as was EBV positivity in CD20-positive tumor-infiltrating lymphocytes (TILs) (74.5% vs. 4.6%, p <0.001). CTNNB1 mutations were not identified in any HCC-IS, while they were present in 24.1% of C-HCC (p = 0.016). Inhibitory and stimulatory immune modulators were expressed at similar levels in HCC-IS and EBV-positive C-HCC. Global hypermethylation, and expression of PD-1 and PD-L1 in TILs, and PD-L1 in tumors, were also associated with HCC-IS (p <0.001), whereas human leukocyte antigen type did not differ according to HCC type or EBV positivity. HCC-IS was an independent factor for favorable recurrence-free survival (adjusted hazard ratio [aHR] 0.23; p = 0.002). However, a subgroup of tumors with a high density of EBV-positive TILs had poorer recurrence-free (aHR 25.48; p <0.001) and overall (aHR 9.6; p = 0.003) survival, and significant enrichment of CD8 T cell exhaustion signatures (q = 0.0296).

CONCLUSIONS

HCC-IS is a distinct HCC subtype associated with a good prognosis and frequent EBV-positive TILs. However, paradoxically, a high density of EBV-positive TILs in tumors is associated with inferior prognostic outcomes. Patients with HCC-IS could be candidates for immunotherapy.

LAY SUMMARY

Hepatocellular carcinomas with histologic evidence of abundant immune cell infiltration are characterized by frequent activation of Epstein-Barr virus in tumor-infiltrating lymphocytes and less aggressive clinical behavior. However, a high density of Epstein-Barr virus-positive tumor-infiltrating lymphocytes is associated with inferior prognostic outcomes, possibly as a result of immune escape due to significant CD8 T cell exhaustion.

Immunotherapy for hepatocellular carcinoma: Current and future

Hepatocellular carcinoma (HCC) arises on the background of chronic liver disease. Despite the development of effective anti-viral therapeutics HCC is continuing to rise, in part driven by the epidemic of non-alcoholic fatty liver disease. Many patients present with advanced disease out with the criteria for transplant, resection or even locoregional therapy. Currently available therapeutics for HCC are effective in a small minority of individuals. However, there has been a major global interest in immunotherapies for cancer and although HCC has lagged behind other cancers, great opportunities now exist for treating HCC with newer and more sophisticated agents. Whilst checkpoint inhibitors are at the forefront of this revolution, other therapeutics such as inhibitory cytokine blockade, oncolytic viruses, adoptive cellular therapies and vaccines are emerging. Broadly these may be categorized as either boosting existing immune response or stimulating de novo immune response. Although some of these agents have shown promising results as monotherapy in early phase trials it may well be that their future role will be as combination therapy, either in combination with one another or in combination with treatment modalities such as locoregional therapy. Together these agents are likely to generate new and exciting opportunities for treating HCC, which are summarized in this review.

Expression of genes that control core fucosylation in hepatocellular carcinoma: Systematic review

BACKGROUND

Changes in N-linked glycosylation have been observed in the circulation of individuals with hepatocellular carcinoma. In particular, an elevation in the level of core fucosylation has been observed. However, the mechanisms through which core fucose is increased are not well understood. We hypothesized that a review of the literature and related bioinformatic review regarding six genes known to be involved in the attachment of core fucosylation, the synthesis of the fucosylation substrate guanosine diphosphate (GDP)-fucose, or the transport of the substrate into the Golgi might offer mechanistic insight into the regulation of core fucose levels.

AIM

To survey the literature to capture the involvement of genes regulating core N-linked fucosylation in hepatocellular carcinoma

METHODS

The PubMed biomedical literature database was searched for the association of hepatocellular carcinoma and each of the core fucose-related genes and their protein products. We also queried The Cancer Genome Atlas Liver hepatocellular carcinoma (LIHC) dataset for genetic, epigenetic and gene expression changes for the set of six genes using the tools at cBioportal.

RESULTS

A total of 27 citations involving one or more of the core fucosylation-related genes (FPGT, FUK, FUT8, GMDS, SLC35C1, TSTA3) and hepatocellular carcinoma were identified. The same set of gene symbols was used to query the 371 patients with liver cancer in the LIHC dataset to identify the frequency of mRNA over or under expression, as well as non-synonymous mutations, copy number variation and methylation level. Although all six genes trended to more samples displaying over expression relative to under-expression, it was noted that a number of tumor samples had undergone amplification of the genes of the de novo synthesis pathway, GMDS (27 samples) and TSTA3 (78 samples). In contrast, the other four genes had undergone amplification in 2 or fewer samples.

CONCLUSION

Amplification of genes involved in the de novo pathway for generation of GDP-fucose, GMDS and TSTA3, likely contributes to the elevated core fucose observed in hepatocellular carcinoma.