胶质母细胞瘤 TCGA-GBM — 六百一十七例与二百六十二例影像

六百一十七例 · 十个数据类别 · 二百六十二例术前影像 · 十三例复发

来源 https://portal.gdc.cancer.gov/projects/TCGA-GBM发布时间: 2026-10-29最后更新: 2026-10-29 阅读 7
胶质母细胞瘤 TCGA-GBM — 六百一十七例与二百六十二例影像

信息速览

数据集名称胶质母细胞瘤 TCGA-GBM — 六百一十七例与二百六十二例影像
数据类型多组学,病理全切片,原始影像,元数据,文献引用
规模六百一十七例(GDC 项目口径;2008 年 Nature 主论文的多平台整合口径为 206 例)
接入方式https://portal.gdc.cancer.gov/projects/TCGA-GBM
AI 就绪度

TCGA-GBM 是 TCGA 计划最早启动、也最早交卷的癌种项目。2005 年立项,2008 年在 Nature 发表首篇多平台整合论文,把 206 例胶质母细胞瘤的拷贝数、表达、甲基化与部分外显子序列放在同一张图上,用成规模的测量说明「癌症是一组通路被反复打断」。此后它成为 TCGA 的方法学模板,其余癌种基本沿用了这套结构。

它的对象是成人最常见的恶性脑肿瘤。这个病的特点决定了数据形态:肿瘤几乎不做根治性切除,标本来自减压手术,正常对照只有血液;病程短、随访事件密集,但影像随访缺失;组织学诊断在 2016 年被分子分型改写,队列里那批按形态入组的标签需要重新推导。

下面按「它解决什么麻烦、里面有什么、怎么取、边界在哪」四段展开,规模数字取自 GDC 项目端点与 TCIA 集合页,生物学结论取自 2008 年 Nature 与 2013 年 Cell 两篇主论文,核验日期写在附录。一屏速览见下表。

维度 内容
全称 The Cancer Genome Atlas Glioblastoma Multiforme
项目编号 TCGA-GBM(GDC project_id)
发起机构 美国国家癌症研究所(NCI)与国家人类基因组研究所(NHGRI)
立项与主论文 2005 年立项,2008 年 Nature 首篇,2013 年 Cell 重分析
病例规模 六百一十七例(GDC 项目页,2026-10 查询)
文件规模 三万零五百八十七份,合计约 263 TB
数据类别 10 类:简单核苷酸变异、测序读段、生物样本、临床、拷贝数、转录组、甲基化、蛋白组、体细胞结构变异、结构变异
数据子类型 25 种(data_type 口径)
实验策略 11 种:WXS、WGS、RNA-seq、miRNA-seq、甲基化阵列、SNP 阵列、表达阵列、RPPA、诊断切片、组织切片、ATAC-seq
样本构成 原发肿瘤六百零六例、复发肿瘤十三例、血液正常五百一十二例、实体组织正常四十例(按病例计,同一病例可有多类样本)
影像集合 TCIA TCGA-GBM,二百六十二例、五百七十五次检查、五千四百一十二个序列
影像模态 MR、CT、DX(平片),不含 PET
影像体量 73.51 GB(DICOM 原始)
病理切片 诊断切片三百八十九例、组织切片六百零六例,存 GDC
蛋白层 RPPA 二百三十七例
主键 患者 barcode(TCGA-XX-XXXX),影像与组学共用
访问层级 开放层免注册;受控层经 dbGaP phs000178
合并队列 GBMLGG 1,122 例(2016 年 Cell,分析层合并,非独立项目)

§1 它要解决的具体麻烦

1.1 一句话结论:它解决的是「通路被打断在哪一节」的问题

2005 年之前的脑肿瘤研究,主流做法是一个基因一个基因地看。EGFR 扩增、TP53 突变、PTEN 缺失各自都有一批论文,但没人能说清这些事件在同一批肿瘤里如何组合,也无法回答「一个肿瘤平均要被打断几条通路才会恶变」。

TCGA-GBM 的设计就是为这个问题准备的:对同一批标本同时测拷贝数、表达与甲基化,挑一部分做外显子测序,然后把三类事件叠在同一套基因与通路上。它得到的结论不是「某个基因很重要」,而是「三条通路几乎在所有样本里都被打断,只是打断的方式不同」。这个结论后来被写进教科书,也成为泛癌分析的起点。

把这句话展开成可操作的判断,可以拆成三层。第一层,如果你要问「某个分子事件在 GBM 里的整体位置」,它是最省事的一站,因为通路层面的整理已经做完。第二层,如果你要问「影像能不能预测分子亚型」,它也能用,但只有二百六十二例有影像,且亚型标签需要自己推导。第三层,如果你要问「复发时发生了什么」,它答不了,复发样本只有十三例。

你想问的问题 需要什么 够不够
三条核心通路在各样本里如何被打断 拷贝数、突变、表达、甲基化同源 够,这是它的设计目标
分子亚型与影像表型是否相关 同源影像加亚型标签 勉强够,影像二百六十二例,标签需自推
复发与进展的分子机制 配对的原发与复发标本 不够,配对样本极少
治疗前后分子变化 治疗前后的连续取样 不够,采集以术前为主

1.2 三条常见的错误期待

第一条,把它当成「六百多例的完整多组学」。项目确实收录六百一十七例,但完整跑过四类以上测量的只是其中一部分,蛋白层只有二百三十七例,RNA-seq 覆盖二百九十三例。按最大病例数规划多组学融合,会在取交集时发现可用样本远少于预期。

第二条,把它当成「带标注的脑肿瘤影像集」。影像确实有二百六十二例,但项目本身不提供任何病灶勾画,也没有统一的分割标签;影像来自常规诊疗,扫描仪、层厚、序列组合由各站点自定。要拿它做分割,得先解决标注问题。

第三条,把它当成「按当代分类整理好的队列」。它按组织学诊断入组,临床表里没有 IDH 状态,也没有 1p/19q 共缺失。2016 年之后的世界卫生组织分类把这些信息放在诊断的核心位置,用它做当代研究必须自己补这一步。

常见期待 实际情况 后果
「六百多例都是全组学」 各层覆盖从二百三十七到六百零六例不等 取交集后样本量骤减
「影像是标准研究序列」 常规诊疗图像,参数由站点决定 跨站点训练容易学到设备特征
「自带肿瘤分割掩膜」 无官方勾画 分割任务需自建标注
「临床表里有 IDH 状态」 无此字段,需从突变或甲基化推导 分型结论依赖自建规则
「复发标本够做纵向分析」 复发肿瘤样本仅十三例 纵向结论没有统计效力
「和低级别胶质瘤是同一批人」 两个项目病例互斥 合并计数会重复或错配

1.3 它能给到的与给不到的

你的需求 它给的 够不够
通路层面的分子图谱 三条核心通路的完整事件表 够,且是文献标准口径
突变与拷贝数的联合分析 外显子突变、基因级拷贝数、片段级拷贝数 够,覆盖数百例
表观遗传分型 甲基化 beta 值,四百二十三例 够做 G-CIMP 一类的分型
蛋白层验证 RPPA 二百三十七例 够做靶点验证,不够做规模训练
影像配分子标签 二百六十二例术前影像配同源组学 够做探索性研究
病理切片弱监督 诊断切片三百八十九例 够做多实例学习
病灶分割掩膜 无 不够,需借 BraTS 等派生集合
复发与进展分析 十三例复发样本 不够
治疗反应标签 缺统一疗效判定与治疗表 不够

1.4 你的处境与它能帮到的事

你的处境 它能帮你什么 下一步该做什么
想验证一个新亚型在公共队列里是否成立 提供表达、甲基化、拷贝数三套矩阵 先固定分型方法,再与已发表标签比对
手上有脑 MRI 模型,缺分子标签做外部验证 二百六十二例同源影像可做一小折测试 按 barcode 取交集,按站点分层报告
做病理基础模型,需要带分子标签的切片 诊断切片三百八十九例配突变与生存 走 GDC 切片下载,注意 SVS 体积
写综述,需要 GBM 多组学的标准参照 2008 年 Nature 与 2013 年 Cell 是通用引用 同时给项目 DOI 与数据查询日期
教学,需要一条完整的多模态病例 一位患者可串起影像、切片、突变、表达、随访 用 §4.3 的案例做贯穿演示
想做复发机制研究 复发样本仅十三例 改用 GLASS 等纵向队列

1.5 一句话判断

它是一个把「通路被打断在哪一节」问到极致的项目,样本量从来不是它的卖点;用它做机制与参照合适,用它做规模与纵向分析不合适。

1.6 不适合谁:四条边界

边界 具体表现 更适合的替代
追求影像规模与标注 二百六十二例,无勾画,参数不齐 BraTS 系挑战集、UCSF-PDGM
追求复发与纵向 复发样本十三例,无影像随访 GLASS、TCGA 之外的纵向队列
追求当代分子诊断标签 无 IDH 与 1p/19q 字段 自行推导,或改用已整理分型的派生队列
追求治疗反应 无统一疗效判定 QIN-GBM 一类治疗反应集合、临床试验队列

