整合人群证据、网络毒理学、单细胞分析及体外研究阐明多环芳烃与卒中关联的机制

《Ecotoxicology and Environmental Safety》：Integrated population-based evidence, network toxicology, single-cell profiling, and in vitro studies elucidate mechanisms connecting PAHs to stroke

　　

整合人群证据、网络毒理学、单细胞分析及体外研究阐明多环芳烃（PAHs）与卒中关联的机制。摘要：环境暴露于多环芳烃（PAHs）与心血管疾病相关，但导致卒中风险的具体代谢物和细胞类型特异性机制仍不清楚。在此，研究人员应用了一个整合的多层次框架，结合基于人群的暴露分析、网络毒理学、转录组学和单细胞验证，以2-萘酚（2-NAP）为代表代谢物，研究PAH暴露如何可能导致卒中。使用逻辑回归和限制性立方样条分析PAH代谢物与卒中的关系，同时使用加权分位数和（WQS）回归及机器学习模型评估混合效应。网络毒理学结合转录组学和机器学习分析用于识别候选基因，并优先确定了四个假定的核心基因（DPP4、IFIH1、KAT6A和LAMP2）。单细胞RNA测序进一步解析了这些基因的细胞类型特异性表达，突出了内皮细胞和巨噬细胞亚群，并涉及转化生长因子-β（TGF-β）相关信号通路。分子对接提示关键代谢物2-NAP与LAMP2之间存在直接相互作用，并在体外进一步验证。功能实验表明，2-NAP暴露可增加LAMP2表达，诱导内皮DNA损伤，损害小胶质细胞迁移，并增强促炎细胞因子的产生。总之，这些发现表明2-NAP可能是连接PAH暴露与卒中的重要因素，并且LAMP2相关通路可能参与PAH相关的脑血管损伤。本研究通过将流行病学关联与明确的分子通路和脑血管细胞群联系起来，提供了机制性见解，增强了PAH相关卒中风险的生物学合理性。

### 论文解读：多环芳烃暴露与卒中关联的机制：一项整合人群、网络毒理学、单细胞分析及体外研究

**研究背景与问题**

卒中（stroke）是全球主要死亡和长期致残原因，传统风险因素（如高血压、糖尿病、吸烟、血脂异常）已明确，但环境暴露在其发病中的作用仍认识不足。多环芳烃（polycyclic aromatic hydrocarbons, PAHs）是普遍存在的燃烧产物，可通过代谢产生活性氧、诱导DNA损伤及表观遗传改变，促进系统性炎症和内皮功能障碍。流行病学证据提示PAHs暴露与不良心血管结局相关，但具体代谢物和细胞类型特异性机制尚不明确。现有研究多聚焦复合心血管结局，缺乏对脑血管损伤独特生物学过程的解析；PAH暴露常被视为均质混合物，未识别主导卒中风险的代谢物；网络毒理学与转录组学分析常脱离人群暴露评估，且缺乏特定脑血管细胞类型的验证。因此，研究人员开展了一项整合多水平证据的研究，旨在将流行病学关联与分子通路及脑血管细胞群联系起来，为PAH相关卒中风险提供生物学合理性。

**研究内容与结论**

研究人员利用美国国家健康与营养调查（NHANES，2005–2016周期，8205名参与者，257例卒中事件）进行人群暴露分析，结合网络毒理学、转录组学、单细胞RNA测序、分子对接及体外实验，构建了从人群到机制的完整证据链。研究以2-萘酚（2-NAP）作为代表性PAH代谢物，揭示了其与卒中的关联，并优先识别出四个核心基因（DPP4、IFIH1、KAT6A、LAMP2），其中LAMP2在内皮细胞和巨噬细胞中上调，且TGF-β信号通路参与其中。体外实验证实2-NAP暴露可诱导内皮细胞DNA损伤、凋亡、屏障破坏，并损害小胶质细胞迁移、促进促炎因子分泌。该研究首次系统整合多层面证据，支持2-NAP作为PAH相关卒中风险的主要标志物，并为LAMP2相关通路在脑血管损伤中的作用提供了机制性证据。论文发表在《Ecotoxicology and Environmental Safety》。

**主要关键技术方法**

