基于整合机器学习与实验验证的SIK1-NR4A2轴：癫痫的核心驱动机制研究

《Frontiers in Cell and Developmental Biology》：The SIK1-NR4A2 axis: a core driver of epilepsy identified by integrated machine learning and experimental validation

　　

摘要专业翻译 背景： 癫痫（Epilepsy, EP）仍是重大临床挑战，提示除经典兴奋-抑制失衡之外，尚有其他机制参与癫痫发生。肝激酶B1（Liver kinase B1, LKB1）信号通路协调能量稳态与细胞存活，但其在癫痫中的关键下游效应分子 largely 未知。 方法： 研究人员整合了基于公共基因表达汇编（Gene Expression Omnibus, GEO）数据集（GSE143272）的LKB1通路基因机器学习筛选，并结合体内（in vivo）与体外（in vitro）实验验证。采用戊四氮（pentylenetetrazol, PTZ）诱导的癫痫持续状态（status epilepticus, SE）大鼠模型、腺相关病毒（adeno-associated virus, AAV）介导的基因敲低/过表达、原代海马神经元慢病毒转染、脑电图（Electroencephalography, EEG）记录、全细胞膜片钳记录及蛋白质印迹（Western blot）分析，以探究候选基因的功能作用及分子机制。 结果： 机器学习鉴定出盐诱导激酶1（salt-inducible kinase 1, SIK1）为LKB1通路中排名首位的癫痫风险基因。SE大鼠海马中SIK1表达显著上调。AAV介导的SIK1敲低延长了癫痫发作潜伏期并缩短发作持续时间，而SIK1过表达则加重发作严重程度。机制上，通过蛋白质-蛋白质相互作用（protein-protein interaction, PPI）网络分析与中介分析，研究人员鉴定出NR4A2为SIK1的下游效应分子。SIK1负向调控NR4A2表达：SIK1敲低增加NR4A2，而SIK1过表达降低NR4A2。拯救实验表明，SIK1与NR4A2共过表达可逆转SIK1单独所致的促癫痫效应。全细胞膜片钳记录显示，在无镁处理神经元中，SIK1敲低或NR4A2过表达可将动作电位幅度恢复至接近正常水平。 结论： SIK1负向调控NR4A2，恢复NR4A2表达可拮抗SIK1诱导的神经元过度兴奋。这些发现提示，靶向SIK1-NR4A2轴可能成为癫痫的潜在治疗策略。

论文解读文章

一、研究背景与问题提出

癫痫是一种以反复癫痫发作为特征的慢性脑部疾病。传统观点认为，癫痫发作源于脑内神经元兴奋与抑制的失衡（E-I imbalance）。然而，尽管大多数临床抗癫痫药物（ASMs）旨在纠正这种失衡，仍有相当比例的患者应答不佳，提示兴奋-抑制失衡理论可能无法完全解释癫痫发生机制。既往研究已将关注点转向其他病理过程，例如脑能量代谢失衡及随之而来的细胞死亡信号通路激活，在癫痫发作中发挥关键作用。在协调能量稳态与细胞存活的众多信号网络中，肝激酶B1（LKB1）信号通路是重要枢纽。LKB1通过激活AMP活化蛋白激酶（AMPK）感知细胞能量状态，从而调节代谢过程并防止能量耗竭所致的细胞损伤。在癫痫背景下，LKB1信号失调可能加剧能量代谢失衡，例如在锌兴奋毒性模型中，LKB1-AMPK-Bim级联被激活，诱导caspase-3依赖性神经元凋亡。此外，LKB1功能障碍可导致mTORC1信号过度激活，促进神经元异位及癫痫样活动。然而，在癫痫这一特定背景下，LKB1信号通路中促成癫痫的关键中介分子 largely 未知。为此，研究人员采用机器学习方法，利用公开的患者转录组数据系统筛选LKB1通路相关基因中可能促成癫痫风险的基因，并在动物模型中验证排名最前候选基因的表达与功能作用，进而探索其潜在下游分子机制及在调节神经元兴奋性中的功能意义。

