帕金森病多组学新型生物标志物与治疗靶点鉴定：从转录组到药物相互作用

《Frontiers in Aging Neuroscience》：Multi-omics identification of novel biomarkers and therapeutic targets for Parkinson's disease: from transcriptome to drug interaction

　　

帕金森病（Parkinson's disease, PD）是一种进行性神经退行性疾病，治疗选择有限。针刺可能作为一种补充干预手段，但其分子机制尚不清楚。本研究整合批量与单细胞转录组数据集，以鉴定PD相关诊断候选基因，并探索其与针刺相关时间表达模式的关系。研究人员从GSE68719中鉴定PD相关差异表达基因。对GSE178470的时间序列分析刻画了一名PD患者在基线和针刺5次及8次后纵向血液转录组中表达轨迹的上升或下降特征。将PD相关失调与针刺后时间变化之间方向相反的基因取交集，随后使用LASSO、随机森林（Random Forest）、Boruta和XGBoost进行特征选择。在GSE68719中评估诊断性能。进一步采用单细胞转录组分析、分子对接和MPP⁺处理的SH-SY5Y细胞进行探索性表征与验证。结果共鉴定出11个候选基因，其中4个被机器学习方法优先选择。该四基因组合在GSE68719中区分PD患者与健康对照的表现：表观AUC为0.873，敏感性为0.690，特异性为0.864，准确度为0.795，精确度为0.769，F1得分为0.727，Brier得分为0.134。这些指标反映的是PD与对照的分类能力，而非对针刺应答的预测。单细胞分析绘制了人类中脑群体中细胞类型特异性表达图谱。分子对接提供了探索性的预测结合姿态。BAG3和HSPB1在MPP⁺处理的SH-SY5Y细胞中表达上调。讨论认为，这些发现为优先筛选具有探索性针刺相关时间模式的PD候选基因提供了一个整合框架。仍需在更大规模的纵向队列和功能研究中进一步验证。

帕金森病（Parkinson's disease, PD）是一种常见的进行性神经退行性疾病，以黑质致密部多巴胺能神经元选择性丢失和路易小体（Lewy bodies）形成为主要病理特征。随着全球人口老龄化，PD的发病率和 societal burden 持续上升。尽管左旋多巴及其衍生物等多巴胺替代疗法仍是症状治疗的主要手段，但常伴随运动波动和异动症等长期并发症，且无法阻止潜在的神经退行过程。因此，迫切需要不仅有效、安全，而且能够改变疾病进程的治疗策略。针刺作为传统中医（traditional Chinese medicine, TCM）的核心疗法，在包括PD在内的神经系统疾病中受到日益关注。临床和实验研究表明，针刺可能改善PD患者的运动症状、情绪和睡眠障碍，并可能延缓疾病进展。机制研究提示针刺可能通过调节神经递质水平、减轻神经炎症和增强线粒体功能发挥神经保护作用。然而，多数研究聚焦于宏观或局部效应，缺乏全面的分子证据——尤其是转录组水平上针刺如何改变与PD发病机制相关的基因表达谱。与此同时，高通量测序和计算生物学的进展使得复杂疾病中基因表达变化的系统研究成为可能。批量RNA测序和单细胞RNA测序（single-cell RNA sequencing, scRNA-seq）等转录组分析已被广泛用于揭示PD的关键分子通路和治疗靶点。尽管如此，利用这些方法研究针刺治疗的分子应答仍很少见。此外，机器学习方法虽越来越多地用于PD诊断或疾病相关基因组合的筛选，但其在针刺相关转录组假说中的应用仍然有限。这构成了一个重要空白，尤其是在建立治疗应答的稳健生物标志物和理解针刺如何调节PD免疫微环境方面。

针对上述空白，研究人员构建了一个整合性转录组框架，以鉴定PD相关候选基因并探索其与针刺相关时间表达模式的关系。研究人员首先利用Mfuzz算法对一名代表性PD患者在针刺治疗期间（基线、5次和8次针刺后）的纵向血液样本进行了时间过程转录组分析，以描述动态基因表达轨迹。通过整合一个独立的PD病例对照数据集，进一步将候选基因缩小至那些在PD中失调且针刺后呈方向相反时间模式的基因。为优先筛选PD相关诊断候选基因，研究人员采用了四种互补的机器学习方法——随机森林（Random Forest）、LASSO回归、XGBoost和Boruta特征选择。这些基因在区分PD患者与健康对照方面表现出区分能力，提示其作为PD相关诊断候选基因的潜在价值。为揭示这些核心基因的功能相关性，研究人员开展了一系列计算机模拟和体外分析。免疫细胞浸润分析将这些基因与PD病理生理相关成分联系起来；单细胞RNA测序将其表达定位至小胶质细胞、内皮细胞和室管膜细胞。分子对接提供了候选蛋白与已上市多巴胺能药物之间可能的蛋白-配体结合姿态的探索性计算机模拟预测。为验证这些计算衍生靶点的内在病理参与，研究人员使用PD细胞模型进行了体外实验。该生物学验证显示，在MPP⁺处理的细胞模型中BAG3和HSPB1发生应激诱导性上调，支持其应激反应性表达而非证明其在针刺诱导神经保护中的因果作用。总之，该整合计算与实验框架优先筛选出具有探索性针刺相关时间模式的PD候选基因，为未来评估基线或纵向基因表达变化是否与PD针刺临床结局相关提供了可检验的假说。