§2 数据构成与规格

2.1 两棵树:GDC 的患者树与 TCIA 的影像树

GDC 的组织方式是四层:case 是患者,sample 是切下来的某一块组织,aliquot 是分装出的一管核酸,analysis 是跑出来的文件。四层都用同一套 barcode 逐级收窄,前十二位在四层里不变。

TCIA 是另一棵树,只有患者、检查、序列、图像四级,且只在患者级与 GDC 相接:DICOM 的 PatientID 就是 GDC 的 case_submitter_id,两者字符串相等,不需要桥表。这个设计在 TCGA 各癌种里是统一的,GBM 也不例外。

两棵树的分工要说清:组学与临床在 GDC,放射影像在 TCIA,病理切片扫描在 GDC 的文件层而不是 TCIA。做三模态对齐时,正确的顺序是先按患者级把两棵树连起来,再在样本级内部分析切片与组学,不要把影像往某个具体等分样上对。

层级 GDC 侧 TCIA 侧 连接键
患者 case,十二位 barcode Patient,PatientID TCGA-XX-XXXX
样本 sample,含样本类型码 无 不适用
检查 无 Study,含检查日期 不适用
序列 无 Series,含模态与序列号 不适用
文件 分析文件,含 UUID DICOM 实例 不适用
临床 病例表、诊断表、随访表 集合页附带的临床压缩包 患者 barcode

2.2 字段与记录的字典

临床侧的主体是三块:人口学、诊断、随访。人口学十五个字段,诊断十五个字段,随访八个字段。这个规模在 TCGA 各癌种里偏小,原因不是采集粗糙,而是脑肿瘤的临床路径短:从诊断到死亡往往不足两年,中间的可记录节点少。

字段块 代表字段 字段数 说明
人口学 sex_at_birth、race、ethnicity、vital_status 15 年龄以诊断时年龄给出
诊断 primary_diagnosis、morphology、tissue_or_organ_of_origin 15 以 ICD-O-3 编码为主
随访 days_to_follow_up、disease_response、timepoint_category 8 多行,需按病例取最后一行
结局 days_to_death、days_to_last_follow_up 2 分列在人口学与随访表中

样本侧的字段决定了一份数据属于肿瘤还是正常、原发还是复发。这一位字符的差别,直接决定下游所有比较的方向。

字段 取值示例 用途
sample_submitter_id TCGA-06-0131-01A 样本级标识
sample_type Primary Tumor / Recurrent Tumor 区分原发、复发与正常
tissue_type Tumor / Normal 组织来源
specimen_type Solid Tissue 标本形态
preservation_method FFPE / Snap Frozen 影响核酸质量
days_to_collection 整数 相对诊断日的天数

文件侧的关键是四个属性:data_category 给粗分类,data_type 给细分类,experimental_strategy 给测量方法,access 给访问层级。四者组合起来才能准确描述一份文件是什么。

字段 取值示例 用途
file_id UUID 字符串 下载与去重的主键
file_name 形如 uuid.htseq.counts.gz 判断内容与格式
data_category Transcriptome Profiling 粗粒度分类
data_type Gene Expression Quantification 细粒度类型
experimental_strategy RNA-Seq 实验方法
access open / controlled 决定是否需要 dbGaP
file_size 字节数 规划存储
md5sum 三十二位十六进制 校验完整性

2.3 一条记录的解剖

把一位患者摊开,可以看到 barcode 如何逐级收窄。前四段跨模态通用,后四段只在组学侧有意义。

TCGA-06-0131-01A-01R-1849-01
│    │  │    │ │ │ │    │
│    │  │    │ │ │ │    └─ 中心代码
│    │  │    │ │ │ └────── 板号
│    │  │    │ │ └──────── 分析物类型(R = RNA)
│    │  │    │ └────────── 部分与孔位
│    │  │    └──────────── 样本类型(01A = 原发肿瘤)
│    │  └───────────────── 患者序号
│    └──────────────────── 组织来源站点
└───────────────────────── 项目前缀

这个结构带来三个实用推论。第一,跨模态连接只用到前十二位,多截一位就会把交集做小。第二,样本类型码的第两位区分原发与复发,做复发分析时不能只按患者去重,否则会把同一患者的两次手术压成一条。第三,切片属于样本级,一个患者可能有多张切片,聚合到患者级时要先定规则。

标识段 位数 含义 是否跨模态
TCGA 4 项目前缀,固定 是
TSS 2 组织来源站点 是
Participant 4 患者序号 是
Sample 3 样本类型与顺序 否,仅组学与病理
Vial/Portion 2 份数与孔位 否
Analyte 2 分析物类型 否
Plate 4 板号 否
Center 2 测序或分析中心 否
对齐目标 应使用的粒度 不能做的事
影像与突变 患者级十二位 用样本级标识去连影像
影像与切片 患者级先连,再在样本内挑切片 假设切片对应影像里的病灶
表达与甲基化 样本级或等分样级 跨样本类型合并
临床与全部 患者级 把随访多行直接并进主表

2.4 数据类别与实验策略

十个数据类别里,文件数最多的是简单核苷酸变异,其次拷贝数,两者合计超过一半。这个形状与其他 TCGA 项目一致:变异调用会按样本、按调用器、按格式重复落盘,文件数天然膨胀;真正的信息量要看覆盖病例。

数据类别 文件数 覆盖病例 说明
Simple Nucleotide Variation 9,789 601 体细胞突变调用,含多种格式
Copy Number Variation 6,981 600 基因级与片段级拷贝数
Sequencing Reads 3,251 428 比对后的读段文件
Biospecimen 2,920 617 样本与标本元数据
Transcriptome Profiling 1,898 572 表达与小 RNA 定量
Structural Variation 1,847 375 结构变异调用
DNA Methylation 1,350 423 探针级 beta 值
Clinical 1,219 617 病例、诊断、随访表
Somatic Structural Variation 1,089 335 体细胞结构变异
Proteome Profiling 243 237 RPPA 蛋白丰度

实验策略的覆盖面更能说明「同一患者被几种尺子量过」。基因分型阵列覆盖最广,接近六百例;外显子与全基因组测序各覆盖三四百例;表达阵列与 RNA-seq 分属两代平台,覆盖合计五百余例;蛋白层最薄。

实验策略 文件数 覆盖病例 主要产物
Genotyping Array 7,348 599 拷贝数与杂合性缺失
WXS 7,897 401 外显子突变
WGS 4,688 347 全基因组变异与结构变异
RNA-Seq 3,519 293 表达定量
miRNA-Seq 834 273 小 RNA 表达
Expression Array 560 540 芯片表达
Methylation Array 1,350 423 甲基化 beta 值
Diagnostic Slide 860 389 诊断切片扫描
Tissue Slide 1,193 606 组织切片扫描
Reverse Phase Protein Array 243 237 蛋白丰度
ATAC-Seq 9 9 染色质可及性,补测小批

这里有一处最容易读错的地方:表达层有两套平台并存,芯片五百四十例、测序二百九十三例,两者基因标识与归一方式都不同,不能直接拼接。做表达分析前先决定用哪一套,把选择写进方法学,比事后校正更省事。

2.5 与同类资源的字段差异

同样是 GBM 多组学,TCGA-GBM 的长处在通路级整合与影像配对,短处在样本量与当代分型标签。把它与 CPTAC 的蛋白组版本、与影像挑战集放在一起看,分工就清楚了。

维度 TCGA-GBM CPTAC GBM BraTS 系
病例规模 六百一十七 九十九 千例级,逐年扩
主打层 突变、拷贝数、甲基化、表达 蛋白组与磷酸化蛋白组 多参数 MRI 与分割掩膜
影像 二百六十二例,无勾画 有配对影像子集 有,且带体素级标注
分型标签 需自推 IDH 与 1p/19q 部分整理 不提供分子标签
典型用途 通路机制与参照 蛋白层验证 分割与影像模型

§3 核心层 A:项目身份与三条边界

3.1 它在 TCGA 里的坐标

TCGA 是计划,GBM 是其中最早的一个项目。计划层面回答「癌症的共同分子逻辑是什么」,项目层面回答「胶质母细胞瘤这一种病被测量到了什么程度」。两个层面的数字不能互换:计划是三万多种样本量级的集合,项目是六百一十七例。

层级 对象 规模口径 回答的问题
计划 TCGA 33 个癌种、上万患者 跨癌种的分子共性
项目 TCGA-GBM 六百一十七例 单一癌种的多组学全貌
合并分析 GBMLGG 见 §3.3 弥漫性胶质瘤的连续谱
影像子集 TCIA TCGA-GBM 二百六十二例 影像与分子的配对

3.2 与泛癌总库的明确分工

