两样本孟德尔随机化法分析抑郁症与间质性肺病关系的研究

目的 采用孟德尔随机化方法分析抑郁症与间质性肺病(ILD)之间的因果关系，为ILD的预防、治疗及预后评估提供新的思路。方法 利用大样本全基因关联研究汇总数据，选择与抑郁症密切相关的遗传位点作为工具变量，逆方差加权法为主，加权中位数法和MR-Egger回归为辅对数据进行两样本孟德尔随机化分析。以OR值评价抑郁症与ILD之间的因果关系，异质性检验、基因多效性检验和敏感性分析3种方法评估结果的稳定性和可靠性。结果 共纳入37个单核苷酸多态性(SNP)位点作为工具变量，逆方差加权法估计抑郁症患者患ILD的风险(OR)是健康人群的1.21倍(95%CI:1.075～1.361,P=0.002)。加权中位数法同样支持抑郁症和ILD之间存在因果效应(95%CI:1.118～1.564,P=0.001)。逆方差加权法和MR-Egger回归的异质性检验结果表明不存在异质性，MR-Egger回归截距项和MR-PRESSO方法检验表明结果受基因多效性影响的可能性较小，leave-one-out分析未发现非特异性SNP。结论 抑郁症与ILD之间可能存在正向因果关联。     

两样本孟德尔随机化研究高血压与焦虑发生风险的因果关系

目的 采用两样本孟德尔随机化研究方法评估高血压与焦虑发生风险之间潜在的因果关系。方法 研究数据来自全基因组关联研究(genome-wide association study, GWAS)汇总数据集,其中焦虑数据包含34 121例样本,高血压数据包含463 010例样本,选取71个与高血压密切相关的遗传位点作为工具变量。分别使用随机效应逆方差加权、加权中位数和Mr-Egger等方法进行两样本孟德尔随机化分析,以OR值评估高血压与焦虑之间的因果关系。同时采用Mr-PRESSO和Mr-Egger检测水平多效性,“leave-one-out”进行敏感性分析。结果 逆方差加权法分析结果显示高血压与焦虑发生之间呈显著正相关,比值比(odds ratio, OR)为1.276 ［95%可信区间(confidence interval,CI):1.074–1.515;P=0.005］;加权中位数回归结果显示了相似的因果关系 ［OR, 1.305, 95%CI:1.016–1.676),P=0.037］。此外,漏斗图和Mr-Egger(P=0.27)截距均没有观察到水平多效性和异质性。结论 研究表明,高血压与焦虑发生之间可能存在正向因果关系。     

孟德尔随机化研究揭示饮茶与胆石症之间的因果关系

目的 采用双向孟德尔随机化方法评估饮茶与胆石症之间的因果关系。方法 饮茶的遗传工具是从UK Biobank全基因组关联研究（genome-wide association study,GWAS）中447 485名参与者确定的,其中共采用与饮茶相关的32个单核苷酸多态性（single nucleotide polymorphism,SNP）进行两样本孟德尔随机化分析以评估饮茶与胆石症之间的因果关系。胆石症的数据来自Finngen联盟公开提供的35 712名胆石症患者和273 442名对照人群的GWAS数据。本研究以inverse-variance weighted(IVW)为主要方法估计因果关系,并进行敏感性分析以确保结果的稳健性。最后,采用反向孟德尔随机化分析验证是否存在反向关联。结果 饮茶与胆石症存在因果关系。在IVW方法中,饮茶每增加一个标准差（2.85杯/d）,胆石症的发生风险减低26.1%(OR=0.739, 95%CI:0.536～0.990)。在weighted mode中也得到了相似的结果（OR=0.652, 95%CI:0.458～0.928）,但没有证据证明存在反向关...     

两样本孟德尔随机化分析肠道微生物与精神分裂症的因果关系

