Single-cell landscape of the ecosystem in early-relapse hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) has high relapse and low 5-year survival rates. Single-cell profiling in relapsed HCC may aid in the design of effective anticancer therapies, including immunotherapies. We profiled the transcriptomes of ∼17,000 cells from 18 primary or early-relapse HCC cases. Early-relapse tumors have reduced levels of regulatory T cells, increased dendritic cells (DCs), and increased infiltrated CD8+ T cells, compared with primary tumors, in two independent cohorts. Remarkably, CD8+ T cells in recurrent tumors overexpressed KLRB1 (CD161) and displayed an innate-like low cytotoxic state, with low clonal expansion, unlike the classical exhausted state observed in primary HCC. The enrichment of these cells was associated with a worse prognosis. Differential gene expression and interaction analyses revealed potential immune evasion mechanisms in recurrent tumor cells that dampen DC antigen presentation and recruit innate-like CD8+ T cells. Our comprehensive picture of the HCC ecosystem provides deeper insights into immune evasion mechanisms associated with tumor relapse.

Glycogen accumulation and phase separation drives liver tumor initiation.

Glucose consumption is generally increased in tumor cells to support tumor growth. Interestingly, we report that glycogen accumulation is a key initiating oncogenic event during liver malignant transformation. We found that glucose-6-phosphatase (G6PC) catalyzing the last step of glycogenolysis is frequently downregulated to augment glucose storage in pre-malignant cells. Accumulated glycogen undergoes liquid-liquid phase separation, which results in the assembly of the Laforin-Mst1/2 complex and consequently sequesters Hippo kinases Mst1/2 in glycogen liquid droplets to relieve their inhibition on Yap. Moreover, G6PC or another glycogenolysis enzyme-liver glycogen phosphorylase (PYGL) deficiency in both human and mice results in glycogen storage disease along with liver enlargement and tumorigenesis in a Yap-dependent manner. Consistently, elimination of glycogen accumulation abrogates liver growth and cancer incidence, whereas increasing glycogen storage accelerates tumorigenesis. Thus, we concluded that cancer-initiating cells adapt a glycogen storing mode, which blocks Hippo signaling through glycogen phase separation to augment tumor incidence.

Onco-fetal Reprogramming of Endothelial Cells Drives Immunosuppressive Macrophages in Hepatocellular Carcinoma.

We employed scRNA sequencing to extensively characterize the cellular landscape of human liver from development to disease. Analysis of ∼212,000 cells representing human fetal, hepatocellular carcinoma (HCC), and mouse liver revealed remarkable fetal-like reprogramming of the tumor microenvironment. Specifically, the HCC ecosystem displayed features reminiscent of fetal development, including re-emergence of fetal-associated endothelial cells (PLVAP/VEGFR2) and fetal-like (FOLR2) tumor-associated macrophages. In a cross-species comparative analysis, we discovered remarkable similarity between mouse embryonic, fetal-liver, and tumor macrophages. Spatial transcriptomics further revealed a shared onco-fetal ecosystem between fetal liver and HCC. Furthermore, gene regulatory analysis, spatial transcriptomics, and in vitro functional assays implicated VEGF and NOTCH signaling in maintaining onco-fetal ecosystem. Taken together, we report a shared immunosuppressive onco-fetal ecosystem in fetal liver and HCC. Our results unravel a previously unexplored onco-fetal reprogramming of the tumor ecosystem, provide novel targets for therapeutic interventions in HCC, and open avenues for identifying similar paradigms in other cancers and disease.

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.

Phase Separation of a PKA Regulatory Subunit Controls cAMP Compartmentation and Oncogenic Signaling.

