整合转录组学-神经影像分析揭示缺血性卒中后半球间功能解耦的多基因相关特征

《Frontiers in Cellular Neuroscience》：Integrative transcriptomic- neuroimaging analysis reveals polygenic correlates of interhemispheric functional decoupling in ischemic stroke

　　

背景：功能磁共振成像（fMRI）已揭示卒中患者存在异常脑活动模式，但功能同伦（functional homotopy，即双侧同源脑区自发活动的同步性）背后的遗传相关特征仍鲜有表征。本研究探讨卒中患者体素镜像同伦连接（voxel-mirrored homotopic connectivity，VMHC）异常的遗传基础。方法：研究人员分析50例卒中患者与50例健康对照（healthy control，HC）的静息态fMRI数据以量化VMHC；利用艾伦人脑图谱（Allen Human Brain Atlas，AHBA）建立空间转录组-神经影像相关，识别VMHC相关基因；转录组分析联合通路中心功能注释（DAVID）与蛋白质-蛋白质相互作用（protein-protein interaction，PPI）网络建模（STRING v12.0）。结果：卒中患者在直回、颞上回、枕中回、楔叶及右距状沟/左后扣带回（p < 0.05，高斯随机场校正）VMHC显著降低；VMHC改变分别与1198个基因呈正负相关；转录组分析显示突触囊泡运输、线粒体能量代谢及神经炎症相关通路显著富集；PPI映射揭示多层网络，枢纽基因包括BRCA1、CDK9、ACTB及ATP6V1A/F，参与转录调控、细胞骨架动力学与囊泡酸化。结论：本多模态整合研究阐明卒中后VMHC异常的多基因相关特征，表明半球间协调可能依赖功能多样基因簇的协同互作；研究为理解卒中后神经网络重组提供分子框架，并建立功能神经影像表型与基因表达的联系，但所有关联仍为相关性，需机制验证。

研究背景与立项依据

缺血性卒中占卒中总数约80%，常累及基底节区，其后的运动、感觉、认知缺陷不仅源于原发灶，也与远隔网络紊乱有关。静息态功能磁共振成像（resting-state fMRI，rs-fMRI）揭示脑自发活动组织为大规模网络，而功能同伦（functional homotopy）作为内在功能架构的基本属性，可通过体素镜像同伦连接（voxel-mirrored homotopic connectivity，VMHC）量化双侧镜像体素BOLD信号的时间相关性。既往研究提示卒中可破坏默认模式网络与感觉运动网络的半球间同伦，但VMHC异常的分子相关尚不清。转录组-影像融合分析借助艾伦人脑图谱（Allen Human Brain Atlas，AHBA）可提供空间分辨的基因表达谱，为揭示脑功能异常的遗传基础提供途径。鉴于卒中多基因风险已被证实，而卒中后VMHC异常的分子基础仍待阐明，研究人员开展本项多模态整合研究，论文发表于《Frontiers in Cellular Neuroscience》。

主要关键技术方法

研究人员纳入江西省人民医院50例首发急性缺血性卒中患者（病灶位于基底节或半卵圆中心）与50例年龄、性别、教育匹配的健康对照（HC）。采集3T rs-fMRI与结构像，经DPABI预处理后计算VMHC，组间采用两样本t检验（体素p<0.001，高斯随机场簇p<0.05）。影像-转录组关联采用AHBA六例尸检脑表达矩阵（abagen处理），以全脑未校正VMHC对比t值为响应变量，基因表达为预测变量，行偏最小二乘（partial least squares，PLS）回归并施以旋转检验（spin-test）1000次置换与十折交叉验证。获选基因分正/负相关集各1198个，行DAVID v6.8的GO/KEGG富集与STRING v12.0（置信≥0.9）PPI网络构建，枢纽基因取节点度前10%，并在人脑转录组图谱中做时空表达映射。

研究结果

VMHC组间比较：卒中组在双侧直回、颞上回、枕中回、楔叶及右距状沟/左后扣带回VMHC显著低于HC。区域平均VMHC与NIHSS总分偏相关为负（部分r=-0.49，p=0.005），其中右距状沟/左后扣带回与NIHSS及意识、凝视分项最强负相关，颞上回与语言分项负相关，两区域合解释NIHSS方差42%。

VMHC相关基因：PLS提取4个显著潜分量，累积R²=0.42，筛得正、负相关基因各1198个（p<0.001）。

基因富集结果：正相关集富集于硫酸软骨素生物合成、翻译起始、RNA聚合酶II转录延伸、凋亡与蛋白水解等GO项及醛固酮合成分泌、凋亡等KEGG通路；负相关集富集于突触囊泡腔酸化、液泡酸化、线粒体电子传递、周期蛋白依赖激酶调控等GO项及氧化磷酸化、糖胺聚糖生物合成等KEGG通路。

PPI网络与枢纽基因：正网络1142节点、负网络1131节点，均具小世界性（聚类系数高于随机p<0.001）。正集枢纽含BRCA1、CDK9、ERCC6；负集枢纽含ACTB、ATP6V1A、ATP6V1F。时空表达显示这些基因在颞中回（BA21）、顶下小叶（BA40）、额上回（BA8）等缺血易感皮质峰值表达。

讨论部分总结

研究人员指出，皮质VMHC降低发生于基底节/半卵圆深部病灶，符合经丘脑皮质或纹状体皮质通路的远隔失联（diaschisis）。正相关基因指向细胞外基质重塑与转录调控，如硫酸软骨素蛋白聚糖（CSPGs）限制可塑性、RNA聚合酶II功能关联发育性认知障碍；负相关基因指向缺血后ATP耗竭、线粒体Ca²⁺超载、线粒体通透转换孔开放及活性氧增多，V-ATPase介导的突触囊泡酸化受损可削弱半球间信号效率。枢纽基因中BRCA1可能经DNA修复与抗氧化应激补偿连接性，CDK9经RNAPII延伸调控p53转录加重缺血代谢障碍，ERCC6联结转录与修复，ACTB影响eNOS/NO通路，ATP6V1A/F决定囊泡充填与释放。所有关联为相关性、源于健康脑转录基线，且AHBA右半球样本受限致左侧偏倚，横断面急性期设计亦限动态推论。

结论翻译

本多模态整合研究识别与卒中患者半球间功能同伦连接（VMHC）异常相关的多基因特征。颞叶、枕叶与扣带回区VMHC降低；共计2396个基因与VMHC改变相关，富集于突触囊泡运输、线粒体能量代谢与神经炎症；PPI网络凸显BRCA1、CDK9、ACTB及ATP6V1A/F等枢纽基因作为潜在分子相关。研究为宏观连接中断与特定分子通路建立框架，为卒中恢复的机制与治疗研究提供可检验假设。所有基因-影像关联均为相关性，需独立验证。