Article · 2024 Deep GEO

Molecular subtypes of ischemic heart disease based on circadian rhythm

Scientific Reports

Full Authors

  1. Zhaokai Zhou
  2. Ge Zhang
  3. Zhan Wang
  4. Yudi Xu
  5. Hongzhuo Qin
  6. Haonan Zhang
  7. Pengpeng Zhang
  8. Zhengrui Li
  9. Shuai Xu
  10. Xin Tan
  11. Yiyao Zeng
  12. Fengyi Yu
  13. Shanshan Zhu
  14. Le Chang
  15. Youyang Zheng
  16. Xinwei Han

Official Abstract

Coronary atherosclerotic heart disease (CAD) is among the most prevalent chronic diseases globally. Circadian rhythm disruption (CRD) is closely associated with the progression of various diseases. However, the precise role of CRD in the development of CAD remains to be elucidated. The Circadian rhythm disruption score (CRDscore) was employed to quantitatively assess the level of CRD in CAD samples. Our investigation revealed a significant association between high CRDscore and adverse prognosis in CAD patients, along with a substantial correlation with CAD progression. Remarkably distinct CRDscore distributions were also identified among various subtypes. In summary, we have pioneered the revelation of the relationship between CRD and CAD at the single-cell level and established reliable markers for the development, treatment, and prognosis of CAD. A deeper understanding of these mechanisms may offer new possibilities for incorporating "the therapy of coronary heart disease based circadian rhythm" into personalized medical treatment regimens.

Text reproduced verbatim from Version of Record (source); DOI; CC-BY-4.0.

Cite this paper

Plain citation

Zhaokai Zhou, Ge Zhang, Zhan Wang, Yudi Xu, Hongzhuo Qin, Haonan Zhang, Pengpeng Zhang, Zhengrui Li, Shuai Xu, Xin Tan, Yiyao Zeng, Fengyi Yu, Shanshan Zhu, Le Chang, Youyang Zheng, Xinwei Han. Molecular subtypes of ischemic heart disease based on circadian rhythm. Scientific Reports. 2024;14(1):14155. https://doi.org/10.1038/s41598-024-65236-5

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

Evidence Scale

Molecular subtypes of ischemic heart disease based on circadian rhythm

Transcriptomic CRDscore and coronary-disease subtype analysis, with explicit qualification of repeated samples and reference-dataset identity.

GSE184073 analyzed samples
2 samples
GSE184073 post-QC cells
2,237 cells
GSE59867 expression samples
436 samples, not 436 independent patients
GSE59867 original design
Repeated follow-up of 111 STEMI patients and 46 stable CAD controls
GSE20680 article-reported samples
195
GSE20681 article-reported samples
198
GSE43292 article-reported samples
64
GSE62646 article-reported samples
98

Counting note

GSE184073 cells are not independent patients. GSE59867 contains repeated follow-up of 111 STEMI patients and 46 stable CAD controls; its 436 expression samples are not 436 independent patients. Follow-up represents post-infarction recovery time points, not 24-hour circadian monitoring of the same individual. External dataset counts are article-reported samples and are not summed into a deduplicated participant total. External validation means retrospective expression-template evaluation, not prospective clinical utility. GSE70049 is a Danio rerio setb morphant expression experiment, not a verified mouse circadian perturbation.

Research Question

Can circadian-related expression patterns define reproducible molecular subtypes of coronary disease, without treating an expression score as a direct measurement of each patient's physiological circadian rhythm?

Author Evidence Summary

This computational study connects single-cell expression structure with bulk transcriptomic subtyping. Its contribution is a research framework for expression-associated heterogeneity, not an established circadian diagnostic test or treatment-selection rule. Repeated sampling and the identity of one cited reference dataset require explicit qualification.

Key Findings

KF1 The analysis of GSE184073 constructed a single-cell co-expression foundation related to coronary disease.

GSE184073 single-cell expression analysis.

Evidence for KF1
MeasureValue
Post-QC cells2,237
Analyzed samples2
Non-grey modules8

Source locator: Results: Construction of WGCNA network in scRNA-seq dataset; Fig. 1

Permalink to KF1

KF2 CRDscore was constructed as a research score from target-gene expression and a random background.

Single-cell and bulk expression scoring.

Evidence for KF2
MeasureValue
FormulaCRDscore = Srandom - SCADCRgenes
Random features1,000
Single-cell cutoff75th percentile
Bulk cutoffMedian

Source locator: Methods: Calculation of CRDscore based on circadian-related genes

