创伤诱导基因型关联性低炎症与危重症中17型免疫应答核心作用的罕见单核苷酸多态性网络分析

《Frontiers in Systems Biology》：Network analysis of rare single-nucleotide polymorphisms suggests a central role for type 17 immune responses in trauma-induced, genotypically associated hypo-inflammation and critical illness

　　

摘要 引言：严重创伤性损伤导致的危重症是全球发病率和死亡率的主要原因，涉及多个器官系统的功能障碍，其驱动因素至少部分与包括17型免疫在内的失调性炎症有关。研究人员及其他团队已证实，与新型单核苷酸多态性（SNPs）相关的遗传易感性可导致不良的危重症结局。 方法：在本研究中，研究人员旨在定义在危重症创伤患者中评估的罕见SNPs（由于研究队列规模，操作性定义为存在于5%-10%的患者中）之间的相互作用网络。 结果：RAR相关孤儿受体A（RORA），一种对Th17细胞分化至关重要的转录因子，通过其SNP rs4774381被推断为枢纽基因，该SNP与低炎症及基因型关联性器官功能障碍轨迹相关。进一步与Th17反应相关的是，非枢纽基因TM4SF19中的SNP rs11919443定义了基因型关联性全身性IL-17A轨迹，并将低炎症与不良结局相关联。相比之下，在对照组SNPs中未识别出Th17相关枢纽基因。二次分析显示，在该对照组中，循环炎症介质（不包括IL-17A）存在若干统计学显著差异，并且同样将低炎症与不良结局相关联。 讨论：因此，本研究拓展了研究人员先前旨在定义危重症背景下遗传易感性导致失调性炎症和病理生理学的工作，指出了17型免疫应答的关键作用，并强调了低炎症作为危重症中一种病理性免疫表型的地位。

论文解读文章

研究背景、存在问题与研究目的

严重创伤性损伤引发的危重症是全球范围内致死致残的主要病因之一，其核心病理特征是多器官功能障碍综合征（Multiple Organ Dysfunction Syndrome, MODS），而失调性炎症被认为是驱动这一过程的关键因素。既往研究表明，白细胞介素-17A（Interleukin-17A, IL-17A）及17型免疫应答在创伤后系统性炎症的放大和器官功能障碍中扮演早期重要角色。尽管临床采用标准化治疗方案，具有相似损伤特征和人口学资料的危重症创伤患者预后仍呈现高度异质性。研究人员及同行提出，患者间的遗传变异可能是这种异质性的重要来源之一。此前的研究已识别出与亚急性死亡率相关的特定单核苷酸多态性（Single-Nucleotide Polymorphisms, SNPs），这些SNPs与早期分化的、以Th17为中心的炎症轨迹及MODS相关。然而，既往研究多聚焦于常见SNPs（在队列中频率≥20%），对于频率较低（5%-10%）的罕见SNPs及其构成的基因网络在创伤后炎症失调中的作用尚不清楚。因此，本研究旨在利用新型生物信息学工具，系统扫描超过500,000个SNPs，以定义与不良预后相关的罕见和常见SNPs网络，重点识别枢纽基因，并根据枢纽基因位点的基因型对患者进行分层，评估其临床结局和炎症谱，从而深入揭示遗传因素如何通过调控17型免疫影响创伤后炎症反应和危重症预后。

主要关键技术方法

1. 1.
   患者队列与样本来源：研究对象为380例钝挫伤存活患者，来源于匹兹堡大学医学中心长老会医院（UPMC Presbyterian hospital，一级创伤中心）急诊科，入组时间为2004年5月5日至2012年5月2日。排除标准包括孤立性头部损伤、妊娠和穿透性创伤。
2. 2.
   基因组分析：从全血样本中提取DNA，使用Illumina? Human Core Exome-24 v1.1 BeadChip芯片进行SNP基因分型，该芯片可检测551,839个SNPs。数据分析使用Illumina Genome Studio 2.0软件。
3. 3.
   炎症生物标志物检测：使用Luminex? 100 IS分析仪和多重试剂盒测定血浆中多种细胞因子/趋化因子的水平，包括IL-17A、肿瘤坏死因子-α（TNF-α）、粒细胞-巨噬细胞集落刺激因子（GM-CSF）等。同时使用ELISA法检测ST2水平，Griess试剂比色法检测NO₂^-/NO₃^-水平。
4. 4.
   生物信息学分析：运用自主研发的Python算法SNPScanner筛选与炎症介质差异相关的SNPs。利用网络应用程序FLAME（v2.0）进行功能富集分析，识别枢纽基因。采用动态网络分析（Dynamic Network Analysis, DyNA）构建随时间变化的炎症介质相互作用网络。通过斯皮尔曼相关性分析推断Th17细胞亚群。