库内已有 TCGA 泛癌总库词条,它把 33 个癌种放在一个入口下,重点在跨癌种比较与统一临床资源。本篇只做一件事:把 GBM 这一个项目的七层组学、影像子集、临床字段与取数路径讲到能直接开工。

分工可以用三句话记住。泛癌总库回答「胶质母细胞瘤在 33 个癌种里处于什么位置」,本篇回答「这一个项目里每一层数据长什么样」。泛癌总库的典型表格是癌种乘模态的覆盖矩阵,本篇的典型表格是患者乘数据层的覆盖与取数路径。泛癌总库的临床口径统一到四类生存终点,本篇要讲清每个字段的取值与缺失形态。

关注点 泛癌总库 本篇
分析单位 癌种 患者与样本
核心表格 癌种乘模态覆盖 病例乘数据层覆盖
影像处理 一句带过,指向 TCIA 模态、站点、序列、对齐逐项展开
临床处理 统一到派生终点 字段块、缺失形态、版本差异
取数代码 通用 GDC 查询 GBM 专用过滤与交集脚本
适用任务 跨癌种迁移、泛癌分型 脑肿瘤多模态、放射基因组学

3.3 与低级别胶质瘤的边界:项目互斥,分析层合并

这是本篇最需要写清的一处。GDC 里 TCGA-GBM 与 TCGA-LGG 是两个独立项目,前者六百一十七例,后者五百一十六例。按患者 barcode 取交集,结果为零;并集是一千一百三十三例。也就是说,在项目层,两者没有一个病例是重复的。

文献里频繁出现的 GBMLGG 不是第三个项目,而是 2016 年 Cell 那篇弥漫性胶质瘤论文的分析队列:把两个项目的样本合起来,按 IDH 突变与 1p/19q 共缺失重新分层,最终纳入一千一百二十二例。与两个项目并集相差的十一例,是该分析未纳入的样本,官方未逐例说明原因。

口径 例数 层级 能否当作 GBM 规模
TCGA-GBM(GDC 项目) 六百一十七 项目 是
TCGA-LGG(GDC 项目) 五百一十六 项目 否,另一癌种
两者并集 一千一百三十三 项目相加 否,含低级别
GBMLGG(2016 合并分析) 1,122 分析队列 否,含低级别

这个结构有一个临床背景:组织学诊断会在病程中改变,一部分低级别胶质瘤在复发时被判为胶质母细胞瘤。但在 TCGA 的入库规则里,患者只归入一个项目,归属由入组时的诊断决定,因此项目层看不到同一个人出现在两边。要做「低级别升级为高级别」的研究,靠这两个项目的交集是找不到的,必须用真正的前后配对队列。

另有一层标签问题需要说明。2016 年的世界卫生组织分类把胶质母细胞瘤拆成 IDH 野生型与 IDH 突变型两个实体,2021 年的第五版进一步把 IDH 突变型的高级别病例归入星形细胞瘤四级,不再使用胶质母细胞瘤这一名称。TCGA-GBM 按形态入组,队列里因此混有一部分按当代标准不该叫 GBM 的病例。这个比例不影响数据质量,但影响结论的表述:凡是用这批数据做分型,都要说明标签是按形态还是按分子。

3.4 与 QIN-GBM、BraTS、CPTAC-GBM 的边界

库内有多个与胶质母细胞瘤相关的词条,容易混。三者的差别不在病种,在「测的是什么」和「标识能不能连」。

QIN-GBM 是麻省总医院提交的治疗反应研究集合,五十四例新诊断患者,术后治疗前各做两次基线扫描,主打灌注与弥散参数的重复性。它是纯影像集合,标识体系不是 TCGA barcode,与 TCGA-GBM 的病例无法按主键连接,也不提供分子数据。

BraTS 是分割挑战赛系列,提供多参数 MRI 与体素级肿瘤亚区标注。它与 TCGA-GBM 有血缘:早期 BraTS 的部分影像取自 TCGA,TCIA 上也有专门的派生集合记录这批转换后的数据。但两者的用途完全不同,BraTS 给的是标注与评测协议,TCGA-GBM 给的是分子与临床。

CPTAC-GBM 是蛋白组学主导的队列,九十九例,主打深度蛋白组与磷酸化蛋白组。它有自己的样本来源与标识体系,与本项目不能按主键直接连接,只能在同一病种的层面参照。

词条 对象 规模 能给什么 与本篇的关系
QIN-GBM 灌注 MRI 重复性 五十四例 双基线参数稳定性 无主键可连,纯影像
BraTS 系 分割挑战 千例级 体素级亚区标注 部分影像同源,用途不同
CPTAC-GBM 蛋白组 九十九例 深度蛋白与磷酸化 病种相同,队列独立
本篇 多组学项目 六百一十七例 突变、拷贝数、甲基化、表达、蛋白、影像 参照层

一句话区分:要做分割,去 BraTS;要做灌注参数重复性,去 QIN-GBM;要做蛋白层深度验证,去 CPTAC-GBM;要做通路机制与多模态关联,回本篇。

3.5 时间线与版本

这个项目的时间线有四段,容易被压成一段。采集与测量在 2008 年前后完成,2013 年做了一次五百余例的重分析,2016 年被并入弥漫性胶质瘤的合并分析,此后 GDC 持续做 harmonized 重处理。

时间 事件 对使用者的意义
2005 项目立项,成为 TCGA 首个癌种 方法学模板由此确立
2008 Nature 首篇,二百零六例多平台整合 三条核心通路成为标准结论
2013 Cell 重分析,五百余例 突变谱与亚型框架更新
2014 TCIA 影像集合首版 影像可公开获取
2016 并入 GBMLGG 合并分析 分型从形态转向分子
2023 TCIA 第五版更新 修复部分 DICOM 标签
2026 GDC 持续重处理 需记录下载时的数据版本

3.6 许可与使用边界

三层许可并存。GDC 开放层遵循 TCGA 数据使用协议,要求致谢研究网络、不得尝试重新识别个体;受控层走 dbGaP 的数据使用认证;TCIA 影像按该集合页标注的访问政策与引用要求使用,引用必须含集合 DOI。

数据层 许可 关键义务 常见误用
GDC 开放层 TCGA 数据使用协议 致谢研究网络,不重新识别 当作无限制公有领域
GDC 受控层 dbGaP 数据使用认证 限定研究用途,按期报告 转发给未授权合作者
TCIA 影像 集合页标注的访问政策 引用含 DOI,致谢 TCGA 只引论文不引数据 DOI
衍生结果 依原许可 标明来源与版本 宣称对原始数据拥有新版权

3.7 标识体系与命名习惯

文献里对这个项目有多种写法,检索时会漏。GDC 用 TCGA-GBM,TCIA 用同名集合,论文里常写 GBM 或 TCGA glioblastoma,旧数据包用 GBM 作为队列缩写。做系统集成时,建议以 project_id 为主键,其余写法建别名表。

写法 出现场景 是否等价
TCGA-GBM GDC、TCIA、API 是,标准写法
GBM 论文、旧数据包 是,需上下文确认
TCGA glioblastoma 综述、新闻 是,口语化
Glioblastoma Multiforme 项目全称 是,正式名
GBMLGG 合并分析 否,含低级别胶质瘤
HGG 旧文献 否,含三级胶质瘤

§4 核心层 B:七层组学与影像子集

4.1 七层组学的覆盖形状

把覆盖病例按层排开,会看到一条很陡的斜线。拷贝数与突变在六百例上下,表达与甲基化在四百到五百例之间,蛋白掉到二百三十七例,补测的染色质可及性只有九例。

数据层 覆盖病例 相对厚度 典型用途
拷贝数与杂合性缺失 599 最厚 受体酪氨酸激酶扩增、染色体臂级事件
体细胞突变 601 最厚 驱动基因、突变签名、超突变表型
表达定量 572 厚 分型、通路活性、预后模型
甲基化 423 中厚 G-CIMP 表观亚型、MGMT 启动子状态
结构变异 375 中 融合与重排
蛋白 237 薄 通路蛋白验证、药物靶点
染色质可及性 9 极薄 机制探索,无统计效力

这张表的读法是「取交集」。若任务需要同时用到突变、表达与甲基化,可用样本会降到四百例以下;再加上蛋白,只剩两百例出头。多组学融合的样本量由最薄的一层决定,这一条在本项目上体现得尤其明显。

另有一处历史层叠需要处理:表达有两代平台,甲基化也有两代芯片。早期用 27K,后期用 450K,两者探针集不同,合并时只能取交集探针,或者只用其中一代。这个选择会直接改变可分型的样本量。

4.2 影像子集:二百六十二例的技术规格

TCIA 的集合页给出一组可直接引用的数字:二百六十二例受试者、五百七十五次检查、五千四百一十二个序列、四十八万余张图像,体量不足八十 GB,模态标注为 MR、CT 与 DX。检查次数多于受试者,说明部分患者有多次扫描;图像张数远大于序列数,是因为一个序列包含几十到几百层。