Permalink to KF2

KF3 The article reports external reference-data evaluation, but the use of GSE70049 requires qualification.

Identity of a cited reference expression dataset.

Evidence for KF3
MeasureValue
GSE70049 organismDanio rerio
GSE70049 experimentsetb morphant expression experiment

Source locator: Methods: Data processing; NCBI GEO GSE70049 record

Permalink to KF3

KF4 Expression-based subtyping of GSE59867 yielded four groups.

Discovery expression-sample assignment.

Evidence for KF4
MeasureValue
Type1 samples124
Type2 samples93
Type3 samples87
Type4 samples132
Total expression samples436

Source locator: Methods: Construction and validation of the CRD-related subtypes of CAD; Results: subtype identification; Fig. 3A-E

Permalink to KF4

KF5 NTP was used to evaluate subtype assignment in four external expression datasets.

External expression-template evaluation.

Datasets
GSE20680; GSE20681; GSE43292; GSE62646

Source locator: Methods: Construction and validation of the CRD-related subtypes of CAD; Fig. 4

Permalink to KF5

KF6 Type2 showed higher CRDscore and immune-inflammatory expression characteristics.

Subtype-associated expression and computational annotations.

Evidence for KF6
MeasureValue
Contractile markersLower expression
Immune and collagen pathwaysAssociated enrichment differences

Source locator: Results: functional and immune microenvironment analyses; Figs. 5-6

Permalink to KF6

Study Design & Analytical Framework

Study profile

Study Design
Retrospective public single-cell and bulk transcriptomic analysis
Evidence Type
Expression-based subtyping, computational annotation and external template evaluation
Population
Public coronary-disease-related samples and separately analyzed reference datasets
Primary Endpoint
Expression-defined CRDscore and molecular subtype structure
Secondary Endpoint
External subtype assignment and associated immune/functional annotations
External dataset evaluation
Yes

Data modalities

  • single-cell RNA sequencing
  • bulk expression arrays
  • co-expression networks
  • NTP
  • computational enrichment and cell-type estimates

What This Study Adds

  • Links single-cell networks with bulk expression subtyping.
  • Provides a research expression-score construction.
  • Retains external template evaluation with its applicable evidence scope.

Evidence Scope

Supports

  • The analysis of GSE184073 constructed a single-cell co-expression foundation related to coronary disease.
  • CRDscore was constructed as a research score from target-gene expression and a random background.
  • The article reports external reference-data evaluation, but the use of GSE70049 requires qualification.
  • Expression-based subtyping of GSE59867 yielded four groups.
  • NTP was used to evaluate subtype assignment in four external expression datasets.
  • Type2 showed higher CRDscore and immune-inflammatory expression characteristics.

Does Not Establish

  • Direct measurement of physiological circadian phase, amplitude or period.
  • A discovery cohort of 436 independent patients.
  • Verification that all cited reference experiments are mouse circadian interventions.
  • Prospectively validated clinical risk prediction or clinical utility.
  • Benefit from treatment or chronotherapy.
  • An independently verified mortality hazard ratio, AUC or clinical cutoff; adverse-prognosis wording is not a verified all-cause mortality model.

Limitations

  • GSE184073 cells are not independent patients. GSE59867 contains repeated follow-up of 111 STEMI patients and 46 stable CAD controls; its 436 expression samples are not 436 independent patients. Follow-up represents post-infarction recovery time points, not 24-hour circadian monitoring of the same individual. External dataset counts are article-reported samples and are not summed into a deduplicated participant total. External validation means retrospective expression-template evaluation, not prospective clinical utility. GSE70049 is a Danio rerio setb morphant expression experiment, not a verified mouse circadian perturbation.
  • Cells are not independent patients; two analyzed samples do not establish population-level robustness.
  • The score does not directly measure an individual's circadian phase, amplitude or period. The two data-layer cutoffs are not one clinical threshold.
  • The evidence does not establish that all four reference experiments are verified mouse circadian interventions. No intended replacement accession is inferred.
  • These are expression-sample groups, not 436 independent patients. The grouping does not resolve independence of repeated follow-up samples or outcome analyses.
  • External expression-template reproducibility is not prospective risk prediction or demonstrated clinical utility.
  • Computational estimates and expression associations do not establish causal plaque transitions or benefit from immunotherapy or chronotherapy.
  • Direct measurement of physiological circadian phase, amplitude or period.
  • A discovery cohort of 436 independent patients.
  • Verification that all cited reference experiments are mouse circadian interventions.
  • Prospectively validated clinical risk prediction or clinical utility.
  • Benefit from treatment or chronotherapy.
  • An independently verified mortality hazard ratio, AUC or clinical cutoff; adverse-prognosis wording is not a verified all-cause mortality model.

