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

目的 采用两样本孟德尔随机化研究方法分析体力活动与卵巢癌发生风险的因果关联。方法 体力活动和卵巢癌的数据分别来自两项公开发表的最新最大的欧洲人群全基因组关联研究汇总数据。从体力活动汇总数据中筛选出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）与逆方差加权法得到的结果一致。敏感性分析表明结果具有稳健性。结论 未有充分的证据支持体力活动与卵巢癌发生风险之间存在因果关联。     

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

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

血糖特征与结直肠良性肿瘤因果关联的孟德尔随机化研究

目的 探索血糖特征（空腹血糖、空腹胰岛素、糖化血红蛋白、胰岛素抵抗)与结直肠良性肿瘤（结直肠腺瘤、直肠良性肿瘤、结肠良性肿瘤）发病之间的因果关系。方法 血糖特征的工具变量数据来自MAGIC和GLGCL联盟,结直肠良性肿瘤的全基因组数据来源于Finn Gen联盟。研究采用随机效应逆方差加权法(inverse-variance weighted, IVW)、MR-Egger回归、加权中位数法(weighted median method,WMM)三种回归模型对四项血糖特征与结直肠良性肿瘤的因果关系进行孟德尔随机化分析。结果 IVW模型计算结果显示空腹血糖（OR=0.999,95%CI:0.997～1.001）、空腹胰岛素（OR=1.002,95%CI:0.998～1.006）、糖化血红蛋白（OR=1.001,95%CI:0.999～1.004）、胰岛素抵抗（OR=0.999,95%CI:0.997～1.002）与结直肠腺瘤之间无因果关系;空腹血糖（OR=0.909,95%CI:0.684～1.209）、空腹胰岛素（OR=1.220,95%CI:0.706～2.107）、糖化血红蛋白（OR=...     

恶性肿瘤与骨密度因果关系的双向两样本孟德尔随机化研究

目的 运用双向两样本孟德尔随机化的方法，探讨恶性肿瘤和骨质疏松之间的因果关系。方法 从Finngen数据库获取20种恶性肿瘤(前列腺癌、淋巴和造血器官癌、膀胱癌、脑癌、乳腺癌、肺和支气管癌、宫颈癌、结肠癌、肾癌、喉癌、肝和肝内胆管癌、食管癌、卵巢癌、胰腺癌、直肠癌、皮肤癌、小肠癌、胃癌、睾丸癌、甲状腺癌)数据集，从GEFOS数据库获取全身骨密度(TB-BMD)数据集。运用逆方差加权法、加权中位数法和MR-Egger法等评估20种恶性肿瘤与TB-BMD的双向因果关系。采用Cochran Q检验、MR-Egger回归法和Leave-one-out法评价研究结果稳定性和可靠性。结果 在乳腺癌、喉癌与TB-BMD的因果关系方面，观察到具有统计学意义的因果关系，其余18种恶性肿瘤与TB-BMD无显著因果关联。乳腺癌可能降低TB-BMD(β=-0.036,95%CI:-0.056～-0.015,P=0.001),喉癌可能增加TB-BMD(β=0.011,95%CI:0.002～0.021,P=0.022),结果具有统计学意义。加权中位数法和MR-Egger法与逆方差加权法关联方向一致。Cochra...     

lncRNAs与肺腺癌因果关联的两样本孟德尔随机化研究

目的 利用两样本孟德尔随机化方法探究长链非编码RNAs(long-chain non-coding RNAs, lncRNAs)与肺腺癌发病风险之间的因果关联。方法 采用肺腺癌全基因组关联分析数据和eQTL Gene联盟的cis-eQTL数据集,将与肺腺癌密切相关的单核苷酸多态性(single nucleotide polymorphisms, SNPs)作为工具变量,运用逆方差加权法、MR-Egger回归模型、加权中位数分析法、简单模式法和加权模式法5种两样本孟德尔随机化(Mendelian randomization, MR)模型来评估lncRNAs与肺腺癌之间的因果效应,并进行异质性检验、基因多效性检验和敏感性分析来评估结果的可靠性和稳定性。结果 PVT1、LINC00824、Z94721.1与肺腺癌的关联效应值差异有统计学意义(均P<0.05),PVT1(与肺腺癌的效应值OR=0.79,95%CI:0.71～0.89)、LINC00824(与肺腺癌的效应值OR=0.59,95%CI:0.42～0.83)降低了原发性支气管肺癌(肺癌)的发病风险;而Z94721.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不存在异质性。结论 克罗恩病与胰腺癌发病存在正向因果关联。加强克罗恩患者胰腺癌相关的筛查、预防对降低克罗恩患者发生胰腺癌的风险具有重要意义。     

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

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

