自身免疫相关间质性肺病中系统性遗传易感性与肺组织重塑的整合研究

《Genes》：Systemic Genetic Susceptibility and Lung Tissue Remodeling in Autoimmune-Associated Interstitial Lung Disease

　　

背景：自身免疫相关间质性肺病（SAID-ILD）涉及在复杂的肺部炎症和纤维化基础上叠加的系统性免疫易感性。方法：本研究采用多层框架，将91种循环炎症蛋白和1400种血清代谢物的双样本孟德尔随机化（MR）分析与肺组织批量转录组学、外植肺单细胞核RNA测序及支气管肺泡灌洗单细胞RNA测序相结合，以将系统性易感性信号与细胞分辨的肺部生物学联系起来。结果：MR分析优先筛选出多种炎症信号，包括CD244和CX3CL1的风险增加估计值；两者均未达到91种蛋白的BH阈值，且OR=1.5时的检验效能均低于25%。批量转录组荟萃分析确定了一个可重复的疾病背景，以上皮损伤、髓系细胞活化和细胞外基质重塑为主导。细胞分辨整合将CD244主要定位于T细胞/自然杀伤细胞，将CX3CL1定位于内皮细胞、上皮细胞和基质细胞区室，支持将CD48–CD244和CX3CL1–CX3CR1作为候选通讯轴进行评估。CX3CL1–CX3CR1募集程序在疾病髓系细胞中升高（Hedges g = 1.27，q = 3.86 × 10?4），而CD48–CD244显示出受体程序兼容性，CD48在21个NicheNet场景中的20个中位列两个受测配体之首。通讯效应在不同肺区室间存在差异，且两个焦点轴均非疾病相关的最强通路；SPP1–CD44显示出最大的固定效应全局变化，尽管其九配对证据在随机效应敏感性分析中不显著。结论：这些结果区分了系统性遗传易感性与局部通讯强度，并确定CD48–CD244和CX3CL1–CX3CR1为生物学上截然不同、依赖于区室的假说，将免疫细胞状态和募集与SAID-ILD中的纤维化肺环境联系起来。

自身免疫相关间质性肺病中系统性遗传易感性与肺组织重塑的整合研究解读

一、研究背景与科学问题

自身免疫相关间质性肺病（SAID-ILD）是一组与显著发病率和死亡率相关的免疫介导肺部疾病。尽管潜在的系统性疾病各异，其肺部表现常共享上皮损伤、免疫及髓系活化、内皮功能障碍、成纤维细胞活化和细胞外基质重塑等病理特征。这些变化并非孤立发生，而是上皮、免疫、内皮和基质细胞在肺微环境中相互作用，推动组织损伤向纤维化进展。循环炎症蛋白（CIPs）和血清代谢物（SMs）反映可能相关的系统性免疫和代谢变化，孟德尔随机化（MR）可利用遗传工具变量优先筛选与疾病易感性相关的循环因子，减少传统观察性研究的混杂和反向因果问题。然而，MR仅回答部分问题：它无法确定效应作用的解剖部位或疾病阶段，循环蛋白的遗传关联也不能直接证明特定肺配体-受体相互作用介导了该关联。批量转录组学可识别肺组织中可重复的疾病相关状态，单细胞和单核RNA测序可揭示涉及的细胞群体并推断潜在配体-受体相互作用。强局部重塑可与有限的系统性遗传支持并存，循环易感性信号在病变肺中可能仅显示温和或不一致的变化。因此，遗传证据可优先筛选潜在易感因子，批量转录组学可定义组织背景，单细胞分析可识别相关细胞群体和受体反应，整合这些层面时还需考虑肺区室间差异。本研究旨在整合双样本MR、批量转录组荟萃分析、外植肺单核RNA测序和支气管肺泡灌洗（BAL）单细胞RNA测序，将系统性易感性信号与细胞分辨的肺生物学联系起来。

二、研究内容与主要发现

研究人员建立了分层、产生假说的框架，将MR优先筛选的系统性易感性信号映射到可重复的上皮-髓系纤维化肺状态和区室依赖性免疫通讯上。CD48-CD244和CX3CL1-CX3CR1代表互补的免疫状态和募集假说，而SPP1-CD44在九个评估配对中显示出最大的固定效应组织兼容性变化，但随机效应敏感性分析存在不确定性。这些计算层面均不能证明受体特异性肺因果性，需要空间验证、离体实验和体内扰动研究。诊断或治疗效用尚未确立。