二、研究设计与主要技术方法

本研究整合了公开转录组数据挖掘与体内外实验验证。研究人员从基因表达汇编（GEO）下载GSE143272数据集（含51例健康个体与91例癫痫患者外周血样本）作为发现队列，并将GSE33000（158例正常脑组织样本）与GSE63808（129例癫痫脑组织样本）合并作为独立外部验证队列。LKB1相关基因取自基因集富集分析（GSEA）数据库的PID_LKB1_PATHWAY基因集。基于"caret"R包构建随机森林（RF）、支持向量机（SVM）、广义线性模型（GLM）及极端梯度提升（XGB）四种机器学习模型，采用重复5折交叉验证优化超参数。利用"pROC"R包绘制ROC曲线，以"rms"R包构建列线图模型，并通过SHAP方法进行模型可解释性分析。中介分析采用"mediation"R包，以SIK1表达为自变量、候选基因为中介变量、癫痫状态为因变量，通过1000次自举重抽样估计平均因果中介效应（ACME）。实验验证方面，研究人员建立戊四氮（PTZ）诱导的癫痫持续状态（SE）大鼠模型，采用腺相关病毒（AAV）介导的海马内SIK1敲低，进行脑电图（EEG）记录；体外实验中，分离新生大鼠原代海马神经元，采用慢病毒（LV）转染实现SIK1敲低/过表达及NR4A2过表达，以无镁人工脑脊液诱导癫痫样活动，通过全细胞膜片钳记录动作电位，并结合免疫荧光染色与蛋白质印迹（Western blot）检测蛋白表达。

三、研究结果

基于机器学习的癫痫风险关键基因鉴定： 基于LKB1相关基因构建的四种机器学习模型中，SVM模型被鉴定为最优模型（AUC = 0.884）。使用合并的GSE33000和GSE63808数据集进行外部验证，RF模型AUC达0.864，优于SVM（AUC = 0.780），提示RF可能具有更好的泛化能力。基于RF模型，排名前五的关键癫痫风险基因为PRKAB1、CTSD、SMARCD3、MLST8和SIK1，并以此构建列线图模型。研究发现，PRKAB1和CTSD表达降低时癫痫风险增加，而SMARCD3、MLST8和SIK1表达升高时风险增加。

SHAP分析揭示关键变量的预测贡献： SHAP分析显示SIK1在五个关键基因中具有最高SHAP值。SIK1和MLST8的SHAP值随基因表达水平升高而增加，CTSD的SHAP值随表达升高而降低，而SMARCD3和PRKAB1的SHAP值与表达的关系并非纯线性。鉴于SIK1具有最高SHAP值，研究人员聚焦于SIK1进行进一步机制探索。

中介分析鉴定NR4A2为SIK1的潜在中介因子： 通过GeneMANIA构建基因相互作用网络及STRING构建PPI网络，SIK1被预测与多个候选基因/蛋白相互作用。中介分析显示，在所有候选基因中，仅NR4A2显示出显著的中介效应（ACME p < 0.001），提示NR4A2可能作为SIK1在癫痫中的关键中介因子。GEO数据集分析确认了SIK1和NR4A2在对照与癫痫样本中的表达。

SIK1在癫痫发作中发挥关键作用： 免疫荧光分析显示SIK1和NR4A2均在大鼠海马中表达。EEG记录显示对照组大鼠无明显癫痫样放电，而SE组大鼠出现显著棘波、棘慢复合波等癫痫样放电，且发作持续时间超过30分钟。与对照组相比，SE组SIK1表达显著增加，NR4A2和CTSD表达显著降低。AAV介导的海马SIK1敲低后，Western blot证实SIK1表达降低而NR4A2表达升高。PTZ攻击后，shRNA.SIK1组与NC组相比，癫痫发作潜伏期显著延长、发作持续时间显著缩短。