The fidelity of intracellular signaling hinges on the organization of dynamic activity architectures. Spatial compartmentation was first proposed over 30 years ago to explain how diverse G protein-coupled receptors achieve specificity despite converging on a ubiquitous messenger, cyclic adenosine monophosphate (cAMP). However, the mechanisms responsible for spatially constraining this diffusible messenger remain elusive. Here, we reveal that the type I regulatory subunit of cAMP-dependent protein kinase (PKA), RIα, undergoes liquid-liquid phase separation (LLPS) as a function of cAMP signaling to form biomolecular condensates enriched in cAMP and PKA activity, critical for effective cAMP compartmentation. We further show that a PKA fusion oncoprotein associated with an atypical liver cancer potently blocks RIα LLPS and induces aberrant cAMP signaling. Loss of RIα LLPS in normal cells increases cell proliferation and induces cell transformation. Our work reveals LLPS as a principal organizer of signaling compartments and highlights the pathological consequences of dysregulating this activity architecture.

Integrated Proteogenomic Characterization of HBV-Related Hepatocellular Carcinoma.

We performed the first proteogenomic characterization of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC) using paired tumor and adjacent liver tissues from 159 patients. Integrated proteogenomic analyses revealed consistency and discordance among multi-omics, activation status of key signaling pathways, and liver-specific metabolic reprogramming in HBV-related HCC. Proteomic profiling identified three subgroups associated with clinical and molecular attributes including patient survival, tumor thrombus, genetic profile, and the liver-specific proteome. These proteomic subgroups have distinct features in metabolic reprogramming, microenvironment dysregulation, cell proliferation, and potential therapeutics. Two prognostic biomarkers, PYCR2 and ADH1A, related to proteomic subgrouping and involved in HCC metabolic reprogramming, were identified. CTNNB1 and TP53 mutation-associated signaling and metabolic profiles were revealed, among which mutated CTNNB1-associated ALDOA phosphorylation was validated to promote glycolysis and cell proliferation. Our study provides a valuable resource that significantly expands the knowledge of HBV-related HCC and may eventually benefit clinical practice.

Landscape and Dynamics of Single Immune Cells in Hepatocellular Carcinoma.

The immune microenvironment of hepatocellular carcinoma (HCC) is poorly characterized. Combining two single-cell RNA sequencing technologies, we produced transcriptomes of CD45+ immune cells for HCC patients from five immune-relevant sites: tumor, adjacent liver, hepatic lymph node (LN), blood, and ascites. A cluster of LAMP3+ dendritic cells (DCs) appeared to be the mature form of conventional DCs and possessed the potential to migrate from tumors to LNs. LAMP3+ DCs also expressed diverse immune-relevant ligands and exhibited potential to regulate multiple subtypes of lymphocytes. Of the macrophages in tumors that exhibited distinct transcriptional states, tumor-associated macrophages (TAMs) were associated with poor prognosis, and we established the inflammatory role of SLC40A1 and GPNMB in these cells. Further, myeloid and lymphoid cells in ascites were predominantly linked to tumor and blood origins, respectively. The dynamic properties of diverse CD45+ cell types revealed by this study add new dimensions to the immune landscape of HCC.

Comprehensive and Integrative Genomic Characterization of Hepatocellular Carcinoma.

Liver cancer has the second highest worldwide cancer mortality rate and has limited therapeutic options. We analyzed 363 hepatocellular carcinoma (HCC) cases by whole-exome sequencing and DNA copy number analyses, and we analyzed 196 HCC cases by DNA methylation, RNA, miRNA, and proteomic expression also. DNA sequencing and mutation analysis identified significantly mutated genes, including LZTR1, EEF1A1, SF3B1, and SMARCA4. Significant alterations by mutation or downregulation by hypermethylation in genes likely to result in HCC metabolic reprogramming (ALB, APOB, and CPS1) were observed. Integrative molecular HCC subtyping incorporating unsupervised clustering of five data platforms identified three subtypes, one of which was associated with poorer prognosis in three HCC cohorts. Integrated analyses enabled development of a p53 target gene expression signature correlating with poor survival. Potential therapeutic targets for which inhibitors exist include WNT signaling, MDM4, MET, VEGFA, MCL1, IDH1, TERT, and immune checkpoint proteins CTLA-4, PD-1, and PD-L1.

