基因预测的血浆脂质代谢物与免疫细胞表型在慢性阻塞性肺疾病中的关联：一项双样本孟德尔随机化研究及探索性反向与中介分析

《International Journal of Chronic Obstructive Pulmonary Disease》：Genetically Predicted Plasma Lipid Metabolites and Immune-Cell Phenotypes in Chronic Obstructive Pulmonary Disease: A Two-Sample Mendelian Randomization Study with Exploratory Reverse and Mediation Analyses

　　

目的：评估基因预测的血浆脂质代谢物与免疫细胞表型同慢性阻塞性肺疾病（COPD）的关联，并评估候选的免疫介导的脂质-COPD通路。方法：研究人员开展了双样本孟德尔随机化（MR）分析，使用全基因组关联研究（GWAS）汇总统计数据，包括458,440名欧洲血统个体的COPD数据（26,274例病例和432,166名对照）、7,174名芬兰参与者的179种血浆脂质代谢物，以及3,757名撒丁岛参与者的731种免疫细胞表型。在P<1×10^?5水平选择的遗传工具经连锁不平衡（LD）聚集、协调并筛查工具强度。逆方差加权（IVW）为主要估计方法，稳健调整轮廓评分（MR-RAPS）用于敏感性分析。Benjamini–Hochberg（BH）校正分别在完整的脂质和免疫细胞检验族中应用。对13种选定脂质的反向MR、使用东亚COPD结局的跨 ancestry 分析，以及18条先前选定的候选通路的两步MR为探索性分析。估计了间接效应，并分别评估了带条件F统计量的补充多变量MR；BH校正覆盖反向和间接效应族。结果：21种脂质代谢物和48种免疫细胞表型显示与COPD的名义IVW关联（P<0.05），然而，无一达到完整族阈值q<0.05。六种免疫细胞表型达到探索性（q<0.10）阈值。MR-RAPS估计与欧洲IVW估计总体方向一致，但未发现任何完整族（q<0.05）关联。跨 ancestry 结局分析未一致支持聚焦的欧洲信号。18个间接效应均未通过BH校正。对于从磷脂酰胆碱（16:0_20:4）通过CD62L^?髓系树突状细胞绝对计数到COPD的通路，间接效应为0.0037（95% CI ?0.0008至0.0083；P=0.110；BH q=0.138）。两个反向MR关联名义显著，但均未通过13项检验的校正（两者BH q=0.140）。尽管多变量MR在选定候选集内识别出若干直接效应模式，但这些发现未确立中介作用。结论：主要欧洲分析对所评估的血浆脂质代谢物或免疫细胞表型与COPD的多重校正关联提供了有限证据，且无候选间接效应获得统计学支持。因此，反向MR信号应被视为假设生成性，需在独立队列和实验研究中验证。现有发现不支持临床风险预测或诊断应用。