指标 数值 口径
受试者 262 TCIA 集合页
检查 575 Study 级
序列 5,412 Series 级
图像 481,158 实例级
体量 73.51 GB DICOM 原始
模态 MR、CT、DX 磁共振、CT 与平片
版本 第五版,2023-08-07 集合页标注 Complete
DOI 10.7937/K9/TCIA.2016.RNYFUYE9 引用必需

三点必须说明。第一,这些影像是常规诊疗的一部分,不是按研究方案采集的,扫描仪厂商、场强、层厚、序列组合都由各站点决定,集合页对此有明确表述。第二,每位患者的基线影像是术前检查,这一点在集合页的标识说明里写明,因此影像与手术标本在时间上是对应的。第三,DICOM 头里的日期被整体平移了三到十年,同一站点内的相对间隔保留,跨站点的绝对日期不可比。

影像覆盖率约为队列的四成,这个比例在 TCGA 各癌种里属于中上。但「有影像」不等于「影像可用」:一部分序列是定位像或平片,真正可用于肿瘤分析的是增强 T1 与 FLAIR 一类的序列。做队列规划时,要按序列类型而不是按患者数估算可用样本。

影像类型 典型用途 使用注意
增强 T1 强化灶范围、切除程度评估 层厚差异大,需重采样
T2 与 FLAIR 水肿范围、非强化灶 序列是否压脂不一致
CT 钙化、出血、术后改变 多为平扫,软组织对比有限
DX 平片 术后或定位参考 对肿瘤分析价值有限

4.3 案例:一位患者从术前扫描到切片扫描的四段轨迹

把一位患者当作具体对象看,比抽象讨论「多模态能不能连」更有用。这位患者在入组站点完成术前 MRI,图像以 DICOM 存入 TCIA;手术取出的标本进入生物样本库,分装后跑拷贝数阵列、外显子测序、RNA-seq 与甲基化;病理科留一张诊断切片,扫描后以 SVS 形式存在 GDC;临床团队按随访表记录生命状态与天数。四段轨迹共用同一个十二位 barcode。

这条轨迹有三个可验证的接点。第一,影像的 PatientID 与 GDC 的 case_submitter_id 完全相等,直接字符串连接即可。第二,影像日期被平移,而 GDC 的天数相对病理诊断日,两边相减得不到「术前多少天」这个变量,只能在同侧内部使用。第三,切片属于样本级,一个患者可能有多张,聚合到患者级前要先定「取第一张」还是「取面积最大的一张」。

轨迹段 存放位置 标识粒度 可连接性
术前影像 TCIA 患者级 与 GDC 患者级直连
手术标本 GDC Biospecimen 样本级 患者级下钻
分子测量 GDC 各数据类别 样本或等分样 患者级聚合
病理切片 GDC Slide Image 样本级 与分子同源
随访结局 GDC Clinical 患者级 与影像同源

对这个案例最常见的误用,是假设影像里的强化灶与某张切片上的肿瘤是同一块组织。实际上切片来自取样部位,影像覆盖全脑,两者只能在患者级对应,不能逐点对应。做放射基因组学时,这种粒度错配是最常见的错误来源。

4.4 病理切片:一个容易被忽略的形态层

GDC 侧有两类切片扫描:诊断切片覆盖三百八十九例,组织切片覆盖六百零六例,文件合计两千零五十三份。这两类不是一回事:诊断切片是病理科用于出报告的 H&E,组织切片往往来自研究取样的蜡块,可能含不同染色的连续切片。

维度 诊断切片 组织切片
覆盖病例 389 606
来源 病理科存档 研究取样蜡块
典型用途 形态学、弱监督分型 与分子同源的区域定位
粒度 样本级 样本级

这个层的价值在于「与分子同源」。切片来自被测量过的那块组织,病人级标签可以直接挂上去,不需要额外的配对工作。代价是标签仍停留在患者级,切片级的分级、坏死比例、微血管增生这些形态学指标需要自己标。

4.5 临床随访的可用性

随访表只有八个字段,主体是随访天数与疾病状态。这个结构足以算总生存,但不足以算无进展生存,因为进展日期没有系统记录。

终点 需要的字段 是否可直接计算 备注
总生存 vital_status、days_to_death、days_to_last_follow_up 是 缺失最少,最常用
疾病特异生存 死因相关字段 部分 死因归类需人工核对
无进展生存 进展日期 否 只有疾病状态,无日期
无病间期 复发日期 否 同上
治疗反应 统一疗效判定 否 无统一判定标准

脑肿瘤的生存曲线下降很快,事件数在这类队列里通常不是瓶颈,这也是它常被用作生存模型基准的原因之一。但要注意删失的形态:失访比例、随访长度在各个入组站点之间不一致,做多因素模型时应当把站点作为分层变量。

4.6 标注与派生资源

原始项目不提供病灶勾画、分割掩膜或影像组学特征。这些东西以「分析结果」的形式挂在 TCIA 的相关资源里,或由第三方论文发布。使用前要区分官方集合与社区派生,两者的许可与质量责任不同。

资源 类型 与原始集合的关系 许可提示
TCGA-GBM-Radiogenomics 分析结果 影像特征与分子关联 随资源页声明
TCGA-GBM-QI-Radiogenomics 分析结果 影像组学特征表 随资源页声明
BraTS-TCGA-GBM 派生集合 转换为挑战赛格式的影像 随资源页声明
GBM-MR-NER-Outcomes 分析结果 命名实体识别与结局 随资源页声明
Pan-Cancer-Nuclei-Seg 分析结果 跨癌种细胞核分割 随资源页声明

派生资源省掉的是重复标注,引入的风险是切分不一致。把派生特征与原始 DICOM 拼接前,必须核对患者清单与序列编号,否则会把训练集的信息泄漏到测试集。这一条在影像任务上比在组学任务上更容易出错,因为序列标识没有 barcode 那样的自解释结构。

4.7 可复现性与版本

可复现性取决于三件事是否被记录:下载时对应的 GDC 数据版本、使用的参考基因组与管线、影像来自 TCIA 的哪一版。数据生成在 2013 年前后冻结,但 harmonized 重处理仍在继续,同一份文件在不同年份可能有不同的 UUID。影像侧同样如此,第五版修的是 DICOM 标签,序列数与患者数没变,但文件内容变了。

复现要素 应记录的内容 缺失后果
数据版本 GDC 数据发布号与查询日期 文件 UUID 漂移,无法复算
参考基因组 GRCh38 或 GRCh37 坐标不一致,变异无法对齐
管线 比对与变异调用器的版本 表达与突变调用不可比
影像版本 TCIA 集合版本与 DOI 文件内容变化,切分不可复现
临床表版本 下载日期与表名 随访更新导致终点变化

§5 核心层 C:怎么拿到、怎么接上

5.1 三条通道与各自的门槛

取数有三条通道:GDC 网页、GDC 接口、TCIA 下载器。网页适合小批量挑选,接口适合可重复的脚本,TCIA 下载器适合整包拉取。三条通道的认证要求不同,规划时要分开算。

通道 适用场景 认证 批量能力
GDC 门户网页 挑选少量文件、看元数据 开放层免注册 弱
GDC 接口 可重复查询、生成清单 开放层免注册 强
gdc-client 按清单批量下载 受控层需 token 强,支持断点
TCIA 网页 浏览集合与引用信息 无 弱
TCIA Data Retriever 整包或按序列下载 无 强
cBioPortal 直接取整理好的矩阵 无 中,口径固定

5.2 用 GDC 接口拿项目概览

项目端点返回病例数、文件数、数据类别与实验策略的汇总。这一段查询是所有后续规划的依据,建议在方法学里写明查询日期。

curl -s "https://api.gdc.cancer.gov/projects/TCGA-GBM?expand=summary,summary.experimental_strategies,summary.data_categories" \
  | python3 -c "import sys,json; d=json.load(sys.stdin)['data']['summary']; print(d['case_count'], d['file_count'], round(d['file_size']/1e12,1))"

返回三个数字:病例数、文件数、总字节数换算成 TB。注意总量包含受控层,不等于开放层可下载的体量;开放层文件约占全部文件的一半不到,主要是因为受控的读段与胚系变异占了大头。

端点 返回内容 典型用途
/projects/TCGA-GBM 项目元数据与汇总 立项、写方法学
/projects/TCGA-GBM?expand=summary 类别与策略汇总 规划存储与算力
/cases 病例级字段 建临床表
/samples 样本级字段 区分原发、复发与正常
/files 文件级字段 生成下载清单
/files/ 单文件元数据 校验与排错

5.3 按数据类别筛文件

文件端点的过滤器是嵌套 JSON,结构固定:最外层是 and,每个条件是一个 in。把项目、数据类别、访问层级三个条件叠起来,就能得到一份干净清单。

import requests

