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.

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.

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.

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.

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.

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.

Ectosomal PKM2 Promotes HCC by Inducing Macrophage Differentiation and Remodeling the Tumor Microenvironment.

Tumor-derived extracellular vesicles are important mediators of cell-to-cell communication during tumorigenesis. Here, we demonstrated that hepatocellular carcinoma (HCC)-derived ectosomes remodel the tumor microenvironment to facilitate HCC progression in an ectosomal PKM2-dependent manner. HCC-derived ectosomal PKM2 induced not only metabolic reprogramming in monocytes but also STAT3 phosphorylation in the nucleus to upregulate differentiation-associated transcription factors, leading to monocyte-to-macrophage differentiation and tumor microenvironment remodeling. In HCC cells, sumoylation of PKM2 induced its plasma membrane targeting and subsequent ectosomal excretion via interactions with ARRDC1. The PKM2-ARRDC1 association in HCC was reinforced by macrophage-secreted cytokines/chemokines in a CCL1-CCR8 axis-dependent manner, further facilitating PKM2 excretion from HCC cells to form a feedforward regulatory loop for tumorigenesis. In the clinic, ectosomal PKM2 was clearly detected in the plasma of HCC patients. This study highlights a mechanism by which ectosomal PKM2 remodels the tumor microenvironment and reveals ectosomal PKM2 as a potential diagnostic marker for HCC.

Two Hippo signaling modules orchestrate liver size and tumorigenesis.

The Hippo pathway is a central regulator of organ size and tumorigenesis and is commonly depicted as a kinase cascade, with an increasing number of regulatory and adaptor proteins linked to its regulation over recent years. Here, we propose that two Hippo signaling modules, MST1/2-SAV1-WWC1-3 (HPO1) and MAP4K1-7-NF2 (HPO2), together regulate the activity of LATS1/2 kinases and YAP/TAZ transcriptional co-activators. In mouse livers, the genetic inactivation of either HPO1 or HPO2 module results in partial activation of YAP/TAZ, bile duct hyperplasia, and hepatocellular carcinoma (HCC). On the contrary, inactivation of both HPO1 and HPO2 modules results in full activation of YAP/TAZ, rapid development of intrahepatic cholangiocarcinoma (iCCA), and early lethality. Interestingly, HPO1 has a predominant role in regulating organ size. HPO1 inactivation causes a homogenous YAP/TAZ activation and cell proliferation across the whole liver, resulting in a proportional and rapid increase in liver size. Thus, this study has reconstructed the order of the Hippo signaling network and suggests that LATS1/2 and YAP/TAZ activities are finetuned by HPO1 and HPO2 modules to cause different cell fates, organ size changes, and tumorigenesis trajectories.

Genomic origin, fragmentomics, and transcriptional properties of long cell-free DNA molecules in human plasma.

