《Molecular Psychiatry》：Transcriptome-informed brain cartography of polygenic risk and association with brain structure in major psychiatric disorders
编辑推荐：
本综述创新性地提出了一种整合成像转录组学的框架，将全转录组关联研究（TWAS）与大脑转录组图谱结合，预测了七种主要精神障碍（ADHD、ASD、AN、BD、MDD、OCD、SCZ）的宏观脑结构异常。研究表明，在MDD、SCZ和ADHD中存在显著的转录组-解剖学空间关联，提示区域转录易感性可能在这些疾病的遗传风险解剖表达中起不同程度作用。通路富集分析揭示了MDD和SCZ以免疫相关过程为主导，而ADHD则与神经发育通路相关。这些发现强调了成像转录组学在连接精神疾病遗传结构与系统水平脑变化之间的潜力与局限。
本文介绍了一项开创性研究，旨在填补遗传风险、基因表达与大脑宏观结构变化之间关联机制的空白。研究团队引入并测试了一个新颖的成像转录组学框架，旨在绘制多基因转录风险到宏观大脑表型的空间图谱。
核心方法：构建GEDAR图谱
该研究整合了来自精神病学基因组学联盟（PGC）的七种主要精神障碍的全基因组关联研究（GWAS）汇总统计数据，包括注意缺陷/多动障碍（ADHD）、自闭症谱系障碍（ASD）、神经性厌食症（AN）、双相情感障碍（BD）、重度抑郁障碍（MDD）、强迫症（OCD）和精神分裂症（SCZ）。研究首先应用S-MultiXcan方法，利用来自基因型-组织表达（GTEx）数据库中11个脑组织的预测模型，对GWAS数据进行全转录组关联研究（TWAS），以识别与每种障碍相关的、遗传预测的差异表达基因。随后，将TWAS推导出的基因权重与来自艾伦人脑图谱（AHBA）的空间转录组数据相结合，生成了一种名为“基于基因表达的障碍关联风险”（GEDAR）的脑图谱。GEDAR图谱量化了与每种障碍多基因结构相关的遗传预测转录失调是如何在成人大脑的皮层和皮层下区域分布的。计算公式为：GEDARj= ∑iwi· Gei,j，其中wi是基因i的TWAS Z分数权重，Gei,j是基因i在脑区j的表达量。
接下来，研究者将计算得到的GEDAR图谱与来自ENIGMA联盟的大规模病例对照荟萃分析的结构性脑图（即皮层厚度和皮层下体积的科恩d值效应图）进行比较。他们计算了区域GEDAR分数与ENIGMA效应大小之间的斯皮尔曼等级相关性，并使用球形旋转置换法（n=1000次）来校正皮层数据的空间自相关，以评估遗传推断的转录易感性与观察到的脑结构改变之间的空间对应关系。研究还探讨了这种转录组-解剖学关联的强度是否与双生子研究得出的遗传力估计值存在系统性的变化，并对显示显著GEDAR-ENIGMA关联的障碍进行了通路富集分析，以确定相关的生物过程。
主要发现：转录与解剖的时空交会
TWAS分析预测了数量相近的差异表达基因（约2759-2768个）。各障碍表现出特征性的U型分布模式，反映了上调（正Z分数）和下调（负Z分数）基因的平衡分布。其中，AN的SPINK8基因下调最显著，OCD以BRD2下调为特征，MDD的NEGR1上调最显著，而SCZ则显示出涉及NOTCH4上调以及BRD2、AGER、HLA-DMB、HLA-DMA、CCHCR1等多个基因下调的显著聚类。
空间相关性分析揭示了具有障碍和脑区特异性的关联模式：
- 1.
合并上下调基因：在皮层区域，MDD在所有基因纳入阈值（前10%、5%、1%）下均显示显著正相关（例如，前10%阈值下ρ=0.529，pspin=0.005），BD在前5%阈值下显示显著正相关（ρ=0.456，pspin=0.01）。在皮层下区域，ADHD在前10%阈值下显示极强的正相关（ρ=0.919，pspin=0.002），SCZ在前10%和5%阈值下显示显著正相关。
- 2.
仅上调基因：在皮层区域，仅MDD在前1%阈值下显示显著正相关；在皮层下区域，仅ADHD在前1%阈值下显示显著负相关。
- 3.
仅下调基因：在皮层区域，ADHD和MDD在多个阈值下显示显著负相关；在皮层下区域，SCZ在前10%和5%阈值下显示显著负相关。
交叉验证与生物通路解析
作为敏感性分析，研究还使用了来自死后脑样本的实证差异表达基因（基于Sadeghi等人的荟萃分析）构建GEDAR图谱。结果发现，仅MDD在皮层和皮层下区域均显示出一致的显著负相关，而ASD、BD和SCZ则未发现显著关联。这进一步支持了MDD中转录组变化与脑结构改变在空间上紧密相连的结论。
通路富集分析为理解这些关联的生物学基础提供了线索：
- •
MDD：下调基因富集了35条通路，其中许多与免疫应答相关，特别是主要组织相容性复合体（MHC）相关通路（如蛋白质复合物结合、抗原处理等），凸显了免疫失调在MDD病理生理学及其相关脑结构改变中的核心作用。
- •
SCZ：下调基因富集了53条通路，同样强烈指向免疫活动（MHC复合体、抗原结合、适应性免疫应答），同时也涉及认知相关的通路（如MAPK/ERK介导的突触可塑性）。
- •
ADHD：下调基因富集的通路与神经发育过程高度相关，包括“神经系统发育的调节”、“神经发生调节”等。在前5%阈值下，连环蛋白（catenin）复合体通路显著富集，该复合体在细胞粘附和Wnt信号传导中起关键作用，对建立和维持神经元连接至关重要。
超越遗传力：对异质性的深入探讨
一个有趣的发现是，GEDAR-ENIGMA空间关联的强度并未简单地与障碍的遗传力估计值成正比。尽管SCZ和ADHD的遗传力较高（约75-80%），但MDD（遗传力约35-45%）却表现出最稳健、最广泛的关联。分析甚至显示，在某些情况下（如皮层区域仅上调基因前1%阈值），遗传力与空间关联强度呈显著负相关。这表明，除了遗传贡献大小外，发育时间、环境相互作用、转录后调控、表观遗传修饰以及疾病特异性分子结构等多种因素，共同调节了遗传风险通过转录变化影响大脑解剖结构的程度。
研究意义与未来方向
本研究提出的GEDAR框架，类似于多基因风险评分（PRS）的空间延伸版本，为在系统水平上绘制多基因转录风险图谱提供了一种新颖且可推广的方法。它首次在多个精神障碍中系统评估了遗传预测的转录失调与观察到的脑结构变化之间的空间对应关系。研究发现强调了在MDD、SCZ和ADHD中，区域性脑结构改变部分反映了遗传风险基因表达失调的空间模式，并且不同障碍涉及主导的生物学通路存在差异（如MDD/SCZ的免疫通路 vs ADHD的神经发育通路）。
当然，研究也存在一些局限性，包括：GEDAR基于预测而非实证的转录变化；AHBA数据来自少量神经典型成人供体，代表性和发育普适性有限；依赖汇总级数据，无法评估个体异质性；使用成年期转录组参考模型，可能无法捕捉早期发病障碍的关键发育期转录模式；TWAS分析平均了所有脑组织面板，而非针对特定脑区进行优化。
总之，这项研究为理解精神疾病中遗传风险如何通过分子水平的转录调控，最终表现为宏观脑结构异常提供了初步证据和有力工具。未来，整合单细胞转录组学、发育时空表达图谱、多组学数据以及个体水平匹配的影像与转录组数据，将有助于构建更全面、更动态的“基因组-转录组-脑解剖-行为”疾病模型，最终推动针对特定生物学通路的精准治疗策略的开发。