研究结果

3.1 生物信息学分析定义了钝挫伤患者罕见SNPs的相互作用网络，并提示Th17分化与系统性低炎症

研究人员将SNPScanner筛选出的与两个高显著性差异炎症生物标志物相关的12个体细胞罕见SNPs输入FLAME进行分析，识别出一个主要枢纽基因和三个次要枢纽基因。其中，位于RORA基因上的SNP rs4774381被鉴定为主要枢纽基因，RORA是Th17细胞分化的关键转录因子。

3.1.1 SNP rs4774381提示主要枢纽基因RORA参与差异生物标志物浓度和创伤后结局

按rs4774381位点纯合基因型（AA vs. BB）对患者分层后，研究发现两组患者在循环IL-17A、IFN-γ、IL-1RA等多种炎症介质上存在统计学显著差异。尽管rs4774381^AA患者表现出更显著的系统性炎症升高，但其Marshall MOD评分却显著低于rs4774381^BB患者，提示可能存在病理性低炎症状态。DyNA结果显示，rs4774381^AA患者的早期炎症网络连接稀疏，而rs4774381^BB患者则高度连接，两者呈现相位偏移式的振荡模式。

3.1.2 RORA SNP联合Th17差异介质分析提示非致病性Th17细胞上调与创伤后低炎症相关

通过相关性分析推断Th17细胞亚群，发现仅在rs4774381^AA患者中存在IL-17A与抗炎因子IL-10的显著正相关，表明非致病性（non-pathogenic）Th17细胞占优势。两组均未观察到代表致病性Th17细胞的IL-17A与GM-CSF的相关性。

3.1.3 其他主要枢纽基因RPP38和UBE2D2的SNPs对创伤患者预后的分层程度低于RORA SNP rs4774381

位于RPP38基因的SNP rs7918717被识别为第二个主要枢纽基因，其纯合子亚组在GM-CSF水平和MODS轨迹上存在差异。位于UBE2D2基因的SNP rs261537被识别为第三个主要枢纽基因，其纯合子亚组在IFN-γ和GM-CSF水平上存在差异，但无显著临床结局差异。

3.2 TM4SF19基因中的rs11919443与多种系统性炎症介质（包括IL-17A）的不同轨迹及低炎症相关，但与不同的Th17细胞亚群无关

位于非枢纽基因TM4SF19的SNP rs11919443，其纯合子亚组（AA vs. BB）在循环IL-17A、IL-6、IL-8等多种炎症介质上存在显著差异。rs11919443^AA患者显示出更长的ICU住院时间和机械通气天数。DyNA分析显示，rs11919443^BB患者的炎症网络连接几乎为零，呈现明显的低炎症表型。与RORA SNP不同，该位点的基因型分组并未导致显著的Th17亚群推断差异。

3.3 对对照SNP组的分析——即算法定义为基因型关联性临床结局差异但无循环炎症介质差异——提示SLC2A1的作用

研究人员选择了一组在炎症介质上无差异但在临床结局上有差异的罕见SNPs作为对照。FLAME分析将其中的SLC2A1（编码葡萄糖转运蛋白GLUT1）识别为唯一枢纽基因。按rs841852位点基因型分层后，BB纯合子患者表现出更差的临床结局（更长的ICU住院时间、总住院天数和更高的MOD评分）。DyNA再次揭示了BB纯合子患者的低炎症表型（零网络连接）。

讨论与结论

本研究通过整合新型计算生物学方法，系统分析了创伤后危重症患者中罕见SNPs的网络，有力地证实并拓展了遗传易感性在创伤后失调性炎症和病理生理学中的关键作用。研究首次明确指出，RAR相关孤儿受体A（RORA）转录因子作为枢纽基因，其SNP rs4774381与Th17分化及随后的系统性低炎症密切相关。这一发现将先前关于IL-17A和致病性Th17细胞驱动过度炎症导致不良预后的认识，延伸至另一个维度：即遗传决定的非致病性Th17细胞上调和IL-10升高，可能构成一种病理性低炎症表型，同样与不良预后相关。此外，研究还发现了TM4SF19基因的rs11919443和SLC2A1/GLUT1基因的rs841852与低炎症和不良临床结局的关联，进一步丰富了遗传因素影响创伤后免疫反应的复杂性。研究结论指出，这项工作突出了17型免疫应答在危重症中的核心作用，并强调了低炎症作为一种重要的病理性免疫表型。这为未来基于基因型和免疫表型的精准医疗策略，如针对特定炎症通路（如IL-17A）或代谢通路（如葡萄糖转运）进行干预，提供了新的理论基础和潜在靶点。