舟山市噪声作业工人高频听力损失的影响因素分析

目的分析舟山市噪声作业工人高频听力损失情况及影响因素,为防制噪声职业病提供依据。方法通过作业场所职业病危害申报与备案管理系统和中国疾病预防控制信息系统职业病与职业卫生信息监测系统收集2015—2019年舟山市噪声职业病危害申报企业资料和噪声职业健康监护资料,分析噪声作业工人的高频听力损失现状,采用多因素Logistic回归模型分析噪声作业工人高频听力损失的影响因素。结果 2015—2019年舟山市共申报噪声职业病危害企业946家,接触噪声危害28 451人,进行健康检查25 066人,体检率为88.10%。检出高频听力损失5 520人,检出率为22.02%。多因素Logistic回归分析显示,地区(嵊泗县:OR=5.474,95%CI:4.090～7.576;定海区:OR=6.747,95%CI:5.224～7.172;岱山县:OR=6.515,95%CI:5.459～7.347)、男性(OR=2.498,95%CI:2.271～2.747)、年龄(30岁～:OR=1.773, 95%CI:1.537～2.044; 40岁～:OR=3.248, 95%CI:2.813～3.749; 50岁～:OR=6.908,95%CI:5.955～8.013;60～74岁:OR=11.860,95%CI:9.892～14.220)、工龄(10年～,OR=1.118,95%CI:1.043～1.199;20年～,OR=1.218,95%CI:1.079～1.376;30～46年,OR=1.475,95%CI:1.260～1.727)、制造业(OR=1.773,95%CI:1.491～2.017)、企业经济类型(集体经济,OR=2.180,95%CI:1.238～3.839;私有经济,OR=3.534,95%CI:1.999～6.250;港澳台经济,OR=1.747,95%CI:1.012～3.275)、小型企业(OR=1.364,95%CI:1.040～1.769)和噪声接触情况(接触噪声+苯:OR=1.418,95%CI:1.094～1.837;噪声+矽尘:OR=2.547,95%CI:1.506～4.308)是接噪工人高频听力损失的影响因素。结论舟山市噪声作业工人听力损害情况较为严峻,应加强年龄较大、工龄较长的男性作业工人、小型私有制造企业的噪声危害管理,关注噪声作业人员的职业健康监护。     

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

目的 采用两样本孟德尔随机化研究方法评估高血压与焦虑发生风险之间潜在的因果关系。方法 研究数据来自全基因组关联研究(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)截距均没有观察到水平多效性和异质性。结论 研究表明,高血压与焦虑发生之间可能存在正向因果关系。     

肠道微生物与失眠之间的因果关系——一项两样本孟德尔随机化分析

目的 利用两样本孟德尔随机化方法探讨肠道微生物与失眠之间的因果关系。方法 使用了来自Mi BioGen联盟进行的最大的全基因组关联研究荟萃分析的肠道微生物群的汇总统计数据。失眠的汇总统计数据来自全基因组关联研究GWAS：失眠测量。采用逆方差加权、MR-Egger、加权中位数、简单众数法和加权众数法来研究肠道微生物群与失眠之间的因果关系。通过一系列敏感性分析来评估结果的稳定性。结果 MR分析显示柔膜菌纲（β=-0.86, OR=0.42,95%CI:0.25～0.71）和未知科（β=-0.86, OR=0.42,95%CI:0.25～0.71）对失眠具有保护作用。此外,本研究还发现了其他14个类群与失眠之间存在潜在因果关系（P<0.05）,分别为草酸杆菌科（F=170.153,P=0.034)、马文布赖恩特菌属（F=162.702,P=0.006)、瘤胃球菌UCG002属（F=298.323,P=0.046)、罕见小球菌属（F=198.442,P=0.045)、厌氧棒状菌属（F=222.074,P=0.032)、伽玛变形菌纲（F=91.816,P=0.006)、巴斯德菌目（F=23...     

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

目的 探究血脂六项指标（甘油三酯、总胆固醇、高密度脂蛋白、低密度脂蛋白、载脂蛋白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与冠心病之间存在因果关系。     

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

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

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

