Search

Patients with hepatocellular carcinoma (HCC) have poor long-term survival following curative resection because of the high rate of tumor early recurrence. Little is known about the trajectory of genomic evolution from primary to early-recurrent HCC. In this study, we performed whole-genome sequencing (WGS) on 40 pairs of primary and early-recurrent hepatitis B virus (HBV)-related HCC tumors from patients who received curative resection, and from four patients whose primary and recurrent tumor were extensively sampled. We identified two recurrence patterns: de novo recurrence (18/40), which developed genetically independently of the primary tumor and carried different HCC drivers, and ancestral recurrence (22/40), which was clonally related to the primary tumor and progressed more rapidly than de novo recurrence. We found that the recurrence location was predictive of the recurrence pattern: distant recurrence tended to display the de novo pattern, whereas local recurrence tended to display the ancestral pattern. We then uncovered the evolutionary trajectories based on the subclonal architecture, driver-gene mutations, and mutational processes observed in the primary and recurrent tumors. Multi-region WGS demonstrated spatiotemporal heterogeneity and polyclonal, monophyletic dissemination in HCC ancestral recurrence. In addition, we identified recurrence-specific mutations and copy-number gains in BCL9, leading to WNT/β-catenin signaling activation and an immune-excluded tumor microenvironment, which suggests that BCL9 might serve as a new therapeutic target for recurrent HCC. Collectively, our results allow us to view with unprecedented clarity the genomic evolution during HBV-related HCC early recurrence, providing an important molecular foundation for enhanced understanding of HCC with implications for personalized therapy to improve patient survival.

Sarcomatoid hepatocellular carcinoma (sHCC) is a rare type of liver malignancy. Currently, the tumor immune features of sHCC are poorly understood. We recruited 31 patients with resected sHCC for whom tissue samples and complete clinicopathologic and follow-up data were available. To understand the immune infiltration of sHCC, immunohistochemical staining was performed on the resected sHCC samples to compare the expressions of programmed death-1 (PD-1), programmed death-ligand 1 (PD-L1), B7-H3, indoleamine 2,3-dioxygenase (IDO), lymphocyte-activation gene 3 (LAG-3), CD8, FOXP3, and CD68 in tumor and peritumoral tissues. Kaplan-Meier and Cox regression analyses were used to assess the predictive value of immune markers. Sarcomatoid components were characterized with significantly higher expression of PD-L1 and B7-H3 in tumor cells than in conventional HCC components, as well as in peritumoral tissue. Additionally, sarcomatoid components had a higher density of FOXP3+ and LAG-3+ cells and a lower density of CD8+ cells than conventional HCC components or peritumoral tissue. Higher expression of PD-L1 in tumor cells significantly correlated with higher densities of CD8+, PD-1+, and LAG-3+ cells. Increased tumor PD-L1 expression and decreased CD8+ T-cell density were associated with poor overall survival (OS) and disease-free survival (DFS) in patients of sHCC. These findings suggest further characterization on relative mechanism of sHCC immune infiltration may identify therapeutic targets for immunotherapy.

