Circadian rhythm disorders represent an abstract concept lacking standardized quantitative metrics. Existing circadian indicators, including traditional rhythm parameters and a limited set of clock gene or physiological biomarkers, are insufficient to robustly capture steady-state endogenous circadian homeostasis in complex disease contexts, thereby constraining quantitative assessment of circadian disruption and limiting its translational applicability. Chronic circadian rhythm disruption is associated with various diseases, including metabolic disorders and malignancies. However, the mechanisms by which circadian disruption influences tumor microenvironment formation and colorectal cancer progression remain incompletely understood. This study employs systems biology analysis to decipher the molecular characteristics of circadian rhythm disruption in colorectal cancer progression. We analyzed single-cell RNA sequencing data from 13 CRC tissue samples and 12 normal mucosal samples, combined with 3733 samples from 34 public batch RNA, microarray, and single-cell RNA sequencing cohorts. We developed and validated the ClockProCRC system, which detects and quantifies intrinsic circadian misalignment in CRC. The ClockProCRC score elucidates how circadian misalignment drives CRC progression trajectories, shapes clinical phenotypes, regulates disease manifestations, and reshapes the tumor microenvironment. SYNE1 gene was identified as a key mediator of circadian misalignment, promoting tumorigenesis by driving epithelial-like phenotypic conversion and demonstrating therapeutic potential in colorectal cancer management. This study establishes a foundation for integrating rhythmic information into clinical practice and advances circadian biology research in the field of CRC.
Introduction Genome‐wide association studies (GWAS) have identified over 200 germline risk loci for colorectal cancer (CRC), yet the causal variants and genes behind most GWAS signals remain unknown and the link between inherited risk and tumor outcome is largely unexplored. Connecting germline single nucleotide polymorphisms (SNPs) to gene expression through expression quantitative trait loci (eQTL) and to clinical outcome is needed to interpret this inherited risk. Methods We combined CRC GWAS summary data (73,149 cases and 112,467 controls of European ancestry) with GTEx v8 eQTL from colon sigmoid, colon transverse, small intestine terminal ileum, and whole blood, integrating causal transcriptome‐wide association study (cTWAS) with SuSiE fine‐mapping, Bayesian colocalization, MAGMA SNP‐to‐gene analysis, and AlphaGenome variant‐effect prediction. Prioritized proteins were assessed by western blot. Four experimentally selected regulatory variants were genotyped in HCT116, SW480, RKO, and NCM460 cells; genotype‐protein associations were tested across cell‐line means, and cis‐regulatory effects were assessed by allele‐specific expression (ASE) and reference‐versus‐alternate dual‐luciferase assays. Prognostic relevance was evaluated in TCGA colorectal tumors. Results cTWAS identified seven genes with posterior inclusion probability (PIP) > 0.50, and colocalization across 78 gene‐tissue pairs revealed 20 associations with PP.H4 > 0.80. Four genes showed convergent evidence: SMAD9 (PIP = 0.916, PP.H4 = 0.981), MAP3K2 (PIP = 0.762, PP.H4 = 0.827), FADS1 (PIP = 0.632, PP.H4 = 0.942), and ACTR1B (PIP = 0.566, PP.H4 = 0.994). All four proteins were reduced in CRC cells. Alt‐allele dosage was inversely associated with SMAD9 (Pearson r = −0.985, BH‐adjusted p = 0.029) and FADS1 protein abundance (r = −0.995, BH‐adjusted p = 0.020), but not with MAP3K2 or ACTR1B. Reporter and ASE assays detected the clearest allele‐specific effects at SMAD9 and FADS1, whereas MAP3K2 and ACTR1B were null in the tested systems. Lower SMAD9 expression was nominally associated with poorer survival (log − rank p = 0.022 to 0.050), but these associations did not survive Benjamini–Hochberg correction across 12 tests. Conclusions Integrating statistical genetics, regulatory prediction, and locus‐directed experiments prioritized SMAD9, MAP3K2, FADS1, and ACTR1B as CRC susceptibility genes. Functional evidence was strongest and most directionally coherent for SMAD9, demonstrated allele‐specific but context‐dependent regulation at FADS1, and placed experimental bounds on the proposed MAP3K2 and ACTR1B mechanisms.
Chengguang Hu, Guang Yang, Han Xiong et al.· Human Mutation· 0 citations
We use cookies to run the site and, with your consent, for analytics and to show ads.
See our Cookie Policy.
