帕金森病候选生物标志物识别的探索性多组学(multi?omics)分析

《Neuroscience》：Exploratory multi?omics analysis for candidate biomarker identification in Parkinson’s disease

　　

帕金森病（Parkinson's disease, PD）是一种缺乏公认临床生物标志物的进行性神经退行性疾病。研究人员对来自河南大学淮河医院PD患者及对照者小规模队列的外周血转录组和蛋白质组数据进行探索性多组学分析，经转录组测序获得124个差异表达基因（differentially expressed genes, DEGs）和28个差异表达蛋白，并进行蛋白?蛋白互作（protein?protein interaction, PPI）分析及基因集变异分析（gene set variation analysis, GSVA）。将结果与来源于Gene Expression Omnibus（GEO）数据集3的2,964个DEGs取交集，得到28个基因组成的Panel。研究人员通过随机森林（random forest, RF）、支持向量机（support vector machine, SVM）及主成分分析（principal component analysis, PCA）对潜在标志物进行优先级排序，并在外部GEO数据、临床血清样本及1?甲基?4?苯基?1,2,3,6?四氢吡啶（1?methyl?4?phenyl?1,2,3,6?tetrahydropyridine, MPTP）诱导的PD小鼠模型中予以验证。PPI网络分析突出C?X?C基序趋化因子配体1（C?X?C motif chemokine ligand 1, CXCL1）、C?X?C趋化因子受体4型（C?X?C chemokine receptor type 4, CXCR4）、S100钙结合蛋白A12（S100 calcium?binding protein A12, S100A12）及C?X?C基序趋化因子配体11（C?X?C motif chemokine ligand 11, CXCL11）为核心枢纽基因（hub genes）；GSVA显示心肌收缩及氧化磷酸化（oxidative phosphorylation, OXPHOS）通路上调。C?C基序趋化因子配体4（C?C motif chemokine ligand 4, CCL4）、CKLF样MARVEL跨膜结构域包含蛋白2（CKLF like MARVEL transmembrane domain containing 2, CMTM2）、 Tudor结构域包含蛋白6（tudor domain containing 6, TDRD6）、钾电压门控通道亚家族S成员1（potassium voltage?gated channel subfamily S member 1, KCNS1）及RAS致癌基因家族成员RAB15（member RAS oncogene family, RAB15）被鉴定为候选标志物，且在临床样本及小鼠模型中表达改变一致。免疫浸润分析显示B细胞、CD8+T细胞及自然杀伤（natural killer, NK）细胞存在差异；LY?303511和HC?toxin被预测为候选化合物。这五个基因构成PD初步候选血液生物标志物，但仍需在更大规模独立队列中进一步验证以实现临床转化。

论文解读：《Neuroscience》——帕金森病候选生物标志物识别的探索性多组学分析

一、研究背景与意义

帕金森病（Parkinson's disease, PD）是一种复杂的多系统神经退行性疾病，核心病理改变为中脑黑质致密部（substantia nigra pars compacta, SNpc）多巴胺能神经元进行性丢失及纹状体（corpus striatum, CPU）突触前末梢变性，导致多巴胺耗竭。目前确诊依赖神经病理学检查，临床缺乏高灵敏度与特异度的早期筛查生物标志物。现有外周血候选标志物如α?突触核蛋白（α?synuclein, α?Syn）、神经丝轻链（neurofilament light chain, NfL）及炎症介质等单指标表现不稳定，且既往组学研究多局限于单一组学层面（转录组或蛋白质组单独分析），难以全面捕捉分子调控变化。尽管外周血与脑组织全局表达谱不同，但PD核心致病通路（免疫激活、炎症、氧化应激、趋化因子信号）在中枢与外周血中均可见一致性失调，使外周血分子特征可作为脑部病理的替代指示物。为此，研究人员开展外周血转录组与蛋白质组联合的多组学整合分析并结合机器学习，旨在克服单组学局限，筛选稳定、可重复的PD血液候选生物标志物，并为后续临床转化提供依据。

二、主要关键技术方法概述