三、关键技术方法

研究整合了多层分析方法。双样本MR分析使用91种CIPs和1400种SMs的遗传工具变量，工具选择阈值p < 1.0 × 10^?5，蛋白分析额外使用全基因组显著变异（p < 5.0 × 10^?8），经LD聚类（r² < 0.001，10,000 kb）保留独立变异，排除次要等位基因频率≤0.01的变异，使用TwoSampleMR协调暴露和结局等位基因。批量转录组分析对五个公共数据集（GSE231693、GSE199152、GSE166036、GSE81292、GSE76808）进行疾病-对照对比，使用limma拟合，方向性Stouffer荟萃分析整合队列水平p值。细胞分辨分析处理了GSE303858/GSE303859 BAL数据集和GSE302151外植肺数据集，使用Seurat处理、DoubletFinder去双细胞，按供体和宽泛细胞类型聚合伪批量谱，TMMwsp归一化CPM，edgeR准似然模型进行差异表达。配体-受体兼容性以配体CPM×受体CPM计算，NicheNet用于受体程序分析，AutoDock Vina用于CX3CR1结构对接探索。

四、研究结果

3.1 遗传分析优先筛选SAID-ILD易感性中的CIP信号：91种CIPs中有10种与SAID-ILD易感性存在名义MR关联，但无一达到91蛋白BH阈值。名义信号包括ADA、IL24、TRANCE/TNFSF11和IL-13的保护性关联，以及CXCL10、CXCL9、CX3CL1、CD244、PD-L1和uPA的风险增加关联。CD244和CX3CL1工具集分别含33和30个SNP，平均F统计量为34.15和24.84。CXCL10是唯一达到严格阈值并满足91蛋白BH阈值的结果（p = 3.21 × 10^?4，q = 0.0292）。OR=1.5时CD244和CX3CL1的估计效能分别为23.5%和17.2%，80%效能可检测OR分别为2.51和3.09。代谢层面，9种SMs显示与SAID-ILD易感性一致的MR关联，包括γ-谷氨酰苏氨酸、3-羟基丁酰甘氨酸、X-12101、X-12714、X-25419、AMP-苯丙氨酸比值、α-生育酚-硫酸盐比值、苯丙氨酸-磷酸盐比值和α-酮戊二酸-脯氨酸比值。配对MR突出三个候选炎症-代谢联系，但中介MR不支持显著的CIPs-SMs-SAID-ILD链。

3.2 肺转录组揭示可重复的上皮-髓系纤维化状态：五数据集批量整合包含82,091条队列水平差异表达记录，覆盖28,997个基因。最复现的信号反映上皮损伤和细胞外基质重塑。POSTN和MMP7位列最强荟萃分析信号，其上调伴随KRT17、COL1A1、COL3A1和SPP1的复现性增加及AGER降低。CD244和CX3CL1在五数据集整合中下调，CD274/PD-L1在组织聚焦重算中达到FDR支持但贡献队列较少。CXCL10显示meta-p = 0.452，CXCR3显示meta-p = 0.0221。三队列RNA-seq肺分析中CD244降低（p = 1.73 × 10^?8），CX3CL1方向混合且名义阳性汇总（p = 0.0349）。代谢物相关基因程序方面，OXCT1和BDH1一致增加，HMGCS2降低；GOT1、GOT2和PSAT1增加，但SLC7A5和复合苯丙氨酸模块未获一致支持。

3.3 CD48–CD244和CX3CL1–CX3CR1的细胞分辨定位：GSE302151严格质控保留9563个核，去除51个预测双细胞后余9512个核，注释保留八种标签。CD244在四个数据集的T/NK细胞中表达最高，CD48分布于T/NK、B/浆细胞、髓系和肥大细胞区室。CX3CL1在肺组织数据集的内皮细胞和BAL的上皮细胞中表达最高，CX3CR1主要检测于T/NK和髓系细胞。跨数据集signed-Z分析识别CX3CL1在内皮（Z = 3.24，q = 0.00642）和髓系（Z = 2.56，q = 0.0277）区室的差异。CX3CL1–CX3CR1募集评分在髓系细胞（Hedges g = 1.27，q = 3.86 × 10^?4）和成纤维细胞/基质细胞（g = 0.90，q = 0.0167）中增加。纤维化-ECM评分在B/浆细胞、内皮、髓系和成纤维细胞/基质区室增加（q < 0.05）。

