Article · 2025 Deep GEO

Integrated-omics analysis defines subtypes of hepatocellular carcinoma based on circadian rhythm

Journal of Integrative Medicine

Full Authors

  1. Xiao-jie Li
  2. Le Chang
  3. Yang Mi
  4. Ge Zhang
  5. Shan-shan Zhu
  6. Yue-xiao Zhang
  7. Hao-yu Wang
  8. Yi-shuang Lu
  9. Ye-xuan Ping
  10. Peng-yuan Zheng
  11. Xia Xue

Cite this paper

Plain citation

Xiao-jie Li, Le Chang, Yang Mi, Ge Zhang, Shan-shan Zhu, Yue-xiao Zhang, Hao-yu Wang, Yi-shuang Lu, Ye-xuan Ping, Peng-yuan Zheng, Xia Xue. Integrated-omics analysis defines subtypes of hepatocellular carcinoma based on circadian rhythm. Journal of Integrative Medicine. 2025;23(4):445–456. https://doi.org/10.1016/j.joim.2025.06.003

Author-controlled citation; verify against the publisher version of record.

Evidence Scale

Integrated-omics analysis defines subtypes of hepatocellular carcinoma based on circadian rhythm

HCCcrds-related subtype analysis separates expression, retrospective outcome associations and external template assessment.

TCGA-LIHC
369 samples in the reported analysis
GSE149614 analyzed cells
11,383 cells
External datasets
5 datasets with different analytical roles

Counting note

Cell counts, candidate genes, HCC cohort samples and separate circadian reference experiments are different units. No complete unique-participant total is inferred. The five external datasets have different analytical uses: the four Fig. 5A and four Fig. 5B sets overlap but are not identical and must not be added as eight independent cohorts. External validation means retrospective cross-dataset expression/classification evaluation, not prospective clinical utility or direct physiological rhythm measurement in HCC patients.

Research Question

Can a circadian-related expression framework distinguish HCC subtypes with reproducible molecular features and retrospective outcome associations?

Author Evidence Summary

This study provides a transcriptomic framework for HCC heterogeneity. The relevant evidence is subtype association and cross-dataset evaluation, not demonstrated clinical benefit from score-guided treatment. Reference circadian experiments and HCC patient cohorts serve different analytical roles.

Key Findings

KF1 Single-cell modules supplied the HCCcrds gene set.

GSE149614 single-cell module analysis.

Evidence for KF1
MeasureValue
Analyzed cells11,383
Genes from modules M2/M4/M5562

Source locator: Author-deposited full text: Results 3.1; Fig. 2

Permalink to KF1

KF2 The score was evaluated in separate clock, feeding and sleep-reference data.

Circadian reference experiments analyzed separately from HCC cohorts.

Reference datasets
GSE10045; GSE150381; GSE39445

Source locator: Author-deposited full text: Results 3.2; Fig. 3

Permalink to KF2

KF3 TCGA-LIHC analysis identified CS-H, CS-L and CS-M, with poorer outcomes associated with CS-H.

Retrospective TCGA-LIHC subtype and outcome analysis.

Evidence for KF3
MeasureValue
TCGA-LIHC n369
SubtypesCS-H; CS-L; CS-M

Source locator: Author-deposited full text: Results 3.3-3.4; Fig. 4

Permalink to KF3

KF4 External subtype and survival analyses used overlapping but different cohort sets.

External analyses reported in Fig. 5A and Fig. 5B.

Evidence for KF4
MeasureValue
Fig. 5A datasetsICGC_LIRI; GSE14520; GSE54236; GSE104580
Fig. 5B datasetsICGC_LIRI; GSE14520; GSE54236; GSE27150

Source locator: Author-deposited full text: Results 3.5; Fig. 5 legend

Permalink to KF4

KF5 Subtypes differed in proliferation/metabolic programs and inferred immune characteristics.

Expression-based functional and immune analyses.

Analyses
Expression enrichment; CIBERSORT; IPS

Source locator: Author-deposited full text: Results 3.6-3.7; Fig. 6

Permalink to KF5

KF6 NNMT and MARCKSL1 were among subtype-associated candidate genes.

Candidate-gene expression and network prioritization.

Evidence
Expression correlation and network prioritization

