产品介绍
双样本孟德尔随机化(Two-Sample Mendelian Randomization)是孟德尔随机化的一种常用设计类型,核心含义是分别从两个独立的人群样本中获取 “暴露相关遗传变异” 和 “结局相关遗传变异” 的汇总数据,再通过统计方法关联两者,推断暴露与结局的因果关系。
应用领域
1、复杂疾病病因学研究
2、生物标志物的临床转化验证
3、生活方式与健康结局的因果推断
4、药物靶点的优先级排序与安全性评估
5、跨种族 / 人群的因果关联验证
技术路线
分析内容
基础因果效应分析
核心模型结果解读:优先关注逆方差加权法(IVW) 的效应量(β 值或 OR 值)、95% 置信区间(CI)及 P 值,判断暴露与结局是否存在统计学显著的因果关联。若置信区间不包含 1(OR 值)或 0(β 值)且 P<0.05,则提示存在显著因果效应。
补充模型结果对比:结合MR-Egger 回归、加权中位数法、加权模式法的结果,验证因果效应的一致性。若不同模型结果方向一致且均显著,说明因果关联稳健;若结果存在差异,需进一步探究异质性或多效性来源。
工具变量有效性检验
关联性检验:计算每个工具变量(SNP)的F 统计量,要求 F>10,排除弱工具变量偏倚。若存在多个弱工具变量,需考虑剔除或采用对弱工具变量更稳健的模型(如加权中位数法)。
工具变量数量评估:确保纳入的 SNP 数量充足(一般建议≥3),避免因工具变量过少导致结果不稳定。
异质性分析
整体异质性检验:通过Cochran's Q 检验评估工具变量间的效应异质性,若 P<0.05,提示不同 SNP 对结局的效应存在差异,可能存在混杂或水平多效性。
异质性来源定位:可通过森林图直观观察单个 SNP 的效应方向和大小,识别可能的异常 SNP;也可采用留一法逐一剔除 SNP,判断异质性是否由个别异常位点导致。
水平多效性检验
MR-Egger 截距检验:若截距项的 P<0.05,提示存在显著的定向水平多效性,此时 IVW 结果可能存在偏倚,需参考 MR-Egger 回归校正后的效应量。
多效性残差和异常值检验(MR-PRESSO):用于检测并校正异常值 SNP 导致的水平多效性,进一步提升结果可靠性。
敏感性分析结果验证
留一法敏感性分析:查看剔除单个 SNP 后剩余 SNP 的因果效应变化,若效应量和显著性未发生明显波动,说明结果不受个别 SNP 影响。
分层 / 亚组分析:若数据允许,按人群特征(如年龄、性别、种族)分层分析,验证因果效应在不同亚组中的一致性。
样本重叠检验:若暴露和结局的 GWAS 样本存在潜在重叠,需采用MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO) 或校正样本重叠的方法,排除样本重叠导致的偏倚。
结果可视化分析
通过森林图展示单个 SNP 的效应量及合并效应,直观呈现因果关联强度;
通过漏斗图判断是否存在发表偏倚或小样本效应,若漏斗图对称,提示无明显偏倚;
通过散点图展示暴露与结局的遗传效应关联,辅助判断因果效应方向。
因果关联合理性讨论
结合生物学机制,解释暴露影响结局的潜在通路,判断因果关联的生物学合理性;
与传统观察性研究、临床试验结果对比,探讨结果的一致性,若存在差异,需分析研究设计或人群特征的影响。
研究结果(部分)