filters = {
    "op": "and",
    "content": [
        {"op": "in", "content": {"field": "cases.project.project_id", "value": ["TCGA-GBM"]}},
        {"op": "in", "content": {"field": "data_category", "value": ["DNA Methylation"]}},
        {"op": "in", "content": {"field": "access", "value": ["open"]}},
    ],
}
payload = {
    "filters": filters,
    "fields": "file_id,file_name,file_size,md5sum,cases.samples.sample_type",
    "format": "JSON",
    "size": 5000,
}
r = requests.post("https://api.gdc.cancer.gov/files", json=payload, timeout=120)
hits = r.json()["data"]["hits"]
print(len(hits), round(sum(h["file_size"] for h in hits) / 1e9, 1), "GB")

这段代码返回文件级记录,不是甲基化矩阵。要做分型,还得把每个样本的 beta 值读进来按探针合并,或者直接取 GDC 打包的矩阵。前者灵活,后者省事,选择取决于是否需要自定义探针过滤。

过滤字段 常用取值 说明
cases.project.project_id TCGA-GBM 项目限定,必填
data_category DNA Methylation 粗分类
data_type Methylation Beta Value 细分类
experimental_strategy Methylation Array 方法限定
access open / controlled 访问层级
cases.samples.sample_type Primary Tumor 样本类型
analysis.workflow_type 具体管线名 口径固定用

5.4 从 TCIA 取影像

TCIA 是另一套接口。集合名固定为 TCGA-GBM,先用患者端点拿清单,再用序列端点拿每位患者的序列,最后用下载器按清单取 DICOM。整个过程不需要账号,但要接受数据使用条款并按要求引用 DOI。

import requests

base = "https://services.cancerimagingarchive.net/nbia-api/services/v1"
patients = requests.get(f"{base}/getPatient", params={"Collection": "TCGA-GBM", "format": "json"}, timeout=180).json()
print("patients", len(patients))
first = patients[0]["PatientId"]
series = requests.get(f"{base}/getSeries", params={"Collection": "TCGA-GBM", "PatientId": first, "format": "json"}, timeout=180).json()
print(first, len(series), series[0]["Modality"])

如果只需要元数据做统计,这两步就够了;要图像则用官方下载器读取集合清单,或把序列实例标识交给下载器逐条取。五千余个序列意味着大量小文件,建议先按患者打包,再解压到按序列分目录的结构。

端点 参数 返回
getCollectionValues 无 全部集合名
getPatient Collection 患者清单
getStudy Collection、PatientId 检查清单
getSeries Collection、PatientId 序列清单
getImage SeriesInstanceUID 图像清单

5.5 用 barcode 做交集

影像与组学的连接只有一步,但写法决定后续正确性。核心是把两边都截断到十二位,再去重、连接。不要用样本级标识去连影像,也不要假设两边大小写一致。

import pandas as pd

clin = pd.read_csv("TCGA-GBM.clinical.tsv", sep="\t", low_memory=False)
tcia = pd.read_csv("tcia_patients.tsv", sep="\t")

clin["pid"] = clin["case_submitter_id"].astype(str).str.strip().str[:12].str.upper()
tcia["pid"] = tcia["PatientID"].astype(str).str.strip().str[:12].str.upper()

tcia_ids = set(tcia["pid"])
clin["has_image"] = clin["pid"].isin(tcia_ids)
merged = clin.merge(tcia[["pid"]].drop_duplicates(), on="pid", how="inner")

print("clinical", clin["pid"].nunique(), "imaging", len(tcia_ids), "overlap", merged["pid"].nunique())
print("primary", (clin["sample_type"] == "Primary Tumor").sum())

连接之后先看三件事:交集是否与官方集合页的受试者数接近、有没有重复的 pid、带影像与不带影像的两组在诊断与生存上是否可比。第三件事最容易被跳过,却是审稿人最先问的问题。§4.3 的案例给出的四段轨迹,正是这一步在单例上的展开。

检查项 期望 异常时的处理
交集数量 接近二百六十二 检查截断位数与大小写
pid 唯一性 无重复 去重后再连接
样本类型 以原发肿瘤为主 排除复发与正常组织
组间可比性 诊断与生存分布接近 做倾向评分或分层
影像日期 早于手术日期 核对日期平移规则

5.6 常见取数故障与口径问题

取数阶段的问题大多不是网络,而是口径。最常见的是把受控层文件写进开放清单,其次是混用新旧版本的文件标识,第三是拿样本级标识去连患者级影像。这些问题在数据量小时不易发现,等到跑完模型才暴露。

现象 原因 处理
清单里有文件下载 403 含受控层文件 按 access 字段过滤或申请 dbGaP
同一患者多份表达文件 不同等分样或不同平台 固定平台与管线再合并
表达矩阵基因数不一致 新旧注释混用 统一到同一注释版本
影像找不到对应患者 barcode 截断错误 统一到十二位并大写
切片与分子对不上 样本级与患者级混用 先患者级连接再下钻
生存天数出现负值 事件日期早于诊断日 按数据字典核查
甲基化样本数对不上 27K 与 450K 混算 按平台分别统计

5.7 存储与算力预算

按开放层全量、影像全量、切片全量三块估算。切片是隐形的大头,单张 SVS 从几十兆到数百兆,六百余例组织切片足以把磁盘吃满;影像虽然只有七十几 GB,但五千余个序列解压后的小文件数量对文件系统是另一种压力。

数据块 体量级 是否需要 GPU 备注
开放层组学矩阵 数十 GB 否 表达、甲基化、拷贝数为主
影像 DICOM 七十余 GB 否 解压后小文件众多
诊断与组织切片 数百 GB 至 TB 级 推理时需要 建议流式读取
受控层读段 TB 级 否 需 dbGaP 与机构存储
中间特征 数十至数百 GB 视方法 建议落盘复用

§6 与库内词条的关系

6.1 关系图谱

在这个库里,TCGA-GBM 处在「计划—项目—派生」三层结构的中间位置。向上是 TCGA 泛癌总库,向下是各类以它为基础再加工的派生集合,横向是脑肿瘤方向的其他队列与挑战集。

库内对象 关系 分工
TCGA 泛癌总库 上级计划 33 癌种的跨癌种比较与通用取数
TCGA-GBM 本篇 单项目的组学、影像、切片与临床细节
TCGA-LGG 姊妹项目 低级别胶质瘤,病例与本篇互斥
GBMLGG 合并分析 分析层合并 并列两项目后的合并队列,非独立项目
QIN-GBM 相邻队列 灌注 MRI 重复性,纯影像
BraTS 系 下游派生 分割挑战,部分影像与本篇同源
CPTAC-GBM 同病种对照 深度蛋白组,队列独立
UCSF-PDGM 影像对照 带分子标签的术前 MRI,规模相当
关系边 方向 判断依据
计划到项目 包含 项目编号以 TCGA- 开头
项目到合并分析 被合并 合并队列标注来源为 GBM 与 LGG
项目到派生 被加工 派生集合标注原始集合名
项目到同类 并列 无包含关系,仅病种相同

6.2 与 QIN-GBM 的边界:同名不同物

库内刚发布的 QIN-GBM 与本篇都叫胶质母细胞瘤,容易混。两者的差别是全方位的:QIN-GBM 是五十四例新诊断患者的灌注 MRI 重复性研究,每位患者术前做两次基线扫描,用来回答「同一个参数隔几天再测,值变不变」;本篇是六百一十七例的多组学项目,影像只是其中一层。

不可连接这一点要写清:QIN-GBM 的标识不是 TCGA barcode,两边的患者无法按主键对上,也没有共同的分子数据。把它们合并会得到一个人为的「更大队列」,实际上只是两批不相干的患者。

维度 本篇 QIN-GBM
对象 多组学项目 影像重复性集合
规模 六百一十七例,影像二百六十二例 五十四例
序列 常规 MRI、CT、平片 七类灌注与结构序列
分子数据 有 无
主键 TCGA barcode 独立的影像标识
典型任务 多模态关联、分型 参数可重复性

6.3 与 BraTS 的边界:同源不同用

BraTS 是分割挑战赛系列,给的是多参数 MRI 与体素级肿瘤亚区标注,配套固定的训练测试划分与评测指标。它与本篇有真实的血缘:早期版本的影像部分取自 TCGA,TCIA 上也有专门的派生集合记录这批转换后的文件。

但两者的用途正好错开。BraTS 的强项是标注与评测协议,弱项是分子与临床;本篇正相反。做分割不要在本篇上从零标,直接用 BraTS;要做影像与分子的关联,BraTS 没有标签,必须回本篇。若两者都用,先核对患者清单的交集,避免同一患者在训练集与测试集之间泄漏。

维度 本篇 BraTS 系
影像来源 常规诊疗,多站点 部分源自本篇,另有新增
标注 无 体素级亚区标注
评测协议 无 固定划分与指标
分子标签 有 无
适用任务 关联与分型 分割与泛化评测

6.4 检索路径

从需求出发找入口,比从数据集名出发更省时间。下面这张表把常见需求映射到先看哪个词条,再看哪个补充资源。