Recent studies have revealed an unexplored population of long cell-free DNA (cfDNA) molecules in human plasma using long-read sequencing technologies. However, the biological properties of long cfDNA molecules (>500 bp) remain largely unknown. To this end, we have investigated the origins of long cfDNA molecules from different genomic elements. Analysis of plasma cfDNA using long-read sequencing reveals an uneven distribution of long molecules from across the genome. Long cfDNA molecules show overrepresentation in euchromatic regions of the genome, in sharp contrast to short DNA molecules. We observe a stronger relationship between the abundance of long molecules and mRNA gene expression levels, compared with short molecules (Pearson's r = 0.71 vs. -0.14). Moreover, long and short molecules show distinct fragmentation patterns surrounding CpG sites. Leveraging the cleavage preferences surrounding CpG sites, the combined cleavage ratios of long and short molecules can differentiate patients with hepatocellular carcinoma (HCC) from non-HCC subjects (AUC = 0.87). We also investigated knockout mice in which selected nuclease genes had been inactivated in comparison with wild-type mice. The proportion of long molecules originating from transcription start sites are lower in Dffb-deficient mice but higher in Dnase1l3-deficient mice compared with that of wild-type mice. This work thus provides new insights into the biological properties and potential clinical applications of long cfDNA molecules.

Ectopic lymphoid structures function as microniches for tumor progenitor cells in hepatocellular carcinoma.

Ectopic lymphoid-like structures (ELSs) are often observed in cancer, yet their function is obscure. Although ELSs signify good prognosis in certain malignancies, we found that hepatic ELSs indicated poor prognosis for hepatocellular carcinoma (HCC). We studied an HCC mouse model that displayed abundant ELSs and found that they constituted immunopathological microniches wherein malignant hepatocyte progenitor cells appeared and thrived in a complex cellular and cytokine milieu until gaining self-sufficiency. The egress of progenitor cells and tumor formation were associated with the autocrine production of cytokines previously provided by the niche. ELSs developed via cooperation between the innate immune system and adaptive immune system, an event facilitated by activation of the transcription factor NF-κB and abolished by depletion of T cells. Such aberrant immunological foci might represent new targets for cancer therapy.

DNAJB1-PRKACA fusion protein-regulated LINC00473 promotes tumor growth and alters mitochondrial fitness in fibrolamellar carcinoma.

Fibrolamellar carcinoma (FLC) is a rare liver cancer that disproportionately affects adolescents and young adults. Currently, no standard of care is available and there remains a dire need for new therapeutics. Most patients harbor the fusion oncogene DNAJB1-PRKACA (DP fusion), but clinical inhibitors are not yet developed and it is critical to identify downstream mediators of FLC pathogenesis. Here, we identify long noncoding RNA LINC00473 among the most highly upregulated genes in FLC tumors and determine that it is strongly suppressed by RNAi-mediated inhibition of the DP fusion in FLC tumor epithelial cells. We show by loss- and gain-of-function studies that LINC00473 suppresses apoptosis, increases the expression of FLC marker genes, and promotes FLC growth in cell-based and in vivo disease models. Mechanistically, LINC00473 plays an important role in promoting glycolysis and altering mitochondrial activity. Specifically, LINC00473 knockdown leads to increased spare respiratory capacity, which indicates mitochondrial fitness. Overall, we propose that LINC00473 could be a viable target for this devastating disease.

Liver cancer development driven by the AP-1/c-Jun~Fra-2 dimer through c-Myc.

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related death. HCC incidence is on the rise, while treatment options remain limited. Thus, a better understanding of the molecular pathways involved in HCC development has become a priority to guide future therapies. While previous studies implicated the Activator Protein-1 (AP-1) (Fos/Jun) transcription factor family members c-Fos and c-Jun in HCC formation, the contribution of Fos-related antigens (Fra-) 1 and 2 is unknown. Here, we show that hepatocyte-restricted expression of a single chain c-Jun~Fra-2 protein, which functionally mimics the c-Jun/Fra-2 AP-1 dimer, results in spontaneous HCC formation in c-Jun~Fra-2hep mice. Several hallmarks of human HCC, such as cell cycle dysregulation and the expression of HCC markers are observed in liver tumors arising in c-Jun~Fra-2hep mice. Tumorigenesis occurs in the context of mild inflammation, low-grade fibrosis, and Pparγ-driven dyslipidemia. Subsequent analyses revealed increased expression of c-Myc, evidently under direct regulation by AP-1 through a conserved distal 3' enhancer. Importantly, c-Jun~Fra-2-induced tumors revert upon switching off transgene expression, suggesting oncogene addiction to the c-Jun~Fra-2 transgene. Tumors escaping reversion maintained c-Myc and c-Myc target gene expression, likely due to increased c-Fos. Interfering with c-Myc in established tumors using the Bromodomain and Extra-Terminal motif inhibitor JQ-1 diminished liver tumor growth in c-Jun~Fra-2 mutant mice. Thus, our data establish c-Jun~Fra-2hep mice as a model to study liver tumorigenesis and identify the c-Jun/Fra-2-Myc interaction as a potential target to improve HCC patient stratification and/or therapy.

Study of prognostic splicing factors in cancer using machine learning approaches.

Splicing factors (SFs) are the major RNA-binding proteins (RBPs) and key molecules that regulate the splicing of mRNA molecules through binding to mRNAs. The expression of splicing factors is frequently deregulated in different cancer types, causing the generation of oncogenic proteins involved in cancer hallmarks. In this study, we investigated the genes that encode RNA-binding proteins and identified potential splicing factors that contribute to the aberrant splicing applying a random forest classification model. The result suggested 56 splicing factors were related to the prognosis of 13 cancers, two SF complexes in liver hepatocellular carcinoma, and one SF complex in esophageal carcinoma. Further systematic bioinformatics studies on these cancer prognostic splicing factors and their related alternative splicing events revealed the potential regulations in a cancer-specific manner. Our analysis found high ILF2-ILF3 expression correlates with poor prognosis in LIHC through alternative splicing. These findings emphasize the importance of SFs as potential indicators for prognosis or targets for therapeutic interventions. Their roles in cancer exhibit complexity and are contingent upon the specific context in which they operate. This recognition further underscores the need for a comprehensive understanding and exploration of the role of SFs in different types of cancer, paving the way for their potential utilization in prognostic assessments and the development of targeted therapies.