IBS和IBD的区别和计算方法介绍

2024-09-08 11:20

本文主要是介绍IBS和IBD的区别和计算方法介绍,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

大家好,我是邓飞。

今天介绍一下IBS和IBD的区别:

IBS(肠易激综合症)和IBD(炎症性肠病)是两种不同的消化系统疾病,主要区别如下:

    • IBS(Irritable Bowel Syndrome):是一种功能性肠道疾病,主要表现为腹痛、腹胀、腹泻或便秘,症状通常与饮食、压力和心理因素相关,没有明显的器质性病变。

    • IBD(Inflammatory Bowel Disease):是一组炎症性疾病,包括克罗恩病(Crohn's disease)和溃疡性结肠炎(Ulcerative Colitis),其特征是免疫系统的异常反应导致肠道慢性炎症,症状通常包括腹痛、腹泻、体重下降和疲劳。

错了,重新来。

IBS(Identical By State,状态同源):IBS指的是两个个体在某一基因座上拥有相同的等位基因,但这些等位基因不一定来源于同一祖先。换句话说,IBS是基于基因型相似性的度量,而不考虑等位基因的具体来源。说人话:IBS是状态同源,是根据基因型相似性进行计算的,不考虑等位基因的来源;

IBD(Identity By Descent,血缘同源:IBD指的是两个个体中共有的等位基因来源于同一祖先。这些等位基因在遗传过程中未经重组而直接传递给后代,因此能够反映个体间的亲缘关系。说人话:IBD是血缘同源,等位基因遗传过程中传递给后代,片段是不变的,需要考虑等位基因的来源。

其实,很久之前,我就写了一篇介绍IBS和IBD的文章:论体型决定性格以及逗比潜质的遗传性分析

文中说:根据SNP或者SSR计算的是IBS矩阵,是不全面的,其实,也可以计算IBD,下面有介绍。

系谱计算的是IBD,同胞,半同胞,表兄弟,都是有亲戚关系的,他们性格,长相相似,是有依据的,因为他们由IBD计算的血缘相似性较高。根据SNP或者SSR计算的是IBS,虽然两个人没有亲戚,但依旧有很高的相似性(比如长相,性格,体型,爱好等),这说明两者的IBS计算的相似性较高。

邓飞2013,公众号:育种数据分析之放飞自我论体型决定性格以及逗比潜质的遗传性分析

说了一大堆概念,代码在哪里???

如何计算IBD?

IBD计算,可以用系谱记录进行计算,计算亲缘关系A矩阵如何利用系谱计算近交系数和亲缘关系系数,怎么推断呢,比如0.25是半同胞,0.5是全同胞或者亲子关系,写到这里,我突然想到,人类的亲疏是根据IBD划分的,亲子关系和同父母的兄弟姐妹的IBD都是0.5,应该是人类中最亲的关系了。

如何计算IBS?

如何计算IBS矩阵呢,其实,Excel也可以计算,比如:ID1和ID2之间的IBS相关系数,比较每个SNP的交叉数,比如SNP1,分别是TT和TT,那么就是2;如果是TT和TA,就是1;如果是TT和AA,就是0。统计一下总和,然后计算一下比值,比如10个SNP总数是20,统计的总和为9,那么IBS为0.45。

把上面的Excel表格的SNP数据,转为plink格式,方法见链接:Excel的SNP数据如何变为plink格式,运行代码:

plink --file file --ibs-matrix

结果如下:

$ cat plink.mibs1 0.45 0.850.45 1 0.40.85 0.4 1

可以看到,ID1 VS ID2为0.45,ID1和ID3是0.85,和Excel计算的结果一样一样滴。

为啥SNP数据还可以计算IBD?

因为官方文档写了呀:https://www.cog-genomics.org/plink/1.9/ibd

PLINK中使用 PI_HAT 值来推定IBD的值。该方法基于**隐马尔科夫模型 hidden Markov model (HMM)**,通过矩估计(method-of-moments)来计算 IBD=1, 2或0 的概率。PLINK1.9中提供了`--genome`的选项,以计算 PI_HAT。

不是直接计算的IBD,而是通过IBS进行推断的,经验值也是0.25是半同胞,0.5是全同胞。

同样的道理,如果是SNP数据计算的G矩阵(如何构建G矩阵-基因组亲缘关系矩阵(Genomic relationships matrix))也应该算是IBD矩阵,和系谱计算的亲缘关系A矩阵一样都属于IBD矩阵。

最后,这个图是几个意思?

这个是表型相似,按道理来说是IBS相似,不可能是IBD矩阵,因为不可能存在亲缘关系。。。感觉这个图让我对我写的一大堆IBS和IBD的内容暴击了1000倍。

这篇关于IBS和IBD的区别和计算方法介绍的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1147969

相关文章

Before和BeforeClass的区别及说明

《Before和BeforeClass的区别及说明》:本文主要介绍Before和BeforeClass的区别及说明,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录Before和BeforeClass的区别一个简单的例子当运行这个测试类时总结Before和Befor

MybatisPlus service接口功能介绍

《MybatisPlusservice接口功能介绍》:本文主要介绍MybatisPlusservice接口功能介绍,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友... 目录Service接口基本用法进阶用法总结:Lambda方法Service接口基本用法MyBATisP

MySQL复杂SQL之多表联查/子查询详细介绍(最新整理)

《MySQL复杂SQL之多表联查/子查询详细介绍(最新整理)》掌握多表联查(INNERJOIN,LEFTJOIN,RIGHTJOIN,FULLJOIN)和子查询(标量、列、行、表子查询、相关/非相关、... 目录第一部分:多表联查 (JOIN Operations)1. 连接的类型 (JOIN Types)

Android学习总结之Java和kotlin区别超详细分析

《Android学习总结之Java和kotlin区别超详细分析》Java和Kotlin都是用于Android开发的编程语言,它们各自具有独特的特点和优势,:本文主要介绍Android学习总结之Ja... 目录一、空安全机制真题 1:Kotlin 如何解决 Java 的 NullPointerExceptio

java中BigDecimal里面的subtract函数介绍及实现方法

《java中BigDecimal里面的subtract函数介绍及实现方法》在Java中实现减法操作需要根据数据类型选择不同方法,主要分为数值型减法和字符串减法两种场景,本文给大家介绍java中BigD... 目录Java中BigDecimal里面的subtract函数的意思?一、数值型减法(高精度计算)1.

Pytorch介绍与安装过程

《Pytorch介绍与安装过程》PyTorch因其直观的设计、卓越的灵活性以及强大的动态计算图功能,迅速在学术界和工业界获得了广泛认可,成为当前深度学习研究和开发的主流工具之一,本文给大家介绍Pyto... 目录1、Pytorch介绍1.1、核心理念1.2、核心组件与功能1.3、适用场景与优势总结1.4、优

Java实现本地缓存的常用方案介绍

《Java实现本地缓存的常用方案介绍》本地缓存的代表技术主要有HashMap,GuavaCache,Caffeine和Encahche,这篇文章主要来和大家聊聊java利用这些技术分别实现本地缓存的方... 目录本地缓存实现方式HashMapConcurrentHashMapGuava CacheCaffe

Linux中的more 和 less区别对比分析

《Linux中的more和less区别对比分析》在Linux/Unix系统中,more和less都是用于分页查看文本文件的命令,但less是more的增强版,功能更强大,:本文主要介绍Linu... 目录1. 基础功能对比2. 常用操作对比less 的操作3. 实际使用示例4. 为什么推荐 less?5.

Spring Security介绍及配置实现代码

《SpringSecurity介绍及配置实现代码》SpringSecurity是一个功能强大的Java安全框架,它提供了全面的安全认证(Authentication)和授权(Authorizatio... 目录简介Spring Security配置配置实现代码简介Spring Security是一个功能强

Java 关键字transient与注解@Transient的区别用途解析

《Java关键字transient与注解@Transient的区别用途解析》在Java中,transient是一个关键字,用于声明一个字段不会被序列化,这篇文章给大家介绍了Java关键字transi... 在Java中,transient 是一个关键字,用于声明一个字段不会被序列化。当一个对象被序列化时,被