Key Points Question What is the prevalence of KRAS variant subtypes and their association with survival and recurrence in patients with surgically treated intrahepatic cholangiocarcinoma (ICC)? Findings In this cohort study including 1024 patients, a total of 14 different subtypes of KRAS somatic variants affecting 127 patients with ICC (12.4%) were identified, including G12D (43.3%), G12V (19.7%), G12C (7.1%), and G13D (6.3%). G12 KRAS variants but not non-G12 KRAS variants were independently associated with worse overall and disease-free survival, and the G12V KRAS variant was the strongest prognostic determinant for the worst overall and disease-free survival. Meaning This cohort study characterized the distribution of KRAS variant subtypes in a large cohort of patients with ICC and showed an association with patient outcome., This cohort study explores the prognostic association of KRAS variant subtypes with survival and recurrence in patients with intrahepatic cholangiocarcinoma., Importance KRAS variants are associated with tumor progression; however, the prevalence of KRAS variant subtypes and their association with survival and recurrence in patients with intrahepatic cholangiocarcinoma (ICC) after curative resection are largely unknown. Objective To explore the prognostic association of KRAS variant subtypes with survival and recurrence in patients with ICC. Design, Setting, and Participants In this cohort study, patients who underwent curative resection for ICC from January 2009 through December 2016 at a single hospital in China were recruited, and whole-exome sequencing, targeted sequencing, and Sanger sequencing were performed to identify KRAS variants. Kaplan-Meier and log-rank tests were used to compare overall survival (OS) and disease-free survival (DFS). Univariate and multivariate analyses were performed using the Cox proportional hazards regression model. Data were analyzed from April 2020 to January 2021. Interventions Hepatectomy in patients with ICC. Main Outcomes and Measures The association of KRAS variant subtypes with OS and DFS. Results Of 1024 included patients with ICC, 621 (60.6%) were male, and the mean (SD) age was 59.2 (10.2) years. A total of 14 different subtypes of KRAS somatic variants affecting 127 patients (12.4%) were identified. G12D was the most frequent allele in this cohort, accounting for 55 of 127 identified KRAS variants (43.3%), followed by G12V (25 [19.7%]), G12C (9 [7.1%]), and G13D (8 [6.3%]). Compared with patients with wild-type KRAS, patients with variant KRAS were more likely to have high levels of carbohydrate antigen 19-9 (92 of 127 [72.4%] vs 546 of 897 [60.9%]; P = .01) and γ-glutamyltransferase (72 of 127 [56.7%] vs 420 of 897 [46.8%]; P = .04). Multivariable analysis revealed that G12 KRAS variants but not non-G12 KRAS variants were independently associated with worse OS (hazard ratio [HR], 1.69; 95% CI, 1.31-2.18; P

Background & aims Intrahepatic cholangiocarcinoma (ICC) is aggressive and has high rates of relapse, conferring poor long-term survival after curative resection. Little is known about the genomic evolution that occurs during ICC relapse. Approach & results We conducted whole-exome sequencing (WES) of 30 paired primary and relapsed tumors from 10 patients with ICC who received curative resection. We sought to identify frequently altered genes, infer tumor subclonal architectures, and track genomic evolution from primary to relapsed tumors. We examined functional effects and the mechanism of action of SLIT2, a gene specifically mutated in relapsed tumors, on tumor growth and metastasis and the tumor microenvironment. Our results indicated that relapsed ICCs were genetically derived from intrahepatic dissemination of primary tumors. But, they acquired additional mutations while maintaining most drivers, such as TP53 and IDH1. Multiregion sequencing suggested polyclonal seeding of ICC dissemination. Four of 10 relapsed ICCs acquired SLIT2 mutations that were not present in the corresponding primary tumors. Validation in an expanded sample revealed SLIT2 mutations in 2.3% (1/44) of primary ICCs and 29.5% (13/44) of relapsed ICCs. Biofunctional investigations revealed that inactivating mutation of SLIT2 resulted in activation of PI3K-Akt signaling in ICC cells, directly enhanced neutrophil chemotaxis, mediated tumor-associated neutrophil infiltration, and contributed to ICC growth and metastasis. Conclusions We characterized genomic evolution during ICC relapse and identified SLIT2 as a driver of tumor dissemination and tumor-associated neutrophil infiltration.

Background & Aims Early recurrence of hepatocellular carcinoma (HCC) after curative resection is common. However, the association between genetic mechanisms and early HCC recurrence, especially in Chinese patients, remains largely unknown. Methods We performed whole-genome sequencing (49 cases), whole-exome sequencing (18 cases), and deep targeted sequencing (115 cases) on 182 primary HCC samples. Focusing on WNK2, we used Sanger sequencing and qPCR to evaluate all the coding exons and copy numbers of that gene in an additional 554 HCC samples. We also explored the functional effect and mechanism of WNK2 on tumor growth and metastasis. Results We identified 5 genes (WNK2, RUNX1T1, CTNNB1, TSC1, and TP53) harboring somatic mutations that correlated with early tumor recurrence after curative resection in 182 primary HCC samples. Focusing on WNK2, the overall somatic mutation and copy number loss occurred in 5.3% (39/736) and 27.2% (200/736), respectively, of the total 736 HCC samples. Both types of variation were associated with lower WNK2 protein levels, higher rates of early tumor recurrence, and shorter overall survival. Biofunctional investigations revealed a tumor-suppressor role of WNK2: its inactivation led to ERK1/2 signaling activation in HCC cells, tumor-associated macrophage infiltration, and tumor growth and metastasis. Conclusions Our results delineate genomic events that characterize Chinese HCCs and identify WNK2 as a driver of early HCC recurrence after curative resection. Lay summary We applied next-generation sequencing and conducted an in-depth genomic analysis of hepatocellular carcinomas from a Chinese patient cohort. The results delineate the genomic events that characterize hepatocellular carcinomas in Chinese patients and identify WNK2 as a driver associated with early tumor recurrence after curative resection.

Background Long noncoding RNAs (lncRNAs) are functionally associated with cancer development and progression. Although gene copy number variation (CNV) is common in hepatocellular carcinoma (HCC), it is not known how CNV in lncRNAs affects HCC progression and recurrence. We aimed to identify a CNV‐related lncRNA involved in HCC progression and recurrence and illustrate its underlying mechanisms and prognostic value. Methods We analyzed the whole genome sequencing (WGS) data of matched cancerous and noncancerous liver samples from 49 patients with HCC to identify lncRNAs with CNV. The results were validated in another cohort of 238 paired HCC and nontumor samples by TaqMan copy number assay. We preformed Kaplan‐Meier analysis and log‐rank test to identify lncRNA CNV with prognostic value. We conducted loss‐ and gain‐of‐function studies to explore the biological functions of LINC01133 in vitro and in vivo. The competing endogenous RNAs (ceRNAs) mechanism was clarified by microRNA sequencing (miR‐seq), quantitative real‐time PCR (qRT‐PCR), western blot, and dual‐luciferase reporter assays. We confirmed the binding mechanism between lncRNA and protein by RNA pull‐down, RNA immunoprecipitation, qRT‐PCR, and western blot analyses. Results Genomic copy numbers of LINC01133 were increased in HCC, which were positively related with the elevated expression of LINC01133. Increased copy number of LINC01133 predicted the poor prognosis in HCC patients. LINC01133 overexpression in HCC cells promoted proliferation and aggressive phenotypes in vitro, and facilitated tumor growth and lung metastasis in vivo, whereas LINC01133 knockdown had the opposite effects. LINC01133 sponged miR‐199a‐5p, resulting in enhanced expression of SNAI1, which induced epithelial‐to‐mesenchymal transition (EMT) in HCC cells. In addition, LINC01133 interacted with Annexin A2 (ANXA2) to activate the ANXA2/STAT3 signaling pathway. Conclusions LINC01133 promotes HCC progression by sponging miR‐199a‐5p and interacting with ANXA2. LINC01133 CNV gain is predictive of poor prognosis in patients with HCC., Genomic copy numbers of LINC01133 were increased in HCC, which were positively related with its elevated expression, and LINC01133 CNV gain predicted the poor prognosis in HCC patients. LINC01133 sponged miR‐199a‐5p, resulting in enhanced expression of SNAI1, which induced EMT in HCC cells. LINC01133 interacted with ANXA2 to activate ANXA2/STAT3 signaling pathway.

BackgroundTumor-associated neutrophils (TANs) and macrophages (TAMs) can each influence cancer growth and metastasis, but their combined effects in intrahepatic cholangiocarcinoma (ICC) remain unclear.MethodsWe explored the distributions of TANs and TAMs in patient-derived ICC samples by multiplex immunofluorescent staining and tested their separate and combined effects on ICC in vitro and in vivo. We then investigated the mechanistic basis of the effects using PCR array, western blot analysis and ELISA experiments. Finally, we validated our results in a tissue microarray composed of primary tumor tissues from 359 patients with ICC.ResultsThe spatial distributions of TANs and TAMs were correlated with each other in patient-derived ICC samples. Interaction between TANs and TAMs enhanced the proliferation and invasion abilities of ICC cells in vitro and tumor progression in a mouse xenograft model of ICC. TANs and TAMs produced higher levels of oncostatin M and interleukin-11, respectively, in co-culture than in monoculture. Both of those cytokines activated STAT3 signaling in ICC cells. Knockdown of STAT3 abolished the protumor effect of TANs and TAMs on ICC. In tumor samples from patients with ICC, increased TAN and TAM levels were correlated with elevated p-STAT3 expression. All three of those factors were independent predictors of patient outcomes.ConclusionsTANs and TAMs interact to promote ICC progression by activating STAT3.

Recent studies suggested that the immune microenvironment and mutational landscape are associated with the response to immune-based therapy in several types of cancer. The roles of those factors in Chinese HCC remain largely unknown. In this study, we obtained 182 FFPE samples of HCC cohort that were previously subjected to NGS (49 WGS, 18 WES, and 115 targeted sequencing). We performed immunohistochemistry to detect CD3, CD4, CD8, CD57, Foxp3, CD68, CD66b, and PD-L1 expression in the samples. We identified diverse associations between the mutational landscape and the immune microenvironment in the HCC samples. High mutational burden and an aristolochic acid-dominated mutational signature were both correlated with elevated tumoral PD-L1 expression and CD3+ T-cell infiltration and high numbers of CD68+ TAMs and CD66b+ TANs. CD4+ and CD8+ T cells exhibited lower infiltration levels in tumors with mutations in AXIN1/CTNNB1 and in tumors with aflatoxin-dominant mutational signatures. Moreover, tumors with TP53 mutations had less CD8+ T-cell infiltration and more Foxp3+ Treg-cell infiltration than those without TP53 mutations. Kaplan–Meier survival analysis revealed that the presence of CD8+, Foxp3+, CD66b+, or CD68+ immune cells; tumoral PD-L1 expression alone; or the presence of CD8+ or Foxp3+ cells combined with TP53 mutation were predictive of recurrence and poor overall survival after curative resection. In conclusion, the association between the mutational landscape and the immune microenvironment warrants further analysis to determine its impact on patient outcomes to guide personalized immune-based therapy for Chinese patients with HCC.

Background and aims There is growing evidence that single-stranded, circular RNA (circRNA) plays a key role in the development of certain cancers, including hepatocellular carcinoma (HCC). It is less clear, however, what role circRNA plays in HCC metastasis. Approach and results In this study, through circRNA sequencing, we identified a circRNA: circASAP1 (a circRNA derived from exons 2 and 3 of the ASAP1 gene, hsa_circ_0085616), which is associated with pulmonary metastasis after curative resection in patients with HCC. CircASAP1 was overexpressed in HCC cell lines with high metastatic potential and in metastatic HCCs. In vitro, circASAP1 promoted cell proliferation, colony formation, migration, and invasion, and in vivo, it enhanced tumor growth and pulmonary metastasis. Mechanism studies showed that circASAP1 acts as a competing endogenous RNA for microRNA 326 (miR-326) and microRNA 532-5p (miR-532-5p), both of which are tumor suppressors in HCC. We found that mitogen-activated protein kinase (MAPK) 1 and colony stimulating factor (CSF)-1 were direct common targets for microRNA 326 (miR-326) and microRNA 532-5p (miR-532-5p), which were regulated by circASAP1. CircASAP1 promotes HCC cell proliferation and invasion by regulating miR-326/miR-532-5p-MAPK1 signaling and, furthermore, mediates tumor-associated macrophage infiltration by regulating the miR-326/miR-532-5p-CSF-1 pathway. Clinical HCC samples exhibited a positive correlation between circASAP1 expression and levels of CSF-1, MAPK1, and CD68+ tumor-associated macrophages, all of which were predictive of patient outcomes. Conclusion We identified circASAP1 as a key regulator of HCC metastasis that acts on miR-326/miR-532-5p-MAPK1/CSF-1 signaling and serves as a prognostic predictor in patients with HCC.

Plasmacytoid dendritic cells (pDCs) are present in various primary and metastatic human neoplasms; however, their clinical significance in hepatocellular carcinoma (HCC) is unclear. In this study, we investigated the distribution, prognostic value, and potential function of pDCs in HCC patients undergoing curative resection. We performed immunohistochemical analyses of whole tumor sections from 224 patients to assess the expression of BDCA2, CD3, CD4, CD8, Foxp3, granzyme B, IL-17, and CD34. The findings were validated using tissue microarrays from another two independent cohorts totaling 841 HCC patients undergoing curative resection. Our results demonstrated that high numbers of BDCA2+ pDCs within tumors correlated with high alpha-fetoprotein levels, greater vascular invasion, advanced tumor-node-metastasis stage, shorter overall survival, and a higher recurrence rate. However, patient outcomes were not associated with pDCs in peritumoral stromal or nontumor tissues. Furthermore, an increase in intratumoral pDCs was associated with increased intratumoral infiltration of Foxp3+ regulatory T cells and IL-17-producing cells and correlated with tumor vascular density. Univariate and multivariate analyses revealed that the presence of intratumoral pDCs alone or in combination with regulatory T and/or IL-17-producing cells was an independent predictor of time to recurrence and overall survival. In conclusion, our study demonstrated that intratumoral infiltration by pDCs is a novel indicator for poor prognosis in patients with HCC, possibly through the induction of an immune tolerogenic and inflammatory tumor microenvironment comprising regulatory T and IL-17-producing cells. An assessment of the combination of these cells represents a superior predictor of patient outcome.