《Frontiers in Immunology》：Transcriptomic analysis reveals immune dysregulation and identifies key genes in ICU patients with severe ARDS
编辑推荐：
目的： 急性呼吸窘迫综合征（ARDS）以严重的免疫失调为特征，但其分子决定因素仍不明确。本研究旨在描绘ICU中ARDS患者的免疫失衡图谱，并通过体外实验验证候选枢纽基因的表达及潜在功能相关性。
方法： 合并批量转录组数据集，并利用差异表达、加权基因共表达网络分
目的： 急性呼吸窘迫综合征（ARDS）以严重的免疫失调为特征，但其分子决定因素仍不明确。本研究旨在描绘ICU中ARDS患者的免疫失衡图谱，并通过体外实验验证候选枢纽基因的表达及潜在功能相关性。
方法： 合并批量转录组数据集，并利用差异表达、加权基因共表达网络分析（WGCNA）和机器学习方法进行分析。评估功能富集和细胞去卷积，随后进行单细胞转录组验证。进行了脂多糖（LPS）诱导的THP-1细胞体外实验以确认候选基因表达。
结果： 联合分析强调了免疫相关通路，并揭示了先天性和适应性免疫亚群的显著改变。两个组蛋白相关基因H2BC4和H2BC12成为候选枢纽基因，它们在髓系群体中优先表达，并在炎症刺激下可诱导上调。
结论： 本研究为ARDS免疫发病机制提供了新见解，并确定了潜在的分子靶点，可为未来的诊断和治疗策略提供参考。
研究背景与意义
急性呼吸窘迫综合征（ARDS）是重症监护病房（ICU）常见的危重症，尽管支持治疗有所进步，其死亡率仍高达30%至40%。ARDS的发病机制涉及肺泡上皮损伤、肺内皮屏障破坏及失控的炎症反应，导致呼吸衰竭。然而，其精确的分子机制仍未完全阐明，缺乏有效的药物治疗手段。越来越多的证据表明，免疫失调在ARDS的发生发展中起核心作用，包括中性粒细胞和单核细胞等先天免疫细胞的过度活化，以及适应性免疫反应的受损。尽管已有研究尝试表征ARDS的免疫相关特征，但对驱动这种免疫失衡的分子决定因素仍缺乏全面理解。随着高通量测序技术的发展，转录组学分析已成为系统解析疾病相关分子网络的强大工具。因此，本研究旨在通过探索转录组学分析结合网络和机器学习方法，描绘ICU中ARDS患者的免疫失衡图谱，确定具有强烈免疫学相关性的枢纽基因，并通过单细胞测序和实验验证其表达模式，从而为ARDS的免疫发病机制和潜在治疗策略提供新见解。该研究成果发表在《Frontiers in Immunology》。
关键技术方法
研究人员从基因表达综合（GEO）数据库获取了ARDS转录组数据，包括GSE200847（41例ARDS患者和5例对照）、GSE163426（32例ARDS患者和5例对照）以及单细胞数据集GSE200848。通过差异表达分析筛选基因，利用加权基因共表达网络分析（WGCNA）构建基因模块，并应用最小绝对收缩和选择算子（LASSO）回归与支持向量机-递归特征消除（SVM-RFE）两种机器学习算法识别枢纽基因。采用CIBERSORT算法进行免疫细胞浸润分析。利用单细胞RNA测序数据（Seurat包处理）验证枢纽基因的特定细胞群表达。体外实验则使用LPS诱导的THP-1细胞模型，通过酶联免疫吸附试验（ELISA）和实时荧光定量PCR（qRT-PCR）进行验证。
研究结果
ARDS中差异表达基因的鉴定
通过对GSE200847和GSE163426两个独立队列的差异表达分析，分别鉴定出2628个和1569个差异表达基因（DEGs）。值得注意的是，组蛋白基因H2BC12（HIST1H2BK）和H2BC4（HIST1H2BC）在两个数据集中均位于上调区域。
批次效应评估与校正
在合并数据集前，主成分分析（PCA）显示存在批次效应。经ComBat函数校正后，不同队列样本在PCA空间中的重叠显著增加，有效消除了批次效应。
WGCNA识别ARDS相关模块
对合并数据集进行WGCNA分析，设定软阈值β为10以构建无标度拓扑结构。共识别出28个模块，其中darkseagreen4模块与ARDS临床状态相关性最强，包含245个基因。将该模块基因与先前鉴定的DEGs取交集，得到159个重叠基因。
机器学习识别稳健的候选枢纽基因
对上述159个重叠基因进行机器学习分析。LASSO回归结合留一法交叉验证（LOOCV）筛选出21个基因，SVM-RFE分析筛选出14个基因。两者取交集后，最终确定H2BC4和H2BC12为候选枢纽基因，并在两个数据集中验证了其在ARDS组的差异表达模式。
功能富集分析揭示ARDS相关关键生物学过程
对WGCNA获得的darkseagreen4模块中的245个基因进行富集分析。结果显示这些基因主要富集于核糖体/翻译相关结构与功能、先天免疫反应、蛋白质转运、氧化磷酸化和蛋白酶体途径。具体包括细胞质大核糖体亚基、细胞质翻译、RNA结合、CCR5趋化因子受体结合以及氧化磷酸化等GO和KEGG条目。
免疫浸润分析揭示ARDS中独特的免疫改变
基于CIBERSORT的细胞去卷积分析显示，ARDS与对照组相比，免疫细胞组成存在显著差异。具体表现为初始B细胞减少，静息CD4记忆T细胞显著减少而活化的CD4记忆T细胞增加，γδ T细胞轻度增加。髓系细胞中，M0巨噬细胞减少，M1巨噬细胞显著增加，静息树突状细胞增多，嗜酸性粒细胞显著减少。
单细胞转录组分析验证枢纽基因的特异性细胞表达
利用GSE200848数据集（包含8名ARDS患者的18717个高质量细胞）进行单细胞分析。UMAP聚类区分了多种免疫和非免疫细胞群。表达图谱显示，H2BC4和H2BC12主要在髓系来源的群体（包括中性粒细胞、单核细胞和巨噬细胞）中表达，而在T细胞、B细胞和NK细胞等淋巴亚群中表达极少。
LPS诱导THP-1细胞模型中候选基因的实验验证
在LPS诱导的THP-1细胞炎症模型中，RT-qPCR分析显示H2BC4和H2BC12在对照组中表达极低或无法检测到，而在LPS刺激后显著上调。同时，ELISA检测证实LPS处理后细胞上清中IL-6和TNF-α的分泌水平较对照组显著增加。
讨论与结论
本研究通过整合批量转录组学、网络分析、机器学习和单细胞转录组学，系统地表征了ICU中ARDS患者的免疫失衡图谱，并通过体外实验进行了验证。研究确定了两个组蛋白相关基因H2BC4和H2BC12作为ARDS相关的候选枢纽基因。这些基因在髓系群体中优先表达，并在炎症刺激下可诱导上调，为ARDS免疫失调的分子基础提供了新见解。与以往主要关注细胞因子和传统炎症介质的研究不同，本研究将组蛋白相关基因确定为潜在贡献者。功能富集分析进一步揭示，ARDS相关基因主要参与蛋白质合成、能量代谢和先天免疫过程相关的通路，表明ARDS中激活的免疫细胞经历了显著的代谢重编程和蛋白质周转增强。免疫浸润分析证实了向髓系主导的免疫反应转变，反映了促炎状态的过渡。单细胞分析提供了细胞水平的分辨率，证实了H2BC4和H2BC12在髓系细胞中的特异性表达。体外实验支持了其在炎症反应中的潜在作用。这些发现提示，H2BC4和H2BC12可能通过调节染色质可及性和转录活性，促进促炎转录程序的激活，从而加剧ARDS的炎症环境。综上所述，本研究为ARDS免疫发病机制提供了新视角，并确定了潜在的分子靶点，有助于开发基于免疫的诊断和治疗策略，强调了系统表征免疫失调对于未来ARDS精准医疗的重要性。