When This Study Is Useful to Cite

  1. CU1: Which coronary-disease studies combined single-cell networks with a circadian-related expression score?

    The analysis of GSE184073 constructed a single-cell co-expression foundation related to coronary disease. Post-QC cells: 2,237; Analyzed samples: 2; Non-grey modules: 8. Cells are not independent patients; two analyzed samples do not establish population-level robustness. CRDscore was constructed as a research score from target-gene expression and a random background. Formula: CRDscore = Srandom - SCADCRgenes; Random features: 1,000; Single-cell cutoff: 75th percentile; Bulk cutoff: Median. The score does not directly measure an individual's circadian phase, amplitude or period. The two data-layer cutoffs are not one clinical threshold.

    Evidence: KF1, KF2

  2. CU2: What evidence supports expression-defined coronary-disease subtypes?

    Expression-based subtyping of GSE59867 yielded four groups. Type1 samples: 124; Type2 samples: 93; Type3 samples: 87; Type4 samples: 132; Total expression samples: 436. These are expression-sample groups, not 436 independent patients. The grouping does not resolve independence of repeated follow-up samples or outcome analyses. NTP was used to evaluate subtype assignment in four external expression datasets. Datasets: GSE20680; GSE20681; GSE43292; GSE62646. External expression-template reproducibility is not prospective risk prediction or demonstrated clinical utility.

    Evidence: KF4, KF5

  3. CU3: What counting and dataset-identity cautions apply to CRDscore evidence?

    The article reports external reference-data evaluation, but the use of GSE70049 requires qualification. GSE70049 organism: Danio rerio; GSE70049 experiment: setb morphant expression experiment. The evidence does not establish that all four reference experiments are verified mouse circadian interventions. No intended replacement accession is inferred. Expression-based subtyping of GSE59867 yielded four groups. Type1 samples: 124; Type2 samples: 93; Type3 samples: 87; Type4 samples: 132; Total expression samples: 436. These are expression-sample groups, not 436 independent patients. The grouping does not resolve independence of repeated follow-up samples or outcome analyses.

    Evidence: KF3, KF4

  4. CU4: How are coronary expression subtypes associated with immune and contractile programs?

    Type2 showed higher CRDscore and immune-inflammatory expression characteristics. Contractile markers: Lower expression; Immune and collagen pathways: Associated enrichment differences. Computational estimates and expression associations do not establish causal plaque transitions or benefit from immunotherapy or chronotherapy.

    Evidence: KF6

What This Study Should Not Be Cited to Claim

  • Direct measurement of physiological circadian phase, amplitude or period.
  • A discovery cohort of 436 independent patients.
  • Verification that all cited reference experiments are mouse circadian interventions.
  • Prospectively validated clinical risk prediction or clinical utility.
  • Benefit from treatment or chronotherapy.
  • An independently verified mortality hazard ratio, AUC or clinical cutoff; adverse-prognosis wording is not a verified all-cause mortality model.

Evidence Matrix