目的 系统性红斑狼疮(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）方法探究生活方式、饮食、睡眠及社会经济地位等因素与支气管哮喘的因果关系。方法 研究主要使用逆方差加权法(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截距检测结果不存在多...     

身体活动水平与精神类疾病的孟德尔随机化研究

目的 应用孟德尔随机化(Mendelian randomization, MR)探索身体活动水平与十种精神类疾病的因果效应。方法 本研究利用身体活动和精神类疾病的全基因组关联分析汇总数据构建两样本MR模型。选择与暴露显著相关并相互独立的遗传变异作为工具变量,应用MR多效性残差和与异质性检测和Radial MR异质性检测方法剔除多效性工具变量。同时应用F统计量,Q统计量,PRESSO global检验及MR-Egger截距项检验等方法评价工具变量的有效性。应用基于修正权重的逆方差加权法作为主要MR模型,并应用极大似然估计法、加权中位数估计法、MR稳健调整轮廓得分法进行敏感性分析。对于与身体活动水平存在显著关联的精神类疾病,进行反向MR分析,探索精神类疾病对身体活动水平的因果效应。结果 MR分析结果显示身体活动水平升高对注意力缺陷多动障碍(Attention deficit hyperactivity disorder, ADHD;OR:1.956;95%CI:1.605～2.384;P=2.45×10^-11)及精神分裂症(Schizophrenia, SCZ;OR:...     

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

目的 验证高密度脂蛋白胆固醇与肝癌发生之间是否存在因果关联。方法 本研究基于公开的全基因组关联分析（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个工具变量均...     

2019年上海某区噪声暴露工人高频听力损失及其影响因素分析

目的 了解上海市某区在岗噪声作业人员中的高频听力损失，分析其影响因素。方法 选取2019年在上海市某区噪声暴露在岗工人为调查对象，通过调查问卷调查噪声暴露工人噪声暴露情况及听力检测数据分析高频听力损失发生情况，采用多因素logistic回归模型分析高频听力损失的相关影响因素。结果 2019年该区在岗噪声职业健康检查人数为13 796人，高频听力损失720人，检出率为5.23%。男性(OR=2.387)、年龄(30～39岁：OR=3.027;40～49岁：OR=7.765;≥50岁：OR=10.292)、工龄(10～19 a:OR=2.085;≥20 a:OR=2.060)、吸烟(OR=2.902);企业规模(小型和微型：OR=2.394)、经济类型(私营：OR=1.581)、行业类型(其他制造业：OR=1.600;通用设备制造业：OR=1.896;金属制品业：OR=2.477;造船业：OR=3.058)、与噪声联合暴露(OR=1.508)是噪声暴露工人高频听力损失的影响因素。结论 该区噪声作业高频听力损失情况不容乐观，男性、年龄≥30岁、工龄≥10 a、吸烟、小微、私营、造船业、与噪声...     

基于两样本孟德尔随机化方法分析血液代谢物与慢性肾病的因果关系

目的 利用孟德尔随机化方法,探究486种血液代谢物与慢性肾病的因果关系.方法 血液代谢物和慢性肾病全基因组关联分析汇总数据均来自于网络公开数据库,主要采用逆方差加权分析法(inverse-variance weighting,IVW),探究每种血液代谢物与慢性肾病之间的因果关系,并以MR-Egger回归、加权中位数法、...     

噪声性听力损失影响因素Logistic分析

目的 探讨噪声性听力损失的影响因素。 方法 选择某市铁路运输设备制造行业全部噪声作业工人519人作为调查对象,调查问卷收集工人一般情况、个人习惯、既往史等资料,收集2012~2015年职业健康监护资料,进行非条件二分类Logistic回归分析。 结果 与噪声性听力损失相关的因素有性别(OR=0.312,95%CI=0.161~0.606)、年龄(50~58岁:OR=15.008,95%CI=5.757~29.815,40~49岁:OR=7.137,95%CI=4.082~112.479;30~39岁:OR=3.440;95%CI=1.940~6.100)、传染病史(OR=2.066,95%CI=1.272~3.355)、现时身体质量指数(BMI)(OR=1.084,95%CI=1.014~1.159)、工作岗位噪声等效声级(LEX,8 h)(OR=1.028,95%CI=1.003~1.055)。 结论 噪声性听力损失影响因素有性别、年龄、传染病史、现时BMI、等效声级。噪声性听力损失是多因素共同作用结果,应积极采取预防措施。