PM2.5多污染物混合暴露加剧帕金森病神经毒性的机制研究：基于INPP5D/SHIP1-PI3K-AKT轴的整合多组学分析

《Journal of Xenobiotics》：Neurotoxicity of PM2.5-Bound Multi-Pollutant Mixtures in Parkinson’s Disease: Exploring the INPP5D/SHIP1-PI3K-AKT Axis via an Integrative Multi-Omics Study

　　

越来越多的流行病学证据将环境PM2.5暴露与帕金森病（Parkinson's disease, PD）风险升高联系起来，然而在真实共暴露场景下PM2.5多污染物混合物所引发的神经毒性仍未得到充分阐明。本研究整合了两样本孟德尔随机化（Mendelian randomization, MR）、面向混合物的网络毒理学、多组学生物信息学、机器学习及体外细胞验证，系统揭示了混合PM2.5污染物诱导PD发病的因果分子网络。基于公开全基因组关联研究（genome-wide association study, GWAS）数据集的两样本MR分析提示，PM2.5暴露与PD易感性升高之间存在潜在的正向遗传关联。研究人员将PM2.5组分靶点、PD风险基因及差异转录本取交集，获得195个共享基因，其中INPP5D（编码SHIP1蛋白）通过LASSO和随机森林算法被鉴定为核心基因。分子对接模拟进一步预测，金属-有机污染物复合物与SHIP1的结合亲和力可能强于单一有机污染物。通路富集筛选出PI3K-AKT级联作为核心调控通路。在SH-SY5Y细胞中，NIST SRM 1648a（城市颗粒物标准物质）处理上调SHIP1表达及磷酸化水平，抑制PI3K-AKT活性，并促进神经炎症和细胞凋亡，而INPP5D敲低可逆转上述病理改变。总体而言，本研究的整合多组学及体外实验证据提出一个合理的机制假说：NIST SRM 1648a模拟的PM2.5相关污染物混合物可能通过SHIP1/PI3K-AKT信号轴加剧PD样神经元损伤。这一多维度证据为混合污染物诱导的神经退行性变提供了分析范式，并突出INPP5D作为PM2.5相关PD的潜在干预靶点。

一、研究背景与科学问题

帕金森病（Parkinson's disease, PD）是全球第二大常见神经退行性疾病，仅次于阿尔茨海默病。随着全球人口老龄化加速和环境暴露因素日益复杂，PD患病率显著上升，已成为目前增长最快的神经系统疾病。流行病学预测显示，到2050年全球PD患者将超过2520万，而中国作为人口最多的国家届时将占全球病例的近一半。既往研究表明，PD发病具有高度异质性，仅5%–15%的患者存在明确的遗传病因或可识别的遗传风险因素，提示环境因素可能在散发性PD发生发展中发挥更为关键的作用。在众多环境风险因素中，细颗粒物（PM_2.5，空气动力学直径≤2.5 μm）作为空气污染的主要成分，已被流行病学研究证实与PD风险增加相关。例如，一项覆盖中国珠江三角洲地区近15万参与者的前瞻性队列研究发现，居住在PM_2.5年均浓度超过50 μg/m³地区的个体，其PD发病风险比低于36 μg/m³地区的个体高1.5倍。

从致病机制角度看，PM_2.5含有可穿越血脑屏障的超细颗粒，可直接浸润中枢神经系统，诱导慢性神经炎症、氧化应激损伤及小胶质细胞异常激活，这些病理过程正是PD发生发展的核心驱动因素。然而，在真实共暴露场景下，PM_2.5多污染物混合物的神经毒性机制仍不明确。因此，系统阐明PM_2.5暴露与PD风险之间的因果关联及潜在分子机制，不仅有助于明确不同空气污染源对PD发病的影响权重，也为制定针对性环境干预策略和PD防治方案提供重要的理论依据。本研究旨在系统表征PM_2.5来源的多污染物混合物在PD发病中的神经毒性网络，并探索INPP5D（编码SHIP1蛋白）作为介导潜在联合毒性损伤的核心枢纽的作用。值得注意的是，INPP5D是神经炎症的关键调控因子，已被证实参与包括阿尔茨海默病在内的神经退行性疾病，但其在PD中尤其是在PM_2.5暴露下的具体机制仍鲜有研究。

二、研究内容与主要结论

研究人员整合两样本孟德尔随机化（MR）、网络毒理学、多组学生物信息学、分子对接、基因集富集分析及体外细胞实验等多种技术，系统揭示了PM_2.5组分通过调控INPP5D介导PD发生发展的分子网络。研究结果提示，INPP5D作为核心信号分子，在PM_2.5诱导的神经元损伤中发挥关键中介作用。该工作为环境污染物诱导神经退行性疾病的分子网络提供了新的理论见解，同时为探究PD环境病因及制定靶向干预策略奠定了基础。