SIK1通过抑制NR4A2加重癫痫发作： 原代海马神经元慢病毒转染实验证实，SIK1敲低显著增加NR4A2表达，而SIK1过表达显著降低NR4A2表达，确认SIK1与NR4A2之间的负向调控关系。拯救实验中，OE.NR4A2单独转染增加NR4A2表达而不影响SIK1水平，OE.SIK1与OE.NR4A2共转染成功过表达两个基因。全细胞膜片钳记录显示，无镁处理显著增加NC神经元动作电位的频率和幅度。SIK1敲低显著降低放电幅度至接近NC对照水平，但放电频率仍升高；SIK1过表达进一步增加放电频率和幅度，导致最严重的神经元过度兴奋；NR4A2单独过表达将放电幅度恢复至接近正常水平，与SIK1敲低效果相似；SIK1与NR4A2共过表达逆转了SIK1单独过表达的加重效应，将放电幅度正常化至与OE.NR4A2组相似水平，但放电频率仍增加。

四、讨论与结论

讨论部分指出，机器学习筛选鉴定SIK1为与癫痫风险相关的排名靠前候选基因，SIK1作为AMPK相关激酶，响应细胞能量状态和cAMP依赖性信号变化，参与神经元兴奋性和应激反应。SHAP分析显示SIK1表达水平与其SHAP值正相关。PTZ诱导的SE大鼠模型中SIK1表达显著上调，与机器学习预测一致。AAV介导的SIK1敲低延长发作潜伏期并缩短发作持续时间，过表达则加重发作严重程度，体外实验也一致地显示SIK1敲低降低动作电位幅度、过表达增加频率和幅度，确立SIK1为癫痫发作的关键功能驱动因子。通过基因相互作用网络和PPI分析鉴定NR4A2为SIK1的潜在相互作用伙伴，中介分析证实NR4A2显著中介SIK1对癫痫风险的影响。Western blot证实SIK1负向调控NR4A2表达。膜片钳记录显示SIK1敲低或NR4A2过表达单独即可将动作电位幅度恢复至接近正常水平，而SIK1与NR4A2共过表达可逆转SIK1单独作用。值得注意的是，幅度可被有效恢复而频率仍升高，提示幅度和频率可能受不同分子机制调控。讨论还指出，GSE143272外周血数据集中SIK1表达未显著升高，可能反映外周血转录特征与脑组织分子改变的差异，而PTZ大鼠海马中SIK1蛋白显著上调，提示SIK1上调在脑组织中可能更为显著。NR4A2（又称Nurr1）是参与神经元发育和功能的核受体转录因子，近期临床遗传学研究在神经发育障碍和癫痫患者中鉴定出NR4A2的de novo错义和功能缺失变异，系统综述提出NR4A2为发育性和癫痫性脑病的一线靶基因，另一研究证明NR4A2上调可保护创伤后癫痫，这些临床和实验发现与本研究观察到的癫痫脑中NR4A2表达降低及恢复NR4A2表达可拮抗SIK1诱导的神经元过度兴奋高度一致。研究局限性包括：仅使用PTZ诱导的急性癫痫模型，未检查KA或毛果芸香碱等其他化学惊厥模型；SIK1负向调控NR4A2的精确分子机制尚未阐明；未检查PTZ模型中SIK1 mRNA水平或内源性基因突变；外部验证队列因脑区、组织采集方法和平台差异应视为初步验证；动作电位幅度和频率的差异调控机制尚不清楚。

研究结论翻译： 总之，通过整合机器学习筛选与体内和体外实验验证，本研究将SIK1-NR4A2轴鉴定为癫痫中一个此前未被认识的调控因子。研究人员证明SIK1在癫痫中显著上调，并负向调控NR4A2——一种近期临床遗传学研究提示与神经发育障碍和癫痫相关的核受体转录因子。机制上，SIK1敲低或NR4A2过表达可保护免受癫痫样神经元放电，而恢复NR4A2表达可逆转SIK1的促癫痫效应。这些发现提示SIK1-NR4A2轴可能代表癫痫中的关键通路，并为靶向该轴作为新型治疗策略提供了潜在依据，但仍需在慢性模型和机制研究中进一步验证。