3.4 肺和BAL区室中的疾病相关配体-受体兼容性：BAL T/NK至T/NK CD48–CD244兼容性Hedges g = 0.41（p = 1.00）。GSE122960内皮至髓系CX3CL1–CX3CR1通路显示未达名义显著的正效应（g = 1.45，p = 0.0525）。无发送-接收通路在数据集内、轴内BH校正后存活。九个评估配对中，SPP1–CD44固定效应疾病-对照变化最大（g = 1.58，p = 6.84 × 10^?7，九配对q = 6.15 × 10^?6），CD48–CD244为?0.68（p = 0.0182），CX3CL1–CX3CR1为0.20（p = 0.470）。REML敏感性分析中SPP1–CD44估计为1.56（p = 0.0291，q = 0.160），CX3CL1–CX3CR1显示显著异质性（I² = 62.7%）。

3.5 NicheNet和受体状态分析识别疾病相关受体程序：CD48在21个CD48–CD244场景中的20个中按Pearson活性排名第一，CX3CL1在21个场景中的3个中排名第一。受体状态基因集富集分析在受体排序接收区室中识别趋化、干扰素/细胞因子反应、细胞毒性/抗原呈递、淋巴细胞或髓系活化和细胞粘附程序。BAL中CD244疾病-对照log2(CPM+1)差异为0.11（p = 0.876），CX3CR1髓系差异为1.32（p = 3.10 × 10^?4）。五模块留一队列外验证AUC为0.730（两焦点轴模块）、0.639（三组织背景模块）和0.742（组合五模块）。CMap/L1000反向签名筛选在150基因阈值下两焦点轴模块最大重叠覆盖率达1.00。CX3CL1–CX3CR1结构分析使用PDB 7XBX实验结构，含19个天然样界面接触对。

五、讨论与结论

研究区分了系统性遗传易感性与局部通讯强度。肺转录组发现的最一致变化反映上皮损伤、髓系活化和纤维化细胞外基质重塑，而多个优先CIPs显示较弱、多变或区室特异性转录变化，这具有生物学合理性：全组织表达受细胞组成和细胞状态强烈影响，而系统性遗传关联可能反映在已建立的纤维化肺组织中无法一致检测的易感性相关过程。CD48–CD244和CX3CL1–CX3CR1显示不同模式：CD244主要表达于T/NK细胞，CD48分布更广，CD48在NicheNet活性中排名第一但无细胞状态模块达到FDR显著性，更符合区室依赖性免疫相互作用和受体状态假说；CX3CL1主要表达于肺组织内皮和BAL上皮，CX3CR1集中于髓系和T/NK细胞，疾病相关变化和FDR支持的募集程序支持其在免疫细胞定位和募集中的潜在作用。SPP1–CD44虽显示最大固定效应变化但随机效应敏感性不显著，且局部表达兼容性大变化不建立遗传支持的因果通路。研究局限性包括MR仅优先筛选循环易感性信号而不直接建立肺受体通路效应、公共数据集异质性、严格质控后部分局部分析仅剩1-2个疾病供体、配体-受体表达产物测量转录兼容性而非分泌、物理接触、受体占据或信号传导，以及对接计算仅为计算推测。研究结论为：该研究建立了分层、产生假说的框架，将MR优先筛选的系统性易感性信号映射到可重复的上皮-髓系纤维化肺状态和区室依赖性免疫通讯上。CD48-CD244和CX3CL1-CX3CR1代表互补的免疫状态和募集假说，SPP1-CD44在九个评估配对中显示最大固定效应组织兼容性变化但随机效应敏感性存在不确定性。这些计算层面均不能证明受体特异性肺因果性，需要空间验证、离体实验和体内扰动研究，诊断或治疗效用尚未确立。