解读共表达模块时要追问的三件事

第一，模块由什么驱动？ 组织样本里，免疫浸润、基质比例、红细胞污染都能形成“漂亮”的模块。可以对照细胞类型标志基因列表，或用去卷积方法估计细胞比例后再看模块是否仍存在。Oldham 等（2008）在人脑转录组中发现，许多共表达模块恰好对应神经元、少突胶质细胞、星形胶质细胞和小胶质细胞等细胞类型——这既说明共表达能揭示组织的细胞组成，也提醒我们“模块 = 调控程序”的解读需要先排除细胞组成。

第二，模块稳定吗？ 换一批样本、换 $\beta$、换最小模块大小，模块是否大体不变？WGCNA 提供模块保留性（module preservation）统计量，用于检验一个模块在独立数据集中是否重现，详见模块保守性与可重复性检验。

第三，hub 是“驱动者”吗？ 共表达 hub 是与模块整体相关最强的基因，常常是该模块最好的标志物，但不一定是上游调控者。转录因子表达量往往较低、变化幅度小，反而不容易成为共表达 hub。代码示例的第三个单元演示了这一点：模拟中模块所有基因都只是同一潜在因子的读数，hub 只是噪声最小的读数。