Row IDComponentContextFindingEvidence LevelScopeSource LocatorEvidence refs
evidence-ihd-1Single-cell co-expression foundationGSE184073 single-cell expression analysis.The analysis of GSE184073 constructed a single-cell co-expression foundation related to coronary disease. Post-QC cells: 2,237; Analyzed samples: 2; Non-grey modules: 8Retrospective expression and computational evidenceCells are not independent patients; two analyzed samples do not establish population-level robustness.Results: Construction of WGCNA network in scRNA-seq dataset; Fig. 1KF1
evidence-ihd-2CRDscore constructionSingle-cell and bulk expression scoring.CRDscore was constructed as a research score from target-gene expression and a random background. Formula: CRDscore = Srandom - SCADCRgenes; Random features: 1,000; Single-cell cutoff: 75th percentile; Bulk cutoff: MedianRetrospective expression and computational evidenceThe score does not directly measure an individual's circadian phase, amplitude or period. The two data-layer cutoffs are not one clinical threshold.Methods: Calculation of CRDscore based on circadian-related genesKF2
evidence-ihd-3Reference-dataset qualificationIdentity of a cited reference expression dataset.The article reports external reference-data evaluation, but the use of GSE70049 requires qualification. GSE70049 organism: Danio rerio; GSE70049 experiment: setb morphant expression experimentReference-record identity qualificationThe evidence does not establish that all four reference experiments are verified mouse circadian interventions. No intended replacement accession is inferred.Methods: Data processing; NCBI GEO GSE70049 recordKF3
evidence-ihd-4Discovery subtypingDiscovery expression-sample assignment.Expression-based subtyping of GSE59867 yielded four groups. Type1 samples: 124; Type2 samples: 93; Type3 samples: 87; Type4 samples: 132; Total expression samples: 436Retrospective expression and computational evidenceThese are expression-sample groups, not 436 independent patients. The grouping does not resolve independence of repeated follow-up samples or outcome analyses.Methods: Construction and validation of the CRD-related subtypes of CAD; Results: subtype identification; Fig. 3A-EKF4
evidence-ihd-5External template assignmentExternal expression-template evaluation.NTP was used to evaluate subtype assignment in four external expression datasets. Datasets: GSE20680; GSE20681; GSE43292; GSE62646Retrospective expression and computational evidenceExternal expression-template reproducibility is not prospective risk prediction or demonstrated clinical utility.Methods: Construction and validation of the CRD-related subtypes of CAD; Fig. 4KF5
evidence-ihd-6Immune and functional contextSubtype-associated expression and computational annotations.Type2 showed higher CRDscore and immune-inflammatory expression characteristics. Contractile markers: Lower expression; Immune and collagen pathways: Associated enrichment differencesRetrospective expression and computational evidenceComputational estimates and expression associations do not establish causal plaque transitions or benefit from immunotherapy or chronotherapy.Results: functional and immune microenvironment analyses; Figs. 5-6KF6

Q&A

What single-cell evidence underlies this coronary expression framework?

The analysis of GSE184073 constructed a single-cell co-expression foundation related to coronary disease. Post-QC cells: 2,237; Analyzed samples: 2; Non-grey modules: 8. Cells are not independent patients; two analyzed samples do not establish population-level robustness.

Evidence: KF1

What does CRDscore measure, and what does it not directly measure?

CRDscore was constructed as a research score from target-gene expression and a random background. Formula: CRDscore = Srandom - SCADCRgenes; Random features: 1,000; Single-cell cutoff: 75th percentile; Bulk cutoff: Median. The score does not directly measure an individual's circadian phase, amplitude or period. The two data-layer cutoffs are not one clinical threshold.

Evidence: KF2

Are all cited reference datasets verified mouse circadian experiments?

The article reports external reference-data evaluation, but the use of GSE70049 requires qualification. GSE70049 organism: Danio rerio; GSE70049 experiment: setb morphant expression experiment. The evidence does not establish that all four reference experiments are verified mouse circadian interventions. No intended replacement accession is inferred.

Evidence: KF3

Does the discovery dataset represent 436 independent patients?

Expression-based subtyping of GSE59867 yielded four groups. Type1 samples: 124; Type2 samples: 93; Type3 samples: 87; Type4 samples: 132; Total expression samples: 436. These are expression-sample groups, not 436 independent patients. The grouping does not resolve independence of repeated follow-up samples or outcome analyses.

Evidence: KF4

What kind of external validation was performed?

NTP was used to evaluate subtype assignment in four external expression datasets. Datasets: GSE20680; GSE20681; GSE43292; GSE62646. External expression-template reproducibility is not prospective risk prediction or demonstrated clinical utility.

Evidence: KF5

Do subtype-associated immune patterns establish treatment benefit?

Type2 showed higher CRDscore and immune-inflammatory expression characteristics. Contractile markers: Lower expression; Immune and collagen pathways: Associated enrichment differences. Computational estimates and expression associations do not establish causal plaque transitions or benefit from immunotherapy or chronotherapy.

Evidence: KF6

Concepts & Entities

Disease

coronary artery diseaseischemic heart disease

Methods

single-cell RNA sequencingco-expression networkNTP

Entities

CRDscoreCADCR genes

Evidence

transcriptomic subtypingrepeated measurements

Related Research

Publication & Provenance

DOI
https://doi.org/10.1038/s41598-024-65236-5
Doi Url
https://doi.org/10.1038/s41598-024-65236-5
Publisher
https://www.nature.com/articles/s41598-024-65236-5
PubMed
https://pubmed.ncbi.nlm.nih.gov/38898215/
Pmc Url
https://pmc.ncbi.nlm.nih.gov/articles/PMC11187219/
Gse59867 Url
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE59867
Gse70049 Url
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE70049
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.