三、主要技术方法

本研究采用的主要关键技术方法包括：基于公开GWAS数据的两样本孟德尔随机化分析（暴露数据来自ukb-b-10817，423,796名欧洲参与者；结局数据来自finn-b-G6_PARKINSON_INCLAVO，含2,167例病例和216,625例对照）；基于OMIM、GeneCards和CTD数据库的PD风险基因筛选；基于CTD数据库的六种PM_2.5有害组分（苯、苯并[a]芘、甲醛、铅、砷、镉）人类靶基因鉴定；基于GSE49036队列的差异表达基因分析；蛋白-蛋白相互作用（PPI）网络构建；LASSO回归与随机森林（RF）机器学习算法筛选核心基因；基于PDB蛋白结构（6XY7）的AutoDock分子对接模拟；基于GSE157783数据集的单细胞转录组分析；以及SH-SY5Y细胞体外验证实验。

三、研究结果

3.1 孟德尔随机化分析PM_2.5暴露与PD风险的关联

基于公开GWAS汇总数据，逆方差加权（IVW）MR分析提示，基因预测的较高PM_2.5暴露水平与PD风险升高之间存在临界显著的正向因果关联（OR = 1.888，95% CI 1.005–3.521，p = 0.048）。所有工具变量F统计量均大于10，排除弱工具变量偏倚。多种敏感性分析方法（MR-Egger、加权中位数、加权众数、简单众数）虽未达统计学显著性，但效应方向均与IVW结果一致。Cochran-Q检验未检测到显著异质性（p = 0.746），MR-Egger截距检验（p = 0.889）和MR-PRESSO全局检验（p = 0.077）均未检测到显著水平多效性，支持因果估计的稳健性。研究提供PM_2.5与PD之间存在正向因果关联的提示性证据。

2.2 筛选PM_2.5暴露相关PD靶基因

研究人员从GeneCards、OMIM和CTD数据库获得PD相关风险基因集，对六种PM_2.5组分的靶基因集、已知PD风险基因集和PD队列差异表达基因进行交集分析，获得195个共享靶基因。基于这些基因构建的PPI网络呈现高度互联的拓扑结构。通过LASSO回归和随机森林算法筛选，两种方法取交集获得两个候选基因：SLC9A1和INPP5D。INPP5D（即SHIP1）作为磷酸肌醇磷酸酶家族成员，已被报道参与免疫调节、神经元存活、突触功能和神经炎症，其失调与神经退行性疾病的炎症微环境破坏和神经元损伤密切相关。在PD样本中INPP5D表达显著上调。鉴于INPP5D在神经炎症和神经元损伤中的明确作用及其在PD中的表达升高，研究人员优先选择INPP5D进行后续机制研究。

3.3 分子对接：PM_2.5组分与SHIP1的结合模式

分子对接结果显示，苯并[a]芘与SHIP1蛋白结合较强，而苯或甲醛单独结合较弱。当有机化合物与金属离子共存时，所有复合物与蛋白的结合均显著增强。苯并[a]芘主要通过疏水作用及与蛋白芳香族残基形成T型π–π堆积嵌入SHIP1催化域；苯主要依赖范德华力；甲醛则主要通过氢键和范德华力相互作用。金属离子与蛋白形成配位键：Pb²⁺与两个侧链羧基形成三个配位键，As³⁺与一个羧基和一个羟基形成两个配位键，Cd²⁺与一个羧基形成两个配位键。金属-有机复合物对SHIP1的结合稳定性显著增强，为真实环境共暴露下PM_2.5污染物混合物加剧神经元损伤提供了直接化学机制解释。

4. 单细胞转录组特征与INPP5D在神经元中的表达

对PD患者和健康对照的单细胞转录组数据进行UMAP降维分析，发现不同组别细胞在UMAP空间中呈现部分重叠和区域分离，表明两组间细胞组成存在差异。各细胞类型形成清晰独立的聚类。在总神经元群体中，INPP5D在PD组与对照组间无显著差异，提示PD黑质中INPP5D整体信号的升高并非源于神经元转录上调，而可能来自组织细胞组成比例的转变。进一步分析显示，兴奋性神经元中INPP5D基线表达最高，PD易感神经元（DaNs）也存在可检测的INPP5D表达。鉴于这些PD易感神经元具有基础INPP5D表达且是PD中受损的主要细胞群体，后续采用神经元样细胞建立体外PD损伤模型，探索内源性SHIP1蛋白激活是否独立于INPP5D转录上调而介导神经元损伤。

5. 功能富集分析与PI3K-AKT通路特征