研究人员采用以下关键技术：1）人群分析：基于NHANES队列（2005–2016），使用多变量逻辑回归、限制性立方样条（RCS）评估单一PAH代谢物与卒中的剂量-反应关系，采用加权分位数和（WQS）回归及极端梯度提升模型（XGBoost）结合Shapley加法解释（SHAP）评估PAH混合物效应并排序贡献；2）网络毒理学：利用CTD、ChEMBL、SwissTargetPrediction、STITCH数据库预测2-NAP靶点，从GeneCards、TTD、OMIM收集卒中相关基因，取交集后构建蛋白-蛋白相互作用（PPI）网络并进行功能富集；3）机器学习：基于GSE16561转录组数据，使用LASSO逻辑回归和支持向量机递归特征消除（SVM-RFE）筛选核心基因；4）单细胞分析：分析GSE174574单细胞RNA测序数据，使用Seurat聚类、CellChat分析细胞通讯、Monocle2进行轨迹推断；5）分子对接：使用AutoDock Vina评估2-NAP及其代谢物与核心蛋白的结合；6）体外实验：使用bEnd.3（脑微血管内皮细胞）和BV2（小胶质细胞）细胞系，进行CCK-8、流式细胞术、γH2A.X蛋白印迹、Transwell迁移、ELISA、TEER及FITC-dextran通透性等检测。

**研究结果**

**3.1 流行病学识别2-NAP为与卒中相关的主要PAH**：通过多变量逻辑回归，发现1-NAP、2-NAP、2-羟基芴（2-FLU）浓度升高与卒中风险增加相关；RCS分析显示正向暴露-反应模式；WQS回归和XGBoost-SHAP分析一致将2-NAP识别为混合物中的主要贡献者。

**3.2 网络毒理学揭示连接2-NAP暴露与卒中的候选靶点**：从2-NAP靶点与卒中基因交集获得205个候选基因，PPI网络显示密集连接模块，功能富集分析涉及炎症信号、内皮调节和血管重塑通路。

**3.3 机器学习优先排序四个核心基因和免疫相关特征**：LASSO和SVM-RFE共同筛选出DPP4、IFIH1、KAT6A和LAMP2。免疫浸润分析显示卒中样本中巨噬细胞和中性粒细胞比例升高，T细胞和NK细胞减少；核心基因与髓系细胞及炎症/免疫调节通路显著相关，其中LAMP2关联最强。

**3.4 单细胞分析突出与PAH相关核心基因表达相关的内皮和巨噬细胞亚群**：单细胞RNA测序解析出19个细胞簇，TGF-β信号在内皮细胞和巨噬细胞中高度富集。核心基因在内皮细胞和巨噬细胞中差异表达，LAMP2在两种细胞中均上调。细胞通讯分析显示，高表达核心基因的内皮细胞（2-NAP_high_En）与髓系细胞、中性粒细胞和巨噬细胞相互作用强度显著高于低表达组。

**3.5 轨迹分析揭示卒中偏向的内皮状态转变**：Monocle2轨迹分析显示，从对照到卒中存在梯度转变，卒中细胞富集于疾病偏向的终末分支，提示疾病相关信号在轨迹后期达到峰值。

**3.6 分子对接支持2-NAP与核心靶点的直接相互作用**：2-NAP及其葡萄糖醛酸和硫酸结合物与四个核心蛋白均显示有利的预测结合亲和力。

**3.7 体外实验支持2-NAP相关的内皮损伤和小胶质细胞功能障碍**：2-NAP处理bEnd.3细胞后，诱导DNA损伤（γH2A.X升高）、LAMP2表达上调、凋亡率升高（最高达81.9%），并降低紧密连接蛋白claudin-5表达、降低跨内皮电阻（TEER）并增加FITC-右旋糖酐通透性，表明内皮屏障破坏。在BV2小胶质细胞中，2-NAP降低迁移能力，增加促炎细胞因子（IL-6、IL-1β、TNF-α）分泌。

**讨论与结论**

讨论部分指出，该研究构建了从流行病学、毒理学、转录组学到功能实验的证据链，系统阐明了PAH暴露（以2-NAP为代表）与卒中的关联。核心基因中，LAMP2在多个层面（人群、组织、单细胞、体外）表现出一致性，其上调与内皮DNA损伤、凋亡和屏障破坏平行，提示自噬-溶酶体轴在PAH驱动的炎症压力下参与损伤清除与免疫重塑。TGF-β信号通路在内皮-巨噬细胞相互通讯中发挥关键作用，可能通过内皮-间质转化（EndoMT）影响血脑屏障完整性。结论部分翻译如下：整合人群分析、网络毒理学、单细胞分析和实验验证，研究结果支持2-NAP作为卒中风险的主要PAH标志物，并阐明了LAMP2在卒中内皮细胞和巨噬细胞中的特异性及其可能作用通路。该研究为理解PAH相关卒中风险提供了机制性锚点，并提出了可用于风险分层和靶向干预的生物标志物。