研究人员采用四个数据集开展分析：(1) Dataset 1：河南大学淮河医院收集特发性PD患者（n=4）与正常对照（n=4）外周血行转录组（RNA?seq）测序；(2) Dataset 2：同批样本行蛋白质组检测获28个差异蛋白；(3) Dataset 3：GEO数据库中脑组织的2,964个PD相关差异表达基因（DEGs）用于交叉验证；(4) Dataset 4：外部独立GEO数据集验证候选标志物。分析流程包括：DEGs筛选（Dataset 1）、差异蛋白鉴定（Dataset 2）、蛋白?蛋白互作（PPI）网络构建（STRING数据库+Cytoscape筛选hub基因）、基因集变异分析（GSVA）富集通路、Dataset 1 DEGs与Dataset 3 DEGs取交集得28基因Panel、随机森林（RF）、支持向量机（SVM）及主成分分析（PCA）对候选基因优先级排序、外部数据集与临床血清样本及MPTP（20 mg/kg/天腹腔注射连续15天建立慢性PD小鼠模型）诱导小鼠模型验证、CIBERSORT评估免疫细胞浸润、Connectivity Map（CMap）预测候选小分子化合物。

三、研究结果

Differential gene analysis of dataset 1（数据集1差异基因分析）

对Dataset 1（4例PD vs 4例对照）进行差异表达分析（预设阈值），筛选出124个差异表达基因（DEGs），火山图与热图显示PD患者与对照组基因表达模式显著分离，初步锁定PD外周血转录组改变靶标（Target Group 1）。

Hub protein screening using dataset 2（基于数据集2的枢纽蛋白筛选）

整合Dataset 2获得的28个差异表达蛋白与124个DEGs，构建PPI网络。拓扑分析确定CXCL1、CXCR4、S100A12及CXCL11为高连接度的hub基因/蛋白，提示趋化因子信号及钙结合炎症相关分子在PD外周血中异常活化。

GSVA indicated upregulated cardiac muscle contraction and oxidative phosphorylation (OXPHOS)（GSVA显示心肌收缩与氧化磷酸化上调）

将Dataset 1的DEGs与Dataset 3（脑组织来源2,964个DEGs）取交集获得28基因Panel。对该Panel行GSVA通路富集，发现心肌收缩通路及氧化磷酸化（OXPHOS）通路在PD中显著上调，反映PD存在系统性能量代谢与线粒体功能相关紊乱。

Candidate marker prioritization and validation（候选标志物优先级排序与验证）

以28基因Panel为输入，分别用随机森林（RF）、支持向量机（SVM）及PCA进行特征重要性评估，综合筛选出CCL4、CMTM2、TDRD6、KCNS1及RAB15五个基因为最优候选标志物。该五基因在外部GEO验证队列、临床PD患者血清样本及MPTP诱导PD小鼠模型外周血中均呈现与发现队列一致的表达变化趋势（上调或下调），证实其跨平台与跨物种的可重复性。

Immune infiltration analysis（免疫浸润分析）

应用CIBERSORT算法评估22种免疫细胞亚群比例，结果显示PD患者外周血中B细胞、CD8⁺T细胞及自然杀伤（NK）细胞比例与对照组存在显著差异，支持PD存在外周免疫微环境重塑。

Candidate compound prediction（候选化合物预测）

基于差异表达基因表达谱，通过CMap数据库反向匹配，预测LY?303511（磷脂酰肌醇3?激酶/蛋白激酶B即PI3K/Akt通路抑制剂）和HC?toxin（组蛋白去乙酰化酶抑制剂）为可能逆转PD分子特征的候选小分子化合物，为后续药物重定位研究提供线索。

四、讨论与结论翻译总结

研究人员通过整合PD外周血转录组与蛋白质组数据，并与脑组织公共DEGs交叉验证，结合机器学习筛选出CCL4、CMTM2、TDRD6、KCNS1及RAB15五个基因组成PD初步血液候选生物标志物Panel；PPI分析提示CXCL1、CXCR4、S100A12、CXCL11为关键hub分子，GSVA揭示OXPHOS及心肌收缩通路异常，免疫浸润分析证实B细胞、CD8⁺T及NK细胞比例改变，CMap预测LY?303511和HC?toxin为潜在干预化合物。结论为：此五基因代表PD初步候选血液生物标志物，但其诊断效能与临床适用性尚需在大样本独立队列中进一步验证方可实现临床转化。