Clustered Mutation Signatures Reveal that Error-Prone DNA Repair Targets Mutations to Active Genes.

Many processes can cause the same nucleotide change in a genome, making the identification of the mechanisms causing mutations a difficult challenge. Here, we show that clustered mutations provide a more precise fingerprint of mutagenic processes. Of nine clustered mutation signatures identified from >1,000 tumor genomes, three relate to variable APOBEC activity and three are associated with tobacco smoking. An additional signature matches the spectrum of translesion DNA polymerase eta (POLH). In lymphoid cells, these mutations target promoters, consistent with AID-initiated somatic hypermutation. In solid tumors, however, they are associated with UV exposure and alcohol consumption and target the H3K36me3 chromatin of active genes in a mismatch repair (MMR)-dependent manner. These regions normally have a low mutation rate because error-free MMR also targets H3K36me3 chromatin. Carcinogens and error-prone repair therefore redistribute mutations to the more important regions of the genome, contributing a substantial mutation load in many tumors, including driver mutations.

OXCT1 functions as a succinyltransferase, contributing to hepatocellular carcinoma via succinylating LACTB.

Metabolic reprogramming is an important feature of cancers that has been closely linked to post-translational protein modification (PTM). Lysine succinylation is a recently identified PTM involved in regulating protein functions, whereas its regulatory mechanism and possible roles in tumor progression remain unclear. Here, we show that OXCT1, an enzyme catalyzing ketone body oxidation, functions as a lysine succinyltransferase to contribute to tumor progression. Mechanistically, we find that OXCT1 functions as a succinyltransferase, with residue G424 essential for this activity. We also identified serine beta-lactamase-like protein (LACTB) as a main target of OXCT1-mediated succinylation. Extensive succinylation of LACTB K284 inhibits its proteolytic activity, resulting in increased mitochondrial membrane potential and respiration, ultimately leading to hepatocellular carcinoma (HCC) progression. In summary, this study establishes lysine succinyltransferase function of OXCT1 and highlights a link between HCC prognosis and LACTB K284 succinylation, suggesting a potentially valuable biomarker and therapeutic target for further development.

Multi-organ Mapping of Cancer Risk.

Cancers are distributed unevenly across the body, but the importance of cell intrinsic factors such as stem cell function in determining organ cancer risk is unknown. Therefore, we used Cre-recombination of conditional lineage tracing, oncogene, and tumor suppressor alleles to define populations of stem and non-stem cells in mouse organs and test their life-long susceptibility to tumorigenesis. We show that tumor incidence is determined by the life-long generative capacity of mutated cells. This relationship held true in the presence of multiple genotypes and regardless of developmental stage, strongly supporting the notion that stem cells dictate organ cancer risk. Using the liver as a model system, we further show that damage-induced activation of stem cell function markedly increases cancer risk. Therefore, we propose that a combination of stem cell mutagenesis and extrinsic factors that enhance the proliferation of these cell populations, creates a "perfect storm" that ultimately determines organ cancer risk. VIDEO ABSTRACT.

Titration of RAS alters senescent state and influences tumour initiation.

Oncogenic RAS-induced senescence (OIS) is an autonomous tumour suppressor mechanism associated with premalignancy1,2. Achieving this phenotype typically requires a high level of oncogenic stress, yet the phenotype provoked by lower oncogenic dosage remains unclear. Here we develop oncogenic RAS dose-escalation models in vitro and in vivo, revealing a RAS dose-driven non-linear continuum of downstream phenotypes. In a hepatocyte OIS model in vivo, ectopic expression of NRAS(G12V) does not induce tumours, in part owing to OIS-driven immune clearance3. Single-cell RNA sequencing analyses reveal distinct hepatocyte clusters with typical OIS or progenitor-like features, corresponding to high and intermediate levels of NRAS(G12V), respectively. When titred down, NRAS(G12V)-expressing hepatocytes become immune resistant and develop tumours. Time-series monitoring at single-cell resolution identifies two distinct tumour types: early-onset aggressive undifferentiated and late-onset differentiated hepatocellular carcinoma. The molecular signature of each mouse tumour type is associated with different progenitor features and enriched in distinct human hepatocellular carcinoma subclasses. Our results define the oncogenic dosage-driven OIS spectrum, reconciling the senescence and tumour initiation phenotypes in early tumorigenesis.

