当前位置: 首页 > news >正文

<span class=“js_title_inner“>结直肠癌五年生存率为何骤降?bulk+单细胞+空间转录组整合揭秘肿瘤微环境“变脸”全过程(附全套复现代码)</span>

为什么结直肠癌的五年生存率在早期接近90%,而一旦进展到IV期,这个数字会断崖式下跌到仅14%左右?

2023年4月,Gastroenterology Report 杂志发表了中山大学附属第六医院兰平教授团队的研究,该研究利用单细胞和空间转录组技术,系统揭示了结直肠癌在进展过程中肿瘤微环境异质性的动态演变规律。

今天我们就来拆解一下这篇文章:Dynamic heterogeneity of colorectal cancer during progression revealed clinical risk-associated cell types and regulations in single-cell resolution and spatial context

研究概述

本研究整合了8个单细胞转录组(scRNA-seq)、3个空间转录组(ST-seq)以及5个批量转录组(bulk RNA-seq)数据集,旨在分析结直肠癌(CRC)从I期进展到IV期过程中,肿瘤细胞、基质细胞和免疫细胞的比例、功能及谱系分化的动态变化。研究鉴定出了与临床风险高度相关的细胞亚型,并发现LIF–LIFR信号通路是晚期结直肠癌的关键调节中心。

实验设计

研究人员收集并处理了来自公共数据库的CRC多组学数据,利用Milo算法检测细胞群丰度在分期进展中的差异,通过Palantir推断分化轨迹,并结合scMetabolism评估细胞代谢状态。随后,利用空间转录组数据验证了这些细胞亚型的空间定位和共定位关系,最后通过80例自主队列样本的qRT-PCR和免疫组化实验,对核心标志物进行了临床验证。

研究结果

Figure 1:研究构建了涵盖结直肠癌I-IV期的单细胞转录图谱,鉴定了包括上皮、基质和免疫细胞在内的八大类细胞。


Figure 2:TM4SF1+恶性上皮细胞在IV期患者中显著增加,这类细胞具有最强的上皮-间质转化能力,且与患者生存期缩短直接相关。


Figure 3:CXCL12+炎症型相关成纤维细胞在IV期肿瘤中富集,其受JUN转录因子调控的特性可能参与了肿瘤纤维化的进展。


Figure 4:CD4+ Treg和Trm细胞在晚期组织中浸润增加,且这些细胞在末端状态表现出极高的代谢水平和TCR信号激活。


Figure 5:CD8+ T细胞主要分化为增殖型和耗竭型两种末端状态,其耗竭程度与细胞的代谢活性呈正相关。


Figure 6:三级淋巴结构(TLS)增强了肿瘤区域的免疫浸润,同时在IV期肿瘤中观察到明显的IgA到IgG的抗体类别转换。


Figure 7:C1QC+巨噬细胞和肥大细胞在IV期显著上调,其代谢重编程过程可能促进了肿瘤的转移。


Figure 8:LIF–LIFR信号通路在晚期肿瘤细胞间高度活跃,临床验证证实LIF高表达是预后不良的独立预测指标。

数据分析

生信分析

本研究涉及scRNA-seq、ST-seq以及bulk RNA-seq三种组学技术。

分析流程

分析首先对单细胞数据进行质控、降维聚类及Harmony批次效应校正;接着运用Milo进行细胞丰度的分期差异分析,Palantir推断发育轨迹,scMetabolism计算代谢通路评分。通过CellChat模拟配体-受体相互作用,结合cell2location对空间转录组点位进行细胞比例反褶积,明确细胞在组织中的物理位置;最后利用CIBERSORTx将单细胞特征映射到bulk转录组中,通过生存分析确认细胞亚型的临床意义。

统计分析

组间显著性差异采用ANOVA分析及Tukey’s多重比较检验;生存曲线分析使用log-rank检验;变量间的相关性通过Pearson相关系数评估,所有统计检验均以P < 0.05为显著性标准。

总结

研究意义

