首页 > 其他分享 >Nature Comm. | CoPheScan:一种考虑连锁不平衡的全表型组关联分析

Nature Comm. | CoPheScan:一种考虑连锁不平衡的全表型组关联分析

时间:2024-09-18 13:51:31浏览次数:13  
标签:LD Nature 变异 关联 表型 Comm CoPheScan 因果

分享一篇最近发表在NC的一篇文章:CoPheScan: phenome-wide association studies accounting for linkage disequilibrium。文章介绍了一种新的贝叶斯方法CoPheScan(Coloc adapted Phenome-wide Scan),用于在考虑连锁不平衡(LD)的情况下进行表型范围关联研究(Phenome-wide association studies, PheWAS)。CoPheScan通过模拟研究展示了其在控制假阳性率方面相比传统方法的优势,并在真实数据集上的应用中展示了其性能。

Nature Comm. | CoPheScan:一种考虑连锁不平衡的全表型组关联分析_数据

研究背景

PheWAS允许研究者探索单个遗传变异与多种表型之间的关联,这对于药物再利用和发现潜在的不良效应具有重要意义。然而,PheWAS受限于无法区分由于LD产生的混杂与真正的多效性。CoPheScan旨在解决这一问题,通过同时处理LD混杂和探索因果关联。

方法

CoPheScan基于coloc方法的改编,用于在已知某个变异对某个表型有因果影响的情况下,对其进行表型范围扫描以测试与其他表型/性状的因果关联。该方法考虑了两个性状在基因组区域内的遗传关联模式,并评估它们是否可能在该区域内共享一个因果变异。CoPheScan可以利用外部协变量数据,如表型之间的遗传相关性,并根据不同的LD信息可用性和分析者是否愿意对测试基因组区域内的因果变异数量做出假设,以不同的方式运行。

Nature Comm. | CoPheScan:一种考虑连锁不平衡的全表型组关联分析_数据_02

模拟研究

作者通过模拟研究展示了CoPheScan的性能。模拟数据包括具有零个、一个或两个因果变异的性状,以对应CoPheScan的三种假设。研究发现,CoPheScan在控制假阳性率方面表现优于传统方法,尽管在存在多个因果变异时敏感性略有下降。

真实数据应用

研究者使用CoPheScan对来自FinnGen疾病端点的136个疾病因果变异、527个血浆蛋白的505个pQTL变异以及3586个蛋白质截断变异(PTV)进行了PheWAS。在2275个英国生物银行(UK Biobank)的疾病表型中,CoPheScan识别了已知多效性基因如APOE的复杂性,并为TGM3在皮肤癌中的新因果角色提供了证据。

结果

CoPheScan在模拟数据上的表现显示,与传统方法相比,它在控制假阳性率方面有显著优势。在真实数据集上,CoPheScan能够识别与已知疾病相关的遗传变异,并揭示了新的可能的因果关联,例如TGM3与皮肤癌的关联。

讨论

CoPheScan通过考虑查询变异集合中的信息,不仅关注给定变异的小p值,还关注与其邻居相比的小p值,并估计由于信息而应向上加权的程度。CoPheScan的一个优势是能够利用SuSIE检测位点上的等位基因异质性。然而,CoPheScan也有其局限性,例如对LD信息的依赖,以及目前仅适用于单祖先研究。

数据和代码

Nature Comm. | CoPheScan:一种考虑连锁不平衡的全表型组关联分析_数据_03

作者:生物信息与育种,请关注同名微信公众号:生物信息与育种。

标签:LD,Nature,变异,关联,表型,Comm,CoPheScan,因果
From: https://blog.51cto.com/u_15668923/12045819

