bsseq是如何进行差异甲基化分析

发布时间:2021-12-28 16:26:32 作者:柒染
来源:亿速云 阅读:165

bsseq是如何进行差异甲基化分析,很多新手对此不是很清楚,为了帮助大家解决这个难题,下面小编将为大家详细讲解,有这方面需求的人可以来学习下,希望你能有所收获。

bsseq 主要用来分析WGBS的数据, 安装过程如下

source(“http://bioconductor.org/biocLite.R“)
biocLite(“bsseq”)

bsseq的分析主要包括以下4步:

  1. 读取原始数据

  2. BSmooth

  3. t-test检验

  4. DMR

1. 读取原始数据

bsseq要求的原始数据格式如下:

bsseq是如何进行差异甲基化分析

共6列数据,制表符分隔,每一行代表一个甲基化位点,前5列很好理解,描述甲基化位点的染色体位置和类别,默认情况下bbseq用于分析CpG类型的甲基化位点。当然其他类型的数据,比如CHG, CHH也支持,但是需要调整参数。Cov代表覆盖到这个位点的reads数,M代表其中发生了甲基化的reads数目。

每个样本一个这样的原始数据,用来表示该样本methylation calling的结果,这样的数据我们从bismark的结果中也可以得到。当原始数据准备好之后,首选需要读取所有样本的原始数据,然后导入到R中,生成一个bbseq定义的对象。在bbseq安装的路径下,提供了一个名为get_BS.chr22.R的脚本,展示了如何从读取所有样本原始数据的过程。

代码如下

bsseq是如何进行差异甲基化分析

这里以mc_imr90_r1_22mc_imr90_r2_22两个样本的原始数据为例,详细展示了读取过程。我们只需要根据自己的数据,适当修改上述代码就可以了。主要注意sampleNamespData数据就可以了。

2 Smooth

已测试数据BS.chr22为例,smooth的过程如下

bsseq是如何进行差异甲基化分析

在实际分析中,由于甲基化位点很多,所以这一步时间特别久,为了提高速度,可以添加mc.cores参数,这个参数指定了CPU个数,用于并行执行。

BS.chr22.1 <- BSmooth(BS.chr22, mc.cores = 2, verbose = TRUE)

3. T-test

在分析之前,有必要过滤掉覆盖度较低的甲基化位点。通常保留在所有样本中覆盖度大于2的甲基化位点,但是也可以修改这个条件。过滤之后,直接通过BSmooth.tstat进行分析

下面的代码基于bsseqData包中的数据,这个数据包含了6个样本,分为normalcancer两组

bsseq是如何进行差异甲基化分析

group1 指定属于treatment组的样本,group2指定属于control组的样本。

4. DMR

通过dmrFinder 函数进行差异甲基化分析, 代码如下:

bsseq是如何进行差异甲基化分析

cutoff 指定DMR的阈值,这个阈值根据t-test的结果进行调整。subset对差异甲基化的结果进行筛选,筛选包含甲基化位点个数大于3而且meanDiff 大于0.1的甲基化区域。

看完上述内容是否对您有帮助呢?如果还想对相关知识有进一步的了解或阅读更多相关文章,请关注亿速云行业资讯频道,感谢您对亿速云的支持。

推荐阅读:
  1. 如何进行不同系统里同一Customizing activity的显示差异分析
  2. 如何进行企业使用ERP和MES集成应用差异分析

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

上一篇:Borg和Kubernetes有什么区别

下一篇:methylKit是进行差异甲基化分析

相关阅读

您好,登录后才能下订单哦!

密码登录
登录注册
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》