ChAMP 分析甲基化芯片数据-差异分析下篇

对于甲基化芯片的差异分析,除了有探针水平的差异分析,还有差异甲基化区域DMR分析。

差异甲基化区域的示意图如下:

该图片来自Bumphunter的文献,图中绿色矩形代表的就是一个差异甲基化区域。

A 图代表的是甲基化位点的在cancernormal 两个group的分布, 每组包含8个生物学重复;B 图代表的是DMR 检测的原理,在正负0.1的两条红线是自定义的差异阈值,而黑色的线是根据CpG位点差异程度拟合出来的线,只有当差异超出了阈值,也就是说黑色的线在两条红色线定义的区域之外的时候,才认为是1个候选的DMR。 从图上来看,每个DMR可以看作是从红色线条定义的区域凸出来的部分,叫做bump。

ChAMP 中,通过champ.DMR函数进行差异甲基化区域分析

用法示例

myDMR <- champ.DMR()

champ.DMR 集成了3种差异甲基化区域分析算法:

  1. Bumphunter

  2. DMRcate

  3. ProbeLasso


默认使用的是Bumphunter算法,三种算法对应的参数会有所不同,具体可以查看函数的帮助文档。

DMR分析完成之后,为了控制假阳性率,都会有一定的过滤手段。在ChAMP中,通过下列两个条件对结果进行过滤

minProbes

DMR区域包括了许多的CpG位点,每个Cp位点对应1个探针,这个参数指定包含的探针的最少个数,默认为7,如果一个DMR覆盖的探针数目少于7个,则不会输出。

adjPvalDmr

这个参数就是我们最常用的校正之后的p值,默认参数为0.05。

如果你的操作系统支持图形界面,可以运行DMR.GUI()命令,在浏览器中交互式的查看结果

DMRtable

差异甲基化区域的染色体区域和p值等基本信息

DMR heatmap

DMR 甲基化水平分布图

每个DMR区间两组样本甲基化水平分布图

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

推荐阅读更多精彩内容