一、Independent significant SNPs 和 Lead SNPs 定义

Independent significant SNPs: P 值小于 5e-8 的 SNP 位点,且 SNP 位点之间的 r2<0.6(注意,这里的 r2 指的是 LD 的 r2 值,不是相关性哦,别误解了,0.6 也是我设定的,不同文献给定的 r2 值不一样,主流设定是 0.6 )

Lead SNPs: Independent significant SNPs 中 r2 值小于 0.1 的 SNP(同样的,这里的 0.1 也是我设定的,不同文献阈值不一样,主流设定是 0.1 )

如果文字描述还不清楚的话,请看图片版的:

在 locus1 中,SNP1, SNP2, SNP3都达到了显著水平(P< 5e-8),SNP1 和 SNP2 的 r2=0.8, SNP2 和 SNP3 的 r2=0.4, SNP1 和 SNP3 的 r2=0.15。

假定我们设定 Independent significant SNPs 的阈值为:P<5e-8 且 r2<0.6,那么locus1 只有 SNP1 和 SNP3 是 Independent significant SNPs。

假定我们设定 Lead SNPs 的阈值为:r2 < 0.1,那么只有SNP3是Lead SNPs

看明白其中的逻辑了吗,也就是说,Lead SNPs 的数量一定是小于等于 Independent significant SNPs 的数量,因为 Lead SNPs 是在 Independent significant SNPs 的基础上进行了二次clump。

以上,就是 Independent significant SNPs 和 Lead SNPs 的定义。

二、FUMA 鉴定 Independent SNPs 和 Lead SNPs

理解了 Independent significant SNPs 和 Lead SNPs 的定义后,我们就可以借助工具鉴定 Independent significant SNPs 和 Lead SNPs。

目前挺多工具支持 Independent significant SNPs 和 Lead SNPs 的识别,比如PLINK。

我来推荐一个比较友好的网页工具 FUMA(https://fuma.ctglab.nl/snp2gene) 。

进入网址:https://fuma.ctglab.nl/snp2gene 后,
只需要输入 GWAS summary 文件,就可以得到 Independent significant SNPs 和 Lead SNPs。

步骤一:上传 GWAS summary文件,把 GWAS summary 文件的表头信息(比如染色体,位置,rs号等)补充在下图中

步骤二:设定 Independent significant SNPs 和 Lead SNPs 的阈值

步骤三:点击 submit job, 然后就可以等结果了

三、输出文件

以下是 FUMA 网站给出的 Independent significant SNPs 和 Lead SNPs 结果:

使用 FUMA 鉴定 Independent SNPs 和 Lead SNPs相关推荐

  1. 本周最新文献速递20220423

    本周最新文献速递20220423 一.精细解读文献 一 文献题目: Nonlinear control of transcription through enhancer–promoter inter ...

  2. 玉米关联群体:155、368、527 自交系群体简介(Yan Jianbing,2010、2013)

    本文将以此介绍 3 个关联分析资源群体,群体大小分别为 155.368.527.155 和 527 两个群体交集为 141,使用的 SNP 芯片一致,可以认为 527 是 155 群体的拓展.368 ...

  3. 高通量测序技术的原理及各平台优势和实践应用的分析

    高通量测序技术的原理及各平台优势和实践应用的分析 2020.9.01 2060 随着人类基因组计划(human genome project )在2003年顺利完成,基因组测序技术取得了长足的进步,这 ...

  4. Identifying genetic relatives without compromising privacy

    文章 He D, Furlotte N A, Hormozdiari F, et al. Identifying genetic relatives without compromising priv ...

  5. 全基因组重测序揭示了野生大豆的局部适应和分化的特征

    近日,凌恩生物客户南京农业大学在<Evolutionary Applications>发表题为"Whole-genome resequencing reveals signatu ...

  6. U-Boot 之一 零基础编译 U-Boot 过程详解 及 编译后的使用说明

      在之前的博文 Linux 之八 完整嵌入式 Linux 环境介绍及搭建过程详解 中我们说了要一步步搭建整个嵌入式 Linux 运行环境,今天就开始编译 U-Boot.我所使用的硬件平台及整个要搭建 ...

  7. STM32MP157开发板调试笔记

    作者:zzssdd2 E-mail:zzssdd2@foxmail.com 一.前言 ​ 最近在ODYSSEY-STM32MP157C板子上移植arm-trusted-firmware.u-boot以 ...

  8. GWAS相关名词解释及基础知识储备[长期更新]

    转载自GWAS综述(生信文献阅读俱乐部精选) 相关名词解释 Genome-wide association studies(关联分析):扫描遗传标记,通常是单核苷酸多态性(SNPs),使用统计学相关的 ...

  9. uboot usb设备驱动

    文章目录 DTS中定义 usb driver bind fastboot fastboot命令处理 fastboot触发 介绍usb device设备驱动 uboot在init_dm中会扫描dtb和代 ...

  10. 单核苷酸多态性(single nucleotide polymorphism,SNP)

    单核苷酸多态性(single nucleotide polymorphism,SNP),主要是指在基因组水平上由单个核苷酸的变异所引起的DNA序列多态性.它是人类可遗传的变异中最常见的一种.占所有已知 ...

最新文章

  1. 含有min()函数的栈,各种操作时间复杂度为O(1)
  2. iOS开发 Xcode8中遇到的问题及改动
  3. windows10下配置环境变量
  4. 为什么(12)式,km不能直接相乘?而要让域k先乘一个代数A里面的单位元,再作用在群M上呢?...
  5. php与java安全之争
  6. 服务器如何返回429状态,Instagram远程服务器返回一个错误:(429)UNKNOWN STATUS CODE...
  7. nginx设置 二级域名 指定端口
  8. QUIC - 低时延互联网传输层协议
  9. 爱奇艺开源轻量级插件化方案 Neptune
  10. oracle 01775,Oracle出现ORA-01775: 同义词的问题
  11. Corosync与Pacemaker实现nfs作为共享存储的MySQL高可用
  12. WEB应用中的信息泄漏以及攻击方法
  13. 【CentOS 7笔记43】,防火墙和iptables filter表#
  14. iis 中php的mime类型,IIS添加MIME扩展类型及常用的MIME类型列表
  15. 用好binutils之gporf
  16. Spring拦截器和过滤器的区别
  17. 单元测试 CI/CD
  18. 计算机在英语教学中的应用课题,浅谈信息技术在英语教学中的应用
  19. 百度网盘提速法,不用插件,不用安装脚本
  20. 106短信平台多少钱一条比较合理?

热门文章

  1. 风靡健身圈的生酮饮食居然有这么多好处,受教了
  2. MuMu模拟器的安装
  3. SAM-BA 2.14 reconfiguration----DIY你的sam-ba
  4. excel 取消合并单元格
  5. 关于编译错误ambiguous call of overridden pre R14 auto-imported BIF get/1
  6. thuwc2019滚粗记
  7. 接口测试 requests的身份认证方式
  8. python中if缩进规则_Python if条件判断
  9. if的用法及常见问题
  10. CSS:border的属性