目的 利用两样本孟德尔随机化方法探究271种肠道微生物与精神分裂症的因果关系。方法 采用公共数据库中FINRISK 2002研究的肠道微生物GWAS数据(n=5 959)和精神疾病全基因组研究联盟的精神分裂症GWAS数据(n=35 802),根据预设的阈值提取与271种肠道微生物相对丰度显著相关的独立遗传位点作为工具变量。主要采用逆方差加权法进行分析,根据效应指标优势比(OR)和95%CI评估结果。使用留一法、异质性检验、水平基因多效性检验、MR Steiger方向性检验来验证结果的稳定性和可靠性。结果 放线菌纲(Actinobacteria)丰度的升高,可以增加精神分裂症的患病风险(OR=1.24,95%CI:1.10～1.40)。留一法分析显示结果稳定,不存在对结果有强影响的工具变量,该因果关系方向为正向,且可以剔除异质性和水平基因多效性对因果效应估计产生的影响。结论 放线菌纲(Actinobacteria)会增加精神分裂症的发生风险。     

基于孟德尔随机化方法探究可干预危险因素与支气管哮喘之间的因果关系

目的 运用两样本孟德尔随机化（Mendelian randomization,MR）方法探究生活方式、饮食、睡眠及社会经济地位等因素与支气管哮喘的因果关系。方法 研究主要使用逆方差加权法(IVW)进行两样本MR分析,同时采用留一法、异质性检验、水平基因多效性检验来验证结果的稳定性和可靠性。结果 IVW模型结果显示吸烟（OR=1.17,95%CI:1.08～1.28）、饮酒（OR=1.22,95%CI:1.05～1.41）、久坐（OR=1.78,95%CI:1.38～2.31）、睡眠时长（OR=2.29,95%CI:1.38～2.31）、体重指数（OR=1.41,95%CI:1.29～1.53）与支气管哮喘患病风险增加有关。新鲜水果摄入（OR=0.50,95%CI:0.34～0.73）、干果摄入（OR=0.54,95%CI:0.34～0.73）、谷类食物摄入（OR=0.62,95%CI:0.41～0.94）、教育程度（OR=0.57,95%CI:0.34～0.73）、社会活动（OR=0.26,95%CI:0.34～0.73）会降低支气管哮喘的发病风险。ME-Egger截距检测结果不存在多...     

酒精依赖和肌萎缩侧索硬化因果关系的孟德尔随机化分析

目的采用孟德尔随机化方法探索饮酒行为和肌萎缩侧索硬化发病风险的因果关系。方法利用大样本全基因组关联研究数据,确定与酒精依赖密切相关的遗传位点作为工具变量,通过逆方差加权法、MR-Egger回归法等两样本孟德尔随机化方法,以比值比(OR值)作为结局指标对因果关系进行评价。结果本研究共纳入476例男性酒精依赖患者和1358例男性对照,IVW方法估计的OR值为1.12(95%可信区间:1.02-1.24,P=0.021),MR-Egger回归结果表明遗传多效性不会对结果造成偏倚(截距=-0.036,P=0.755)。结论本研究采用两样本孟德尔随机化的方法,表明在男性中严重的酒精依赖会增加肌萎缩侧索硬化的风险。     

两样本孟德尔随机化：肠道菌群与偏头痛

目的 利用两样本孟德尔随机化方法探究211种肠道菌群与偏头痛(migraine)的因果关系。方法 采用公共数据库MiBioGen研究的肠道菌群GWAS数据(n=13 266)和芬兰数据库(The FinnGen Biobank)的偏头痛GWAS数据(n=280 567),根据预设的阈值提取与211种肠道菌群相对丰度显著相关的独立遗传位点作为工具变量(Instrumental variables, IVs)。主要采用逆方差加权法进行分析，根据效应指标优势比(Odds ratio,OR)和95%置信区间(Confidence interval,CI)评估结果。使用留一法、异质性检验、水平基因多效性检验来验证结果的稳定性和可靠性。结果 双歧杆菌目(Bifidobacteriales)、双歧杆菌科(Bifidobacteriaceae)、双歧杆菌属(Bifidobacterium)丰度的升高，可以降低偏头痛的患病风险，其优势比与置信区间分别为(OR=0.77,95%CI:0.70～0.83)、(OR=0.77,95%CI:0.70～0.83)和(OR=0.82,95%CI:0.74～0.91)...     