Source locator: Author-deposited full text: Results 3.8; Table 2; final PubMed abstract

Permalink to KF6

Study Design & Analytical Framework

Study profile

Study Design
Public single-cell and bulk transcriptomic integration with retrospective subtype evaluation
Evidence Type
Human molecular association and cross-dataset classification analysis
Population
HCC expression cohorts, with separate circadian reference experiments
Primary Endpoint
HCCcrds-associated molecular subtype structure
Secondary Endpoint
Retrospective outcome associations and computational biological annotations
External dataset evaluation
Yes

Data modalities

  • single-cell RNA sequencing
  • bulk transcriptomics
  • co-expression analysis
  • LRA clustering/NTP
  • immune deconvolution

What This Study Adds

  • Provides an HCC expression-subtyping framework.
  • Distinguishes external subtype assignment from outcome analysis.
  • Identifies candidate markers and functional annotations within the stated evidence limits.

Evidence Scope

Supports

  • Single-cell modules supplied the HCCcrds gene set.
  • The score was evaluated in separate clock, feeding and sleep-reference data.
  • TCGA-LIHC analysis identified CS-H, CS-L and CS-M, with poorer outcomes associated with CS-H.
  • External subtype and survival analyses used overlapping but different cohort sets.
  • Subtypes differed in proliferation/metabolic programs and inferred immune characteristics.
  • NNMT and MARCKSL1 were among subtype-associated candidate genes.

Does Not Establish

  • Prospective diagnostic or prognostic clinical utility.
  • Direct individual measurement of physiological circadian phase or amplitude.
  • Circadian disruption as a demonstrated cause of cancer.
  • Clinical benefit from score-guided treatment.
  • IPS as measured response in patients receiving immunotherapy.
  • Therapeutic validation of the candidate genes.
  • Completed page-by-page equivalence checking between the author-deposited text and final typeset version.
  • Unprovided AUCs, hazard ratios, participant totals or additional omics modalities.

Limitations

  • Cell counts, candidate genes, HCC cohort samples and separate circadian reference experiments are different units. No complete unique-participant total is inferred. The five external datasets have different analytical uses: the four Fig. 5A and four Fig. 5B sets overlap but are not identical and must not be added as eight independent cohorts. External validation means retrospective cross-dataset expression/classification evaluation, not prospective clinical utility or direct physiological rhythm measurement in HCC patients.
  • Cells and candidate genes are not patient counts or demonstrated causal regulators.
  • Reference-experiment evaluation does not mean HCC patients underwent longitudinal physiological circadian measurements.
  • Association does not establish that circadian disruption causes liver cancer or constitute a validated clinical decision rule.
  • Do not claim that all five cohorts validated survival, or add the two four-cohort analyses as eight independent cohorts.
  • IPS and immune deconvolution are not an observed immunotherapy-response trial.
  • This evidence synthesis does not claim independent causal intervention or therapeutic validation of these genes.
  • Prospective diagnostic or prognostic clinical utility.
  • Direct individual measurement of physiological circadian phase or amplitude.
  • Circadian disruption as a demonstrated cause of cancer.
  • Clinical benefit from score-guided treatment.
  • IPS as measured response in patients receiving immunotherapy.
  • Therapeutic validation of the candidate genes.
  • Completed page-by-page equivalence checking between the author-deposited text and final typeset version.
  • Unprovided AUCs, hazard ratios, participant totals or additional omics modalities.
  • Bibliographic identity was checked against the final PubMed/publisher record. Detailed evidence locators refer to the author-deposited full text; final typeset figure-level equivalence was not independently re-established.