Deep whole-genome analysis of 494 hepatocellular carcinomas.

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

Acetyl-CoA, protein acetylation, and liver cancer.

Using multi-omics approaches, Park et al. show that reduced cellular acetyl-CoA and protein hypoacetylation promote liver cancer growth and dedifferentiation.

Pseudouridine synthases modify human pre-mRNA co-transcriptionally and affect pre-mRNA processing.

Pseudouridine is a modified nucleotide that is prevalent in human mRNAs and is dynamically regulated. Here, we investigate when in their life cycle mRNAs become pseudouridylated to illuminate the potential regulatory functions of endogenous mRNA pseudouridylation. Using single-nucleotide resolution pseudouridine profiling on chromatin-associated RNA from human cells, we identified pseudouridines in nascent pre-mRNA at locations associated with alternatively spliced regions, enriched near splice sites, and overlapping hundreds of binding sites for RNA-binding proteins. In vitro splicing assays establish a direct effect of individual endogenous pre-mRNA pseudouridines on splicing efficiency. We validate hundreds of pre-mRNA sites as direct targets of distinct pseudouridine synthases and show that PUS1, PUS7, and RPUSD4-three pre-mRNA-modifying pseudouridine synthases with tissue-specific expression-control widespread changes in alternative pre-mRNA splicing and 3' end processing. Our results establish a vast potential for cotranscriptional pre-mRNA pseudouridylation to regulate human gene expression via alternative pre-mRNA processing.

Transcription factors TEAD2 and E2A globally repress acetyl-CoA synthesis to promote tumorigenesis.

Acetyl-coenzyme A (acetyl-CoA) plays an important role in metabolism, gene expression, signaling, and other cellular processes via transfer of its acetyl group to proteins and metabolites. However, the synthesis and usage of acetyl-CoA in disease states such as cancer are poorly characterized. Here, we investigated global acetyl-CoA synthesis and protein acetylation in a mouse model and patient samples of hepatocellular carcinoma (HCC). Unexpectedly, we found that acetyl-CoA levels are decreased in HCC due to transcriptional downregulation of all six acetyl-CoA biosynthesis pathways. This led to hypo-acetylation specifically of non-histone proteins, including many enzymes in metabolic pathways. Importantly, repression of acetyl-CoA synthesis promoted oncogenic dedifferentiation and proliferation. Mechanistically, acetyl-CoA synthesis was repressed by the transcription factors TEAD2 and E2A, previously unknown to control acetyl-CoA synthesis. Knockdown of TEAD2 and E2A restored acetyl-CoA levels and inhibited tumor growth. Our findings causally link transcriptional reprogramming of acetyl-CoA metabolism, dedifferentiation, and cancer.

EPDR1 promotes PD-L1 expression and tumor immune evasion by inhibiting TRIM21-dependent ubiquitylation of IkappaB kinase-ß.

While immune checkpoint blockade (ICB) has shown promise for clinical cancer therapy, its efficacy has only been observed in a limited subset of patients and the underlying mechanisms regulating innate and acquired resistance to ICB of tumor cells remain poorly understood. Here, we identified ependymin-related protein 1 (EPDR1) as an important tumor-intrinsic regulator of PD-L1 expression and tumor immune evasion. Aberrant expression of EPDR1 in hepatocellular carcinoma is associated with immunosuppression. Mechanistically, EPDR1 binds to E3 ligase TRIM21 and disrupts its interaction with IkappaB kinase-b, suppressing its ubiquitylation and autophagosomal degradation and enhancing NF-κB-mediated transcriptional activation of PD-L1. Further, we validated through a mouse liver cancer model that EPDR1 mediates exhaustion of CD8+ T cells and promotes tumor progression. In addition, we observed a positive correlation between EPDR1 and PD-L1 expression in both human and mouse liver cancer samples. Collectively, our study reveals a previously unappreciated role of EPDR1 in orchestrating tumor immune evasion and cancer progression.

