商务服务
bismark 识别甲基化位点
2024-11-19 03:07

在中,根据甲基化的所处的上下文环境,分成以下3类;

bismark 识别甲基化位点

  1. CpG

  2. CHG

  3. CHH

代表磷酸二酯键,指的是甲基化的C的下游是1个碱基。代表除了碱基之外的其他碱基,即中的任意一种,代表甲基化的C下游的2个碱基是和, 表示甲基化的C下游的两个碱基都是。

bismark 识别甲基化位点

比对完之后,会生成1个bam 文件。使用命令可以从bam 文件中识别到甲基化的C,命令如下

bismark_methylation_extractor —comprehensive test/test_data_bismark_bt2.bam

只有1个参数,这个bam 文件是比对生成的bam文件,每个样本一个bam文件。

默认情况下,软件会自动根据两个因素生成结果文件

  1. 甲基化的C的类型就是前面提到的, , 3种类型

  2. 比对情况包括比对到四条链上, , , 4种情况所以会生成 3 X 4 = 12 个文件,对于链特异性文库来说,会生成3 X 2 = 6 个文件,这6个文件内容是类似的,都是记录了甲基化的C的染色体位置。

选项的作用就是在生成最终文件时,只考虑3种甲基化类型,将所有的比对情况进行合并,这样最终只会生成3个文件.

CpG_context_test_data_bismark_bt2.txtCHG_context_test_data_bismark_bt2.txtCHH_context_test_data_bismark_bt2.txt

以为例,内容如下:

Bismark methylation extractor version v0.19.0SRR15024317_length=86   -       1       57798691        zSRR15024319_length=86    +       2       10166600        ZSRR15024331_length=86  +       11      77736289        ZSRR15024338_length=86  +       3       197272186       Z

共5列,第一列为比对上的序列ID,第二列为基因组的正负链信息,第三列为染色体编号,第四列染色体上的位置,第5列为甲基化的C的状态。

不同字母表示不同的甲基化C:

X 代表CHG中甲基化的Cx  代笔CHG中非甲基化的CH 代表CHH中甲基化的Ch  代表CHH中非甲基化的CZ  代表CpG中甲基化的Cz  代表CpG中非甲基化的CU 代表其他情况的甲基化C(CN或者CHN)u  代表其他情况的非甲基化C (CN或者CHN)

对于, 采用字母的大小写来表征甲基化状态;对于, 采用字母的大小写来表征甲基化状态;对于, 采用字母 的大小写来表征甲基化状态。

上面的文件是methylation calling 最直接的证据,但是对于甲基化水平的定量来说,缺少了相关信息。运行时,除了生成上述文件之外,还会有下列3个文件

test_data_bismark_bt2_splitting_report.txttest_data_bismark_bt2.M-bias.txttest_data_bismark_bt2.M-bias_R1.png

记录了该样本甲基化的汇总信息

Final Cytosine Methylation ReportTotal number of C’s analysed:    40348Total methylated C’s in CpG context:    1365Total methylated C’s in CHG context:    21Total methylated C’s in CHH context:    103Total C to T conversions in CpG context:    678Total C to T conversions in CHG context:    10076Total C to T conversions in CHH context:    28105C methylated in CpG context:    66.8%C methylated in CHG context:    0.2%C methylated in CHH context:    0.4%

定义了每一个甲基化位点的详细信息,就是我们定量常用的beta 值部分文件内容如下

CpG contextposition        count methylated        count unmethylated      % methylation   coverage1       42      13      76.36   552       31      9       77.50   40

bismark 识别甲基化位点

双坐标轴图,左侧的纵轴代表甲基化比例,右侧的纵轴代表甲基化的数目,横坐标代表测序读长。

    以上就是本篇文章【bismark 识别甲基化位点】的全部内容了,欢迎阅览 ! 文章地址:http://www78564.xrbh.cn/news/27555.html 
     文章      相关文章      动态      同类文章      热门文章      栏目首页      网站地图      返回首页 迅博思语移动站 http://www78564.xrbh.cn/mobile/ , 查看更多   
最新文章
迅雷发布2025年第二季度财报:总营收1.04亿美元 同比增长30.6%
8月14日,迅雷发布了截至2025年6月30日的第二季度未经审计财报。据财报显示,本季度迅雷总营收为1.04亿美元,同比增长30.6%。总
原创赵露思风波,带货成就:娱乐圈的两面人生
近期,娱乐圈炙手可热的焦点无疑是赵露思,她身上接连发生的事件引发了舆论的广泛关注与热烈讨论。事件的起点是赵露思与她的经纪
定做出口彩膜编织袋生产企业
彩膜编织袋,适合于包装危化品(危险品)饲料和食品添加剂,可双面彩印膜及柔性印刷。 适合 10kg、20kg、25kg、40kg、50kg等规格
乐山到北京天坛公园景点, 2025必看_实测攻略_秘籍曝光_省钱避坑指南,3天玩转核心价值,解锁诱人数字!
遇到交通绕路、门票贵、人多到爆?90%的人都做错了,别急。今天手把手教你从乐山出发3天玩转北京天坛公园,省下上千块还能避开人
代理DINALON B1X25 G30V04129
供应 PA6 DINALON B1X25 G30 W 生产商: Grupo Repol供应 PA6 DINALON B1X25 G30V0 生产商: Grupo Repol供应 PA6 DINALON B1X25 G
从常德到厦门曾厝埯,5天4夜怎么玩才不踩坑?(避坑指南+省钱攻略)
曾厝垵文创村:这里汇聚了众多文艺小店、咖啡馆和特色民宿,是体验厦门慢生活的好去处。胡里山炮台:作为历史遗迹这里不仅能感受到
太原市献血条例
太原市献血条例(2014年8月29日太原市第十三届人民代表大会常务委员会第二十三次会议通过2014年11月28日山西省第十二届人民代表大
智慧识破诈骗陷进,园区建行多方联动织密“防护网”
近日,建行园区斜塘支行突然接到一通紧急电话,客户A某语气焦急,称其在附近网点办理业务,但因名下银行卡账户交易受限,需开户
既实惠又好的普洱茶,选购指南,技巧让你轻松找到好茶
既实惠又好的普洱茶,选购指南,技巧让你轻松找到好茶你是不是也这样?买普洱茶不是被坑了就是花了冤枉钱,别急我懂你!许多人说
2003年公斤飞天茅台酒目前回收价免费鉴定
2003年公斤飞天茅台酒目前回收价免费鉴定2003年10月,出厂价提高23%,价格为268元,市场零售价在320左右;2004年2005年,出厂价