When This Study Is Useful to Cite

  1. CU1: Which HCC studies derived circadian-related subtypes from transcriptomic data?

    Single-cell modules supplied the HCCcrds gene set. Analyzed cells: 11,383; Genes from modules M2/M4/M5: 562. Cells and candidate genes are not patient counts or demonstrated causal regulators. TCGA-LIHC analysis identified CS-H, CS-L and CS-M, with poorer outcomes associated with CS-H. TCGA-LIHC n: 369; Subtypes: CS-H; CS-L; CS-M. Association does not establish that circadian disruption causes liver cancer or constitute a validated clinical decision rule.

    Evidence: KF1, KF3

  2. CU2: How was a tumor-associated expression score evaluated in external circadian reference data?

    The score was evaluated in separate clock, feeding and sleep-reference data. Reference datasets: GSE10045; GSE150381; GSE39445. Reference-experiment evaluation does not mean HCC patients underwent longitudinal physiological circadian measurements.

    Evidence: KF2

  3. CU3: What cross-cohort evidence supports HCC circadian-related subtype classification?

    External subtype and survival analyses used overlapping but different cohort sets. Fig. 5A datasets: ICGC_LIRI; GSE14520; GSE54236; GSE104580; Fig. 5B datasets: ICGC_LIRI; GSE14520; GSE54236; GSE27150. Do not claim that all five cohorts validated survival, or add the two four-cohort analyses as eight independent cohorts.

    Evidence: KF4

  4. CU4: Which immune and candidate-marker associations characterize HCCcrds subtypes?

    Subtypes differed in proliferation/metabolic programs and inferred immune characteristics. Analyses: Expression enrichment; CIBERSORT; IPS. IPS and immune deconvolution are not an observed immunotherapy-response trial. NNMT and MARCKSL1 were among subtype-associated candidate genes. Evidence: Expression correlation and network prioritization. This evidence synthesis does not claim independent causal intervention or therapeutic validation of these genes.

    Evidence: KF5, KF6

What This Study Should Not Be Cited to Claim

  • Prospective diagnostic or prognostic clinical utility.
  • Direct individual measurement of physiological circadian phase or amplitude.
  • Circadian disruption as a demonstrated cause of cancer.
  • Clinical benefit from score-guided treatment.
  • IPS as measured response in patients receiving immunotherapy.
  • Therapeutic validation of the candidate genes.
  • Completed page-by-page equivalence checking between the author-deposited text and final typeset version.
  • Unprovided AUCs, hazard ratios, participant totals or additional omics modalities.

Evidence Matrix

Row IDComponentContextFindingEvidence LevelScopeSource LocatorEvidence refs
evidence-hcc-1Score derivationGSE149614 single-cell module analysis.Single-cell modules supplied the HCCcrds gene set. Analyzed cells: 11,383; Genes from modules M2/M4/M5: 562Retrospective expression and computational evidenceCells and candidate genes are not patient counts or demonstrated causal regulators.Author-deposited full text: Results 3.1; Fig. 2KF1
evidence-hcc-2Circadian reference evaluationCircadian reference experiments analyzed separately from HCC cohorts.The score was evaluated in separate clock, feeding and sleep-reference data. Reference datasets: GSE10045; GSE150381; GSE39445Retrospective expression and computational evidenceReference-experiment evaluation does not mean HCC patients underwent longitudinal physiological circadian measurements.Author-deposited full text: Results 3.2; Fig. 3KF2
evidence-hcc-3Clinical subtype associationRetrospective TCGA-LIHC subtype and outcome analysis.TCGA-LIHC analysis identified CS-H, CS-L and CS-M, with poorer outcomes associated with CS-H. TCGA-LIHC n: 369; Subtypes: CS-H; CS-L; CS-MRetrospective expression and computational evidenceAssociation does not establish that circadian disruption causes liver cancer or constitute a validated clinical decision rule.Author-deposited full text: Results 3.3-3.4; Fig. 4KF3
evidence-hcc-4Different external evaluation setsExternal analyses reported in Fig. 5A and Fig. 5B.External subtype and survival analyses used overlapping but different cohort sets. Fig. 5A datasets: ICGC_LIRI; GSE14520; GSE54236; GSE104580; Fig. 5B datasets: ICGC_LIRI; GSE14520; GSE54236; GSE27150Retrospective expression and computational evidenceDo not claim that all five cohorts validated survival, or add the two four-cohort analyses as eight independent cohorts.Author-deposited full text: Results 3.5; Fig. 5 legendKF4
evidence-hcc-5Functional and immune associationsExpression-based functional and immune analyses.Subtypes differed in proliferation/metabolic programs and inferred immune characteristics. Analyses: Expression enrichment; CIBERSORT; IPSRetrospective expression and computational evidenceIPS and immune deconvolution are not an observed immunotherapy-response trial.Author-deposited full text: Results 3.6-3.7; Fig. 6KF5
evidence-hcc-6Candidate marker interpretationCandidate-gene expression and network prioritization.NNMT and MARCKSL1 were among subtype-associated candidate genes. Evidence: Expression correlation and network prioritizationRetrospective expression and computational evidenceThis evidence synthesis does not claim independent causal intervention or therapeutic validation of these genes.Author-deposited full text: Results 3.8; Table 2; final PubMed abstractKF6