A non-canonical repressor function of JUN restrains YAP activity and liver cancer growth.

Yes-associated protein (YAP) and its homolog, transcriptional coactivator with PDZ-binding motif (TAZ), are the main transcriptional downstream effectors of the Hippo pathway. Decreased Hippo pathway activity leads to nuclear translocation of YAP/TAZ where they interact with TEAD transcription factors to induce target gene expression. Unrestrained YAP/TAZ activity can lead to excessive growth and tumor formation in a short time, underscoring the evolutionary need for tight control of these two transcriptional coactivators. Here, we report that the AP-1 component JUN acts as specific repressor of YAP/TAZ at joint target sites to decrease YAP/TAZ activity. This function of JUN is independent of its heterodimeric AP-1 partner FOS and the canonical AP-1 function. Since expression of JUN is itself induced by YAP/TAZ, our work identifies a JUN-dependent negative feedback loop that buffers YAP/TAZ activity at joint genomic sites. This negative feedback loop gets disrupted in liver cancer to unlock the full oncogenic potential of YAP/TAZ. Our results thus demonstrate an additional layer of control for the interplay of YAP/TAZ and AP-1.

Endogenous retrotransposition activates oncogenic pathways in hepatocellular carcinoma.

LINE-1 (L1) retrotransposons are mobile genetic elements comprising ~17% of the human genome. New L1 insertions can profoundly alter gene function and cause disease, though their significance in cancer remains unclear. Here, we applied enhanced retrotransposon capture sequencing (RC-seq) to 19 hepatocellular carcinoma (HCC) genomes and elucidated two archetypal L1-mediated mechanisms enabling tumorigenesis. In the first example, 4/19 (21.1%) donors presented germline retrotransposition events in the tumor suppressor mutated in colorectal cancers (MCC). MCC expression was ablated in each case, enabling oncogenic ß-catenin/Wnt signaling. In the second example, suppression of tumorigenicity 18 (ST18) was activated by a tumor-specific L1 insertion. Experimental assays confirmed that the L1 interrupted a negative feedback loop by blocking ST18 repression of its enhancer. ST18 was also frequently amplified in HCC nodules from Mdr2(-/-) mice, supporting its assignment as a candidate liver oncogene. These proof-of-principle results substantiate L1-mediated retrotransposition as an important etiological factor in HCC.

Identification of liver cancer progenitors whose malignant progression depends on autocrine IL-6 signaling.

Hepatocellular carcinoma (HCC) is a slowly developing malignancy postulated to evolve from premalignant lesions in chronically damaged livers. However, it was never established that premalignant lesions actually contain tumor progenitors that give rise to cancer. Here, we describe isolation and characterization of HCC progenitor cells (HcPCs) from different mouse HCC models. Unlike fully malignant HCC, HcPCs give rise to cancer only when introduced into a liver undergoing chronic damage and compensatory proliferation. Although HcPCs exhibit a similar transcriptomic profile to bipotential hepatobiliary progenitors, the latter do not give rise to tumors. Cells resembling HcPCs reside within dysplastic lesions that appear several months before HCC nodules. Unlike early hepatocarcinogenesis, which depends on paracrine IL-6 production by inflammatory cells, due to upregulation of LIN28 expression, HcPCs had acquired autocrine IL-6 signaling that stimulates their in vivo growth and malignant progression. This may be a general mechanism that drives other IL-6-producing malignancies.