你的问题 先看 再看
胶质母细胞瘤的多组学全貌 本篇 §2、§3 TCGA 泛癌总库的癌种章节
低级别胶质瘤的数据 TCGA-LGG 本篇 §3.3 的合并口径
脑肿瘤分割的训练数据 BraTS 系 本篇 §4.6 的派生集合
灌注参数的重复性 QIN-GBM 本篇 §6.2 的边界说明
蛋白层的深度验证 CPTAC-GBM 本篇 §4.1 的覆盖形状
带分子标签的术前 MRI 本篇 §4.2 UCSF-PDGM 作外部验证
生存模型的公共基准 本篇 §4.5 泛癌总库的统一终点资源

6.5 组合使用的三种典型配方

配方一,影像预训练加分子微调:先用 BraTS 训练一个脑肿瘤编码器,再用本篇二百六十二例做小样本微调,预测甲基化亚型或生存。配方二,多组学发现加影像验证:先在本篇找表达或甲基化签名,再看签名与影像特征的关联,最后用 UCSF-PDGM 一类独立队列验证。配方三,通路机制加蛋白确认:以本篇的通路事件为骨架,用 CPTAC-GBM 的蛋白与磷酸化数据验证哪一层真正传导到了蛋白。

配方 主数据 辅助数据 注意
影像预训练加微调 本篇影像子集 BraTS 站点差异会导致负迁移
组学发现加影像验证 本篇组学 UCSF-PDGM 验证队列的扫描协议不同
通路机制加蛋白确认 本篇通路事件 CPTAC-GBM 两边样本不完全重叠

§7 质量评估与局限性

7.0 阅读本节前的一个前置校正

评价这个项目时,最容易发生的误判是把「计划的首创性」当成「项目的规格」。它是 TCGA 的第一个癌种,方法学意义上的地位无可替代;但落到可用数据上,蛋白层只有二百三十七例,复发样本只有十三例,影像没有标注。这些数字与「首个多组学项目」的光环不是一回事。

第二个前置校正是把「七层组学」当成「七层都厚」。层数多解决的是「能不能问跨层的问题」,不解决「每个问题有多少样本」。第三,把「按组织学入组」当成「按分子分型整理好」,2016 年之后的分类标准需要自己补。带着这三个校正再看评分,分数才有意义。

7.1 五维评分与 DAIMS 对照

评分采用五分制,每一维给一个整数并写明理由。分数不是排名,是让人快速判断哪一维需要额外投入。DAIMS 的三个描述层在这里对应为:描述层看文档与元数据,结构层看标识体系与层级,访问层看许可与取数通道。

维度 分数 理由
规模 3 病例数够做单癌种分析,蛋白与复发样本偏少
结构 5 barcode 体系严谨,跨模态一次字符串连接即可
配对 4 影像、切片、分子、随访同源,但影像仅覆盖四成
可复现 4 冻结加版本化,重处理导致标识漂移
合规 4 三层许可清楚,受控层申请周期长
文档 5 两篇主论文加数据字典,闭环完整
DAIMS 层 覆盖内容 缺口
描述层 项目页、数据字典、主论文、集合页 分子分型标签无官方字段
结构层 患者、样本、等分样、文件四级 影像树与组学树仅在患者级相接
访问层 开放层与受控层双通道 受控层审批周期以月计

7.2 分数怎么读

规模给 3 分,不是因为它小,而是与任务不匹配:做通路机制绰绰有余,做影像深度学习不够。结构给 5 分,是十二位 barcode 让跨模态连接变成一次字符串比较,省掉大量数据工程。配对给 4 分,扣在影像只覆盖四成,且没有统一协议。

可复现给 4 分,扣在 harmonized 重处理上:数据生成冻结了,但文件标识仍在变。合规给 4 分,扣在受控层的审批周期,开放层本身很顺。文档给 5 分,因为两篇主论文、数据字典与集合页形成了完整闭环,任何字段都能查到定义。

分数段 含义 建议动作
5 业内标杆,可直接依赖 按官方文档执行即可
4 有小坑,不影响结论 记录版本,做一次校验
3 需额外投入才能用 预留清洗与对齐时间
2 只能作为补充 不作为主数据
1 不建议使用 换数据集

7.3 已知的、可枚举的偏差

偏差一,选择偏差。入组要求标本量足够完成多平台测量,且患者需能耐受手术,因此体能状态差、肿瘤位置不适合切除的病例被系统性排除,队列的预后好于真实人群。

偏差二,站点偏差。影像来自多家影像来源站点,组学来自数十个组织来源站点,扫描仪与建库批次叠加。§4.3 的案例里那条「一位患者四段轨迹」,在跨站点时并不成立:不同站点的序列参数差异足以让模型学到设备而不是病变。

偏差三,时间偏差。采集集中在 2005 到 2013 年,治疗以放疗加替莫唑胺为标准,与今天的电场治疗、免疫治疗背景不同,用它训练的预后模型会高估风险。偏差四,标签偏差。按形态入组,含一部分按当代标准应归为 IDH 突变型星形细胞瘤的病例。

偏差 来源 影响方向 缓解方式
选择偏差 手术与标本门槛 预后偏好 按体能状态与切除程度分层
站点偏差 多中心采集 模型学到站点特征 留一站点验证
批次偏差 不同年份、不同平台 表达与甲基化漂移 分位数归一或批次校正
时间偏差 2005–2013 采集 治疗背景不同 不与当代队列直接比
标签偏差 组织学入组 分型结论偏移 自行推导 IDH 与 1p/19q

7.4 外部验证矩阵

一个数据集的可信度,很大程度上取决于它被外部重复验证过多少次。下面把常见结论与已有证据类型列出来,方便判断哪些可以直接引用。

结论 内部证据 外部验证 可信度
三条核心通路在多数样本中被改变 2008 年多平台整合 多个独立队列复现 高
四种表达亚型的存在 表达聚类 多篇独立分析与 meta 分析 高
MGMT 启动子甲基化与治疗结局相关 甲基化加临床 临床试验与多项回顾研究 高
甲基化分型可复现 IDH 与 1p/19q 分层 2016 年合并分析 独立队列与临床检测一致 高
影像特征可预测分子亚型 影像子集小样本 复现有限 中
影像组学可直接预测生存 单中心或少数站点 缺统一协议 低到中

7.5 它为何能持续:因为它成了脑肿瘤研究的坐标系

数据集的生命周期通常止于论文发表,这个项目不一样。数据在 2013 年前后冻结,引用量仍在增长,原因是它变成了一个参照系:新方法要在它上面证明自己不比旧方法差,新队列要拿它做外部验证,教科书要用它讲「癌症是多条通路被打断」。

稳定的代价是它无法回答新问题。单细胞、空间转录组、治疗前后配对、影像随访、复发机制,这些当代问题它都答不了。它的持续不是因为内容常新,而是因为它把「同一批肿瘤的多把尺子」这件事一次做完了,后来者要么重复这件事,要么在它的坐标系里定位自己。

§8 生态位与历史

8.1 历史沿革

这个项目的历史可以分成四段:立项与首篇、重分析与亚型、合并与分型转向、影像与派生的二次开发。每一段的重心不同,对今天使用者的意义也不同。

阶段 时间 关键动作 留下的资产
立项与首篇 2005–2008 多平台整合、通路分析 三条核心通路、方法学模板
重分析与亚型 2010–2013 五百余例重分析、表达亚型 突变谱、亚型框架
合并与分型转向 2016 与低级别胶质瘤合并 GBMLGG 队列、甲基化分型
影像与派生 2014–2023 影像集合发布与更新 二百六十二例术前影像
持续复用 2016– 作为基准被反复调用 参照系地位

8.2 首篇论文的三条结论

2008 年那篇论文的结论可以压成三条,每一条都改变了后来的研究方式。

第一条,三条核心通路。受体酪氨酸激酶信号、p53 通路、RB 通路,在这批肿瘤里几乎都被某种事件打断,且同一条通路内的事件互斥。这条结论把「找单个驱动基因」的范式改成了「看通路整体状态」。

第二条,新基因的发现。那篇论文报告了此前未被认识的 PIK3R1 频繁突变,此后它成为磷脂酰肌醇三激酶通路研究的标准引用之一。同一批数据也重新定位了 ERBB2、NF1、TP53 在这类肿瘤里的角色。

第三条,表观与治疗的连接。论文把甲基化、突变与临床治疗信息叠在一起,发现接受过治疗的肿瘤里,MGMT 启动子甲基化与错配修复缺陷导致的超突变表型相关。这一条把表观遗传、治疗史与突变负荷三者连成了一条线。

结论 证据类型 后续影响
三条核心通路被打断 拷贝数、突变、表达联合 通路级分析成为标准
PIK3R1 频繁突变 外显子测序 通路研究的标准引用
MGMT 甲基化与超突变相关 甲基化加治疗史 表观与治疗的连接范式

