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

目的 采用孟德尔随机化方法分析抑郁症与间质性肺病(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之间可能存在正向因果关联。     

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

目的 采用两样本孟德尔随机化方法研究肾功能对骨质疏松症的因果关系。方法 对全基因关联研究(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)。结论 肾功能降低对骨质疏松症的发生可能存在因果关联。     

弓形虫感染与精神分裂症的两样本孟德尔随机化研究

目的 采用两样本孟德尔随机化方法探究弓形虫感染和精神分裂症之间的因果关系.方法 利用汇总的大样本GWAS数据提取与弓形虫血清抗体密切相关的遗传位点作为工具变量,分别运用MR-Egger回归、加权中位数和逆方差加权法进行孟德尔随机化分析,以OR值及95％CI评价弓形虫感染与精神分裂症之间是否存在关联.采用Egger-in...     

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

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

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

目的 采用两样本孟德尔随机化研究方法分析体力活动与卵巢癌发生风险的因果关联。方法 体力活动和卵巢癌的数据分别来自两项公开发表的最新最大的欧洲人群全基因组关联研究汇总数据。从体力活动汇总数据中筛选出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-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数据库中的大样本汇总数据，以单核苷酸多态性作为工具变量，在两样本孟德尔随机化模型中，运用5种孟德尔随机化方法评估16个可干预风险因素与肺结核发生之间因果关联。统计结果报告为优势比(Odds ratio,OR)及其95%置信区间(Confidence intervals,CI)。使用留一法、水平多效性检验进行敏感性分析。结果 目前吸烟情况(OR=39.24, 95%CI:4.12～373.63)会增加结核病发生的风险，受教育时间(OR=0.60, 95%CI:0.37～0.96)及转铁蛋白饱和度(OR=0.69, 95%CI:0.50～0.96)与TB发生风险降低之间也存在因果关系。留一法显示结果稳健，MR-Egger截距值未显示水平多效性的存在。结论 基因预测目前吸烟情况、受教育时间以及转铁蛋白饱和度与肺结核的患病风险存在因果关系。     

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

目的 本文通过两样本孟德尔随机化的方法探究克罗恩病与胰腺癌发病之间的因果关系。方法 从全基因组关联分析研究数据中筛出与胰腺癌相关的单核苷酸多态性作为基因工具变量,用逆方差加权分析方法、加权中值方法和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不存在异质性。结论 克罗恩病与胰腺癌发病存在正向因果关联。加强克罗恩患者胰腺癌相关的筛查、预防对降低克罗恩患者发生胰腺癌的风险具有重要意义。     

基于孟德尔随机化的血脂六项与冠心病因果关系研究

目的 探究血脂六项指标（甘油三酯、总胆固醇、高密度脂蛋白、低密度脂蛋白、载脂蛋白A、载脂蛋白B）与冠心病之间的因果关系，为血脂六项指标与冠心病发病风险关联提供遗传学证据支持。方法 研究采用MR-Egger回归、加权中位数法（weight median estimator, WME）、随机效应逆方差加权法（inverse-variance weighted, IVW）、加权模型4种回归模型对血脂六项指标与冠心病的因果关系进行孟德尔随机化分析。结果 IVW模型显示甘油三酯（OR=1.402）、总胆固醇（OR=1.414）、低密度脂蛋白（OR=1.867）、载脂蛋白B（OR=1.738）为冠心病的危险因素，高密度脂蛋白（OR=0.791）为冠心病的保护因素，IVW模型则显示载脂蛋白A为有统计学意义的影响因素；MR-Egger回归截距项P值均>0.05，即筛选出的SNP不存在基因多效性，因此孟德尔随机化在本研究中为因果推断的有效方法；异质性检验结果P值均<0.001，故应重点关注随机效应IVW模型。结论 利用两样本孟德尔随机化方法排除混杂因素和反向因果关联后，得到无偏估计的结果，据此可以确定甘油三酯、总胆固醇、高密度脂蛋白、低密度脂蛋白、载脂蛋白B与冠心病之间存在因果关系。     

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

目的采用孟德尔随机化方法探索调整体质指数的臀围与冠心病患病风险的因果关系。方法利用大样本全基因组关联研究数据,从人体测量特征基因研究(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。结论臀围增加会降低冠心病的患病风险。     

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

目的 本文通过孟德尔随机化（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患者卵巢癌相关指标的筛查,预防患者卵巢...     

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

目的 基于孟德尔随机化方法探究尿酸（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的发生风险。     

孟德尔随机化探究类风湿性关节炎与冠状动脉疾病因果关系以及步行速度的介导作用

目的 类风湿性关节炎(rheumatoid arthritis, RA)患者的冠状动脉疾病(coronary artery disease, CAD)发生风险会升高。由于关节疼痛/僵硬或其他关节问题, RA患者的身体素质通常会有所降低,比如步行速度下降。但是,RA、步行速度和CAD这3者之间的因果关联目前尚不清楚。本研究旨在采用网络孟德尔随机化方法研究RA与CAD因果关系及步行速度在该因果关系中的介导作用。方法 公共数据库下载RA、CAD和步行速度的全基因组关联研究（genome - wide association study, GWAS）统计结果数据,进行两样本的孟德尔随机化分析。RA数据共有103 688例样本,CAD数据有547 261例样本,步行速度有452 264例样本。选取与RA和步行速度相关的独立SNP作为工具变量。使用逆方差加权法与加权中值法进行因果关系分析,并采用MR - Egger和留一法进行模型评价。结果 分析结果显示,RA与CAD的总体因果效应是正相关(OR = 1.03, 95% CI:1.01~1.04, P = 3.90×10 - 5),RA与正常步行速度负相关(β = - 0.02, 95% CI: - 0.04～- 0.01, P = 3.15 × 10 - 4)。步行速度与CAD有显著的负因果效应(OR = 0.841, 95% CI: 0.81～0.88, P = 3.43×10 - 15),中介分析显示,步行速度对RA与CAD的关联效应的贡献比为15.8%。结论 总之,本研究结果表明RA会增加CAD的发生风险且步行速度可能在此因果关系中起到重要的介导作用。     

Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator

Developments in genome‐wide association studies and the increasing availability of summary genetic association data have made application of Mendelian randomization relatively straightforward. However, obtaining reliable results from a Mendelian randomization investigation remains problematic, as the conventional inverse‐variance weighted method only gives consistent estimates if all of the genetic variants in the analysis are valid instrumental variables. We present a novel weighted median estimator for combining data on multiple genetic variants into a single causal estimate. This estimator is consistent even when up to 50% of the information comes from invalid instrumental variables. In a simulation analysis, it is shown to have better finite‐sample Type 1 error rates than the inverse‐variance weighted method, and is complementary to the recently proposed MR‐Egger (Mendelian randomization‐Egger) regression method. In analyses of the causal effects of low‐density lipoprotein cholesterol and high‐density lipoprotein cholesterol on coronary artery disease risk, the inverse‐variance weighted method suggests a causal effect of both lipid fractions, whereas the weighted median and MR‐Egger regression methods suggest a null effect of high‐density lipoprotein cholesterol that corresponds with the experimental evidence. Both median‐based and MR‐Egger regression methods should be considered as sensitivity analyses for Mendelian randomization investigations with multiple genetic variants.     

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

目的 系统性红斑狼疮(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...     

孟德尔随机化探究2型糖尿病与急性心肌梗死的性别特异性因果关系

目的 观察性研究结果显示2型糖尿病（type 2 diabetes, T2D）是急性心肌梗死（acute myocardial infarction, AMI）发生的风险因素，且这种风险存在一定的性别差异。然而观察性研究常受到混杂因素的影响，明确两者的因果关系有助于AMI的防控。本研究旨在利用孟德尔随机化(Mendelian randomization, MR)模型揭示T2D与AMI之间的因果效应及性别特异性。方法 T2D和AMI的全基因组关联研究(genome - wide association study, GWAS)结果汇总数据来源于欧洲。T2D数据共包含890 002例样本，其中男性425 613例样本，女性464 389例样本。AMI数据一共355 246例样本。利用T2D显著相关且独立的单核苷酸多态性(single nucleotide polymorphism, SNPs)作为工具变量，在总体及不同性别中利用两样本MR方法评估T2D与AMI的因果关联。以逆方差加权法(inverse variance weighted, IVW)作为主要的分析方法，使用Cochran Q检验检测性别异质性。结果 分析结果显示，整体人群中T2D对AMI的因果效应呈正相关(OR = 1.06; 95%CI: 1.04 ~ 1.09)。男性人群与总体人群的结果相似(OR = 1.03; 95%CI: 1.02 ~ 1.04)，但是在女性人群中并没有观察到T2D与AMI的因果关系(OR = 1.01; 95%CI: 1.00 ~ 1.02)。异质性分析显示男性和女性之间的结果确实存在显著差异(P = 0.004)。结论 本研究结果表明T2D会增加AMI的患病风险，且这种效应在男性群体中更加明显，具有性别特异性。