关键词:
淋巴瘤,非霍奇金
乙型肝炎病毒
孟德尔随机化分析
因果关系
摘要:
目的采用双样本双向孟德尔随机化(MR)分析探讨非霍奇金淋巴瘤(NHL)和慢性乙型肝炎病毒(HBV)感染的相关性。方法NHL的遗传变异数据来自芬兰数据库(FinnGen)联盟2021年公开的全基因组关联分析(GWAS)数据集,包括1088例NHL患者和299952例对照受试者。慢性HBV感染的GWAS数据集来自2021年发表的GWAS分析,包括145例NHL患者和351740例对照受试者。先将NHL作为暴露因素,与NHL显著相关的单核苷酸多态性(SNP)作为工具变量,慢性HBV感染作为结局变量,采用逆方差加权(IVW)法进行双样本MR分析。再将慢性HBV感染作为暴露因素,与慢性HBV感染显著相关的SNP作为工具变量,NHL作为结局变量,进行反向双样本MR分析。IVW法利用各工具变量的方差倒数作为权重进行拟合,对SNP进行比值法逐一测算并加权回归,从而获得总体估计值。利用MR⁃Egger回归法、加权中位数(WME)法作为IVW法的补充。在敏感性分析中,采用留一法敏感性分析评估单个SNP的影响;通过Cochran Q检验分析所选工具变量的异质性;利用MR⁃Egger回归法衡量工具变量平均水平多效性的大小,计算出方向性P值;使用MR⁃多效性残差和与异质性(MR⁃PRESSO)Global Test排除可能的水平多效性离群值,减少偏倚。结果在留一法敏感性分析中剔除对因果关联影响较大的SNP。在正向MR分析中,工具变量为10个与NHL相关的SNP;IVW法表明NHL与慢性HBV感染不存在因果关联(OR=0.979,95%CI:0.925~1.036,P=0.465);MR⁃Egger回归法(OR=0.992,95%CI:0.926~1.062,P=0.825)和WME法(OR=0.992,95%CI:0.934~1.055,P=0.805)作为IVW法的补充,获得了一致的结果。在敏感性分析中,Cochran Q检验结果显示工具变量间不存在异质性(IVW法P=0.271,MR⁃Egger回归法P=0.239);在MR⁃Egger回归法(截距为-0.01,P=0.778)和MR⁃PRESSO Global Tes(t P>0.05)中未发现水平多效性,提示结果稳健。在反向MR分析中,工具变量为8个与慢性HBV感染相关的SNP。IVW法(OR=1.117,95%CI:0.942~1.324,P=0.202)亦未发现慢性HBV感染与NHL存在显著的因果关系,MR⁃Egger回归法(OR=0.777,95%CI:0.450~1.343,P=0.401)与WME法(OR=1.120,95%CI:0.887~1.415,P=0.351)也获得了类似的风险估计。敏感性分析同样提示结果具有一致性和可靠性,Cochran Q检验显示工具变量间不存在异质性(IVW法P=0.775,MR⁃Egger回归法P=0.903);MR⁃Egger回归法(截距为0.102,P=0.548)和MR⁃PRESSO Global Tes(t P>0.05)未发现水平多效性。结论MR分析表明NHL与慢性HBV感染之间不存在因果关系。