《Brain Research》：Integration of multi-omics Mendelian randomization, neuroimaging mediation, and functional validation identifies PYGL as a neuroprotective causal gene for stroke prognosis
编辑推荐：
卒中仍是全球长期残疾的主要原因，即使考虑了初始严重程度、年龄和合并症等已确立的临床预测因素后，患者的功能恢复仍存在巨大差异。这种无法解释的异质性激发了对卒中预后遗传架构日益增长的兴趣。全基因组关联研究（GWAS）已鉴定出与改良Rankin量表（mRS）评估的功
卒中仍是全球长期残疾的主要原因，即使考虑了初始严重程度、年龄和合并症等已确立的临床预测因素后，患者的功能恢复仍存在巨大差异。这种无法解释的异质性激发了对卒中预后遗传架构日益增长的兴趣。全基因组关联研究（GWAS）已鉴定出与改良Rankin量表（mRS）评估的功能结局相关的遗传变异，包括GISCOME联盟报告的LOC105372028附近的rs1842681多态性。然而，将这些位点与恢复轨迹联系起来的生物学机制仍不清楚。一个根本性挑战是，大多数GWAS信号位于非编码区域，而这些变异影响卒中结局所通过哪些基因在很大程度上是未知的。弥合这一鸿沟需要能够从遗传关联数据中优先排序因果基因的分析框架。基于汇总数据的孟德尔随机化（SMR）通过将GWAS汇总统计与表达数量性状位点（eQTL）数据整合，以鉴定基因表达水平通过多效性或因果性与复杂性状关联的基因，从而应对了这一挑战。与通常评估单一暴露的常规孟德尔随机化（MR）方法不同，SMR利用异质性依赖工具（HEIDI）检验以区分真正的因果关联与由连锁不平衡（LD）驱动的虚假信号。对于脑相关性状，近期BrainMeta资源的可用性——包含来自2,865个欧洲血统皮层样本的eQTL数据——为专门绘制中枢神经系统内部基因表达的遗传调控图谱提供了机遇。值得注意的是，此前尚无研究将（i）来自人脑单核RNA测序的细胞类型特异性eQTL数据、（ii）跨多个脑区的组织类型特异性eQTL验证、（iii）结合脑影像衍生表型（IDP）的两步MR中介分析以及（iv）体外功能验证整合在单一分析框架中。研究人员的"一体化"多组学MR设计通过对来自遗传工具、细胞类型特异性基因调控、脑结构影像和实验验证的证据进行三角互证，弥补了这一空白。此外，PYGL此前未曾被任何全基因组或全转录组方法认定与卒中预后相关的神经保护基因。在本研究中，研究人员采用了一套全面的分析流程来鉴定卒中预后的因果基因。研究人员首先进行了SMR分析，将卒中结局GWAS数据与来自BrainMeta的脑eQTL整合，使用HEIDI检验排除异质性驱动的关联。通过严格显著性阈值的候选基因随后进行贝叶斯共定位分析，以评估GWAS和eQTL信号是否共享一个共同的因果变异。为剖析解剖和细胞特异性，研究人员将分析扩展到跨多个脑区的组织类型MR，以及检验神经元、少突胶质细胞、星形胶质细胞和小胶质细胞中表达的细胞类型特异性MR。研究人员进一步进行了中介分析，以探索脑影像衍生表型是否连接基因表达与功能结局。最后，研究人员在海马神经元细胞中利用氧糖剥夺/再灌注（OGD/R）实验对优先排序的神经保护基因PYGL进行了功能验证。
**论文解读：整合多组学孟德尔随机化、神经影像中介与功能验证鉴定 PYGL 作为卒中预后的神经保护因果基因**
**一、研究背景与问题**
卒中（Stroke）是全球范围内导致长期残疾的首要原因之一。即使在考虑了初始卒中严重程度、年龄和合并症等已确立的临床预测因素后，患者的功能恢复仍表现出巨大的个体间差异。这种无法由已知临床变量解释的异质性，促使研究人员日益关注卒中预后的遗传架构。全基因组关联研究（GWAS）已鉴定出与改良Rankin量表（mRS，一种衡量卒中后残疾和依赖程度的标准化评分工具）评估的功能结局相关的遗传变异，例如GISCOME（Genetics of Ischemic Stroke Functional Outcome）联盟报告的位于LOC105372028附近的rs1842681多态性。然而，将这些遗传位点与恢复轨迹联系起来的生物学机制仍基本未知。一个根本性挑战在于，绝大多数GWAS信号位于基因组的非编码区域，而这些遗传变异究竟通过调控哪些基因来影响卒中结局，在很大程度上仍未明确。弥合从遗传关联到因果基因这一鸿沟，需要能够从海量遗传数据中优先鉴定因果基因的分析框架。
在此背景下，基于汇总数据的孟德尔随机化（SMR）方法应运而生。SMR通过整合GWAS汇总统计数据与表达数量性状位点（eQTL，即影响基因表达水平的遗传变异）数据，鉴定基因表达水平是否通过多效性或因果途径与复杂性状相关联。相较于通常仅评估单一暴露因素的传统MR方法，SMR的核心优势在于其利用HEIDI（Heterogeneity in Dependent Instruments，异质性依赖工具）检验来区分真正的因果关联与由连锁不平衡（LD）驱动的虚假信号，从而显著提高了因果推断的可靠性。对于脑相关性状，BrainMeta联盟近期发布的资源——涵盖2,865例欧洲血统皮层样本的eQTL数据——为在中枢神经系统内部绘制基因表达遗传调控图谱提供了独特契机。
值得注意的是，在研究人员开展此项研究之前，尚没有研究在同一分析框架中同时整合以下四个层面的证据：（1）来自人脑单核RNA测序的细胞类型特异性eQTL数据；（2）跨多个脑区的组织类型特异性eQTL验证；（3）结合脑影像衍生表型（IDP）的两步MR中介分析；（4）体外细胞功能实验验证。研究人员设计的"一体化"多组学MR研究正是为了填补这一方法学空白，通过对来自遗传工具变量、细胞类型特异性基因调控、脑结构影像以及实验验证的多维证据进行三角互证（Triangulation），全面评估卒中后功能恢复的因果性可成药靶点。尤其值得关注的是，PYGL基因此前未曾被任何全基因组或全转录组方法鉴定为卒中预后的神经保护基因。
**二、研究内容与意义**
研究人员首先利用SMR分析整合卒中结局（mRS）的GWAS汇总数据与BrainMeta的脑eQTL数据，并应用HEIDI检验排除由异质性驱动的虚假关联。通过严格统计学阈值的候选基因随后接受贝叶斯共定位分析（Colocalization Analysis），以评估GWAS信号与eQTL信号是否共享同一个因果变异。为进一步解析基因效应的解剖学与细胞学特异性，研究人员将分析扩展至基于GTEx（Genotype-Tissue Expression，基因型-组织表达项目）数据的跨多脑区组织类型MR，以及基于Bryois等人单核RNA测序数据的细胞类型特异性MR，系统检验了基因在兴奋性神经元、抑制性神经元、少突胶质细胞、少突胶质前体细胞（OPC）、星形胶质细胞和小胶质细胞中表达水平的因果效应。此外，研究人员在SMR和MR分析中观察到基因对mRS的直接因果效应后，进一步利用英国生物银行（UK Biobank）的大样本GWAS数据进行了两步MR中介分析，以探究脑影像衍生表型（IDP）是否在基因表达与功能结局之间充当中介变量。最后，研究人员在小鼠海马神经元细胞系HT22中，通过氧糖剥夺/再灌注（OGD/R，一种体外模拟缺血性卒中的经典模型）实验，对优先鉴定的神经保护基因PYGL进行了功能层面的实验验证。
**三、主要关键技术方法**
本研究的技术方法体系涵盖以下核心环节：（1）暴露数据来源：卒中预后GWAS汇总数据来自GISCOME联盟，包含6,021例欧洲血统缺血性卒中患者，主要结局为卒中后约3个月的mRS评分（二分类：0–2分为结局良好，3–6分为结局不良）；（2）基因表达工具变量：脑eQTL数据来自BrainMeta联盟（2,865例皮层样本），组织特异性eQTL来自GTEx V10，细胞类型特异性eQTL来自Bryois等人对192例人死后脑样本的单核RNA测序研究，涵盖前额叶皮层、颞叶皮层和白质；（3）因果推断方法：采用SMR分析联合HEIDI检验，并使用coloc R包进行贝叶斯共定位分析以评估GWAS与eQTL共享因果变异的证据强度；（4）中介分析：脑影像衍生表型（IDP）的GWAS汇总数据来自英国生物银行，采用两步MR中介分析方法估计中介效应；（5）功能验证：利用HT22小鼠海马神经元细胞系构建OGD/R模型，通过质粒过表达和小干扰RNA（siRNA）敲低PYGL，采用TUNEL染色和流式细胞术评估细胞凋亡。
**四、研究结果**
**（1）SMR分析鉴定卒中结局的推定因果基因**
SMR分析共鉴定出6个通过显著性阈值（P
SMR < 0.001）且无异质性证据（P
HEIDI > 0.05）的推定因果基因。其中，PYGL的SMR效应估计为负值（b
SMR = ?0.39, P
SMR = 2.06 × 10
?4），提示PYGL表达水平越高，卒中不良结局风险越低，即具有保护效应。同样，GSTA1（b
SMR = ?0.19, P
SMR = 7.84 × 10
?4）和PTGIS（b
SMR = ?0.55, P
SMR = 9.15 × 10
?4）也被鉴定为保护性因素。相反，NFATC1表现出正向因果效应（b
SMR = 0.78, P
SMR = 7.39 × 10
?4），提示NFATC1表达升高与不良预后风险增加相关。此外，SULT1B1和非编码转录本AL358334.2也被鉴定为显著结果。
**（2）共定位分析支持基因表达与卒中结局之间的共享因果变异**
贝叶斯共定位分析显示，NFATC1具有最高的共享因果变异概率（PP.H4 = 0.901），PYGL（PP.H4 = 0.859）、PTGIS（PP.H4 = 0.837）和GSTA1（PP.H4 = 0.802）也均显示出强有力的共定位证据。在所有情况下，共定位概率（PP.H4）均大幅超过独立因果变异概率（PP.H3 < 0.09），强烈提示观察到的关联是由同一遗传信号驱动，而非LD区域内的不同变异所致。
**（3）组织类型特异性MR分析**
基于GTEx的组织类型特异性MR分析（图2）揭示，PYGL在多个脑区——包括小脑半球、小脑、皮层和额叶皮层——的高表达与不良结局风险的显著降低相关（比值比OR范围为0.577至0.74，均P < 0.005）。PYGL在壳核中的表达在主要二分类mRS分析中也显示出保护效应（OR = 0.379, P = 0.0001），但该估计值在有序mRS敏感性分析中未保持显著。对于GSTA1，额叶皮层中表达升高与较低的不良预后风险相关（OR = 0.835, P = 0.0038）。相比之下，PTGIS和NFATC1在大多数脑区未显示出显著关联。
**（4）细胞类型特异性MR分析**
细胞类型特异性MR分析（图3）进一步揭示了基因效应的细胞层面特异性。NFATC1在少突胶质细胞中的表达升高与不良预后风险显著增加相关（OR = 1.69, P = 1.0 × 10
?4），而其在天文胶质细胞、小胶质细胞和OPC中的表达未见显著关联。值得注意的是，PYGL在兴奋性神经元中的表达与不良结局风险显著降低相关（OR = 0.62, P = 0.012），明确了其保护效应定位于兴奋性神经元。PTGIS在兴奋性神经元、抑制性神经元和OPC中均未观察到显著关联。
**（5）中介分析结果**
中介分析（图4）揭示了连接基因表达与功能结局的脑影像中介通路。对于PTGIS，右侧矢状层白质纤维束完整性（IDP dMRI ProbtrackX 12 ptr 1）部分中介了基因表达对mRS的效应，间接效应显著（β = ?0.0324, P = 0.015），直接效应仍然显著（β = ?0.455, P = 0.0045）。对于GSTA1，多个影像特征显示出显著的中介效应，包括功能连接（ICA100 cdqc 849）、白质微结构（TBSS ICVF 右侧矢状层）以及胼胝体完整性（TBSS L3 膝部）。对于PYGL，BA45区皮层厚度（BA-exvivo lh thickness BA45）显著中介了基因表达与mRS的关联（间接效应β = 0.0089, P = 0.011），直接效应仍显著（β = ?0.364, P = 0.0051）。此外，左侧尾状核体积（IDP T1 FAST ROIs L caudate）的间接效应显著（β = ?0.0205, P = 0.036），但直接效应不显著（β = ?0.074, P = 0.088），提示该通路以完全中介为主。
**（6）PYGL在神经元保护中的功能验证**
在HT22细胞OGD/R模型中，免疫荧光染色证实PYGL蛋白主要定位于神经元内，与神经元标志物NeuN呈现明确共定位。OGD/R刺激显著增加了细胞凋亡率。PYGL过表达（OE-PYGL）显著减轻了OGD/R诱导的凋亡，表现为TUNEL阳性细胞比例降低及流式细胞术检测的凋亡率下降。相反，通过siRNA沉默PYGL（si-PYGL）则加剧了细胞损伤，导致凋亡细胞比例进一步升高。这些实验发现为MR结果提供了功能性证据，确认PYGL是一种能够减轻缺血性神经元损伤的神经保护因子。
**五、讨论与结论**
本研究通过整合多组学MR分析、神经影像中介分析和体外功能实验，系统地鉴定并验证了与卒中预后（mRS）相关的因果基因。SMR和共定位分析稳健地优先排序了四个关键基因：PYGL、GSTA1和PTGIS为保护性因子，NFATC1为不良功能结局的风险因子。研究人员将分析扩展至组织和细胞类型特异性层面，并描绘了中间脑影像表型。
细胞类型特异性MR鉴定出PYGL的保护效应主要由其在兴奋性神经元中的表达所驱动，这一发现得到了体外功能实验的证实——PYGL过表达显著减轻了HT22细胞中OGD/R诱导的凋亡。PYGL编码糖原磷酸化酶的肝脏亚型，催化糖原分解的限速步骤。尽管PYGL被经典归类为肝脏亚型，但研究人员的免疫荧光数据证实了神经元中PYGL蛋白的表达。脑糖原是缺血应激期间的关键能量储备，糖原分解受损会加剧再灌注后的神经元损伤。兴奋性神经元因代谢需求高且特别依赖糖原来源的葡萄糖，因而对缺血性损伤尤为脆弱。中介分析进一步揭示，PYGL的保护效应部分通过保留左侧尾状核体积和皮层厚度来中介，提示遗传决定的高PYGL表达可能赋予兴奋性神经元对抗缺血损伤的代谢弹性，从而保护关键脑结构并促进更好的临床恢复。
GSTA1被鉴定为保护性因子，与谷胱甘肽代谢在脑内氧化应激反应中的确立作用一致。组织特异性分析表明，GSTA1在额叶皮层的表达与改善的预后显著相关。中介分析证明，GSTA1部分通过调节白质微结构和功能连接发挥其保护效应。PTGIS催化前列腺素H
2向具有血管舒张、抗炎和细胞保护特性的前列环素（PGI
2）的转化。中介分析揭示了白质纤维束完整性中介了PTGIS表达与功能结局的关联。
相反，NFATC1被鉴定为不良卒中结局的显著风险因子。NFATC1是一种对钙依赖性基因调控至关重要的转录因子。细胞类型特异性分析揭示NFATC1的有害效应特异于少突胶质细胞，提示少突胶质细胞中NFATC1信号的失调可能损害卒中后白质修复。少突胶质细胞负责轴突的髓鞘形成，其功能障碍导致白质损伤和不良恢复。共定位分析显示NFATC1具有最高的共享因果变异概率（PP.H4 = 0.901），强烈支持NFATC1表达与卒中结局之间的共享因果变异。
**研究结论部分**：本研究为卒中预后背后的因果基因调控网络提供了全面的证据。PYGL被鉴定为一种神经保护因子，并在细胞模型中得到了验证，这使其被提名为促进缺血性卒中后恢复的有前景的治疗靶点。组织和细胞类型特异性分析突显了在解释遗传关联时考虑解剖和细胞背景的重要性。中介分析建立了通过神经影像表型将基因表达与功能结局联系起来的机制通路，为未来的干预研究提供了潜在的生物标志物。