KEGG通路富集分析显示，195个共享靶基因主要富集于神经退行性疾病、细胞应激与生存调控、免疫与炎症相关通路，其中PI3K-AKT信号通路富集得分最高，提示其作为连接PM_2.5暴露与PD病理的核心通路。GO功能富集分析显示，靶基因显著富集于细胞应激反应和细胞命运调控相关的生物过程，核心条目包括细胞凋亡、氧化应激和炎症反应。GSEA分析显示PD组PI3K-AKT信号通路显著下调（NES = ?1.527153468，FDR = 0.0005284），表明该通路活性整体受抑与PD发病进展密切相关。INPP5D编码的SHIP1蛋白特异性水解第二信使磷脂酰肌醇(3,4,5)-三磷酸[PI(3,4,5)P₃]的5'-磷酸基团，生成磷脂酰肌醇(3,4)-二磷酸[PI(3,4)P₂]，从而负向调控PI3K-AKT信号通路的激活。INPP5D活性增加可直接导致磷酸化AKT（p-AKT）水平降低，进而影响细胞代谢和存活相关下游信号转导。

6. SRM 1648a暴露激活SHIP1并抑制PI3K-AKT信号轴促进体外PD模型细胞凋亡

采用NIST SRM 1648a城市颗粒物标准物质模拟真实复杂PM_2.5暴露条件。经RA诱导后，分化SH-SY5Y细胞呈现细长形态，并上调神经元特异性蛋白MAP2表达。MPP⁺和SRM 1648a均剂量依赖性地降低细胞活力。Western blot结果显示，MPP⁺和SRM 1648a处理组中SHIP1及其磷酸化形式（p-SHIP1）表达均升高，联合暴露进一步增强。p-PI3K表达显著下调，p-AKT和总AKT表达均降低，SRM 1648a刺激进一步加剧这些指标的下调。在联合暴露下，促炎因子IL-6和凋亡相关蛋白cleaved Caspase-3（c-Caspase-3）表达显著升高。结果表明，SHIP1的高表达和激活是SRM 1648a加剧PD病理损伤的关键机制，涉及抑制PI3K-AKT通路活性和上调炎症及凋亡相关蛋白。

7. INPP5D敲低通过缓解PI3K-AKT通路抑制保护神经元样细胞

通过siRNA介导的INPP5D敲低实验，INPP5D蛋白水平降低22%。与si-NC对照组相比，si-NC + MPP⁺ + 1648a组的结果与前述一致。而INPP5D基因敲低（si-INPP5D + MPP⁺ + 1648a）则逆转了这些病理表型：p-PI3K和p-AKT表达水平恢复，IL-6和cleaved Caspase-3表达降至与基线组相当的水平。虽然总PI3K水平在所有条件下保持不变，但INPP5D敲低显著恢复了总AKT表达。结果表明，INPP5D在SRM 1648a加剧PD相关病理中发挥关键致病作用：在SRM 1648a与MPP⁺共暴露时，该基因通过抑制PI3K-AKT生存通路驱动神经炎症和凋亡反应，而单独敲低时则未产生这些病理效应。

四、讨论与结论

本研究整合MR、网络毒理学、机器学习、分子模拟和细胞实验等多种技术，系统阐明了PM_2.5组分通过调控INPP5D介导PD发生发展的分子机制。两样本MR提供了PM_2.5暴露与PD风险关联的提示性遗传流行病学证据，与既往多项观察性研究一致。通过多数据库交叉分析，识别出PM_2.5有害组分靶基因与PD相关基因的重叠，构建了PPI网络，表明PD发生发展受多基因调控网络驱动。LASSO与RF联合应用提高了核心中介基因筛选的可靠性，确定INPP5D为核心候选基因。分子对接揭示金属-有机复合物与SHIP1结合增强，提示真实环境混合暴露可加剧神经毒性。单细胞分析表明INPP5D在神经元中表达，但整体信号升高可能源于组织细胞组成变化。功能富集和GSEA证实PI3K-AKT通路为连接PM_2.5暴露与PD的核心通路，且该通路在PD中显著下调。体外实验证实SRM 1648a通过激活SHIP1、抑制PI3K-AKT通路、上调IL-6和c-Caspase-3形成病理级联，而INPP5D敲低可逆转这些损伤。

研究结论：本研究构建了整合两样本孟德尔随机化、面向混合物的网络毒理学、多组学生物信息学和体外功能验证的证据链，提示PM_2.5多污染物混合物暴露是PD的潜在风险因素，并确定INPP5D/SHIP1作为联合毒性神经元损伤的核心介导因子。机制上，PM_2.5混合物中有机与金属毒物的联合作用可能上调SHIP1表达和磷酸化，抑制神经保护性PI3K-AKT信号级联，并引发持续性神经炎症和多巴胺能细胞凋亡。INPP5D敲低可有效逆转所有混合物诱导的病理改变。本研究建立了探索复杂环境污染物混合物驱动神经退行性变的整合分析框架，并强调INPP5D作为长期暴露于多污染物空气混合物人群PD预防的潜在治疗靶点。