8.3 衍生资源的谱系

围绕这个项目长出了一整条衍生链:从原始文件,到表达与甲基化矩阵,到亚型标签,到影像组学特征与分割掩膜,再到预训练模型权重。每一层都加上自己的假设,越往下游,越要核对它与原始集合的对应关系。

层级 产物 依赖 风险
原始层 测序文件、DICOM、SVS 无 体量大,需自行处理
矩阵层 表达、甲基化、拷贝数矩阵 原始层 平台与归一不统一
标签层 亚型、IDH 状态、分级 矩阵层 推导规则各异
特征层 影像组学特征、细胞核特征 原始层 口径不统一
结论层 预后模型、关联分析 以上各层 外推性有限

8.4 同类项目对照

把本篇放进更大的脑肿瘤数据生态里,它的位置很清楚:不是最大的影像队列,不是最新的分型队列,也不是蛋白层最深的队列,但它是唯一把七层组学、病理切片、术前影像与生存放在同一主键上的公共资源。

项目 影像 分子 随访 定位
TCGA-GBM 二百六十二例 七层 生存为主 多模态参照系
TCGA-LGG 有配对影像 七层 生存为主 低级别对照
CPTAC-GBM 有配对子集 蛋白为主 有限 蛋白层验证
QIN-GBM 五十四例灌注 无 有限 参数重复性
BraTS 系 千例级带标注 无 无 分割评测
UCSF-PDGM 数百例 有 有 影像与分子配对

§9 使用指引

9.1 三种角色的路径

影像工程师的路径是「先影像、后标签」:从 TCIA 取 DICOM,按站点分层,先做分割或分类的预训练,再回 GDC 取甲基化亚型与生存做小样本微调。这条路线正是 §4.3 的案例所展开的四段轨迹,只是把单例放大成队列。

病理工程师的路径是「先切片、后生存」:从 GDC 取 SVS,做切片级自监督,再聚合到患者级做生存回归。生物信息学者的路径是「先矩阵、后关联」:直接取表达、甲基化、拷贝数矩阵,做分型与通路分析,需要影像时再按 barcode 连接。

角色 第一步 第二步 第三步
影像工程师 取 DICOM 与站点表 重采样与分层 小样本微调
病理工程师 取 SVS 与切片表 切片自监督 患者级聚合
生物信息学者 取表达与甲基化矩阵 分型与差异分析 按 barcode 关联影像
临床研究者 取临床表与随访 定义终点 多因素建模
角色 常见误区 纠正
影像工程师 把有影像的与没影像的混在一起训练 先按有无影像分层,再按站点分组
病理工程师 用切片级标签直接做患者级预测 先聚合到患者级
生物信息学者 把芯片与测序表达直接合并 固定一套平台
临床研究者 把随访多行当独立样本 每例只保留最后一行

9.2 一条最短的可用流程

从零到一张可解释的图,最短路径是:查项目概览、拉临床表、筛表达文件、按 barcode 与影像取交集、做一次分组比较。下面把前四步串起来,假设网络可达且本地已装 requests 与 pandas。

mkdir -p tcga-gbm/data && cd tcga-gbm/data
curl -s "https://api.gdc.cancer.gov/projects/TCGA-GBM?expand=summary" -o project.json
curl -s -X POST "https://api.gdc.cancer.gov/files" \
  -H "Content-Type: application/json" \
  -d '{"filters":{"op":"and","content":[{"op":"in","content":{"field":"cases.project.project_id","value":["TCGA-GBM"]}},{"op":"in","content":{"field":"data_category","value":["Transcriptome Profiling"]}},{"op":"in","content":{"field":"access","value":["open"]}}]},"fields":"file_id,file_name,file_size","format":"JSON","size":5000}' \
  -o transcriptome_files.json
python3 -c "import json; d=json.load(open('transcriptome_files.json')); print(len(d['data']['hits']))"

产出是清单,不是矩阵。下一步用 gdc-client 按清单下载,或直接取打包矩阵。无论哪条路,都建议先落一份患者清单与文件清单,再开始下载,这样中断可以续,换机器可以复现。

步骤 命令或端点 产出
项目概览 /projects/TCGA-GBM 规模与类别
文件清单 /files 待下载列表
批量下载 gdc-client download -m 本地文件
影像清单 TCIA getPatient 患者与序列
交集 pandas merge 分析队列
质控 缺失与重复检查 可用样本表

9.3 引用方式

引用要分三层:数据引用、论文引用、资源引用。只引论文不引数据 DOI,会让数据生产者的贡献被隐没;只引数据不引论文,会让使用者缺少生物学背景。影像侧要按集合页要求给出 DOI 与版本,组学侧要致谢研究网络的原句。

引用对象 应写的内容 常见错误
数据(组学) 项目编号、GDC、访问日期 只写 TCGA 不写项目
数据(影像) 集合名、版本、TCIA、DOI 漏版本号
首篇论文 Nature 2008,DOI 引成综述
重分析论文 Cell 2013,DOI 与首篇混引
合并分析 Cell 2016,DOI 当作原始数据引用
致谢 研究网络原句 改写致谢措辞

9.4 合规清单

使用前过一遍这张表,能避免大部分返工。尤其是受控层的转发限制与重新识别禁令,违反的后果不在技术上,在合规上。

检查项 要求 状态
影像引用 含 DOI 与版本号 待核
组学致谢 含研究网络原句 待核
受控数据 有 dbGaP 批准号 待核
合作者 均在授权名单内 待核
重新识别 未尝试 待核
二次分发 未分发原始受控文件 待核
衍生结果 标明来源与版本 待核

§10 AI 使用声明

项目 声明
数据来源 GDC 项目端点、TCIA 集合页、两篇主论文摘要,查询日期 2026-10-29
生成方式 人工撰写,AI 辅助整理结构与措辞
数字口径 只采用官方页面与论文摘要,未做二次统计
未核实项 各数据层的精确交叉覆盖、影像按序列类型的病例数
分布统计 不提供,避免把探索性统计写成规格
医学建议 无,本文不构成诊疗建议
更新机制 官方版本变化时修订,修订日期见页脚
责任边界 使用后果由使用者承担,引用以官方为准

常见误读(避坑清单)

常见说法 实际 证据
「TCGA-GBM 没有影像」 TCIA 有二百六十二例术前影像 TCIA 集合页受试者数
「GBM 与 LGG 有病例重叠」 两个项目按 barcode 交集为零 GDC 病例清单取交集
「GBMLGG 是第三个项目」 是 2016 年的分析队列 Cell 2016 摘要与 cBioPortal 条目
「六百多例都能做多组学」 蛋白层仅二百三十七例 GDC 策略汇总
「自带肿瘤分割掩膜」 无官方勾画 集合页与 GDC 文件类型
「临床表里有 IDH 状态」 无此字段,需自推 GDC 诊断字段构成
「复发样本够做纵向分析」 复发肿瘤样本仅十三例 GDC 样本类型统计
「病理切片存放在 TCIA」 切片在 GDC,TCIA 只有放射影像 两个门户的分工
「barcode 可以任意截断」 患者级十二位,样本级十六位 官方 barcode 定义
「影像是按研究方案采集的」 常规诊疗图像,参数由站点自定 集合页说明

FAQ

Q1:TCGA-GBM 和 TCGA-LGG 是什么关系,会不会重复计数?
在 GDC 里这是两个互斥的项目,按患者 barcode 取交集为零,并集是一千一百三十三例。文献里的 GBMLGG 是 2016 年把两者合并后的分析队列,规模略小于两者之和,既不是 GBM 的规模,也不是 LGG 的规模。

Q2:它有影像吗?
有。TCIA 的 TCGA-GBM 集合收录二百六十二例术前影像,含 MR、CT 与平片,共五百七十五次检查、五千四百一十二个序列,体量不足八十 GB。影像不在 GDC,GDC 侧只存病理切片扫描。

Q3:为什么 GDC 写六百一十七例,论文写二百零六例?
2008 年论文报告的是完成多平台整合的 206 例,其中 91 例做了外显子测序;GDC 的项目计数包含所有入库病例,含诊断未报告与仅有部分测量的记录。写方法学要说清用的是哪个口径。

Q4:组学数据能不能直接下载?
突变调用、表达矩阵、甲基化 beta 值、拷贝数、RPPA、临床表与切片扫描属开放层,注册账号即可下载;原始读段与胚系变异属受控层,需经 dbGaP 申请。

Q5:能不能直接拿到 IDH 野生型与突变型标签?
不能。项目按组织学诊断入组,临床表没有 IDH 字段,需要自己从突变调用或甲基化谱推导,并把规则写进方法学。

Q6:为什么同一患者有多份表达文件?
可能来自不同等分样、不同平台或不同管线。取数时应固定平台与工作流,并按样本类型筛选原发肿瘤,避免把多份文件当成多个患者。

