双向孟德尔随机化分析总胆固醇升高风险与胆石症发生风险的关系.docx
胆道疾病DOI:10.12449/JCH240322双向孟德尔随机化分析总胆固醇升高风险与胆石症发生风险的关系赵伟伟L2,杜晓旭戈宏读,1内蒙古民族大学临床医学院，内蒙古通辽0280002兴平市汤坊镇卫生院，陕西兴平7131003内蒙古民族大学医学院，内蒙古通辽028000三f言储:物颍，gehongyan_1999(ORQD:00-0003-3691-2805)摘要：目的运用双向孟德尔随机化评估总胆固醇升高风险与胆石症发生风险之间的关联方法总胆固醇与胆石症相关的单核甘酸多态性数据来自OPengwas公共雌库，对基于全基因组相关联研究(GWAS)的所有，匚总数据进行二次数据分析，分别选择与总胆固醇或者胆石症密切关联的遗传位点作为暴露或者结局变量，用Egger回归法、加权中位数法、IVW随机效应模型、IVW固定效应模型等方法进行双向孟德尔随机化分析，以。国直评价总胆固醇和胆石症之间的因果关系。结果以总胆固醇为暴露、胆石症为结局，去除异质性前总胆固醇-胆石症的总体OR=O.91z95%C:0.85-0.97;去除异质性后总胆固醇-胆石症的总体O8=093,95%C.0.89-0.97.以fl旦石症为暴露、总胆固醇为结局，去除异质性前胆石症-总胆固醇的总体OR=O.20,95%C:0.06-0.65;去除异质性后胆石症-总胆固醇的总体OR=Q.28,95%C:0.10-0.83.基因预测的总胆固醇和胆石症之间存在双向因果关系.结论本研究证实了总胆固醇与胆石症之间的双向因果关系，随着遗传上有关能使总胆固醇水平升高的等位基因的增加，发生胆石症的风睑降低反之，随着遗传上有关能使胆石症发病增高的等位基因的增加，总胆固醇水平升高的风险降低。关键词：总胆固醇；胆石症；双向孟德尔随机化分析Associationbetweentheriskofincreaseintotalcholesterolandtheriskofcholelithiasis:AbidirectionalMendelianrandomizationstudyZHAOWeiwei''2,DUXiaoxur,GEHongyarP.(1.CollegeofClinicalMedicine,InnerMongoliaMinzuUniversity,Tongliao,InnerMongolia028000,China;2.TangfangTownshipHospital,XingpinglShaanxi713100,China;3.MedicalCollege,InnerMongoliaMinzuUniversity,Tongliao,InnerMongoIia028000,China)Correspondingauthor:GEHongyan,gehongyan1999(3)(ORCID:0000-0003-3691-2805)Abstract:ObjectiveToinvestigatetheassociationbetweentheriskofincreaseintotalcholesterol(TC)andtheriskofCholelitbiasisbyusingbidirectionalMendelianrandomization(MR).MethodsTheopengwaspublicdatabasewasusedtoobtainthesinglenucleotidePOIymOfPhiSmdataassociatedwithTCandcholelithiasis,andasecondarydataanalysiswasperformedforallsummarydataofgenome-wideassociationstudies.TbegeneticlocicloselyassociatedwithTCorcholelithiasiswereselectedasexposureoroutcomevariables,andthebidirectionalMRanalysiswasperformedusingthemethodssuchasEggerregressi,Weightedmedian,IVWrandomeffectsmodel,andIVWfixedeffectsmodel,withoddsratio(OR)valuesforealuatingthecausalrelationshipbetweenTCandcholelithiasis.ResultsWithTCastheexposureandcholelithiasisastheoutcome,TC-ChoIeIithiasishadanoverallORvaIueof0.91(95%confidenceintervalCl:0.850.97)beforeeliminationofheterogeneityand0.93(95%C:0.890,97)aftereliminationofheterogeneity.WithcholelithiasisastheexposureandTCastheoutcome,TC-ChoIeIithiasishadanoverallORvalueof020(95%C:0.060.65)beforeeliminationofheterogeneityand028(95%C:0.100.83)aftereliminationofheterogeneity.TherewasabidirectionalcausalrelationshipbetweengeneticallypredictedTCandcholelithiasis.CondusionThisstudyconfirmsthebidirectionalcausalrelationshipbetweenTCandcholelithiasis.TheriskofcholelithiasisdecreaseswiththeincreaseinallelesassociatedwiththeelevationofTClevel;onthecontrary,theriskofelevatedTCleveldecreaseswiththeincreaseinallelesassociatedwiththeonsetOfcholelithiasis.Keywords:TotalCholesterol;Giolelithiasis;BidirectionalMendelianRandomizationAnalysis胆石症发病率在中西方国家都比较高，在西方国家一些大城市的医院中，大约有10%的成年人患有胆石症田。目前关于胆石症的风险因素有很多研究，但在这些研究中存在互相矛盾的地方，如有的研究m认为较低的总胆固醇水平是胆石症的危险因素，但有的研究【3】认为较低胆固醇水平是胆石症的保护性因素关于胆石症对总胆固醇水平的影响和因果关系尚未得到证实已有的研究大多依靠传统的观察性研究得出关联结论。传统的观察性研究由于本身的缺陷，受到大量混杂因素的干扰，因此本研究引入孟德尔随机化(MR)来进行因果推断，MR利用等位基因在减数分裂时随机分配到子代的原理模拟随机对照试验，从而规避了不同混杂因素对暴露和对照的影响I4.MR在观察性研究中使用遗传变异作为工具变量来检测和量化因果关系。本研究使用双向孟德尔随机化(bidirectionalMR)来探索总胆固醇水平和胆石症之间的风险关系和因果推断。1资料与方法1.1数据来源与研究设计通过搜索OPengwas(https:/gwas.mrcieu.ac.uk/datasets)数据库中有关总胆固醇水平和胆石症相关的全基因组相关联研究(GWAS)汇总数据，得到两组样本量相对较大和单核首酸多态性(SNP)个数比较多的GWAS数据，具体信息见表12°将其中一组数据中的任意一个GWAS数据和另外一组数据中的任意一个GWAS数据进行组合，得到四种不同的暴露和结局的组合。利用设计的bidirectionalMR模型来评估暴露和结局(总胆固醇和胆石症)之间的风险关系和因果懒M图6StandardMRBidirectionalMRU注：左侧，三SMR;谢bidirectionalMR.Gx,可案露X对结果Y的因果影响.Gyi工具变量，可用于估计展露丫对结果X的因果影响。U:混杂因素.假设1,与暴露显著相关；假设2,与混杂因素不相关；假设3,与结局不相关。图1bidirectionalMR的模型图Figure1ModeldiagramofbidirectionalMendelianrandomization1.2工具变量的选择从全基因组(GWAS)上获得P<5x10-8水平的与暴露相关的遗传变异。该遗传变异必须满足MR分析的三大假设:(1)关联性假设，SNP与暴露强相关；(2)独立性假设，遗传变异与影响"暴露和结局”的混杂因素相独立；(3)排他性假设，遗传变异只能通过暴露对结局发生作用，而不通过其他途径设置连锁不平衡参数(r2<0OO1),即两个SNP之间是完全连锁平衡的(这两个SNP的分配是完全随机的遗传距离设定为10000kb以选择SNP1从而确保其独立性，并排除连锁不平衡(linkagedisequilibrium,LD)对结果触响。表1总胆固醇GWAS数据汇总信息Table 1 SummaryoftheGWASdatafortotalcholesterol数超来源种族来源样本SSNP个数公布时间PMlD或数隔曦全球脂质遗传学协会(GLGC)混合人群18736524469822013年24097068代谢组学研究(met)欧洲人群II5078123218752020年Trait:TotalChoIesteroI-EUOpenGWASproject(mrcieu.ac.uk)表2胆石症GWAS数据汇总信息Table 2 SummaryoftheGWASdataforcholelithiasis数据来源种族来源SsSNP个数公布时间PMID或辘链接英国生物银行(Ukb)欧洲人群337199108945962017年Trait:Diagnoses-mainICDlO:K80ChoIeIithiasis-IEUOpenGWASproject(mrcieu.acUk)芬兰生物银行(finn)欧洲人群214167163804522021年Trait:GioIeIithiasis-IEUOpenGWASproject(mrcieu.acUk)通过GWAS提取满足MR三大假设的SNP,合并暴露与结局的数据集，并删除回文SNP。最后剩余的SNP即为指代暴露的最佳工具变量。1.3 bidirectionalMR方案本研究分为两部分，第一部分分别以表1总胆固醇GWAS数据为暴露、表2胆石症GWAS数据为结局，应用MR分析方法中的Egger回归法、中位数加权法、逆方差加权法随机效应模型和固定效应模型等主要方法来初步计算因果效应值（b值）、OH值及95%C瞥。然后将上述暴露和结局调换，以胆石症GWAS数据为暴露、以总胆固醇水平GWAS为结局进行MR分析。第二部分将第一部分的数据去除异质性后做同样分析。1.4 敏感性分析用MR分析R包对所涉及的SNP进行C。Chran'sQ检验，用于评估个体遗传变异之间的异质性o如果Cochran/sQ检验的P<Q.05,代表基因间存在异质性可能。一般来说当异质性很大（P<0.OOl）的时候，需要使用随机效应模型，当异质性小的时候，使用固定效应模型。检测SNP的水平多效性，当P<0.05代表存在水平多
