整合毒理学分析将大气污染相关急性呼吸窘迫综合征(ARDS)特征谱关联至苯并[a]芘(BaP)诱导的急性肺损伤(ALI)

《Ecotoxicology and Environmental Safety》：Integrated toxicological analysis links air pollution-associated acute respiratory distress syndrome signatures to benzo[a]pyrene-induced acute lung injury

　　

摘要：急性呼吸窘迫综合征(acute respiratory distress syndrome, ARDS)病死率高，且日益增多的证据表明空气污染可能参与其发生发展，但污染物暴露导致严重炎症性肺损伤的分子机制尚不完全清楚。研究人员整合本地临床关联数据、生物信息学分析、分子相互作用预测及体内毒理学验证，探讨污染相关呼吸损伤潜在机制。扬州地区月度ARDS入院人数与大气颗粒物浓度（PM2.5和PM10）呈正相关。分析GSE76293数据集，鉴定出167个大气污染相关靶点与ARDS差异表达基因(differentially expressed genes, DEGs)的交集基因，主要富集于炎症及免疫应答通路。机器学习筛选出三个候选靶点——NLRC4、ETS2和VNN1。鉴于多环芳烃(polycyclic aromatic hydrocarbons, PAHs)是颗粒物代表性毒性组分，选用苯并[a]芘(benzo[a]pyrene, BaP)进行体内验证。分子对接(molecular docking)及分子动力学(molecular dynamics, MD)模拟显示PAHs类化合物与候选蛋白可稳定结合。BaP染毒小鼠出现急性肺损伤(acute lung injury, ALI)，表现为肺组织病理损伤、肺湿/干重比(W/D ratio)升高及支气管肺泡灌洗液(bronchoalveolar lavage fluid, BALF)炎性介质升高，同时肺组织中NLRC4/GSDMD信号通路相关分子及ETS2、VNN1表达上调。上述发现支持颗粒物相关毒物与ARDS相关急性炎症性肺损伤之间存在潜在机制联系。

论文解读：《Integrated toxicological analysis links air pollution-associated acute respiratory distress syndrome signatures to benzo[a]pyrene-induced acute lung injury》发表于《Ecotoxicology and Environmental Safety》

一、研究背景与立题依据

急性呼吸窘迫综合征(acute respiratory distress syndrome, ARDS)是一种以急性起病的低氧血症、双侧非心源性肺水肿及广泛肺组织炎性损伤为特征的临床危重综合征，ICU患者ARDS发生率约10%，住院病死率高达40%，目前缺乏特效药物治疗。除感染、创伤等传统诱因外，环境暴露日益受到关注——短期及长期暴露于空气污染物（尤其是PM_2.5和PM₁₀）与ARDS发生风险增加及预后恶化相关，颗粒物可诱导肺部炎症、升高IL-6、IL-1β、TNF-α等炎性介质。然而，大气污染物（复杂混合物）致ARDS严重炎性肺损伤的具体毒理学分子机制尚未阐明。多环芳烃(polycyclic aromatic hydrocarbons, PAHs)是PM_2.5/PM₁₀中具代表性的有毒有机组分，其与ARDS相关炎性通路的关联尚不明确。为此，Song Lin、He Aifeng、Jiang Wei等研究人员采用临床流行病学关联分析联合生物信息学筛选、网络毒理学、机器学习、分子对接/分子动力学模拟及体内动物实验的整合策略，探究颗粒物相关毒物（以BaP为代表PAHs）致ARDS相关急性肺损伤的潜在分子机制与关键靶点。

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

研究人员收集2020年7月至2025年6月扬州地区ARDS入院病例及同期大气污染物（PM_2.5、PM₁₀、SO₂、CO、NO₂、O₃）月均浓度行Pearson相关性分析；从CTD、Swiss Target Prediction、TargetNet数据库获取5种空气污染物（CO、NO₂、O₃、PAHs、SO₂）相关靶点基因，从GEO数据库下载ARDS数据集GSE76293（12例ARDS患者外周血中性粒细胞/多形核中性粒细胞与12例健康对照）筛选差异表达基因(differentially expressed genes, DEGs)（|log₂FC|≥2，adj. p<0.05），取污染物靶点与ARDS DEGs交集；对交集基因行STRING-PPI网络构建（Cytoscape分析）、GO功能注释与KEGG通路富集（Metascape平台）；采用LASSO回归（十折交叉验证）与随机森林(Random Forest)算法筛选核心基因；利用CIBERSORT评估免疫细胞浸润并与核心基因做Pearson相关分析；对核心蛋白与污染物小分子行AutoDock Vina分子对接，选取NLRC4-PAHs复合物行100 ns GROMACS分子动力学(molecular dynamics, MD)模拟；选用苯并[a]芘(benzo[a]pyrene, BaP)经鼻滴注给予C57BL/6J小鼠连续14天建立急性肺损伤模型，检测肺组织H&E染色病理、肺湿/干重比(W/D)、BALF中IL-6、IL-1β、IL-18、TNF-α（ELISA）、肺组织核心基因mRNA（qRT-PCR）及蛋白（Western Blot、免疫荧光）。

三、研究结果

3.1. Ambient particulate matter was associated with local ARDS admissions

研究人员对扬州地区2020—2025年月均大气污染物浓度与月度ARDS入院人数行Pearson相关分析，结果显示月度ARDS入院数与PM₁₀（R=0.363，p=0.0043）及PM_2.5（R=0.273，p=0.035）呈显著正相关，而与SO₂、CO、NO₂、O₃无统计学显著相关；PM_2.5与PM₁₀、SO₂、CO强正相关，O₃与多数气态污染物弱或负相关。提示本地ARDS住院与颗粒物污染存在关联性，支持颗粒物（含PAHs组分）为ARDS相关严重呼吸损伤的重要环境暴露因素。