相关文章

  • Java调用Apache commons-text求解字符串相似性
    前言    在之前的一篇漂亮国的全球的基地博客中,我们曾经对漂亮国的全球基地进行了一些梳理。博文中使用的数据来源,重点是参考以为博主分享的KML的数据,同时针对其国内的基地部署信息,我们从互联网百科的数据中搜寻到一些。其实拿到这两份数据的时候,是存在一些问题的,比如,KML的......
  • commonJs和ESModule的区别
    1.规范的区别ESM模块的导入使用Import关键字,导出使用export关键字。commonJs导入使用require关键字,导出使用module.export 2.文件名后缀不一样(这一点不用记,针对Node)在node.js,默认将.js后缀文件识别为CJS模块,.cjs也是CJS模块,.mjs文件识别为esm模块。 3.模块加载时机......
  • SQLSTATE[HY000] [2013] Lost connection to MySQL server at 'reading initial commu
    错误信息 SQLSTATE[HY000][2013]LostconnectiontoMySQLserverat'readinginitialcommunicationpacket',systemerror:111 翻译成中文为:在读取初始化数据包时失去到MySQL服务器的连接,系统错误111。通讯包初始化失败,估计是不允许连接访问引起的。解决办法以下是......
  • WPF DataGrid ContextMenu CommandParameter Relative x:Type ContextMenu ,Path=Plac
    //xaml<DataGrid.ContextMenu><ContextMenu><MenuItemHeader="SerializeBinary"Command="{BindingBinSerializeCmd}"CommandParameter="{BindingRelativeSource={Relativ......
  • INFO20003 SQL Requesting Communications
    INFO20003S22024–ASSIGNMENT2v1.41INFO20003Semester2,2024Assignment2:SQLDue:Week8-Sunday15thSeptember2024,5:59pmMelbourneTime.Submission-ViaLMShttps://canvas.lms.unimelb.edu.au/Case:“Slarc”App“Slarc”:SuperLovelyAppfor......
  • COMM 1100 Foundations of Communication
    COMM1100(A11)FOUNDATIONSOFCOMMUNICATIONSTUDIESFall2024COURSEDESCRIPTIONThiscourseoffersacomprehensiveoverviewofwhatitmeanstostudycommunications.Studentswillexploreclassicdefinitionsandmodelsofcommunicationsandtracehowth......
  • PCIe进阶之TL:Common Packet Header Fields & TLPs with Data Payloads Rules
    1TransactionLayerProtocol-PacketDefinitionTLP有四种事务类型:Memory、I/O、Configuration和Messages,两种地址格式:32bit和64bit。构成TLP时,所有标记为Reserved的字段(有时缩写为R)都必须全为0。接收者Rx必须忽略此字段中的值,PCIeSwitch必须对其进行原封不......
  • 脚本 $'\r': command not found异常解决办法
    1、运维人员在日常工作中需要写一些shell脚本监控服务器资源使用情况:2、下面一个监控脚本:echo-e"服务器名称,CPU使用率,内存使用率,磁盘总空间,磁盘使用率,磁盘剩余空间"cpu=$(top-b-n1|grep"^%Cpu"|awk'{printf"%.2f%%",$2+$4}')mem=$(free-m|awk'NR==2{printf......
  • [WPF]XAML中使用IMultiValueConverter实现Command的多参数传参
    对ICommand进行多参数传参问题如何对ICommand传入多个参数?背景最近在做一个WPF的开发,有多个相近的功能写了不同的Command,因为要对应不同的对象。因为是CtrlCV,显得代码有点冗赘不够优雅,但是ICommand又只能接受一个参数。思路使用MultiBinding,对CommandParameter进行绑定,然后......
  • [JOI2018]定期券 (Commuter Pass)
    \(\mathtt{TAG}\):最短路,DP,拓扑排序题意给定一个\(n\)个点\(m\)条边的无向图,边有边权。给定两对点\(s_1,t_1\)和\(s_2,t_2\)。你可以选定\(s_1\)到\(t_1\)的一条最短路径,使得这些边的边权变为\(0\),要求操作之后\(s_2\)到\(t_2\)的最短路长度最小。First.\(......