慢性阻塞性肺疾病（COPD）是一种以持续性气流受限、小气道疾病和肺实质破坏为特征的异质性慢性呼吸系统疾病，仍是全球发病和死亡的主要原因。一项系统综述和建模分析估计，2019年全球30–79岁人群中约有3.919亿COPD患者，在GOLD固定比率定义下患病率为10.3%（95% CI 8.2–12.8），其中约80.5%的病例位于低收入和中等收入国家。人口老龄化、烟草暴露和环境空气污染预计将维持这一负担。尽管当前治疗可缓解症状并减少急性加重，但不能完全逆转已建立的结构性肺损伤。临床和公共卫生负担使得区分疾病相关生物标志物与可能影响COPD易感性的暴露因素变得重要。脂质参与肺表面活性物质组成、膜结构和炎症反应。观察性研究报告了COPD个体循环脂质谱和肺表面活性物质磷脂的改变。既往以脂质为重点的MR研究也探讨了脂质特征与COPD易感性之间的关联。然而，观察性脂质谱可能受吸烟、疾病严重程度、合并症或治疗的影响，而MR发现取决于工具变量有效性、统计功效和多重检验的适当控制。此外，单个分子脂质种类可能具有不同的生物学特性，不应假定其在更广泛的脂质类别内表现一致。COPD还以复杂的先天性和适应性免疫失调为特征。巨噬细胞、树突状细胞和淋巴细胞群体参与持续性炎症、组织损伤和急性加重。临床研究报告了调节性T细胞、B细胞活化、CD8⁺ T细胞和黏膜免疫群体的变化，以及循环免疫细胞和基因表达模式与COPD急性加重的关联。MR研究进一步评估了免疫细胞表型与COPD的关系。然而，免疫表型高度相关且具有情境依赖性，观察性关联无法确定这些细胞特征是先于COPD、跟随疾病活动还是反映共享的上游生物学。孟德尔随机化（MR）使用遗传变异作为工具变量，在工具相关性、独立于混杂因素以及工具到结局不存在替代通路等假设下评估暴露-结局关系。Cao等人使用中介MR框架检查了731种免疫细胞表型和1,400种血浆代谢特征，重点关注循环代谢物是否介导免疫细胞特征与COPD之间的关联。他们的研究为研究COPD中的代谢-免疫关系提供了重要的遗传流行病学框架。然而，互补方向——基因预测的分子脂质种类是否可能通过免疫细胞表型与COPD相关——仍较少被表征。血浆代谢组学面板涵盖化学上异质的化合物，包括氨基酸、类固醇、胆汁酸、外源性物质、代谢物比率和未完全表征的特征。以脂质为重点的分析可能提供关于分子磷脂酰胆碱、磷脂酰乙醇胺和三酰甘油种类的更具体信息，这些种类可能与膜组成、肺表面活性物质和炎症信号传导相关。此外，高维MR发现的解释需要适当控制多重检验、评估弱工具不确定性，并明确量化间接效应的不确定性。当对GWAS样本量相对适中的分子特征使用宽松的工具选择阈值时，这些考虑尤为重要。在既往工作基础上，研究人员使用基于GWAS汇总数据的双样本MR框架评估了179种分子表征的血浆脂质种类、731种免疫细胞表型和COPD。研究人员在完整族Benjamini–Hochberg（BH）校正下使用逆方差加权（IVW）评估脂质-COPD和免疫细胞-COPD关联，并使用稳健调整轮廓评分方法评估其稳健性。选定的反向关联作为探索性分析进行检查，并使用东亚COPD结局评估跨 ancestry 可迁移性。候选脂质-免疫-COPD通路使用两步MR和带条件工具强度诊断的多变量MR进一步研究。因此，本研究设计为既往代谢-免疫MR研究的以脂质为重点的扩展，在保守证据框架下解决互补的生物学方向，而非复制相同的中介模型。研究目标是确定在考虑多重检验、弱工具不确定性和间接效应估计不确定性后，种类特异性脂质关联和候选免疫介导通路是否仍获支持。研究人员开展了双向双样本MR设计，基于公开可用的GWAS汇总统计数据。主要分析分别评估179种血浆脂质代谢物和731种免疫细胞表型与COPD的关系。IVW作为主要估计方法，稳健调整轮廓评分MR（MR-RAPS）用于评估估计对弱工具不确定性的稳健性。反向MR、跨 ancestry 结局分析以及18条先前识别的脂质-免疫-COPD通路的评估被视为探索性。研究方案和分析计划未预注册。主要欧洲分析对所评估的血浆脂质种类或免疫细胞表型与COPD的稳健关联提供了有限证据。名义脂质关联和探索性免疫细胞信号识别了进一步研究的候选者，但所评估的脂质-免疫-COPD间接效应均未获统计学支持。这些发现通过区分生物学上合理的候选信号与多重校正的中介证据，完善了现有证据。需要更大的暴露GWAS、独立的 ancestry 匹配复制和实验研究来确定这些候选者中是否有任何参与COPD发病机制。研究人员为开展研究用到的主要关键技术方法包括：基于公开GWAS汇总数据的双样本孟德尔随机化设计；逆方差加权（IVW）作为主要估计方法；稳健调整轮廓评分（MR-RAPS）用于弱工具不确定性敏感性分析；Benjamini–Hochberg（BH）多重检验校正分别在完整脂质和免疫细胞检验族中应用；连锁不平衡（LD）聚集和协调进行工具选择；两步MR估计间接效应；多变量MR（MVMR）评估互调直接效应并计算条件F统计量；探索性反向MR；使用东亚COPD GWAS的跨 ancestry 可迁移性分析。样本队列来源包括：欧洲血统COPD GWAS（458,440人，26,274例病例和432,166名对照）、芬兰GeneRISK队列脂质GWAS（7,174人）、撒丁岛免疫细胞表型GWAS（3,757人）、东亚BioBank Japan COPD GWAS（4,017例病例和162,653名对照）。研究结果方面：血浆脂质代谢物与COPD的关联：在所有179种脂质代谢物中，21种IVW关联名义P<0.05（6种正向和15种反向），但无一达到BH q<0.05（最小q=0.114）。所有15种反向名义关联涉及磷脂酰胆碱（PC；12种）或磷脂酰乙醇胺（PE；3种）。最强的名义信号为PC（O-16:0_20:4）（OR 1.057，95% CI 1.024–1.091；P<0.001；q=0.114）。这些为探索性关联，非FDR支持的发现。免疫细胞表型与COPD的关联：在731种免疫细胞表型中，48种具有名义IVW P<0.05，无一在完整族BH校正后达到q<0.05。六种表型IVW q值低于0.10，包括调节性T细胞（Treg）、自然杀伤（NK）细胞和效应记忆（EM）T细胞特征：CD4 on CD39⁺ activated Treg（q=0.081）、Naive CD4⁺ %CD4⁺（q=0.081）、HLA-DR⁺ NK AC（q=0.081）、Activated and resting Treg % CD4 Treg（q=0.081）、Secreting Treg % CD4 Treg（q=0.081）和CD3 on EM CD8br（q=0.097）。这些q<0.10信号作为探索性排名描述性报告，而非确定性发现。免疫细胞表型对脂质-COPD关联的探索性中介：在18条候选脂质-免疫-COPD通路中，17条间接效应95% CI包含零。剩余通路PC（18:1_18:2）→ Naive CD4⁺ %CD4⁺ → COPD的间接估计为0.0042（95% CI 0.0001–0.0084；P=0.046），但未通过BH校正（q=0.138）。18个间接效应均未达到q<0.05。对于PC（16:0_20:4）→ CD62L^?髓系DC AC → COPD，间接估计为0.0037（95% CI ?0.0008至0.0083；P=0.110；q=0.138），无统计学支持的中介证据。补充MVMR分析评估相同18条候选通路，每条通路保留28–44个独立SNP。18条通路中14条标准条件F统计量≥10，相同14条通路使用异质性稳健条件F统计量达到该基准。四条通路至少一个条件F<10，因此被认为易受条件弱工具偏倚影响。对于PC（16:0_20:4）和CD62L^?髓系DC AC，条件F统计量分别为49.7和31.1。PC（16:0_20:4）与COPD的互调直接关联为OR 1.031（95% CI 0.996–1.066；P=0.089），相应免疫直接关联为OR 1.028（95% CI 0.996–1.060；P=0.098），均无统计学显著性。探索性反向MR分析：基因预测的COPD与PC（18:2_20:1）和PC（O-16:0_20:4）显示名义关联（P=0.019和P=0.021），但均未通过13项选定检验的BH校正（两者q=0.140）。稳健性与跨 ancestry 分析总结：MR-RAPS成功估计所有179种脂质和731种免疫细胞暴露。在所有910种暴露中，IVW和MR-RAPS效应估计强相关（Pearson r=0.87），792/910保留相同效应方向。东亚COPD敏感性分析产生900/910种暴露的可评估IVW估计，保留原始脂质工具的中位52.9%和原始免疫工具的45.0%。与欧洲IVW估计的方向一致性在所有可评估暴露中为494/900，在结果定义的聚焦集中为12/27。在欧洲主要分析中具有名义IVW证据的21种脂质特征中，仅两种与东亚COPD结局名义相关且均显示不一致效应方向；六种欧洲IVW q<0.10的免疫表型均未显示东亚名义IVW证据。讨论部分总结：本研究使用GWAS汇总水平数据评估了179种血浆脂质代谢物和731种免疫细胞表型与COPD的关系。在主要欧洲分析中，21种脂质和48种免疫细胞关联达到名义显著性，但无一在各自完整暴露族内通过Benjamini–Hochberg校正。六种免疫细胞表型达到探索性（q<0.10）阈值，而无脂质种类达到该水平。MR-RAPS对若干聚焦信号产生广泛一致的效应方向，但未识别任何满足完整族（q<0.05）的关联。此外，18个候选间接效应均未通过多重检验校正。总之，这些发现识别了一组有限的假设生成性脂质和免疫细胞特征，但未提供免疫介导的脂质-COPD通路的统计学支持证据。21种名义脂质关联包括6种正向和15种反向名义关联；所有15种反向关联涉及PC（12种）或PE（3种）。PC（O-16:0_20:4）具有最强的正向信号（OR 1.057，95% CI 1.024–1.091；q=0.114）。反向示例包括PC（16:1_18:1）（OR 0.944，95% CI 0.906–0.983；q=0.170）和PE（18:1_18:1）（OR 0.951，95% CI 0.918–0.985；q=0.170）。然而，这些关联均未通过179项脂质检验的完整族校正。单个脂质种类间的不同方向表明，潜在相关的脂质效应不能仅从广泛脂质类别标签推断。脂肪酸组成、链长、饱和度和醚键的差异可能影响膜行为和脂质信号传导。尽管如此，现有结果不支持将PC或PE种类整体归类为有害或保护性。脂质面板还包含许多PC种类，增加了名义关联在该类别内聚集的机会。现有脂质组学研究为检查单个分子种类提供了生物学背景。Agudelo等人报告COPD患者支气管肺泡灌洗液中表面活性物质脂质减少，以及磷脂可用性与肺功能之间的相关性。在FINRISK中，Shoghli等人发现循环PC（14:0_22:6）与 incident COPD之间的反向关联。这些观察支持继续研究分子脂质谱，但不直接验证现有信号。脂质种类、生物隔室、暴露尺度和分析设计在各研究间存在差异。此外，基因预测的循环脂质水平无法识别关联背后的肺部过程，如表面活性物质重塑、膜破坏或炎症信号传导。观察到的效应估计通常不大。其比值比不能直接转化为COPD风险的临床有意义变化，因为脂质GWAS效应尺度不对应标准化临床浓度增量。绝对风险差异也未估计。因此，名义关联应被视为候选优先级信号，而非临床重要脂质效应的证据。免疫细胞发现同样需要表型特异性解释。六种表型达到探索性（q<0.10）阈值，包括与na?ve CD4⁺细胞、调节性T细胞、NK细胞和效应记忆CD8⁺细胞相关的特征。然而，无一达到主要（q<0.05）阈值。若干领先信号涉及Treg相关特征。既往研究描述了COPD中Treg丰度和功能的改变，尽管这些改变的方向和幅度因疾病阶段、吸烟状态、解剖隔室和表型定义而异。这种异质性相关，因为绝对细胞计数、相对比例和标志物表达强度代表不同的生物学特征，不应互换解释。现有数据不支持Treg活性在COPD中统一保护或有害的一般结论。HLA-DR⁺ NK绝对计数显示反向估计，探索性q值为0.081。Guan等人先前报告HLA-DR⁺ NK AC与COPD之间的反向MR关联。方向的一致性支持进一步评估该特定表型，但不确立NK细胞作为类别的保护作用。该关联未达到主要阈值，相应MR-RAPS结果在完整族校正后未获支持。IVW和MR-RAPS估计在方向和幅度上广泛一致。在所有910种暴露中，估计强相关，聚焦展示集中所有27种特征保留相同效应方向。21种名义脂质信号中12种在MR-RAPS下仍名义显著。这种一致性表明若干欧洲估计的方向对MR-RAPS实施的弱工具调整不高度敏感。然而，IVW和MR-RAPS之间的一致性不构成独立复制。两种估计使用相同暴露和结局GWAS以及 largely 相同工具。此外，MR-RAPS未识别任何完整族（q<0.05）关联。四种脂质种类达到探索性MR-RAPS（q=0.093），但这些发现未达到主要阈值，也未推翻IVW分析的证据分类。因此，MR-RAPS支持选定点估计的稳定性，而未将名义关联转化为主要发现。东亚COPD结局分析未一致支持聚焦的欧洲信号。27种聚焦特征中仅12种在欧洲和东亚分析间显示一致效应方向。在21种名义欧洲脂质信号中，两种与东亚结局显示名义关联，但均与欧洲估计方向不一致。六种探索性欧洲免疫细胞信号均未显示东亚名义IVW证据。在完整东亚筛选中，两种脂质特征和一种免疫细胞表型达到IVW（q<0.05），但在东亚MR-RAPS下无一保持低于该阈值。这些东亚-only发现可能反映 ancestry 特异性遗传结构、工具的不完全可迁移性或偶然，而非欧洲关联的复制。欧洲工具保留的中位比例脂质特征仅为52.9%，免疫细胞表型为45.0%。因此，等位基因频率、连锁不平衡结构和SNP-暴露效应的差异可能影响了估计。由于SNP-暴露关联来自欧洲血统暴露GWAS，该分析不是正式的 ancestry 匹配复制。它应被解释为评估欧洲衍生工具是否在东亚COPD结局中产生兼容估计。反向MR产生基因预测COPD易感性与PC（18:2_20:1）和PC（O-16:0_20:4）的名义关联。均未通过13项选定反向检验的BH校正。因此，发现不提供COPD遗传易感性影响这些脂质特征的统计学支持证据。反向分析限于选定候选集，在工具选择、工具强度和统计功效方面与正向分析不同。因此，未能获得校正证据不应解释为COPD对脂质代谢无影响的证据。相反，现有分析未在评估的候选集内检测到稳健的反向关联。18个候选间接效应均未通过多重检验校正。仅从PC（18:1_18:2）通过na?ve CD4⁺细胞占CD4⁺细胞百分比到COPD的通路显示名义间接效应证据。其校正q值为0.138，排除了统计学支持的中介结论。CD62L^?髓系DC AC作为PC（16:0_20:4）与COPD之间的候选中间体特别受关注。髓系细胞参与COPD炎症提供了生物学背景，但估计的间接效应为0.0037（95% CI ?0.0008至0.0083）。此外，两步计算使用单变量中介-结局系数。补充MVMR直接效应估计未提供暴露调整的间接效应。因此，两种分析处理相关但不同的量。当间接效应未获支持时，MVMR直接效应的存在或生物学合理性均不能确立中介作用。若干脂质和免疫细胞直接效应系数在MVMR分析中达到候选集调整阈值。这些结果表明选定特征在受限的18通路集内保留了与COPD的互调关联。然而，这些通路通过原始数据依赖工作流选择，候选集内校正未考虑先前多阶段选择过程。四个模型还包含至少一个条件F统计量低于10的暴露，一条通路在校正后保留残余异质性证据。因此，这些MVMR发现应被视为通路水平敏感性证据，而非独立发现。中介比例未解释，因为无间接效应获多重校正支持。当总效应和间接效应小或估计不精确时，其比率可能不稳定、为负或大于一。报告此类比例可能暗示数据不支持的机制确定性程度。Cao等人使用相同731表型免疫GWAS以及两步和多变量MR检查了免疫-代谢与COPD的关联。他们的1,400特征代谢面板和免疫-代谢物-COPD重点与本研究179种GeneRISK脂质种类和脂质-免疫-COPD候选不同。他们将免疫细胞表型作为暴露，循环代谢物作为候选中介，而研究人员评估互补方向，以免疫细胞表型作为候选中间体。研究使用不同COPD GWAS登录号（分别为GCST90018807和GCST90473707），但不能假定这些结局资源之间的参与者独立性。他们报告的CD62L^?髓系DC AC–1-棕榈酰-GPG（16:0）通路与本研究PC（16:0_20:4）候选不同。他们的代谢面板还包括化学异质化合物，而本研究聚焦于定义明确的脂质种类。本研究通过检查互补生物学方向并应用完整族多重校正、弱工具敏感性分析和条件工具强度诊断，扩展了相关文献。其主要贡献是界定所评估候选的有限统计学支持，而非确立额外脂质-免疫机制。本研究使用明确证据等级系统评估分子脂质种类和免疫细胞表型。完整族BH校正分别应用于所有179项脂质和731项免疫细胞检验，从而避免仅对名义选定结果进行校正。MR-RAPS应用于完整暴露族而非仅选定阳性发现。分析还区分间接效应与互调直接效应，评估MVMR中的条件工具强度，分别校正候选间接效应和反向MR族，并明确将跨 ancestry 分析视为可迁移性评估而非独立复制。主要、探索性和名义证据的透明分类是额外优势。该框架降低了相关估计间一致性被误解为关联确认的风险。应考虑若干统计学局限性。第一，脂质和免疫细胞GWAS样本量相对适中，分别为7,174和3,757名参与者。因此，工具在P<1×10^?5而非传统全基因组阈值选择。尽管所有保留的单变量工具F统计量高于19.5且使用MR-RAPS评估弱工具不确定性，但不能排除赢家诅咒和残余弱工具偏倚。对许多暴露而言，限制于全基因组显著工具的敏感性分析不可行，因为可用变异太少。第二，MR-RAPS可在其模型下容纳SNP-暴露估计不确定性和平衡多效性效应，但通常不能解决相关水平多效性。未实施专门建模相关多效性的方法，如CAUSE。此外，非显著MR-Egger截距或MR-PRESSO检验不能证明水平多效性不存在。第三，18条候选通路源于数据依赖的多阶段选择过程。对这些通路的BH校正仅控制选定集内的多重性，未提供所有可能脂质-免疫组合的全局误差控制。此外，修订两步分析中使用的脂质-免疫成分估计从原始成分结果恢复，而非在新的脂质-免疫筛选中独立重新估计。因此，中介发现应被视为历史候选的重新评估，而非独立通路发现分析。第四，四个MVMR模型至少一个条件F统计量低于10，表明可能存在条件弱工具偏倚。一条通路在校正后还检测到残余异质性。条件F统计量和多效性诊断改善解释但不保证工具有效性。第五，主要分析限于欧洲血统个体。跨 ancestry 分析使用欧洲SNP-暴露估计与东亚COPD结局，工具覆盖不完整。因此不能提供 ancestry 匹配复制。暴露与主要结局GWAS之间的潜在参与者重叠无法量化，主要结局数据集的COPD病例定义和原始协变量调整的完整信息不可用。最后，MR估计反映终生基因预测差异，不直接捕捉与COPD阶段、急性加重、戒烟或治疗相关的代谢或免疫变化。诊断区分度、校准和增量预测价值未评估，现有结果不支持临床筛查、诊断或治疗靶向。无候选通路经实验验证。未达到多重校正阈值不排除小的生物学效应，但在推断此类效应前需要独立且充分把握度的研究。结论：主要多重控制分析对所评估的血浆脂质种类或免疫细胞表型与COPD的稳健关联提供了有限证据。名义脂质关联和探索性免疫细胞信号识别了进一步研究的候选者，但所评估的脂质-免疫-COPD间接效应均未获统计学支持。这些发现通过区分生物学上合理的候选信号与多重校正的中介证据，完善了现有证据。需要更大的暴露GWAS、独立的 ancestry 匹配复制和实验研究来确定这些候选者中是否有任何参与COPD发病机制。