3.2. Transcriptomic analysis identified ARDS-associated gene signatures and key co-expression modules

对GSE76293数据集标准化后筛选ARDS vs 对照DEGs（|log₂FC|≥2，adj. p<0.05），热图与火山图显示两组基因表达谱明显分离。加权基因共表达网络分析(weighted gene co-expression network analysis, WGCNA)确定最佳软阈值幂β=8，识别出若干共表达模块，其中turquoise模块与ARDS表型相关性最强（p<1e-200），模块内基因连接度(module membership, MM)与基因显著性(gene significance, GS)高度相关（r=0.94，p<1e-200），表明该模块包含ARDS关键枢纽基因。

3.3. Pollutant-related candidate genes were enriched in inflammatory and immune-response pathways

从CTD等3个数据库获取5种空气污染物相关基因去重得17995个，与ARDS DEGs取交集获得167个共同基因。STRING构建PPI网络，GO富集主要条目为炎症反应(inflammatory response)、免疫活性调节、淋巴细胞活化调控、MHC II类分子结合等；KEGG富集涉及免疫调节、炎症级联与细胞应激信号通路。表明污染相关ARDS特征基因主要富集于炎症与免疫应答生物学程序。

3.4. Machine learning prioritized NLRC4, ETS2, and VNN1 as candidate targets

对167个交集基因分别构建Random Forest模型（按MeanDecreaseGini排序）和LASSO回归（十折交叉验证选定λ），两算法共同筛选出NLRC4、ETS2、VNN1三个核心候选基因。GSE76293中三者于ARDS组较对照组显著上调（p<0.001）。CIBERSORT免疫浸润分析显示ARDS存在固有/适应性免疫重塑，三基因与不同免疫细胞亚型具特异性相关性，提示其可能参与ARDS免疫微环境重构。

3.5. Molecular docking and dynamics analyses supported stable interactions between PAH-related compounds and candidate proteins

分子对接显示五种污染物小分子与三候选蛋白对接结合能均为负值（有利结合），其中PAHs与NLRC4（-10.20 kcal/mol）、VNN1（-9.80 kcal/mol）、ETS2（-7.40 kcal/mol）结合能显著低于气态污染物（-1.40～-3.90 kcal/mol）。PyMOL展示PAHs在三者蛋白口袋中的结合构象。对NLRC4-PAHs复合物行100 ns MD模拟，均方根偏差(root mean square deviation, RMSD)、均方根涨落(root mean square fluctuation, RMSF)、回转半径(radius of gyration, Rg)及埋藏溶剂可及表面积(solvent-accessible surface area, SASA)表明复合物在模拟期间保持稳定，二维相互作用图显示疏水作用与氢键参与配体结合。ETS2-PAHs与VNN1-PAHs复合物MD也显示稳定行为，结合自由能均为负值，支持PAHs可与三候选靶蛋白自发稳定结合。

3.6. BaP exposure induced acute inflammatory lung injury and was accompanied by increased candidate target expression and NLRC4/GSDMD-related signaling in mice

C57BL/6J小鼠经鼻给予BaP（2 mg/kg，连续14 d）建立模型。H&E染色示BaP组肺泡结构破坏、炎细胞浸润、肺泡隔增厚；肺W/D比值显著升高，提示肺水肿加重；BALF中IL-6、IL-1β、IL-18、TNF-α水平均显著高于对照组。Western Blot显示肺组织ETS2、VNN1蛋白表达上调，qRT-PCR显示NLRC4、ETS2、VNN1 mRNA均上调，免疫荧光示NLRC4与Gasdermin D（GSDMD）荧光信号增强。表明BaP诱导小鼠急性肺损伤伴随三候选靶点表达升高及NLRC4/GSDMD(pyroptosis executor)信号激活。

四、讨论与结论总结（翻译浓缩）

讨论指出本研究整合临床关联、转录组筛选、网络毒理学、机器学习、分子对接/MD及体内验证，鉴定出167个污染-ARDS交集基因富集于炎症免疫通路，机器学习锁定NLRC4（NLRC4炎症小体激活可致GSDMD依赖细胞焦亡及IL-1β/IL-18释放，促肺泡上皮/内皮损伤与肺水肿）、ETS2（转录因子，可通过MAPK/NF-κB调控炎症基因）、VNN1（参与氧化应激与炎症调节）为核心靶点。BaP作为代表性PAH，在体内诱导ALI样改变并上调三靶点与NLRC4/GSDMD信号，计算机模拟支持PAHs与三靶点稳定结合。局限性包括单中心月度聚合数据分析未校正COVID-19混杂及短时暴露滞后效应、BaP不能完全代表大气颗粒物复杂混合物、GEO数据集样本量偏小、未行基因敲除/挽救实验验证因果性及终末麻醉剂乌拉坦潜在肺影响。

研究结论： 本研究结果支持颗粒物相关毒物（以PAHs如BaP为代表）与ARDS相关急性炎症性肺损伤间存在潜在机制联系；PAHs可能通过协调调控NLRC4、ETS2和VNN1的表达及NLRC4/GSDMD相关信号通路参与此过程，其中NLRC4炎症小体-细胞焦亡(pyroptosis)通路尤具进一步研究价值。这些结果为未来ARDS环境污染易感机制研究与候选干预靶点提供了毒理学依据。