Objective Selecting interventions for patients with solitary hepatocellular carcinoma (HCC) remains a challenge. Despite gross classification being proposed as a potential prognostic predictor, its widespread use has been restricted due to inadequate studies with sufficient patient numbers and the lack of established mechanisms. We sought to investigate the prognostic impacts on patients with HCC of different gross subtypes and assess their corresponding molecular landscapes.
Design A prospective cohort of 400 patients who underwent hepatic resection for solitary HCC was reviewed and analysed and gross classification was assessed. Multiomics analyses were performed on tumours and non-tumour tissues from 49 patients to investigate the mechanisms underlying gross classification. Inverse probability of treatment weight (IPTW) was used to control for confounding factors.
Results Overall 3-year survival rates varied significantly among the four gross subtypes (type I: 91%, type II: 80%, type III: 74.6%, type IV: 38.8%). Type IV was found to be independently associated with poor prognosis in both the entire cohort and the IPTW cohort. The four gross subtypes exhibited three distinct transcriptional modules. Particularly, type IV tumours exhibited increased angiogenesis and immune score as well as decreased metabolic pathways, together with highest frequency of TP53 mutations. Patients with type IV HCC may benefit from adjuvant intra-arterial therapy other than the other three subtypes. Accordingly, a modified trichotomous margin morphological gross classification was established.
Conclusion Different gross types of HCC showed significantly different prognosis and molecular characteristics. Gross classification may aid in development of precise individualised diagnosis and treatment strategies for HCC.
- hepatocellular carcinoma
- surgical oncology
Data availability statement
The datasets analysed during the current study are available from the corresponding author on reasonable request.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
WHAT IS ALREADY KNOWN ON THIS TOPIC
The accuracy of various staging systems for solitary hepatocellular carcinoma (HCC) in guiding treatment decisions and predicting prognosis remains limited.
Gross classification can predict the prognosis of patients undergoing hepatic resection for small solitary HCC, but it is not widely used in clinical practice.
A few stemness and fibrotic stroma marker genes are differently expressed in HCCs with each gross subtype.
WHAT THIS STUDY ADDS
The gross classification is validated as an effective factor in predicting prognosis in patients with HCC using a relatively large prospective cohort and the inverse probability of treatment weight method to eliminate potential confounding effects.
Distinct molecular expression patterns, gene mutations and components of tumour microenvironment (TME) are discovered among the four gross subtypes.
Infiltrative type HCC exhibits the most similarities to intrahepatic cholangiocarcinoma in terms of gross appearance, prognosis and downregulated expression profiles.
Only infiltrative type HCC shows response to adjuvant transcatheter arterial chemoembolisation (TACE).
An easy-to-use modified gross classification system (MMC) for HCC based solely on margin morphology is proposed.
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
Gross classification of HCC is supported by both clinical and molecular rationales, making it a comprehensive and reliable method for stratifying patients. A priority assessment for infiltrative type in treatment decision-making is recommended.
Adjuvant TACE should be performed in infiltrative type HCC following hepatectomy.
The components of immune TME in each gross subtype, which can be evaluated by imaging, may provide evidence to choose the optimal immunotherapy-based strategy for advance HCC.
MMC can improving the clinical applicability of gross classification for HCC.
More attention should be paid on macroscopic features of all tumours, instead of focusing solely on microscopic features.
Hepatocellular carcinoma (HCC) is the most common cancer of hepatobiliary system and remains as one of the leading causes of cancer-related death. Although the prognosis of patients with HCC has improved significantly with rapidly evolving in various therapeutic strategies, recurrence remains as the main challenge in the management of HCC. Local therapies, including surgery, radiofrequency ablation (RFA) and transcatheter arterial chemoembolisation (TACE), considered as standard care for early-stage or intermediate HCC are recommended by many guidelines,1–3 but clinical studies comparing these interventions have not reached consistent conclusions.4 The dilemma of choosing interventions persists for HCC, especially for the cases with single nodule. There is still a lack of a non-invasive and simple indicator of prognosis together with aiding in clinical decision-making in HCC. On the other hand, atezolizumab plus bevacizumab was confirmed to provide survival benefits for patients with advanced HCC.5 To date, no promising biomarkers of response to the first-line therapy have been identified.6 7
There is an old adage in pathological diagnosis that goes ‘macroscopy accounts for 70%, microscopy accounts for 30%’. Typically, the gross pathology of HCC is often simplistically equated with the size, number or texture of the lesions. In 1987, Moriyama’s group classified HCC into five gross subtypes based primarily on tumour shape: single nodular type (type I), single nodular type with extranodular growth (type II), contiguous multinodular type (type III), poor demonstrated type (also named as infiltrative type, type IV) and early HCC type.8 Several limited retrospective studies have suggested a correlation between tumour shape and prognosis in patients with HCC,8–11 even imaging features reflecting HCC gross appearance have been proposed to be predictive of outcomes after RFA, TACE, even treated with lenvatinib.12 13 Despite the utility of the classification system for HCC, its clinical implementation has been limited due to diagnostic challenges associated with certain subtypes and a scarcity of cases representing each subtype that would facilitate robust clinical investigation.14 15 Describing the molecular landscape of HCC subtypes can help explain the rationale behind this classification and expedite its optimisation.
Primarily using a prospective cohort of patients undergoing hepatectomy for different gross subtypes HCC, we sought to investigate the role of tumour morphology in oncological survival. Integrative multiomics analyses were employed to portray the molecular landscapes among distinct HCC gross subtypes and optimise the classification system, culminating in the development of a margin morphology classification (MMC) system. The current study can establish a strong foundation for informed precise and effective personalised therapies of patients with HCC, particularly those with a single nodule when imaging can accurately depict the margin morphology of the tumour.
Materials and methods
The gross classification of resected HCC was separately assessed by one surgeon and one pathologist mainly according to the definition raised by Liver Cancer Study Group of Japan.8 16 Details are as follows: type I: single nodule with distinct margin, usually round with complete tumour envelope; type II: single nodule with extranodular growth, no more than three extranodular points; type III: a unifocal lesion composed of confluent multiple nodules, distinct boundaries among the nodules; type IV: infiltrative nodule, with poor demarcated boundary and especially multiple extranodular points. The gross morphology of multilayered tumour specimens before and after formalin fixation were both included to be evaluated, and we also referred to the imaging features at the largest tumour dimension to ensure the credibility of the gross classification (figure 1A,B).
Study population and tumour samples
A prospective cohort of patients who underwent hepatectomy for HCC from 2017 to 2021 at First Hospital of Jilin University was constructed, and was retrospectively analysed. The inclusion criteria were patients: (1) had HCC which was confirmed by pathological examination of resected specimens; (2) had single HCC nodule detected by CT/MRI imaging and (3) no extrahepatic spread. The exclusion criteria were patients who: (1) had HCC nodule >15 cm; (2) with multiple HCC nodules diagnosed by imaging; (3) with ruptured and recurrent HCC; (4) with palliative surgical resection (R1 or R2 resection); (5) with intrahepatic cholangiocarcinoma (ICC) or other malignant tumours diagnosed pathologically; (6) had perioperative death (≤30 days); (7) had inconsistency in determining the gross classification of HCC between surgeons and pathologists and (8) had missing data on essential variables. Another prospective cohort of patients with advanced HCC who received immunotherapy combined with antiangiogenesis targeted therapy (≥2 cycles) was also enrolled. Residual surgical excision specimens for pathological diagnosis were collected and stored in −80℃ freezer. Samples were consecutively collected from 12 to 13 patients for each gross subtype, resulting in a total of 49 patients with solitary HCCs measuring ≤5 cm in diameter. Total RNA and DNA were isolated from each sample for subsequent whole-genome sequencing and expression analysis. Whole-genome sequencing analysis was conducted on 49 paired tissues from the 49 patients. Additionally, expression analysis was performed on a subset of 69 samples including 39 tumour tissues and 30 were non-tumour tissues from 39 patients (figure 1A). Notably, the collected samples did not include any capsular tissues. The non-tumour tissues used for sequencing were identified as free from contamination with tumour tissues by histology.
A detailed description of other methods used in current study can be found in online supplemental methods.
Strong differences in prognosis across the four gross subtypes of HCC
The patient flow chart is presented in figure 1A. Among the 400 patients undergoing hepatic resection for HCC with single nodule, 280 (70.0%) were hepatitis B virus (HBV)-related individuals. With a median follow-up of 25.5 months, 67 of 396 patients (16.9%) died, and 147 (37.4%) developed recurrence of HCC. A total of 52 (13.0%) individuals had type IV nodules, while 118 (29.5%), 129 (32.3%) and 101 (25.3%) had HCCs belonging to type I, type II and type III nodules, respectively. The comparisons of clinicopathological characteristics and operative variables among patients with different gross types of HCC are noted in online supplemental table 1. In type IV group, the proportions of patients with microvascular invasion (88.4%), macrovascular invasion (48.1%) and satellite nodules (40.4%) and HBV infection (92.3%) were significantly higher than patients in other three groups (all p<0.001). Meanwhile, type IV tumours were found with larger tumour size, poorer differentiation and less complete tumour capsule than other three types HCC (all p<0.001). The proportions of patients with HCC with type IV at BCLC C stage (61.5%) were significantly higher than other subtypes (10.2%, 14.7% and 13.9% for types I, II, III respectively; p<0.001). HCC notably, the significant linear-by-linear association demonstrated that a progressive change in type IV compared with types I, II and III (all p<0.001, online supplemental table 1). Other variables were comparable among different gross types (all p>0.05). The patients with type I HCCs had superior overall survival (OS) and recurrence-free survival (RFS) than patients with type II or III, with type IV demonstrating the worst prognosis (figure 1C,D; p<0.0001). Similar trends were also found in subgroup analyses enrolling patients with ≤5 cm HCC or with HBV-related HCC (online supplemental figure 1).
Distinct transcriptomic profiles among the four gross subtypes of HCC
To gain a better understanding of the underlying mechanisms involved in gross classification and to clarify the rationale for molecular-level gross classification, mRNA, lncRNA and protein expression levels were detected using transcriptomics and proteomics methods in tumour (≤5 cm) and non-tumour samples (figure 2, online supplemental figure 2–5). The results showed that HCCs with type III and types IV had more differentially expressed genes (DEGs) than those in type I or type II group (figure 2A,B). Hierarchical clustering revealed that type II and type III HCCs had similar expression patterns, while type I or type IV HCCs had less expressional overlaps with other types (figure 2C). Furthermore, simplified enrichment analysis of DEGs was performed to cluster functional enrichment results. Type II and type III HCCs were characterised by downregulated immune-related pathways, including leucocyte activation and cytokine production. Type IV HCCs exhibited markedly downregulated metabolic pathways and lipid transport, together with upregulated biogenesis and nucleocytoplasmic transport (figure 2D,E, online supplemental figure 6 and 7), which was also demonstrated at the protein level (online supplemental figure 5). The relative gene expression patterns across gross subtypes were calculated by GSVA, while ESTIMATE and xCell were used to evaluate the abundance of immune and stromal cells. The four gross subtypes had three main expression modules, whereas type II and type III HCCs shared similar expression characteristics (figure 2F–I and online supplemental figure 8). HCCs with type IV exhibited significantly higher stromal scores (figure 2J) and a higher median immune score than HCCs with other types (figure 2K). A high presence of CD4+ and CD8+ T effector memory cells was found in type I, and the increased abundance of fibroblasts in type IV was calculated (online supplemental figure 9).
Type IV as an independent risk factor for HCC with special transcriptional characteristics
The worse prognosis remained in type IV HCC compared with other three gross subtypes after using stabilised inverse probability of treatment weight (IPTW) to balance the distribution of baseline characteristics due to possible selection bias among patients (figure 3A–D, online supplemental table 2). Subsequently, the univariate and multivariate Cox-regression analyses of risk factors for OS and RFS were performed for patients undergoing hepatic resection of HCC with type IV and other subtypes. Multivariate analyses demonstrated that type IV was independently associated with poorer OS after hepatectomy (HR 2.50, 95% CI 1.37 to 4.56, p=0.003), as well as poorer RFS (HR 1.65, 95% CI 1.05 to 2.60, p=0.031), respectively (figure 3E,F). Although some common pathways, such as spindle and condensed chromosome-related cellular components, were enriched across all subtypes, type IV HCC exhibited a higher number of gene counts and more invasion-related pathways, including focal adhesion and cell leading edge (figure 3G). Compared with non-tumour tissues, HCC with type IV was showed with enhanced angiogenesis, epithelial mesenchymal transition (EMT) and protein secretion using GSEA, which were shown with opposite expression patterns in other groups (figure 3H–L and online supplemental figure 10). Given that the gross appearance of HCCs with type IV is indistinguishable from ICC and both have shortened patient prognosis, we investigated whether they shared comparable expression patterns. Enrichment analysis revealed that HCCs highly expressed cell cycle associated genes than ICC, while type IV not other three types HCC displayed analogous downregulated pathways to ICC (figure 3M and online supplemental figure 11). Additionally, according to Bagaev’s classification for cancers, only type IV HCC can be identified as having an ‘immune-enriched and fibrotic’ tumour microenvironment (TME), while types II and III are characterised by a depleted TME, and type I tumours have either immune-enriched or depleted TME features (figure 3N).17 Histological examination revealed an increase in infiltrated leukocytes (CD45+ cells) in HCCs with type IV. However, within these infiltrating cells, there appeared to be a predominance of macrophages (CD68+ cells) rather than CD8+ T cells (figure 4). We further confirmed the presence of differences in TME-related factors, including VEGFA, VEGFC, CD34, PDGFRA, PDGFRB, TGFβ1 and HGF across different gross subtypes. Higher expression levels of VEGFA and TGFβ1 were identified in the type IV group compared with other groups. Additionally, the CD34-positive blood vessels in type IV HCCs appeared dilated and longer compared with HCCs with other subtypes (online supplemental figure 12).
Non-tumour tissues of patients with type IV HCC are characterised by immunosuppressive microenvironment
When comparing the microenvironment score between tumour tissues and non-tumour tissues, the differences in type IV HCC are comparatively less pronounced than those observed in other subtypes (online supplemental figure 13A). The dendrogram analysis revealed that 75% of non-tumour tissues were clustered into a single group, while 75% (6/8) of the non-tumour tissues in the type IV group could be grouped into one cluster (online supplemental figure 13B). Notably, we identified 308 downregulated DEGs in type IV non-tumour samples compared with the samples of the other three types, the majority of which were enriched in immune-related pathways, such as antigen processing and presentation, and lymphocyte activation (online supplemental figure 13C,D). Furthermore, our GSEA results, using hallmark gene sets, demonstrated a decrease in inflammatory response, complement, IL2-Stat5 signalling and EMT, along with upregulation of fatty acid metabolism and oxidative phosphorylation (online supplemental figure 13E–J).
Genomic landscape of HCC with different gross classification
Analysis of gene mutations in the cohort of 49 patients revealed TTN (45%), TP53 (43%), MUC16 (22%), CTNNB1 (20%) and TRPA1 (20%) as the top five frequently mutated genes. The observed mutation frequencies of TP53 were significantly higher in type IV HCC samples (66.7%) than those in type I, type II and type III HCC samples (30.8%, 25.0% and 50.0%, respectively) (figure 5A). Notably, TP53 was identified as the main significantly mutated gene across all subtypes except for type I HCC (figure 5B). No significant differences in tumour mutation burden harboured by HCC were detected among different groups (figure 5C). Mutational signatures of each sample are shown in figure 5D,E, and no significant differences in signatures were observed among different subtypes. Three mutational signatures across all samples, namely, COSMIC signature 5, signature 16 and signature 22, were calculated and identified (figure 5F). Previous studies have suggested that tumours exhibiting an aristolochic acid signature (signature 22) have a higher number of infiltrating CD8+ T cells,18 the abundances of immune cells with/without specific signature fails to reach statistical significance (figure 5G). Structure variations were shown in figure 5H.
Altered HCC gross classification associated gene expression contributed by copy number variation
The probabilistic scoring of amplification and deletion alterations in chromosome regions is depicted in figure 6A,B. Notably, no significant enrichment of alterations was detected in any of the gross subtypes. To identify altered genes in each chromosome region, copy number variation (CNV)-mRNA correlation analysis was employed, as shown in figure 6C. Following an overlap with upregulated DEGs of tumours with different gross subtypes, type IV HCCs were found to be enriched with 12 gene amplifications, as illustrated in figure 6D. It is noteworthy that UCK2 and DSN1 were amplified across all HCC subtypes. Additionally, four, three, nine and two genes were identified as deletion genes in types I–IV HCCs, respectively (figure 6E). Furthermore, the correlated expressions of CR1 and FOSB were detected across all gross classifications.
Functional modules and hub genes of HCC associated with different gross subtypes
Using weighted correlation network analysis (WGCNA), a total of 12 gene modules were determined, as illustrated in figure 7A. The MEpink module exhibited a positive correlation with gross subtype (correlation 0.43, p=0.007), microvascular invasion (correlation 0.55, p=0.0003), macrovascular invasion (correlation 0.33, p=0.04), satellite (correlation 0.42, p=0.007), recurrence (correlation 0.39, p=0.02) and death (correlation 0.36, p=0.02). Conversely, the MEred module was negatively associated with microvascular invasion (correlation −0.49, p=0.001), AFP value (correlation −0.47, p=0.003), and death (correlation −0.33, p=0.04) (figure 7B). The findings presented in figure 7C demonstrated that the gene modules significantly related to type IV HCC were MEpink (correlation 0.67, p=3e-06), MEbrown (correlation 0.56, p=2e-04) and MEturquoise (correlation 0.45, p=0.004). Moreover, hierarchical clustering analysis revealed that these three modules exhibited more prominent coexpression patterns among themselves compared with other modules (online supplemental figure 14A). As such, we identified the genes in these three modules together with CNV-mRNA correlated genes as the ‘malignant gene set (MGS)’. GO and KEGG analysis showed that the MEpink module containing 194 genes were enriched in GO pathways including macrophage activation and cell development, and KEGG pathways such as tyrosine metabolism (figure 7D,E and online supplemental figure 14). MEbrown module genes were enriched in TME remodelling associated pathways such as extracellular matrix organisation, regulation of angiogenesis and response to TGF-β (figure 7D,F). The enrichment analysis of MEturquoise module revealed that these genes were primarily associated with immune-related pathways (figure 7D,G). MEred module was characterised as metabolic process related module (figure 7D,H). The genes were included in MEgreenyellow module and enriched into response to metal ion (online supplemental figure 14B). The module eigengene (ME) values of pink, brown and turquoise were higher in type IV HCC than the other groups, while the ME value of red was lower. Turquoise and brown modules had decreased ME values in type II and type III HCCs (figure 7F,G). The hub genes of specific in each gross type identified by the module membership (MM) value and gene significance (GS) value were shown in (figure 7I–L, online supplemental figure 14 and online supplemental table 3).
Prognostic prediction model for HCC based on MGS
To enhance the practicality and generalisability of this gross classification, we integrated MGS and developed a prognostic prediction model for HCC based on data from TCGA-LIHC cohort. LASSO regression analysis was employed, resulting in a development of a HCC prognosis risk score model consisting of 18 genes (termed as ‘HEPAR-18’). Importantly, our model demonstrated excellent performance, with area under the curve (AUC) values of 0.810, 0.764 and 0.780 for predicting survival probability at 1-year, 2-year and 3-year intervals, respectively (figure 8A–C). Each variable in this model can significantly distinguish the survival benefits of individuals (figure 8D and online supplemental figure 15). A subsequent multivariate Cox regression analysis revealed that all 18 genes in the ‘HEPAR-18’ model were independent risk factors for HCC and the coefficient values were shown in figure 8E. The Risk Score were evaluated following the gene expression, subsequently the distribution of Risk Score and corresponding survival event of each individual in the training cohort were described. As demonstrated in figure 8F,G, the patients with high Risk Score were demonstrated with inferior survival benefits. Meanwhile, the AUC values of HEPAR-18 and TNM stage were 0.776 vs 0.674 (figure 8H). The samples at TNM III/IV were presented with higher Risk Score (figure 8I). Same tendencies are robustly confirmed in the validation sets: International Cancer Genome Consortium cohort (figure 8J–M), HBV-related HCC cohort (figure 8N–Q) and the current cohort (figure 8R,S). Furthermore, the samples with larger tumour size and at BCLC intermediate/advanced stage were presented with higher Risk Score using the HBV-related HCC cohort (figure 8T,U).
Gross classification would be a valuable tool guiding HCC therapy
The histology of type IV HCC was distinguished by high stroma content (eg, CAFs, vessels) and increased immune cell infiltration, especially monocytes and macrophages (figure 4, figure 9A), while relatively infiltrated-depleted TME was observed in the HCC with type II or type III, which was consistent with the transcriptomic findings discovered by GSVA and xCell (figure 2E–K). To further explore the clinical value of HCC gross subtypes and their molecular characteristics, we initially compared the efficacy of adjuvant TACE treatment in patients undergoing partial hepatectomy for different gross subtypes of HCC. As successful TACE efficacy requires that the tumour has sufficient nutrient vessels, we found that postoperative TACE improved OS in type IV HCC as expected (475.3 days vs 640.1 days, p=0.065). Given that the potential presence of more severe disease in the TACE group, our analysis confirmed significantly beneficial effects of adjuvant TACE therapy on OS after correcting confounding factors using stabilised IPTW (p=0.039) (figure 9B). The therapeutic effects on HCC were not significant in the Type I/II/III group before or after adjustment by IPTW (figure 9C).
We also encompassed another prospective cohort containing 79 patients with advanced HCC, who received immunotherapy combined with antiangiogenesis targeted therapy as their primary treatment modality. The gross subtype of each tumour was evaluated by two radiologists using enhanced CT/MRI, with blinding to the outcomes of patients. The majority of patients (81%) in this cohort harboured at least one type IV nodule with a higher incidence of macrovascular invasion (59.4% vs 46.7% for other three types). The objective response rate (ORR) in patients with type IV HCC (30%) were seemly slightly higher compared with the patients with type II/III (20%), although shortened survival benefits were still observed in patients with type IV HCC (figure 9D). In a small subset of patients (12 patients) with type IV and non-type IV HCC nodules simultaneously, better response to immune checkpoint blockade therapy combined with antiangiogenesis targeted therapy in type IV than other type was found in some cases (figure 9E). Type III HCCs can be further classified into two subtypes: type IIIA, characterised by confluent multiple nodules without any extranodular parts, typically exhibiting a regular round smooth margin; and type IIIB, characterised by confluent multiple nodules with extranodular parts, usually with irregular margin (online supplemental figure 16A). These two subtypes of type III HCCs have distinct patterns of gene expression and differences in survival outcomes (online supplemental figure 16B–D).
Clinical staging for HCC based on imaging findings is extensively used, such BCLC staging system, which often requires clinicians to make treatment decisions without a histological diagnosis.17–19 However, current guidelines based on these staging systems often lack specific management recommendations for patients with solitary HCC, especially those less than 5 cm in size (BCLC stage 0 /A).2 20 To date, no consistent conclusions have been reached regarding the preference of treatment strategies such as surgical or ablative treatment in these patients.21–23
Gross classification, endorsed by LCGSJ and Korean Liver Cancer Association, shows promise to stratify HCC management, but it has yet to gain widespread adoption in clinical settings.8 24 In this study, we constructed a relatively large prospective cohort of 400 surgical patients with solitary HCC with different gross subtypes. Types II and III HCC were demonstrated to have worse OS and RFS than type I, while type IV HCC presented the worst prognosis. Although similar trends were observed in previous studies, the differences were not found to be statistically significant due to lack of samples with specific subtype of HCC or multiple nodules.10 11 25 Particularly, type IV HCC is very rare in Japanese patient cohorts, who more commonly suffer from HCV.14 In contrast, it appeared to be relatively frequent among Korean and Chinese HCC patient cohorts, which are predominantly associated with HBV infection.10 11 As such, the distribution of HCC gross classification may be influenced by different aetiologies.
Consistent with previous studies, our study confirmed that histopathological features including vascular invasion, micrometastasis or tumour size increased in severity in the order of type I, type II/III to type IV.10 26 Therefore, in non-type I HCC (especially type IV), hepatic resection which can achieve more extensive excision than other curative-intent modalities, would be preferred.14 The differences in patient’s outcomes between type II/III and type I HCCs were no longer significant, those between type IV and the other three gross subtypes remained significant after IPTW adjustment. These findings suggest that type IV HCC may involve undiscovered mechanisms contributing to its aggressive features.
In current study, HCC tumours demonstrated downregulated metabolic pathways, with the most significant alterations observed in type IV. This subtype displayed a similar transcriptional profile to the most malignant molecular subtype S-Pf defined in a previous study.27 We further revealed that lipid metabolism enriched red module was significantly negatively associated with type IV, recurrence and death by WGCNA analysis. To date, only a few genes involved in different metabolism pathways (eg, ALDOB, GSTZ1 and CPT2) have been identified as tumour suppressors in HCC.28–30 Our findings support that reprogramming the intrinsic metabolic capacity of cancer has the potential to provide therapeutic opportunities.31 Deciphering the metabolic features of distinct cellular components within TME would be critical to achieve specific elimination of cancer cells while preserving immune cell activity simultaneously.32 33
Furthermore, we delineated the organisation of TME for each gross subtype HCC spatially and transcriptomically: Type IV HCC is characterised by fibrotic or immune-enriched fibrotic TME, while type II/III has features of depleted TME. Type I HCC is distinguished by depleted or immune-enriched TME. Conversely, the previous study suggests HCC with depleted TME subtype calculated only based on RNA-seq data exhibited inferior prognosis than other TME subtypes,34 which may result from the limitations of predicting components based only on gene expression.
Of note, type IV HCC exhibits similar macroscopic features (greyish white, tough and poorly defined) and postoperative prognosis to those of ICC (median OS: 634 days) in our centre. Type IV HCC presented with the highest frequency of TP53 mutations (66.7%), and with the greatest similarity to ICC in terms of downregulated gene enrichment pathways across the four gross subtypes, possessing molecular features reminiscent of the recently identified ICC-like subtype,35 which is partially identified with significant downregulated bile acid metabolic pathways and TP53 mutation. TP53 mutations can induce the dedifferentiation of mature hepatocytes into progenitor-like cells and potentially contribute to the formation of ICC-like HCC.36 While bile acid is reported to contribute to HCC oncogenesis in the context of Mst1/2 or Sirt5 knockout genomic background,37 38 the role of bile acid metabolism in driving HCC progression warrants further exploration. On the other hand, we found the somatic mutations in type IV HCC are distinct from those observed in ICC. For instance, FGFR2 fusion mutation which was prevalent in ICC is not detected in type IV HCC.39 In future research, it is crucial to prioritise the exploration of the relationship and distinctions between the tumour components in combined HCC-ICC and type IV HCC.
Our analysis first revealed a type IV strongly correlated ‘pink module’, which also exhibited the strongest positive correlation with aggressive behaviours of HCC by WGCNA, suggesting its crucial role in driving HCC progression. Four alcohol dehydrogenase genes ADH1A, ADH1B, ADH1C and ADH6 are included in it, while ADH1A was already shown to be a prognostic marker for HCC and negatively correlated with cell proliferation.27 40 Together with the higher abundance of macrophages detected histologically in type IV HCC, three hub genes of pink module, TLR7, GPR34 and TREM2 are primarily enriched in macrophages (http://www.proteinatlas.org/),41 suggesting the pivotal involvement of macrophages in promoting invasive and metastasis. TLR7 is regarded as a molecular marker for M2-type tumour associated macrophages. The potential therapeutic value of activating TLR7 in HCC and ICC has been investigated in clinical trials (NCT04338685).42 The increase of Trem2+ macrophages after TACE treatment suppresses recruitment of CD8+ T cells to the tumour lesion, potentially contributing to HCC recurrence.43 Further investigation into targeting diverse populations of macrophages in advanced HCC is necessary.44 Collectively, we believe that gross classification may aid in assessing the TME, especially in immunophenotyping TME throughout HCC cancer stages. Moreover, devising strategies specifically targeting cell groups tailored to distinct gross subtypes would be promising after obtaining a deeper understanding of immune clusters in the future.
Previous studies have shown that adjuvant TACE did not improve the prognosis of patients undergoing hepatectomy for large HCC.45 Intriguingly, our study found that adjuvant TACE prolonged the prognosis of patients in type IV subgroup, but not for patients with the other three types when compared with hepatectomy alone. Actually, more occult micrometastases outside the surgical region can be detected by digital subtraction angiography in patients with type IV HCC. In addition, we found the satellite nodules in type IV exhibit a highly vascularised TME histologically, which may enhance their responsiveness to TACE. Collectively, adjuvant TACE should be considered for patients with type IV HCC. However, TACE alone would not be recommended for treating non-type I HCC.13
A previous study showed that in unresectable patients, type 4 HCC (similar to type IV here) were more likely to achieve an overall response to antiangiogenesis drug lenvatinib.46 The majority of patients in our non-resectable HCC cohort were found to have type IV nodules. Despite attaining higher ORR to the combination of antiangiogenesis targeted therapy and immunotherapy in type IV HCCs, their prognosis remained worse compared with those only with other three subtypes of nodules. Including more patients would help to clarify the issue. Immunotherapy together with locoregional treatment (eg, TACE) might be propitious in type IV HCC, which has been shown to be effective for advance HCC in a single-arm phase II study.47
To promote the clinical application of gross classification, we recommend evaluating HCC gross classification based on fresh tissue specimens with a priority assessment for type IV in treatment decision-making (figure 9F). Considering the similarities in pathology, molecular features, and prognosis between type II and type III HCC, we propose a modified gross classification system MMC for HCC based solely on margin morphology. The MMC system includes: MMC-I (smooth type), nodules with smooth near-rounded margins including gross classification type I and type IIIA; MMC-II (extranodular growth type), nodules with extranodular margins comprising ≤50% of the tumour circumference, or ≤3 directions, including type II and type IIIB; and MMC-III (infiltrative type), infiltrative nodules with irregular margins comprising >50% of the tumour circumference, or >3 directions, including type IV (figure 9G–H). In clinical practice, priority should be given to assessing whether the HCC nodule belongs to MMC-III.
In conclusion, our study reveals significant differences in molecular and pathological characteristics as well as prognosis among different gross subtypes of HCC (figure 10). These findings provide a biological basis and clinical rationale for the development of personalised and precise treatment plans for HCC. Furthermore, we propose that gross classification, which can be easily obtained by radiological examinations, can serve as a foundation for refined stratified management of HCC.
Data availability statement
The datasets analysed during the current study are available from the corresponding author on reasonable request.
Patient consent for publication
This study was conducted in accordance with the Declaration of Helsinki and the Ethical Guidelines for Clinical Studies, and approved by the Institutional Review Boards of First Hospital of Jilin University (No: 2020-675-1).
We would like to express our sincere gratitude to Professor Run-Sheng Chen’s team for their valuable assistance and support during the research. We wish to acknowledge the technical assistance provided by Professor Shu-Li Xie and Dr. Ming-Qian Li. We want to appreciate Dr. Meng-Han Gao for language polishing and logical organisation. We would like to thank Jia-Ao Yu, Shu-Xuan Li, Shi-Fei Song, Kai Kou and Xing Lv for their hard work in collecting samples. We thank the Department of Biobank, Division of Clinical Research, First Hospital of Jilin University for constructing tissue microarray. We also would like to appreciate Professor Yan-Hua Wu and Professor Jing Jiang for statistical supports. Some published data used in this study was generated by the TCGA Research Network: https://www.cancer.gov/tcga. A few figures were generated on BioRender.com.
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
- Data supplement 1
- Data supplement 2
- Data supplement 3
- Data supplement 4
- Data supplement 5
- Data supplement 6
- Data supplement 7
- Data supplement 8
- Data supplement 9
- Data supplement 10
- Data supplement 11
- Data supplement 12
- Data supplement 13
- Data supplement 14
- Data supplement 15
- Data supplement 16
- Data supplement 17
- Data supplement 18
- Data supplement 19
- Data supplement 20
ZF and MJ contributed equally.
Contributors Conception: GL, ZF and GW; study design: ZF, MJ, LZ and NW; administrative support: GL, GW, LQ and HZ; data collection and acquisition: ZF, MJ, LZ, CW, HW, PZ, FW, XD, XSun, WQ, MW, XShi, JY, CJ, JZ, WC, BH, KC, YC, WM, YH, YL, GW, TL and JQ; data analysis: ZF, MJ, ML and YH; manuscript preparation: ZF, MJ and ML; critical revision: TY and GL; guarantor: GL; final approval of manuscript: all authors.
Funding This study was supported by the Natural Science Foundation of Jilin Province (No. YDZJ202201ZYTS014) and Natural Science Foundation of China (No. 81602059, 82241223 and U20A20360), Foundation of Health Commission of Jilin Province (No. 2018J052).
Competing interests None declared.
Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.