Mendelian Randomization (MR) 文章学习

有关MR方法的学习请参考:https://www.mubu.com/doc/5Ke3fGOVKrY

今日学习一篇2023年发表于BMC medicine上的文章,题目是:Circulating levels of micronutrients and risk of infections: a Mendelian randomization study。

摘要

背景:微量元素的缺失增加感染的易感性。目前有关微量元素和感染之间的观察性研究及RCT很有限。
方法:Two-sample MR,数据来源为欧洲血统(European ancestry)的独立队列(independent cohorts)。表型(危险因素:3种感染)数据来源为UK Biobank and FinnGen。统计分析为:反方差加权法(IVW,Inverse variance-weighted MR)和一系列敏感性分析。统计显著性P设置为2.08E-3
结果:研究发现血液循环中的铜的水平和胃肠道感染明显相关。一个SD的增加与0.91的OR值相关。这些发现被后续的敏感性分析证实为稳健的。

方法

1.报道按照STROBE-MR的方法。(Skrivankova VW, Richmond RC, Woolf BAR, Davies NM, Swanson SA, VanderWeele TJ, et al Strengthening the Reporting of Observational Studies in Epidemiology using Mendelian Randomisation (STROBE-MR): explanation and elaboration. BMJ. 2021;375:n2233.
2.Flow chart:

Flow chart

3.数据:共14 GWASs数据集:8个暴露因素数据集和6个结局数据集。
为了减少人群分层的偏倚,数据选择了同为欧洲血统(European ancestry)的数据。

  1. 表型(危险因素)GWAS数据的两个搜索平台:GWAS Catalog and PubMed
  2. 表型(危险因素)GWAS数据的临床显著性确定:P ≤ 5E−08
  3. 结局(疾病)GWAS数据的两个搜索平台:UK Biobank 和 FinnGen Release 6。
  4. 使用SAIGE对与遗传相关的性别、年龄及前4个(UK Biobank)或前10个主成分进行调整(FinnGen)。
  5. 汇总统计(summary statistics)的meta分析采用METAL的随机效应模型(random-effects model)(Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of
    genomewide association scans. Bioinformatics. 2010;26(17):2190–1.
    )。
  6. Cochran’s Q statistical test (included in METAL)用于评价两个队列结局研究的基因工具的异质性。
  7. 统计效能: F统计量(后续再添加
  8. MR主分析方法:
    11.1 计算每个SNP的Wald比(即SNP-结果关联/SNP-暴露关联)。若不同SNP对应一个危险因素(微量元素),则选择使用IVW(反向方差加权法)计算效应量(Wald比)。
    11.2 每个遗传预测的SD增加一个单位, 报道的结局关联增加1个OR。
    11.3 MR分析针对UK Biobank、FinnGen Release 6及meta分析的汇总分别进行。
    11.4 名义上的P<0.05。 多重检验的矫正统计量p=0.05/(8 exposures*3 outcomes ) = 0.05/24 = 2.08E−03。
  9. MR敏感性分析方法:
    12.1 水平多效性检测(Cochran’s Q statistical test)
    12.2 评估基因工具的独立性及排他性:PhenoScanner V2工具http://www.phenoscanner.medschl.cam.ac.uk/ (accessed 30 October 2022).)
    12.3 leave-one-out 分析 (>2 SNP): 检验IVW估计值的稳健性。及是否存在特殊的SNP(促进这种关联,也就是多效性)
    12.4 Cu的潜在多效性分析
  10. 非严格标准的基因工具的二次分析
    13.1 纳入低标准的基因工具变量(r 2 < 0.01 and P ≤ 5E−06):可能引入弱相关变量
    13.2 MR-RAPS:纠正引入弱相关变量的偏倚。
  11. Post hoc analyses (事后分析)
    14.1 分析结果中Cu与胃肠道感染之间的关联。
    14.2 两种事后分析:1.提取另外的肠道感染的GWAS数据,进行独立队列的MR分析。2. 分析反向因果。进行遗传预测的Cu血清水平和遗传预测的肠道感染的风险之间的MR分析。
  12. 数据分析工具:所有MR分析用的是R语言中的TwoSampleMR package。METAL则应用于对结局的Meta分析。

想法:MR研究的指标较多,但也有一定的规律,重点在于多读文献,筛选合适的数据集、工具变量,适当的分析及恰当的语言组织能力(telling a good story)。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 220,295评论 6 512
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,928评论 3 396
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 166,682评论 0 357
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,209评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,237评论 6 397
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,965评论 1 308
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,586评论 3 420
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,487评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 46,016评论 1 319
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,136评论 3 340
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,271评论 1 352
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,948评论 5 347
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,619评论 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,139评论 0 23
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,252评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,598评论 3 375
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,267评论 2 358

推荐阅读更多精彩内容