孟德尔随机化研究甲状腺功能减退症与多囊卵巢综合征的因果关系

目的 采用孟德尔随机化(Mendelian randomization, MR)方法研究甲状腺功能减退症(Hypothyroidism)与多囊卵巢综合征(Polycystic ovarian syndrome, PCOS)之间的双向因果关系。方法 根据预设的阈值从甲状腺功能减退症的全基因组关联研究数据(Genome-wide association study, GWAS)中提取出与其强相关的独立单核苷酸多态性位点(Single nucleotide polymorphism, SNP)作为工具变量，主要采用逆方差加权法(inverse-variance weighted, IVW)评估甲状腺功能减退症相关表型与PCOS发生风险的关联；同时采用Cochran Q检验评估SNP的异质性。使用MR-PRESSO检验检测是否存在离群的SNP。MR-Egger截距测试检验SNP的水平多效性。“留一法”敏感性分析检验MR结果是否受单个SNP的影响。结果 共筛选到141个与甲状腺功能减退症相关的SNP,IVW分析法显示甲状腺功能减退症和PCOS之间存在因果关系(OR=16.725,95%CI:1....     

基于两样本孟德尔随机化的克罗恩病与胰腺癌发病的因果关系研究

目的 本文通过两样本孟德尔随机化的方法探究克罗恩病与胰腺癌发病之间的因果关系。方法 从全基因组关联分析研究数据中筛出与胰腺癌相关的单核苷酸多态性作为基因工具变量,用逆方差加权分析方法、加权中值方法和MR-Egger方法三种孟德尔随机化方法研究克罗恩病与胰腺癌发病的因果关联。比值比和95%可信区间作结果,OR值评价克罗恩病与胰腺癌发病的因果关系。结果 筛选到80个相关SNPs。逆方差加权分析法显示克罗恩病是胰腺癌发病的危险因素（OR=1.106;95%CI:1.009～1.213;P<0.05）,MR-Egger效应估计值（OR=1.328;95%CI:1.042～1.693;P<0.05）的统计结果与之类似;加权中值法（OR=1.140;95%CI:0.996～1.305;P>0.05）的统计结果无统计学意义。敏感性分析显示研究结果没有产生偏倚,异质性检验表明SNPs不存在异质性。结论 克罗恩病与胰腺癌发病存在正向因果关联。加强克罗恩患者胰腺癌相关的筛查、预防对降低克罗恩患者发生胰腺癌的风险具有重要意义。     

体力活动与卵巢癌因果关联的孟德尔随机化研究

目的 采用两样本孟德尔随机化研究方法分析体力活动与卵巢癌发生风险的因果关联。方法 体力活动和卵巢癌的数据分别来自两项公开发表的最新最大的欧洲人群全基因组关联研究汇总数据。从体力活动汇总数据中筛选出16个高度相关的单核苷酸多态性位点作为工具变量。利用逆方差加权法、MR-Egger回归和加权中位数估计法评估体力活动与卵巢癌之间的因果效应。另外采用排除回文结构工具变量、排除多效性工具变量和“留一法”进行敏感性分析。结果 逆方差加权法发现遗传预测的体力活动与卵巢癌发生风险无关（OR=0.91, 95%CI:0.66～1.24）,MR-Egger回归（OR=0.61, 95%CI:0.13～2.87）和加权中位数估计法（OR=1.16, 95%CI:0.79～1.69）与逆方差加权法得到的结果一致。敏感性分析表明结果具有稳健性。结论 未有充分的证据支持体力活动与卵巢癌发生风险之间存在因果关联。     

双向孟德尔随机化法探究酮体与认知表现的因果关系

目的 采用双向两样本孟德尔随机化(MR)方法来探究酮体与认知表现之间的因果关系。方法 采用欧洲全基因组关联研究汇总数据中的单核苷酸多态性(SNPs)作为工具变量。将逆方差加权法作为主要方法来评估酮体工具变量与认知表现的因果关系,并进行通路和功能富集分析预测潜在的机制。使用Bonferroni法校正显著性阈值确定显著性(P=0.017)。结果 两样本孟德尔随机化分析显示,乙酰乙酸(β=0.075, 95%CI:0.002～0.148,P=0.045)和β-羟基丁酸(β=0.082, 95%CI:0.039～0.126,P=2.36×10^-4)与认知表现有因果关系,但经Bonferroni法校正后,乙酰乙酸与认知表现的因果关系不再显著。此外,发现丙酮与认知表现之间不存在因果关系(β=0.041, 95%CI:-0.018～0.099,P=0.175)。通路和功能富集分析表明,β-羟基丁酸对认知表现的影响机制可能通过多种方式,如胆固醇代谢和胰岛素抵抗。在反向孟德尔随机化分析中,我们没有发现认知表现对酮体方向上的因果关系。结论 β-羟基丁酸水平与更好的认知表现之间存在因果...     

年龄相关听力损失与阿尔茨海默病相关性的两样本孟德尔随机化分析

目的 本研究通过双样本孟德尔随机化分析探讨年龄相关性听力损失与阿尔茨海默病之间的因果关系。方法 采用逆方差加权、MR-Egger回归、加权中位数、简单模型和加权模型五种方法来评估年龄相关性听力损失与阿尔茨海默病风险之间的关系,采用敏感性分析（多效性、异质性和留一性检验）来评估结果的稳健性。结果 统计结果显示与年龄相关的听力损失和阿尔茨海默病之间未见因果关联（逆方差加权法：OR=1.052 6,95%CI:0.715 5～1.548 5;MR-Egger:OR=1.134 7,95%CI:0.212 3～6.066 0;加权中位数法：OR=0.890 8,95%CI:0.528 1～1.502 5;简单模型法：OR=0.715 7,95%CI:0.250 5～2.044 5;加权模型法：OR=0.747 0,95%CI:0.315 3～1.769 8）,这些结果与阿尔茨海默病代理数据相一致(逆方差加权法：OR=0.956 0,95%CI:0.900 8～1.014 6;MR-Egger:OR=0.988 7,95%CI:0.772 9～1.264 7;加权中位数法：OR=0.948 7,...     

两样本孟德尔随机化研究肾功能对骨质疏松症的因果关系

目的 采用两样本孟德尔随机化方法研究肾功能对骨质疏松症的因果关系。方法 对全基因关联研究(GWAS)的汇总数据进行分析。eGFR和OP的GWAS汇总数据来源于东亚。eGFR数据共包含143 658例样本,OP数据一共212 453例样本。利用与肾功能(以eGFR衡量)密切关联的单核苷酸多态性作为工具变量,分别用逆方差加权法、MR-Egger回归和加权中位数法行两样本孟德尔随机化分析,以OR值评价肾功能对骨质疏松症风险的因果关系。结果 共纳入66个单核苷酸多态性作为工具变量,MR-Egger回归结果表明不存在基因的水平多效性(P=0.449)。逆方差加权法结果显示eGFR的降低可增加患骨质疏松症的风险,对数变换后的eGFR每降低一个标准差,会导致骨质疏松的风险增加24%(OR=1.24, 95%CI:1.05～1.46)。结论 肾功能降低对骨质疏松症的发生可能存在因果关联。     

基于孟德尔随机化方法分析尿酸与慢性肾脏病的因果关系

目的 基于孟德尔随机化方法探究尿酸（uric acid,UA）与慢性肾脏病（chronic kidney disease,CKD）之间的因果关系。方法 以UA与CKD形成的全基因组关联研究数据库为基础,筛选出单核苷酸多态性位点(single nucleotide polymorphisms,SNPs),参数设为P＜5×10-8、连锁不平衡r2＜0.01、kb = 10 000,运用逆方差加权法(inverse - variance weighted,IVW)、MR - Egger回归法、加权中位数法、基于众数的简单估计和基于众数的加权估计5种互补的方法进行孟德尔随机化分析,以OR值及95%CI评价UA与CKD之间的因果关系,利用异质性检验、水平多效性、逐个剔除检验等进行敏感性分析。设UA与痛风（GOUT）之间的因果关系为阳性对照。结果 共筛选24个与UA和CKD均相关的SNPs,IVW结果显示UA与CKD之间具有相关性,但不存在因果关系,UA升高不是CKD发生的危险因素（OR = 0.95,95%CI: 0.78～1.16）。以GOUT作为阳性对照,IVW结果显示UA与GOUT之间具有显著的因果关系（OR = 1.60,95%CI: 1.22～2.10）,随着UA水平的不断升高GOUT的患病率逐渐增加。敏感性分析表明MR结果具有稳健性。结论 UA与CKD之间有相关性,但未有充分的证据表明UA与CKD之间存在因果关系,故单纯降低UA不能减少CKD的发生风险。     

高密度脂蛋白胆固醇与肝癌发病的孟德尔随机化研究

目的 验证高密度脂蛋白胆固醇与肝癌发生之间是否存在因果关联。方法 本研究基于公开的全基因组关联分析（GWAS）数据库，其暴露因素为高密度脂蛋白胆固醇（HDL-C），来源于为期随访12年的韩国基因组与流行病学研究队列（Ko GES），其结局变量为肝癌，来源于为期随访10年的日本生物样本库（BBJ）。分别用IVW、MR-Egger、WME三种方法对其进行两样本孟德尔随机化分析，并利用Cochran Q检验和MR PRESSO进行MR前提假设验证，以及leave-one-out法进行敏感性分析。结果 两样本孟德尔随机化的IVW显示，HDL-C与肝癌发生（OR=0.938, 95%CI:0.779～1.129, P=0.496）无统计学意义。MR-Egger回归法（OR=1.167, 95%CI:0.858～1.587, P=0.329）、加权中位数（OR=0.921, 95%CI:0.722～1.175, P=0.508）分析也有相似的结果。经检验得分析结果中未发现存在多效性偏倚（截距项为-0.019, P=0.089），异质性较小；敏感性分析显示MR分析结果稳健。本次所选的52个工具变量均...     

臀围与冠心病的因果推断——孟德尔随机化分析

目的采用孟德尔随机化方法探索调整体质指数的臀围与冠心病患病风险的因果关系。方法利用大样本全基因组关联研究数据,从人体测量特征基因研究(GIANT)数据库(211 114名欧洲人)中筛选出与臀围密切相关且相互独立的遗传位点,进而在冠心病全基因重复验证与Meta分析合作组(CARDIoGRAM)数据库(86 995名欧洲人)中选出对应的遗传位点,确定工具变量,通过逆方差加权法的两样本孟德尔随机化方法,以OR值反映臀围与冠心病的因果关系。结果共筛选出70个单核苷酸多态性(SNP)作为工具变量,逆方差加权法结果表明,臀围每增加1倍标准差,冠心病的风险降低16.9%(OR=0.831,95%CI 0.730～0.946);MR-Egger回归截距表明遗传多效性不会对结果造成偏倚(截距为-0.0012,P=0.875);"leave-one-out"分析未发现非特异性SNP。结论臀围增加会降低冠心病的患病风险。     

甲状腺功能障碍与肩周炎的因果关系：一项双向、两样本孟德尔随机化研究

目的 流行病学和观察性研究表明，甲状腺功能障碍(包括甲状腺功能减退和甲状腺功能亢进)与肩周炎之间存在关联。然而，它们之间的因果关系仍不清楚。因此，本研究旨在使用双向两样本孟德尔随机化(MR)分析方法探究甲状腺功能障碍与肩周炎之间的因果关系。方法 从公开的全基因组关联研究(GWAS)中提取先天或后天的甲状腺功能减退症(n=86 169)、甲状腺功能亢进(n=173 938)和肩周炎(n=451 099)的数据。通过设置关联强度和消除连锁不平衡，筛选出最后作为工具变量(IVs)的SNP。应用以逆方差加权(IVW)为主要分析方法，加权中位数、MR-Egger、简单模式和加权模式为补充分析方法来评估因果效应。并且进行了异质性和多效性检验。同时使用留一分析评估结果的稳定性。结果 遗传预测甲状腺功能减退与肩周炎风险增加具有因果关系(IVW,OR=1.074,95%CI:1.010～1.143,P=0.023)。相比之下并未在甲状腺功能亢进对肩周炎及肩周炎对甲状腺功能障碍中发现关联。根据敏感性分析，本研究研究不具有水平多效性。结论 本研究证实了甲状腺功能减退会增加肩周炎的风险，在甲状腺功能亢进对肩周...     

系统性红斑狼疮与心力衰竭：两样本孟德尔随机化研究

目的 系统性红斑狼疮(systemic lupus erythematosus, SLE)与心力衰竭之间的关系已受广泛关注,但观察性研究中二者的因果关联是不清楚的,因此本研究旨在使用两样本孟德尔随机化(Mendelian randomization, MR)研究来探究SLE和心力衰竭之间的因果关联。方法 从全基因组关联研究(Genome-wide association study, GWAS)汇总数据中提取与SLE相关的单核苷酸多态性(single nucleotide polymorphisms, SNPs),并选择符合标准的工具变量。然后,主要用逆方差加权法(inverse variance-weighted, IVW)进行两样本孟德尔随机化分析,并用MR-PRESSO、加权中值法、MR-Egger法和最大似然比来评估结果的稳健性。另外采用MR-Egger截距项评估工具变量的水平多效性,Cochran Q检验及留一法用来评估异质性。此外,多变量MR研究被用于进一步验证结果的稳健性。结果 IVW结果显示SLE的遗传易感性与心力衰竭相对风险增加呈正相关(OR=1.024, 95%CI...     

克罗恩病与卵巢癌发病风险的因果关联一项两样本孟德尔随机化方法的探究

目的 本文通过孟德尔随机化（Mendelian randomization,MR）方法,探究克罗恩病（Crohn disease,CD)与卵巢癌（ovarian cancer)两种疾病间发病的因果关系。方法 从克罗恩病的全基因组关联研究数据（genome wide association study,GWAS）中提取出合适的单核苷酸多态性（single nucleotide polymorphism,SNP）作为基因工具变量,利用逆方差加权法、加权中值方法和MR-Egger法探究克罗恩病与卵巢癌发病风险的因果关联。分析结果以比值比和95%可信区间显示,以比值比评价克罗恩病与卵巢癌之间的因果关系。结果 共筛选到110个与克罗恩病相关的SNP,逆方差加权分析法显示克罗恩病和卵巢癌之间存在因果关系(OR=1.11,95%Cl 1.04～1.19),加权中值法(OR=1.14,95%Cl 1.03～1.26)计值进一步肯定CD是卵巢癌发病的危险因素。水平多效性、异质性检验及留一法显示工具变量不会对结果造成偏倚。结论CD与卵巢癌之间存在正向因果关联。加强CD患者卵巢癌相关指标的筛查,预防患者卵巢...     

饱和脂肪酸与心血管疾病的风险：基于孟德尔随机化的研究

目的 本研究旨在探究饱和脂肪酸(Saturated fatty acid, SFA)与心血管疾病之间是否具有因果关联。方法 采用两样本孟德尔随机化(Mendelian randomization, MR)研究,利用与SFA相关的遗传变异作为工具变量,从GWAS数据集中提取冠心病、心房颤动、心力衰竭、缺血性卒中和10种心血管危险因素作为主要结局指标。以随机效应逆方差加权法(IVW)作为主要的分析方法,并通过敏感性分析——异质性检验、水平多效性检验及逐个剔除检验来评估MR结果的可靠性。结果 SFA与冠心病(OR:1.30, 95%CI:1.15～1.46,P=1.46×10^-5)、总胆固醇(OR:2.02, 95%CI:1.86～2.18,P=6.71×10^-67)、低密度脂蛋白胆固醇(OR:1.46, 95%CI:1.39～1.54,P=2.29×10^-46)及甘油三酯(OR:2.22, 95%CI:2.03～2.42,P=4.43×10^-72)之间的关联具有统计学意义,与心房颤动、心力衰竭、缺血性卒...