本研究详尽描绘了结直肠癌进展过程中肿瘤微环境的动态重塑图谱,不仅明确了与转移相关的特定恶性细胞亚群,还解析了免疫抑制环境形成的轨迹。研究确立了TM4SF1和LIF作为评估CRC进展和预后的关键生物标志物,为晚期患者的靶向治疗提供了新的潜在路径。

文章复现

这篇文章的原始数据和生信分析代码都公开了,非常全面。

原始数据分别存放在GEO和10x Genomics官方平台,具体编号如下:

  • • 单细胞数据集:GSE161277, GSE132465, GSE144735, GSE146771, GSE164522, GSE178318, GSE188711, GSE200997

  • • Bulk RNA-seq数据集:TCGA-COAD, TCGA-READ, GSE17536, GSE17537, GSE39582

  • • 空间转录组数据集:GSE188711及10x Genomics官网下载链接

代码仓库地址: https://github.com/KeHax/SingleCellCRC-pipeline


推荐阅读

中国银河生信云平台(UseGalaxy.cn)致力于零代码生信分析。平台拥有海量计算资源、3000 多个生信工具和数十条生信流程,并且为用户提供 200G 免费存储空间。进群交流请先加 usegalaxy 为好友。

最佳Galaxy生信云平台教程:从入门到精通(图文版)

转录组分析流程和工具大全(最强总结)

全网最佳WGCNA分析教程,一键完成

一文搞懂GSEA富集分析

一文详解细菌耐药性生信分析:从下机数据到耐药基因鉴定

一文学会从测序数据到构建系统发育树:超全面的详细步骤与软件指南

平台经典生信流程

推荐课程

我们还为进阶用户提供高质量培训课程,欢迎参加:

RNA-seq数据分析实战 | 2026年第2期,开启你的生信学习之旅

http://www.jsqmd.com/news/335184/

相关文章:

  • 可调谐石墨烯超材料吸收体FDTD仿真模拟 【案例内容】该案例提供了一种可调谐石墨烯超材料吸收体
  • 速进学习!AI应用架构师分享法律文本AI理解系统的优化技巧
  • 电气工程专业代码研究:考虑区域供热网络热惯性的AA-CAES电站与热电联产集成的优化调度模型
  • 网页设计过程中常见的误区有哪些?
  • 小白到专家:大模型学习指南与应用场景解析
  • MATLAB代码:计及碳排放交易及多种需求响应的微网/虚拟电厂日前优化调度 关键词
  • [工业自动化-27]:“机器和生物一样,都能通过‘反馈’来自动调整行为,实现目标。”
  • AI与大模型的本质,程序员入门必看!
  • 德语初学感受日志(未完成模板版本)
  • 锂电池充电器用不对称半桥反激变换器电路仿真 两个管子均可实现ZVS 模型包含开环和电压闭环控制
  • 2026豆包AI推广服务商权威评测:五家主流GEO公司深度对比与选型指南 - 品牌2025
  • 稀疏信号代码详解
  • Python异步编程全解析:从asyncio到FastAPI的性能优化实践
  • React19事件调度的设计思路
  • 基于粒子群算法(PSO)优化BP神经网络权值与阈值的实现
  • STM32_GPIO四种输出模式
  • 基于LangChain构建企业级RAG应用的关键架构设计
  • 学习记录260202
  • C++模板编程:泛型代码的终极武器
  • <span class=“js_title_inner“>揭秘LATS:为何这种Agent设计模式让AI决策能力突飞猛进?</span>
  • Flutter 三端应用实战:OpenHarmony “拾光匣”——在匆忙尘世中,为你收藏一缕微光
  • C++内存管理全攻略
  • 基于卷积神经网络(CNN)的图像融合方法详解
  • SQL Backup Master(文件备份软件)
  • Flutter 三端应用实战:OpenHarmony “微光笔记”——在灵感消逝前,为思想点一盏灯
  • MATLAB中LASSO方法的特征矩阵优化与特征选择实现
  • C++核心三要素:封装、实例化与this
  • Flutter 三端应用实战:OpenHarmony “呼吸之境”——在焦虑洪流中,为你筑一座内心的岛屿
  • Recovery Toolbox for Word(Word修复软件)
  • Recovery Toolbox for PDF(PDF文件修复工具)