Q7:可以用它做外部验证吗?
可以,但要说明验证集性质。它是回顾性、多站点、常规检查的集合,适合验证方法稳健性,不适合宣称临床可用性;验证时应按站点分层并报告站点间差异。

Q8:切片数据有多大,需要什么配置?
单张 SVS 从几十兆到数百兆,六百余例组织切片可达数百 GB 至 TB 级。建议流式读取或按需取块,不要一次性解压到内存。

术语表

术语 解释
GBM Glioblastoma Multiforme,胶质母细胞瘤
LGG Lower Grade Glioma,低级别胶质瘤
GBMLGG 把 GBM 与 LGG 合并的分析队列
G-CIMP 胶质瘤 CpG 岛甲基化表型
IDH 异柠檬酸脱氢酶基因,分型的核心标志
1p/19q 共缺失 少突胶质细胞瘤的分子标志
MGMT 启动子甲基化影响烷化剂疗效的标志
GDC Genomic Data Commons,NCI 的基因组数据门户
TCIA The Cancer Imaging Archive,NCI 支持的影像归档
barcode 逐级编码样本来源的字符串标识
TSS Tissue Source Site,组织来源站点代码
aliquot 从样本中分装出的一份可用于测量的等分样
WXS 全外显子测序
WGS 全基因组测序
RPPA 反相蛋白阵列,直接测量蛋白丰度
MAF Mutation Annotation Format,突变注释格式
CNV 拷贝数变异
GISTIC 拷贝数显著区域分析工具
beta 值 甲基化探针的甲基化比例估计
SVS Aperio 全切片图像格式
DICOM 医学影像通信标准与文件格式
DX 数字平片,本集合中的第三类模态
DSS 疾病特异生存
OS 总生存
DAIMS 数据集资产信息模型,统一描述描述层、结构层与访问层

附录 A:核心数字核对表

项 值 来源 来源层级
病例数 617 GDC 项目端点 官方接口
文件数 30,587 GDC 项目端点 官方接口
总体量 约 263 TB GDC 项目端点 官方接口
开放层文件 13,117 GDC 文件端点 官方接口
数据类别 10 GDC 项目端点 官方接口
数据子类型 25 GDC 文件聚合 官方接口
实验策略 11 GDC 项目端点 官方接口
影像受试者 262 TCIA 集合页 官方页面
影像检查 575 TCIA 集合页 官方页面
影像序列 5,412 TCIA 集合页 官方页面
影像图像 481,158 TCIA 集合页 官方页面
影像体量 73.51 GB TCIA 集合页 官方页面
影像模态 MR、CT、DX TCIA 集合页 官方页面
首篇论文规模 206 Nature 2008 摘要 论文
合并队列 1,122 Cell 2016 摘要 论文
结构项 值 来源
顶层单位 case(患者) GDC 数据字典
患者级标识长度 12 位 官方 barcode 定义
样本级标识长度 16 位 官方 barcode 定义
影像粒度 患者级 TCIA 说明
切片粒度 样本级 GDC 文件结构
跨模态连接方式 字符串相等 两个门户的说明

附录 B:数字的来源分层

层级 数字 说明
官方接口 病例数、文件数、类别与策略汇总、各层覆盖数 2026-10-29 查询 GDC 端点
官方页面 262、575、5,412、481,158 TCIA 集合页第五版
官方页面 2005 立项、2014 影像首版、2023 第五版 项目与集合说明
论文摘要 二百零六例与九十一例测序、三条核心通路 Nature 2008
论文摘要 五百余例、千余例合并 Cell 2013 与 Cell 2016
未核实 影像按序列类型的病例数、各层交叉覆盖 官方未给,未做二次统计

附录 C:完整参考文献

  1. The Cancer Genome Atlas Research Network. Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature, 2008, 455:1061–1068. DOI: 10.1038/nature07385.
  2. Brennan CW, Verhaak RGW, McKenna A, et al. The somatic genomic landscape of glioblastoma. Cell, 2013, 155:462–477. DOI: 10.1016/j.cell.2013.09.034.
  3. Ceccarelli M, Barthel FP, Malta TM, et al. Molecular profiling reveals biologically discrete subsets and pathways of progression in diffuse glioma. Cell, 2016, 164:550–563. DOI: 10.1016/j.cell.2015.12.028.
  4. Verhaak RGW, Hoadley KA, Purdom E, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma. Cancer Cell, 2010, 17:98–110.
  5. Scarpace L, Mikkelsen T, Cha S, et al. The Cancer Genome Atlas Glioblastoma Multiforme Collection (TCGA-GBM), Version 5. The Cancer Imaging Archive, 2016. DOI: 10.7937/K9/TCIA.2016.RNYFUYE9.
  6. Liu J, Lichtenberg T, Hoadley KA, et al. An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics. Cell, 2018.
  7. Genomic Data Commons. GDC Data Portal, TCGA-GBM project page. https://portal.gdc.cancer.gov/projects/TCGA-GBM.
  8. cBioPortal for Cancer Genomics. Glioblastoma Multiforme (TCGA, Firehose Legacy). https://www.cbioportal.org/study/summary?id=gbm_tcga.

附录 D:任务代码

下面的脚本把「取病例、取样本类型、统计可分析队列」三步串起来,用于在开题时快速判断规模。

import requests

FIELDS = "submitter_id,primary_diagnosis,samples.sample_type,demographic.vital_status"
payload = {
    "filters": {"op": "in", "content": {"field": "project.project_id", "value": ["TCGA-GBM"]}},
    "fields": FIELDS,
    "format": "JSON",
    "size": 2000,
}
cases = requests.post("https://api.gdc.cancer.gov/cases", json=payload, timeout=120).json()["data"]["hits"]
primary = sum(1 for c in cases if any(s["sample_type"] == "Primary Tumor" for s in c.get("samples", [])))
recurrent = sum(1 for c in cases if any(s["sample_type"] == "Recurrent Tumor" for s in c.get("samples", [])))
print("cases", len(cases), "primary", primary, "recurrent", recurrent)

输出三个数字,分别回答「项目有多大」「多少例有原发肿瘤样本」「多少例有复发样本」。第三个数字决定了能不能做复发相关分析,通常在最开始就该确认。

输出 含义 用途
cases 项目病例数 队列规模
primary 有原发肿瘤样本的例数 主分析队列
recurrent 有复发样本的例数 决定纵向分析是否可行

附录 E:如果只记五件事

序号 要点
一 它是 TCGA 的首个癌种项目,六百一十七例,价值在通路机制与参照
二 有影像,二百六十二例术前 MR、CT 与平片,存 TCIA 而非 GDC
三 与 TCGA-LGG 病例互斥,GBMLGG 是 2016 年的分析队列,不是第三个项目
四 按组织学入组,IDH 与 1p/19q 需自行推导,复发样本仅十三例
五 开放层免注册,读段与胚系变异走 dbGaP 受控

编辑说明

本篇核验日期为 2026-10-29,规模数字取自 GDC 项目端点与 TCIA 集合页,生物学结论取自 2008 年 Nature、2013 年 Cell 与 2016 年 Cell 三篇摘要。与任务简报的一处出入已按官方页面更正:TCGA-GBM 在 TCIA 有独立的影像集合(二百六十二例),并非无影像。GBM 与 LGG 的「重叠」也按 GDC 病例清单核实为项目层互斥、分析层合并。未做全量下载与逐条统计,分布类数字一律不写;官方页面之间若存在口径差异,正文并列呈现而不做裁定。影像按序列类型的病例数、各数据层的精确交叉覆盖未在官方页面给出,标注为未核实。修订时优先更新附录 A 与附录 B,再同步正文概览表。


相关数据集导航

以下为站内 AI-Ready 数据集百科中与本词条共享多个主题标签的相关数据集,按相关度降序排列:

  • tcga-lgg — 共享标签:基因组学与多组学 / 肿瘤学 / 脑肿瘤 / 多模态
  • upenn-gbm — 共享标签:基因组学与多组学 / 肿瘤学 / 脑肿瘤 / 病理图像
  • soft-tissue-sarcoma — 共享标签:基因组学与多组学 / 肿瘤学 / 病理图像 / 多模态
  • hnsc — 共享标签:基因组学与多组学 / 肿瘤学 / 病理图像 / 多模态
  • hancock-hnscc — 共享标签:肿瘤学 / 病理图像 / 多模态
  • icgc — 共享标签:基因组学与多组学 / 肿瘤学 / 病理图像
  • human-cell-atlas — 共享标签:基因组学与多组学 / 肿瘤学 / 病理图像
  • ivy-gap — 共享标签:肿瘤学 / 脑肿瘤 / 病理图像
  • hubmap — 共享标签:基因组学与多组学 / 肿瘤学 / 病理图像
  • pathseg — 共享标签:肿瘤学 / 病理图像 / 多模态

导航说明:本章节由全站统一标签体系自动计算生成(标签重合度算法),双向可达;点击链接可跳转至对应数据集词条。

返回 AI-Ready 数据集