本文为 浙大二院神经内科 刘恺鸣 教授和放射科 孙建忠 教授团队发表于 The Journal of Headache and Pain 的工作,题为 “Cortical network dysfunction in migraine:linking genes, metabolism and clinical disability”。DOI:10.1186/s10194-025-02182-1。
Abstract
Background:偏头痛是一种涉及遗传易感性和环境因素的多因素神经血管疾病。尽管已有广泛研究,但遗传因素如何导致偏头痛患者大脑功能障碍和易感性的确切机制尚不明确。
Methods:本研究纳入139名偏头痛患者和65名健康对照者。我们基于局部形态测量指标构建了形态学逆散度(MIND)网络,以探究大脑结构改变、认知功能、代谢分布与基因表达谱之间的关系。MIND网络进一步用于预测偏头痛患者的个体功能残疾程度。
Results:转录组分析鉴定出124个与突触信号传导和神经递质转运相关的基因,这些基因主要在内皮细胞、少突胶质细胞以及兴奋性/抑制性神经元中表达。代谢失调可解释5.61%的变异,其中氧葡萄糖指数(OGI)是主要驱动因素。偏头痛患者在与语言、扣带-岛盖和躯体运动系统相关的MIND网络中出现显著紊乱,这些紊乱与认知域(如记忆、疼痛、语言)相关,并能够预测功能残疾(r = 0.23,Pperm = 0.015)。
Conclusion:本研究揭示了偏头痛中宏观大脑结构异常与微观转录模式之间的新联系。这些发现增进了我们对偏头痛神经生物学的理解,并为未来干预提供了潜在的治疗靶点。
Keywords: Migraine, Morphometric inverse divergence network, Transcriptional patterns, Cell type-specific signatures
Neurobiol Dis.:转录组与神经递质与脊髓小脑性共济失调结构和功能空间变异模式的关联性研究
BMC Medicine: 首发精神分裂症患者炎症标志物转录水平与短期脑结构变化的关系
Trends Neurosci.:使用fMRI数据提取大脑功能网络方法综述
Molecular Psychiatry:通过一种新的大脑结构系统方法解决精神分裂症的异质性:个体化结构协方差网络分析
与帕金森病灰质损失相关的扣带回网络显示健康大脑中胆碱能基因的高表达
Commun Biol.:脑连接组梯度在终末期肾病患者中的功能异常及其与临床表型及认知缺陷的关联
Introduction
偏头痛的典型特征为反复发作的搏动性头痛和感觉过敏。皮质扩散性抑制和三叉神经血管系统激活等病理生理机制可部分解释该疾病,但我们对这种复杂脑功能紊乱的认识仍十分有限。近期神经影像研究提供了偏头痛患者灰质改变的有力证据,包括颞叶与额叶皮质的灰质体积变化以及躯体感觉区皮质厚度异常。然而,这些单变量研究方法常聚焦孤立结构指标,难以捕捉偏头痛中复杂的多变量皮质重组模式,导致研究结果时有矛盾。这一局限使得宏观结构改变与其潜在神经生物学基础之间的可能联系难以显现。
连接组图谱研究的最新进展为弥合这一鸿沟提供了新机遇。形态测量逆散度(MIND)方法通过量化皮质厚度、曲率和沟深等多特征在顶点层面的多变量相似性,实现了研究范式的革新。与传统依赖简化区域汇总的形态相似性网络不同,MIND在保持皮质组织丰富空间信息的同时,展现出更优异的生物学效度——其与金标准轴突连接测量的强相关性、与既定细胞构筑边界的一致性,以及与区域基因共表达模式的稳健关联均印证了这一点。在此框架下,较高的MIND值意味着皮质区域间更大的形态相似性,通常反映更同步的皮质组织和更强的结构整合;而较低的MIND值则表征更大的形态差异性和更弱的整合度。
通过艾伦人脑图谱(AHBA)等资源,将此类先进连接组学与分子谱分析相结合已日趋可行。AHBA提供的高分辨率全皮质转录组图谱,使研究者能够探索基因表达的区域差异如何影响(或受制于)宏观网络改变。该方法已为多发性硬化至重度抑郁症等多种神经疾病研究带来新见解,揭示了与网络功能障碍相关的细胞类型特异性转录特征。将MIND的精确结构定位与AHBA分子数据相结合应用于偏头痛研究,有望阐释特定皮质网络为何更易受偏头痛相关改变影响,以及遗传风险变异如何通过特定转录程序影响这种易感性。
基于上述方法的协同优势,本研究采用多尺度方法探究偏头痛皮质网络功能障碍的神经生物学基础。主要目标包括:(1)运用MIND方法在大规模偏头痛患者队列中系统绘制相较于健康对照的皮质结构网络异常模式;(2)通过验证“偏头痛易感区具有独特分子与代谢特征”的假说,探索这些结构改变的生物学基础。为此,我们将结构发现与AHBA皮质基因表达数据整合以识别相关转录特征和细胞类型,同时结合标准正电子发射断层扫描(PET)脑代谢图谱评估能量相关过程的作用;(3)通过检验网络水平变化与认知功能元分析图谱的空间重叠性,并测试个体MIND模式能否预测患者头痛相关功能残疾,确立这些改变的临床相关性。通过建立宏观脑结构与微观生物学及临床结局的联系,本研究旨在为偏头痛病理生理学提供更整合的理解框架。
Methods
Participants
本研究方案已获得浙江大学医学院附属第二医院医学伦理委员会的正式批准。所有参与者在入组前均签署书面知情同意书。患者招募严格遵循《国际头痛疾病分类(第三版)》(ICHD-3)中关于偏头痛的诊断标准。偏头痛患者纳入标准包括:[1]确诊偏头痛;[2]年龄范围18-70岁;[3]偏头痛病史超过一年;[4]每月头痛天数≤14天。所有参与者均适用以下全面排除标准:[1]冠状动脉疾病;[2]糖尿病;[3]高血压;[4]活动性感染性疾病;[5]慢性炎症或自身免疫性疾病;[6]高胆固醇血症;[7]严重全身性疾病;[8]当前吸烟者;[9]妊娠或哺乳期;[10]肥胖(身体质量指数≥30);[11]药物过度使用性头痛或物质成瘾。
健康对照参与者通过定向抽样从当地社区招募。所有对照者均满足与患者参与者相同的排除标准,并通过结构化筛查额外确认完全无头痛疾患或精神疾病。最终研究队列共纳入204人,包括139名偏头痛患者和65名人口统计学特征匹配(年龄与性别匹配)的健康对照者。
Fig.1展示了整合脑成像、临床数据和转录组分析以研究偏头痛皮质形态学测量网络异常的整体研究流程,具体方法学构成如下所示。
Fig. 1 The analytical workflow integrating neuroimaging and transcriptomic data. The neuroimaging pipeline (A) constructed individual 360 × 360 MIND matrices from five morphometric features, with regional MIND values obtained by averaging all connections. Case-control differences were quantified using Cohen’s d and analyzed for cognitive associations, clinical prediction, and metabolic correlates. Transcriptomic processing (B) utilized the abagen toolbox to map Allen Human Brain Atlas expression data to the HCP atlas through probe selection, normalization, and multi-stage registration. Connectome-transcriptome integration (C) employed PLS regression to identify spatial covariation between MIND alterations and gene expression, followed by functional annotation of significant gene sets through Gene Ontology enrichment and cell-type specificity analysis using established markers
Clinical assessment
所有参与者均完成了综合性医学访谈,记录人口统计学特征及家族病史。针对偏头痛患者,我们采用经过验证的偏头痛失能评估(MIDAS)问卷定量评估头痛相关功能障碍。
MRI acquisition and pre-processing
所有患者在扫描时均确认处于无发作(发作间期)状态。所有神经影像数据均使用联影医疗MR790 3.0T扫描仪采集。首先进行常规临床脑部扫描以排除继发性头痛病因。高分辨率解剖图像采用三维T1加权快速扰相梯度回波序列获取,具体参数如下:repetition time = 7.5 ms, echo time = 2.9 ms, field of view = 256 × 240 mm², flip angle = 8°, isotropic voxel size = 1 × 1 × 1 mm³, and 208 contiguous axial slices。
所有结构磁共振数据均通过FreeSurfer 6.0版进行预处理。通过全自动流程从3D T1加权图像重建每位受试者的皮质表面,该流程包括颅骨剥离、运动校正、皮质灰白质分割、坐标转换、大脑半球及皮层下结构分离,并生成灰白质界面与软脑膜表面。
Construction of individual MIND network
我们使用五种标准化顶点形态测量特征构建个体MIND网络:皮质厚度、灰质体积、表面积、平均曲率和脑沟深度。对每位参与者,形态学特征在人类连接组计划图谱定义的360个皮质脑区内进行聚合,形成多变量结构特征分布。
采用K近邻算法计算脑区间成对相似性,通过估计多变量分布间的对称Kullback-Leibler散度实现。所得散度值随后转换为0至1标度的MIND相似性度量,数值越高代表相似性越强。在区域分析中,我们通过计算360个皮质区域各自所有连接权重(未设阈值)的平均值,得到节点级MIND值。该方法在保留相似性测量连续性的同时,提供了结构整合的局部化指标。
Statistical analysis
Case-control analysis of regional MIND network
我们使用半参数广义可加模型(GAMs)分析MIND指标的组间差异,模型构建如下:
Y = β0 + β1 · group + f1(age) + β2 · gender + β3 · TIV,
其中Y代表局部MIND值,group为分类变量(偏头痛组vs.健康对照组),𝑓1(age)是用于解释潜在非线性年龄效应的非参数平滑函数,性别和颅内总容积(TIV)作为协变量。此建模方法可灵活估计年龄相关效应,无需对函数形式施加参数假设。组间比较的效应量采用Cohen’s d量化。为控制所有360个皮层脑区的多重比较,我们应用错误发现率(FDR)校正,显著性阈值设为P < 0.05。在网络水平分析中,我们通过将相同的GAM框架应用于网络平均MIND值,检验了12个典型功能网络内的MIND差异。此互补性分析有助于识别与偏头痛相关的结构连接模式的系统水平改变。
Estimation of regional gene expressions
AHBA数据集包含六名成年捐献者死后大脑的全脑微阵列基因表达数据。我们使用abagen工具箱处理来自AHBA的全脑基因表达数据。预处理流程如下:[1]更新探针到基因的注释;[2]选择在供体间区域变异模式最一致的探针;[3]将样本分配到脑区;[4]使用稳健sigmoid函数对组织样本表达值和基因表达值进行标准化。由于仅有两名供体包含右半球数据,我们的分析仅从左半球提取基因表达值。此外,仅选择在人类蛋白质图谱(HPA)中通过表达水平识别的大脑表达基因进行后续分析,最终获得一个脑区×基因表达矩阵(180 × 7411)。
Transcriptomic analysis
为探究局部MIND改变与底层遗传组织之间的关联,首先使用病例对照比较(偏头痛vs.健康对照组)获得的Cohen’s d值量化偏头痛相关MIND异常的脑区分布图。随后,应用偏最小二乘回归分析d值分布图(左半球180个皮层脑区)与AHBA脑区基因表达谱矩阵之间的空间关联。此流程可识别与偏头痛皮层网络功能障碍模式空间共定位的基因表达特征。统计显著性通过置换检验结合空间自相关校正确定。基因权重通过Bootstrap重采样(1000次迭代)转换为z分数,并分为PLS+和PLS-组(FDR校正,P < 0.05)。
Gene set enrichment analysis
我们使用Metascape对基因本体论(GO)基因集进行功能富集分析。通过计算Kappa检验评分度量任意两个富集条目间的成对相似性。接着,对相似性矩阵进行层次聚类,并将构建的树状图修剪为独立的簇。富集通路的显著性阈值为0.05(FDR校正)。为获取细胞类型特异性标记基因,既往研究对成人大脑进行了全面的转录组和表观基因组单细胞分析。我们提取了25类脑细胞的细胞特异性基因集,主要包括星形胶质细胞、内皮细胞、小胶质细胞、神经元、少突胶质前体细胞和少突胶质细胞。为校正基因集大小,我们计算了每个基因集的富集分数,并与置换检验(n = 10000)估计的分数进行比较,从而得到标准化富集分数(NES)。
Metabolism profiles association
为探究偏头痛皮层异常的空间模式是否与大脑底层代谢格局相关,我们利用来自33名健康对照的PET衍生代谢图谱规标准数据集进行了空间相关性分析。需注意,此为间接推断,因本研究未获取偏头痛患者的PET数据。标准代谢数据包括五个特征:脑血流量、脑氧代谢率、脑葡萄糖代谢率、氧-葡萄糖指数和糖酵解指数。这些图谱来自Vaishnavi等人的研究,其使用[¹⁵O]水和[¹⁸F]氟代脱氧葡萄糖PET对静息态年轻健康成人进行测量。利用neuromaps工具箱将代谢测量数据转换至HCP_MMP图谱。基于HCP_MMP图谱,将五个代谢系统图谱分配至360个区域,得到5 × 360的代谢矩阵。我们通过多元线性回归检验了代谢特征矩阵与偏头痛相关MIND异常(Cohen’s d值)之间的关联。量化了每个代谢特征的相对重要性。
Cognitive association analysis
鉴于本研究队列缺乏专用认知测试,我们利用Neurosynth框架进行功能解码分析,以推断结构发现与认知功能的相关性。该广泛应用的方法基于超过14,000项神经影像研究的元分析数据,生成与特定认知术语稳定相关的脑区图谱。针对23个认知域,我们计算了受累脑区内重叠体素占总体的比例。此外,进行了认知术语与偏头痛相关MIND之间的相关性分析(FDR校正)。
Functional disability prediction
使用MIDAS评分量化偏头痛相关功能障碍。我们利用LIBSVM工具箱实现支持向量回归模型,根据全脑MIND模式预测个体功能障碍评分。通过留一法交叉验证评估模型性能,使用预测值与实际MIDAS评分之间的Pearson相关性评估预测准确性。
Internal validation analysis
为确保研究结果的稳健性,我们使用随机选取的70名偏头痛患者和33名HC子集进行内部重复验证。该验证队列约占总体样本量的一半。应用与发现队列完全相同的处理和分析流程,在重复数据集中计算局部MIND的组间差异。通过计算两个队列得到的Cohen’s d分布图之间的空间相似性,正式评估主要发现的重复性。此定量比较评估了效应量在皮层区域分布的 一致性,为我们观察到的偏头痛MIND改变的可推广性提供了经验证据。
Results
本研究探索队列共纳入139名偏头痛患者及65名健康对照者,验证队列包含70名偏头痛患者与33名健康对照者。人口统计学特征与临床资料汇总于Supplement Table s1。
Global and regional MIND alterations in migraine
与HC相比,偏头痛患者的全局及局部形态相似性网络均出现显著紊乱。两组间全局MIND值差异显著(d = -0.47,P = 0.002)。脑区分析显示,偏头痛患者内侧颞叶皮层的MIND值较高,而额下皮层和背外侧前额叶皮层的MIND值较低(Fig.2A-B及Supplement Table s2)。在网络层面,偏头痛患者的语言网络、扣带-岛盖网络和躯体运动网络的MIND值均降低(Fig.2C)。MIND改变的空间分布在发现队列与验证队列间呈现高度一致性(r(358) = 0.87,Pspin < 0.001)(Supplement Fig.s1)。
Fig.2 Morphometric network alterations in migraine and their cognitive and clinical correlates. A Average MIND spatial distribution patterns in HCs and migraine patients. B Case-control differences in MIND values between migraine patients and HCs expressed as a Cohen’s d map, with statistically significant regions shown below (PFDR < 0.05). C Network-level MIND comparisons showing reduced values in language, cingulo-opercular, and somatomotor networks in migraine patients (asterisks indicate PFDR< 0.05). D Cognitive terms associated with regions showing increased (top) or decreased (bottom) MIND in migraine patients. E Prediction of individual MIDAS scores based on MIND maps using support vector regression. The scatter-plot show positive relationship between actual and predicted MIDAS scores. An asterisk represents significant differences (PFDR < 0.05).
Figure s1. Replicable case-control differences of MIND. (A) Average MIND spatial distribution patterns in HCs and migraine patients. (B) Case-control difference of regional MIND between migraine patients and HCs expressed as a Cohen’s d map. Some cortical regions display statistically significant differences in MIND between migraine patients and HCs (PFDR <0.05). (C) The case-control differences of MIND in replication cohort were positively correlated with the differences in discovery cohort.
Transcriptional and cellular signatures of MIND alterations
为了识别这些结构改变背后的转录特征,我们采用PLS回归分析方法,将左半球180个皮质脑区的MIND异常空间模式与这些相同区域的7411个大脑表达基因的表达谱进行关联。PLS回归分析显示,局部MIND差异与基因表达模式之间存在显著关联(PLS1:解释方差14.9%,Pspin < 0.01;PLS2:解释方差16.8%,Pspin < 0.01)。其中PLS2的表达模式与观测到的MIND差异呈现最强正相关性(r(178)=0.41,Pspin < 0.001,Fig.3A-C)。该分析筛选出124个显著关联基因,包括74个过表达基因(PLS2+)和50个低表达基因(PLS2-)(Fig.3D及Supplementary Table s3)。对这些基因进行功能富集分析发现,它们在突触信号传导和神经递质运输调控相关生物过程中显著富集(Fig.3E)。细胞类型特异性分析进一步表明,这些转录特征主要与内皮细胞、少突胶质细胞、兴奋性神经元和抑制性神经元相关(Fig.3F)。
Fig. 3 Transcriptomic signatures of case-control differences in MIND. A Average migraine-associated MIND alterations show on the cortical surface. B The weighted gene expression profile of PLS2 across cortical regions, represented on the surface of the brain. C Scatter plot displaying the correlationbetween case-control MIND differences (d values) and regional PLS2 scores. D Ranked list of genes with significant PLS2 loadings, showing the top and bottom-ranked genes related to migraine (P < 0.05 after FDR correction). E Gene Ontology (GO) terms enriched in migraine-related MIND alterations, with circle size representing gene count and color intensity reflecting statistical significance (P < 0.05 after FDR correction). F Brain cell-type specificity of gene expression, with normalized enrichment scores (NES) showing the association of specific genes with different cell types. The bolded and underlined terms and cell types are significant after FDR correction (P < 0.05).
Metabolic signatures associated with MIND changes
本研究通过将偏头痛相关MIND改变的空间模式与皮层灰质标准代谢图谱进行关联分析,发现两者存在显著整体关联。多元线性回归模型可解释偏头痛相关MIND改变5.61%的变异度(F(5, 354) = 4.21,Pspin = 0.001)。模型预测的MIND差异值与实际观测值呈显著正相关(r(358) = 0.24,Pspin = 0.002,Fig.4A-B)。其中氧-葡萄糖指数(OGI)贡献度最高(相对贡献率 = 41.0%,PFDR < 0.05),其次为糖酵解指数(GI)(相对贡献率 = 28.0%,PFDR < 0.05),其余代谢特征贡献均未达显著性水平(所有all PFDR > 0.05)。Fig.4A-C展示了代谢特征与MIND改变的空间关联分布,结果显示颞叶及前额叶区域为主要关联脑区。
Fig.4 Cortical MIND network alterations and metabolic associations in migraine. A Fitted MIND difference between migraine patients and HCs with metabolic profiles expressed as a Cohen’s d map displayed on the cortical surface. B Scatter plot showing the correlation between fitted and observed d values with metabolic profiles. C Radar chart depicting the relative contributions of metabolic markers to MIND differences. An asterisk indicates the significant difference after FDR correction (P < 0.05).
Cognitive correlates of MIND abnormalities
偏头痛患者MIND值升高的脑区与“记忆”、“意象”及“情绪”认知域存在关联(Fig.2D),而MIND值降低的脑区则与“动作”、“语言”及“疼痛”认知域相关(Fig.2D)。这些空间关联分析表明,偏头痛患者MIND指标的改变与认知功能障碍存在联系。具体认知术语与偏头痛相关MIND的关联关系详见Table S4。
Prediction of disability by MIND patterns
使用局部MIND值的支持向量回归模型能够显著预测偏头痛患者的MIDAS评分(r(137) = 0.23,Pperm = 0.015)。该模型纳入了存在显著病例-对照组差异脑区的MIND值作为预测特征。通过留一交叉验证对预测准确性进行评估,实际与预测MIDAS评分的对应关系如Fig.2E所示。
Discussion
Structural network abnormalities and symptom heterogeneity
本研究揭示了偏头痛患者大脑皮层重塑的独特模式,为该疾病复杂的临床症状提供了全新见解。我们观察到内侧颞叶——记忆巩固与情感处理的关键脑区——的MIND值显著增高。虽然本研究未直接进行认知或精神评估,但海马与内嗅皮层等内侧颞叶区MIND值的升高现象,为偏头痛患者常见认知与情绪症状提供了神经生物学解释。这些参与记忆巩固与情感处理的关键脑区表现出的超同步结构改变,可能构成发作期相关现象的病理基础,该发现为未来研究此类结构改变与情绪障碍、记忆相关症状等发作期现象的关联提供了可验证的假说。与此形成互补的是,我们在前额叶区域(包括背外侧前额叶皮层和额下回)观察到MIND值降低,该现象与执行控制功能紊乱及下行疼痛调节受损的已知神经基础相吻合。研究提示了一种潜在机制:这些执行与调控网络的结构完整性受损可能削弱患者对传入伤害性信号的“自上而下”调控能力,这既可部分解释难治性偏头痛的形成机制,也能说明患者常见的注意力与认知缺陷。前额叶区域的受累与先前关于扣盖网络在疼痛处理中作用的研究一致,并进一步佐证了偏头痛作为多网络功能障碍疾病的认识。我们的研究继承并拓展了既往偏头痛神经影像学发现。静息态功能磁共振研究持续显示感觉运动网络(SMN)亚区的功能改变——该网络是整合感觉、运动与认知处理的关键枢纽。初级感觉皮层、运动皮层、前运动皮层及辅助运动区等SMN成分的功能改变,已被证实与情感、认知及疼痛处理功能紊乱相关。本研究的结构发现与这些功能观察结果相互印证,共同强化了偏头痛(特别是疼痛相关及感觉运动网络)存在广泛网络功能障碍的证据。语言网络在本研究中显现为受影响的有趣系统,这与近期偏头痛语言相关异常的报告相符。偏头痛发作期常见的失语现象在我们的研究中找到了结构对应——语言相关皮层区域的MIND值改变。值得注意的是,受影响的MIND网络区域关联的认知功能(包括行动、语言、疼痛及情绪处理)与偏头痛患者报告的认知异常谱系存在特定一致性。这些发现在内部验证队列中成功复现,显示发现样本与验证样本的MIND改变具有高度空间相似性,进一步强化了结果的可靠性。特定MIND改变与具体认知域的关联性,解释了为何蒙特利尔认知评估(MoCA)等标准化筛查工具在偏头痛患者中常出现结果不一致的现象。由于认知缺陷并非整体性而是选择性出现——由离散的网络功能障碍驱动——通用筛查工具可能无法捕捉细微病理。这与临床研究显示MoCA客观评分变异度大、通常仅揭示执行功能与记忆等领域选择性缺陷,且受年龄、发作时长等因素显著影响的结果一致。更重要的是,我们的研究突破了将偏头痛传统视为单纯疼痛疾病的局限,将其重新定义为影响多认知域的广泛性网络功能障碍疾病。这一新认知促使我们在临床实践中需常规评估头痛特征之外的认知症状,并考虑针对特定认知改变模式患者实施认知康复等靶向干预。
在神经影像学层面,我们的发现得到近期大规模荟萃分析的支持并与之相互整合。我们观察到的额叶脑回与内侧颞叶异常,与激活似然估计(ALE)荟萃分析识别的脑区高度一致。偏头痛患者额叶脑回与脑岛的聚类结构改变(灰质体积减少)及丘脑功能激活改变也见诸报道。这些脑区在多项研究中的集中关注强化了我们的结论:这些脑区是偏头痛病理的核心节点,而非随机发现。此外,功能荟萃分析反复发现的默认模式网络(DMN)连接异常,在概念层面与我们的MIND改变相契合,共同反映了发作间期内在网络的功能障碍。
Evidence for abnormalities in energy metabolism and mitochondrial function in migraine
神经活动对能量的高需求早已被广泛证实,大量研究表明大脑能量消耗与葡萄糖-氧代谢之间存在紧密耦联关系。我们的研究提供了有力证据,表明代谢失调在偏头痛患者脑结构网络改变中起着根本性作用。偏头痛相关脑结构网络异常与特定代谢模式(尤其是能解释41%变异度的氧-葡萄糖指数)之间的强关联,为偏头痛病理生理机制提供了新见解。这种代谢-结构关联有助于解释偏头痛临床观察中多个长期存在的现象:偏头痛发作期血糖急性升高这一公认现象,可能是大脑应对发作期能量需求增加的代偿机制。葡萄糖作为脑代谢主要供能物质,不仅能生成ATP,还可维持氧化还原平衡并支持关键神经递质与结构成分的合成。在偏头痛发作期间,能量需求增加与潜在葡萄糖缺乏的共同作用可能导致氧化应激,形成恶性循环,进而促成我们在脑结构网络中观察到的结构改变。我们的发现与既往研究一致,即偏头痛症状严重程度与线粒体功能障碍程度及糖酵解活性降低相关。动物研究表明糖酵解抑制剂可缓解偏头痛样症状,这进一步支持了代谢途径调控可能成为可行治疗策略的观点。氧-葡萄糖指数的突出影响提示,针对有氧糖酵解的干预措施(如生酮饮食或线粒体支持疗法)可能有助于阻止网络退化。虽然PET成像仍是金标准,但血清乳酸/丙酮酸比值等更易获取的代谢标志物,可帮助识别最能从代谢疗法中获益的患者。通过证明代谢异常与结构改变相关,我们的研究结果表明早期代谢干预可能对偏头痛治疗有益。这种代谢视角为开发偏头痛神经保护策略提供了新机遇。
Transcriptional signatures: rethinking therapeutic targets
我们的转录组分析揭示了偏头痛相关皮层变化的关键分子特征。MIND改变脑区内突触基因与内皮细胞、少突胶质细胞标志物的共同富集现象超越了单纯的神经元过度兴奋模型,这支持了偏头痛的多细胞疾病机制。这些发现与现有遗传学证据高度契合:全基因组关联研究已确认突触功能与髓鞘完整性相关通路在偏头痛中的作用。内皮细胞特征验证了临床观察中血脑屏障功能障碍的关联,或可解释合并血管性共病患者对CGRP抑制剂反应的差异性。鉴于这类生物制剂主要靶向神经元通路,其在血管共病或白质病变患者中疗效减弱的现象,可能反映了内皮功能障碍的重要作用。我们的转录特征成功捕捉到偏头痛发作期导致内皮损伤与血脑屏障破坏的炎症过程,提示血管靶向治疗可能惠及特定患者亚群。研究中对少突胶质细胞的突出发现具有同等重要性,这与遗传学和组织学研究的多重证据相互印证。全基因组关联研究持续提示髓鞘完整性基因(如MBP基因)与偏头痛易感性相关,近期神经病理学工作更在偏头痛患者三叉神经通路中发现髓鞘板层结构破坏。兴奋性与抑制性神经元特征的同步识别支持皮层兴奋性失衡假说,而炎症标志物则进一步佐证内皮功能障碍机制。这些结果为开发针对偏头痛病理生理中突触、血管及胶质成分的靶向疗法提供了分子框架。我们的遗传学发现与偏头痛遗传学研究进展高度一致。大规模全基因组关联研究(GWAS)荟萃分析极大拓展了已知偏头痛风险位点,其通路富集分析持续揭示共同的生物学机制,包括神经元兴奋性、离子通道、血管因子及突触功能。与我们MIND模式相关的基因组在谷氨酸能信号传递、突触小泡运输等神经生物学通路上呈现显著富集,这为遗传风险提供了功能性神经解剖学关联。通过证明遗传背景相关的功能改变与临床功能障碍(MIDAS评分衡量)存在相关性,我们的研究在群体遗传数据与个体脑网络功能障碍之间构建了桥梁,为临床转化提供了重要洞见。
Significance and limitations
我们的研究结果具有多方面的转化意义,可促进对偏头痛的临床理解和治疗方法的发展。首先,MIND异常定位于内侧颞叶(负责记忆和情绪处理)和前额叶皮层(负责执行控制),这为解释偏头痛发作时常见却令人困惑的临床症状(如认知模糊和情绪不稳定)提供了直接的神经解剖学依据。其次,我们观察到的MIND模式与功能障碍(MIDAS评分)之间的密切关联表明,MIND可作为一种亟需的客观生物标志物。这对于量化超越患者主观报告之外的疾病负担具有重要价值,尤其是在自我报告困难的人群中,如儿科或老年患者。第三,OGI作为关键病理驱动因素的发现,为我们提供了一个引人注目的机制靶点。这一发现为使用生酮饮食或其他旨在稳定大脑能量稳态的非药物干预措施作为潜在治疗途径提供了理论依据。总之,通过将偏头痛的病理生理学基础定位于多尺度网络(基因、代谢、功能)的紊乱,本研究主张临床关注点应实现关键转变:从纯粹的症状管理转向基于病理学的分层模型。本质上,这种方法可能为真正针对生物学的靶向治疗铺平道路。
在解读我们的发现时,需考虑几个局限性。首先,基因表达分析仅限于左半球,并排除了皮层下脑区,这可能导致偏侧化偏差,且未能全面表征皮层下脑区的作用。未来的研究应探讨右半球和皮层下脑区的基因-结构关系,以提供更全面的理解。其次,横断面设计和相对有限的样本量限制了我们建立因果关系的能力。更大规模的纵向研究将有助于阐明MIND异常与分子变化之间的时间动态关系。第三,临床特征主要集中于偏头痛核心症状;纳入更详细的表型数据,包括头痛特征和相关症状,将使结构-临床相关性分析更为细致。第四,我们的代谢和认知分析分别依赖于来自规范性PET数据集和Neurosynth元分析数据库的间接推断,而非直接的患者测量结果。虽然这种"空间关联"方法在生成假设方面非常有力,但我们的发现仅能表明偏头痛影响脑区具有特定的规范性特征;它们并不能确认我们队列的实际代谢或认知状态。未来的研究需要纳入这些全面的患者测量数据,以验证这些发现并在个体患者中建立直接的结构-功能-症状关系。第五,尽管所有患者均在无发作期接受扫描,但MRI扫描相对于偏头痛周期的具体时间并未得到系统控制。其潜在影响可能微乎其微,因为本研究重点关注皮层形态计量学——这是一种结构测量,与更易变的功能性测量相比,其在短时间间隔内相对稳定。因此,虽然缺乏精确的时间数据是一个局限性,但它对我们的结构研究结果的潜在影响可能很小。另一个需要考虑的因素是偏头痛人群内部的异质性。我们的队列包括了有先兆和无先兆的患者。为了保持统计效力,我们将这些组别合并,这一决定得到了敏感性分析的支持,该分析显示我们的主要研究结果高度稳健(r(358) = 0.99 (Pspin < 0.001),Supplemen Fig.S2)。
Figure s2. High consistency of MIND alterations. A scatter plot showing the strong spatial correlation between the case-control differences (Cohen's d values) from the primary analysis (x-axis) and the MwoA-only sensitivity analysis (y-axis). Each point represents one of the 360 cortical regions.
Conclusion
基于常规MRI获取的MIND网络改变对偏头痛相关的临床功能障碍与认知损害具有敏感性,其背后机制涉及有氧糖酵解失调以及突触信号传导、神经递质转运和特定神经血管细胞类型相关的转录特征富集。本研究将偏头痛的宏观皮层形态测量与微观分子机制相连接,为阐明其病理生理学提供了新视角,并识别出潜在的治疗靶点。
点击“阅读原文”可查看更多相关文章。原文及补充材料可联系作者获取。
对神经影像、功能磁共振成像及组学,机器学习相关内容感兴趣者可以扫码入群交流学习:
欢迎微信扫码关注神经影像石大夫
点亮“在看” ,让更多人看到