Q&A

How was the HCCcrds gene framework derived?

Single-cell modules supplied the HCCcrds gene set. Analyzed cells: 11,383; Genes from modules M2/M4/M5: 562. Cells and candidate genes are not patient counts or demonstrated causal regulators.

Evidence: KF1

What role did the circadian reference datasets play?

The score was evaluated in separate clock, feeding and sleep-reference data. Reference datasets: GSE10045; GSE150381; GSE39445. Reference-experiment evaluation does not mean HCC patients underwent longitudinal physiological circadian measurements.

Evidence: KF2

Which HCC subtypes were associated with clinical outcomes?

TCGA-LIHC analysis identified CS-H, CS-L and CS-M, with poorer outcomes associated with CS-H. TCGA-LIHC n: 369; Subtypes: CS-H; CS-L; CS-M. Association does not establish that circadian disruption causes liver cancer or constitute a validated clinical decision rule.

Evidence: KF3

Were external classification and survival analyses performed in identical cohorts?

External subtype and survival analyses used overlapping but different cohort sets. Fig. 5A datasets: ICGC_LIRI; GSE14520; GSE54236; GSE104580; Fig. 5B datasets: ICGC_LIRI; GSE14520; GSE54236; GSE27150. Do not claim that all five cohorts validated survival, or add the two four-cohort analyses as eight independent cohorts.

Evidence: KF4

Do the immune analyses demonstrate immunotherapy benefit?

Subtypes differed in proliferation/metabolic programs and inferred immune characteristics. Analyses: Expression enrichment; CIBERSORT; IPS. IPS and immune deconvolution are not an observed immunotherapy-response trial.

Evidence: KF5

What evidence supports the nominated molecular markers?

NNMT and MARCKSL1 were among subtype-associated candidate genes. Evidence: Expression correlation and network prioritization. This evidence synthesis does not claim independent causal intervention or therapeutic validation of these genes.

Evidence: KF6

Concepts & Entities

Disease

hepatocellular carcinoma

Entities

HCCcrdsNNMTMARCKSL1

Methods

transcriptomic subtypingLRA clusteringNTPimmune deconvolution

Evidence

retrospective outcome association

Related Research

Publication & Provenance

DOI
https://doi.org/10.1016/j.joim.2025.06.003
Doi Url
https://doi.org/10.1016/j.joim.2025.06.003
Publisher
https://www.sciencedirect.com/science/article/abs/pii/S2095496425000834
PubMed
https://pubmed.ncbi.nlm.nih.gov/40581550/
Author Deposited Full Text
https://www.researchgate.net/publication/392697303_Integrated-omics_analysis_defines_subtypes_of_hepatocellular_carcinoma_based_on_circadian_rhythm
Source Version
Bibliographic identity was checked against the final PubMed/publisher record. Detailed evidence locators refer to the author-deposited full text; final typeset figure-level equivalence was not independently re-established.
Review Basis
Evidence synthesis from the cited main text, figure legends and designated primary records. No independent reanalysis of participant-level data. Unavailable supplementary or figure-only details are not reconstructed.
Author-controlled evidence synthesis, distinct from the published abstract and not a replacement for the version of record. Source-specific counting and interpretation limits are retained.

Study-level evidence table

Study-level author-maintained evidence summary. Not participant-level data and not a ready-to-pool meta-analysis dataset. For formal citation, use the original article DOI.

Download study-level evidence CSV (UTF-8)

Formula safety: cells starting with =, +, -, @ (after leading whitespace), or tab/CR/LF are prefixed with an apostrophe. Treat all cells as text; remove only that safety prefix when recovering the original string.

Verify publication identity

These records identify the same original work; database inclusion is not evidence of scientific validity or clinical benefit.