A significant enrichment that is not: spatial nulls, co-expression, and the imaging transcriptomics of EEG alpha-power genetics

技术摘要：“显著富集却并非如此：空间零假设、共表达以及 EEG $\alpha$ 波功率遗传学的影像转录组学”

问题陈述
本研究旨在解决影像转录组学中的一个关键空白：即特定神经生理特征的遗传架构是否在产生该特征的大脑区域内具有空间组织性。具体而言，作者调查了与 EEG $\alpha$ 波功率相关的基因是否在已识别出的 $\alpha$ 节律生成器（41 个皮层区域，见 Tabarelli et al., 2022）中过度表达。

在方法论上，本文针对该领域标准实践中的一个已知局限性提出了挑战。富集分析通常依赖于“旋转测试”（空间置换）来控制空间自相关。然而，这种方法并未考虑到集合内的基因通常是共表达的，这会产生独立于目标特征之外的空间结构。作者假设，在空间零假设下表现显著的结果，可能仅仅是该皮层领地内任意基因集的通用属性，而非目标特征特有的信号。

研究方法
本研究采用了一个多阶段流程，整合了 GWAS 汇总统计数据、区域基因表达以及严格的零模型构建：

1. 数据来源：

   • 遗传学： 来自 ENIGMA-EEG 的汇总统计数据（N $\approx$ 6,000–7,700），涵盖六种表型：中央 $\alpha$ 波功率（主要研究对象）、枕叶 $\alpha$ 波功率、枕叶 $\alpha$ 波峰值频率，以及 $\theta$、$\beta$ 和 $\delta$ 波功率。
   • 转录组学： 来自 Allen 人脑图谱（AHBA）的区域基因表达数据，通过 abagen 进行处理，并映射至 Glasser HCP-MMP1.0 图谱（360 个区域）；同时进行了基于 Schaefer-100 和 Yan-600 分区的敏感性分析。
   • 目标区域： 由 Tabarelli 等人（2022）利用 MEG 源重建定义的 41 个皮层 $\alpha$ 源区域（双侧，82 个分区）。

2. 基于基因的分析：

   • 使用 MAGMA 计算基因水平的关联统计量。
   • 每个区域定义了两种富集评分：连续评分（区域表达量与全基因组基因 z 统计量之间的 Spearman 相关系数）和 Top-100 评分（按 p 值排序的前 100 个基因的平均表达量）。

3. 零模型与验证：

   • 空间零假设（旋转测试）： 在皮层球面上对分区质心进行 10,000 次空间旋转，以保留空间自相关性。
   • 基因集零假设： 使用 10,000 个大小匹配的随机基因集，以测试观察到的富集是特定于 $\alpha$ 相关基因，还是任何共表达基因集的共有属性。
   • 阳性对照： 针对 AHBA 表达矩阵的第一主成分（PC1，一种已知的皮层梯度）测试该流程的敏感性。
   • 替代验证： 通过交叉检查自定义的旋转实现与 brainsmash 替代物的一致性。

关键结果

1. 初步空间显著性： 在仅使用领域标准的旋转测试下，研究初步显示支持该假设。中央 $\alpha$ 波功率在 $\alpha$ 源区域显示出显著富集（连续 $p_{spin} = 0.022$；Top-100 p_{spin} = 0.030}）。这一结果得到了 brainsmash 替代物（$p = 0.018$）和阳性对照（$p_{spin} = 2 \times 10^{-4}$）的证实。
2. 缺乏特异性： 富集现象并非频带特有的。$\theta$、$\beta$ 和 $\delta$ 波功率在 $\alpha$ 源区域也表现出相当或更强的富集。此外，该结果在 $\alpha$ 表型之间未能复现：枕叶 $\alpha$ 波功率和 $\alpha$ 波峰值频率均不显著。
3. 基因集零假设的消解： 当针对大小匹配的随机基因集进行评估时，$\alpha$ Top-100 集合与随机集合无异（$p_{geneset} = 0.33$）。观察到的富集是枕顶额领地内任意基因集的通用属性，而这种属性在仅考虑空间的零假设中是不可见的。
4. 缺乏稳健性： 在 24 单元分析网格中，没有任何测试能通过错误发现率（FDR）校正（最小 $q = 0.127$）。结果对分区的改变并不稳健；它仅在原始 Glasser 空间中具有名义上的显著性，而在 Schaefer-100 和 Yan-600 中则消失或减弱。

核心贡献

• 实质性发现： 本研究提供证据表明，EEG $\alpha$ 波功率的遗传学在 $\alpha$ 节律的皮层生成器中并不存在区域特异性的转录组特征。任何微弱的富集倾向在不同频带间是共享的，这更有可能是由 EEG 功率的通用遗传相关性驱动，而非 $\alpha$ 特有的分子架构。
• 方法论演示： 本工作作为一个“实操演示”，证明了仅纠正空间自相关是必要但不充分的。它表明，一个结果可以同时具备旋转显著性、替代物验证支持以及机制上的合理性，但完全可以被基因集的共表达结构所解释。
• 提出的标准： 作者认为，影像转录组学的富集声明必须在报告空间零假设的同时，也报告基因集零假设，以避免由共表达混淆导致的假阳性。

意义与主张
本文声称，$\alpha$ 相关基因在 $\alpha$ 生成器中的所谓富集，实际上是由于共表达结构导致的统计伪影，而非生物学真相。作者强调，该“阴性结果”是稳健的，因为研究设计本身就旨在抵御过度解读，利用了内部功率匹配对照、阳性对照和多种零模型。

其方法论结论是主要的普适性启示：仅仅依赖空间零假设（如旋转测试）可能会导致将皮层的一般性表达模式错误地归因于特定的特征遗传学。作者总结道，未来对于仅通过空间零假设检验的富集声明，应将其视为未解释的结果，且增加基因集零假设的成本，相对于导致文献充斥着“皮层的属性”而非“基因的属性”这一风险而言，几乎可以忽略不计。

该研究并非声称发现了新的 $\alpha$ 波功率生物学机制，而是通过严谨地测试并驳斥了关于其区域转录组组织结构的特定假设，同时揭示了该领域的一个方法论漏洞。