研究人员使用了多种公开转录组数据集。GSE178470包含一名代表性PD患者接受针刺治疗的纵向血液转录组谱（基线、5次和8次治疗后）。GSE68719包含PD患者与健康对照的差异基因表达谱。独立验证队列来自GSE239454。单细胞RNA测序数据来自人类中脑GSE157783。体外实验采用人神经母细胞瘤SH-SY5Y细胞系（来源于中国典型培养物保藏中心），用0.5 mM 1-甲基-4-苯基吡啶（MPP⁺）处理24小时建立PD细胞模型。关键方法包括：使用limma R包进行差异表达分析；使用Mfuzz进行针刺干预期间的时间表达模式聚类；使用四种机器学习算法（LASSO、随机森林、Boruta、XGBoost）进行特征选择和诊断模型构建；使用单样本基因集富集分析（ssGSEA）评估通路相关性；使用Harmony进行批次校正和UMAP降维分析单细胞数据；使用AutoDock Vina进行分子对接；通过RT-qPCR和Western blotting验证BAG3和HSPB1的表达变化。

结果部分包括以下内容：探索性时间表达模式：对GSE178470中一名针刺应答性PD患者的时间过程转录组应用Mfuzz聚类，鉴定出6个时间表达模式，其中4个簇呈持续上升或下降轨迹；7,878个基因表达逐渐增加，3,631个基因持续下降。功能富集分析显示这些动态表达基因富集于神经系统稳态相关通路，KEGG分析突出“多巴胺能突触”和“神经退行性疾病—多种疾病”通路，GO术语在突触功能和可塑性方面显著富集。由于时间过程RNA-seq缺乏生物学重复，这些结果仅作为假说产生。

鉴定PD相关诊断候选基因：在GSE68719中鉴定出49个上调基因和17个下调基因；与针刺相关时间基因集取交集获得11个候选基因。其中8个基因在PD中低表达而针刺后升高，3个基因（HSPB1、KCNE4、BAG3）在PD中高表达而针刺后降低。四种机器学习算法特征筛选后锁定四基因组合（GPCPD1、HSPB1、KCNE4、BAG3）。基于该组合的逻辑回归模型在GSE68719中表观AUC为0.873（95% CI 0.789–0.958），敏感性0.690，特异性0.864，准确度0.795，精确度0.769，F1得分0.727，Brier得分0.134。5折交叉验证逻辑模型AUC为0.776，XGBoost模型为0.900。这些指标反映PD与健康对照的诊断区分。

四基因模型通路关联与单细胞表达模式：在GSE68719中，四个核心基因在PD患者与健康对照间显著差异表达。ssGSEA显示这些基因与氧化应激（相关系数0.58–0.65）、自噬（0.40–0.73）和突触信号（0.57–0.65）呈中度至强相关，而与多巴胺信号相关较弱（0.01–0.16）。单细胞分析显示PD与对照间室管膜细胞和小胶质细胞比例存在描述性差异；GPCPD1在健康对照的小胶质细胞、室管膜细胞和内皮细胞中表达水平较高，提示其可能参与细胞稳态和神经胶质支持。

探索性分子对接：对GPCPD1、KCNE4、BAG3、HSPB1与左旋多巴、美左旋多巴和乙左旋多巴进行分子对接。GPCPD1与三种药物结合能最低（?7.1、?7.1、?6.9 kcal/mol），BAG3亲和力较弱（?4.4至?4.8 kcal/mol），KCNE4与美左旋多巴和乙左旋多巴中度结合（?4.6、?4.7 kcal/mol），HSPB1仅与乙左旋多巴表现出较强结合（?6.1 kcal/mol）。这些结果为优先验证GPCPD1提供了结构假说。

BAG3和HSPB1在MPP⁺诱导细胞应激模型中上调：在MPP⁺处理的SH-SY5Y细胞中，BAG3和HSPB1的mRNA水平（P=0.0004和P=0.0198）和蛋白水平（P=0.0458和P=0.0286）均显著高于对照组，表明这两种基因对MPP⁺诱导的细胞应激有应答。

讨论部分总结：该研究整合批量转录组、单细胞转录组、机器学习、探索性对接和靶向体外表达分析，优先筛选出PD相关诊断候选基因，并生成关于针刺相关时间表达变化的假说。多层验证增强了结果的稳健性和临床相关性。局限性包括：GSE178470时间过程来自单一代表性患者，Mfuzz结果应视为描述性和假说生成性；四基因模型用于区分PD与对照而非预测针刺应答；特征筛选未完全嵌套于交叉验证中；单细胞、分子对接和MPP⁺实验提供支持性但探索性的证据。结论：该研究提供了一个系统性整合转录组分析框架，优先筛选出具有探索性针刺相关时间模式的PD诊断候选基因。四基因组合可作为未来诊断和机制研究的有用候选集；BAG3和HSPB1的应激反应性表达为其在PD相关实验模型中的进一步研究提供了生物学支持。未来需纳入重复针刺纵向队列、临床结局指标和功能验证，以确定这些候选基因是否在机制上参与针刺对PD的相关效应。