建立基因组索引
调用bowtie2,因为bowtie2支持插入缺失。

bismark_genome_preparation--path_to_bowtie /usr/local/bowtie/<path_to_genome_folder> 参考基因路径--verbose 输出log信息

比对

bismark --bowtie2 -N 0      允许错配数-L 20    seed大小--quiet --un      过滤多处匹配reads--ambiguous   多处匹配reads信息独立记录--sam   输出格式为sam-o  输出目录<genome_folder>-1 read_1.fq-2 read_2.fq

去重复
duplicate_bismark -p <input.sam>

提取甲基化位点

bismark_methylation_extractor-P pair-end--comprehensive    输出CHG CHH CpG的甲基化信息--no-overlap --bedGraph   输出bedGraph文件--counts 每个C上甲基化reads和非甲基化reads的数目--buffer_size 20G--report  一个甲基化summay--cytosine_report  报道全基因组所有CpG--genome_folder  <path_to_reference_genome>input.sam-o output_dir


解读甲基化信息
bismark2report --dir <output_dir> --alignment_repot <report_path>
针对所有样本进行汇总.
bismark2summary sample_bismark_bt2.bam

结果解读
CpG:甲基化C下游是个G碱基。
CHH:甲基化C下游的2个碱基都是H(A、C、T)。
CHG:甲基化的C下游的2个碱基是H和G。

col1 : 比对上的序列ID
col2 : 基因组正负链:+ -
col3 : 染色体编号
col4 : 染色体位置
col4 : 甲基化C的状态(XxHhZzUu)X 代表CHG中甲基化的C
x  代笔CHG中非甲基化的C
H 代表CHH中甲基化的C
h  代表CHH中非甲基化的C
Z  代表CpG中甲基化的C
z  代表CpG中非甲基化的C
U 代表其他情况的甲基化C(CN或者CHN)
u  代表其他情况的非甲基化C (CN或者CHN)


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

col1 : position
col2 : 甲基化count
col3 : 非甲基化count
col4 : 甲基化率(beta)
col5 : coverage

col1 : 染色体编号
col2 : 染色体起始位置
col3 : 染色体终止位置
col4 : 甲基化率 (5/5+6)
col5 : 甲基化数目
col6 : 非甲基化数目

col1 : 染色体编号
col2 : 染色体起始位置
col3 : 正负链信息
col4 : 甲基化碱基数目
col5 : 非甲基化碱基数目
col6 : CG
col7 : CG背景(CG+一个碱基)

DNA甲基化比对:Bismark相关推荐

  1. DNA甲基化与表观遗传学数据挖掘与分析学习会 (10月26-27 上海)

    两天一夜高强度训练,理论与实战相结合 赠送高清视频供学员复习使用 会务背景 表观遗传学调控多种生命活动及疾病发生发展,已成为近年来生命科学领域的研究热点,促使了生物学多领域研究的突破性进展.随着实验和 ...

  2. DNA甲基化数据分析专题

    欢迎关注"生信修炼手册"! DNA 甲基化作为重要的表观遗传学的标记,能够调控基因的表达,在生长发育和疾病相关研究领域都有着重要意义.测定甲基化的手段有很多,芯片作为一种成熟的手段 ...

  3. 表观遗传学及DNA甲基化数据分析

    关于举办"表观遗传学及DNA甲基化数据分析"专题班通知 各有关单位: 表观遗传学是研究基因的核苷酸序列不发生改变的情况下,基因表达的可遗传的变化的一门遗传学分支学科.随着实验技术的 ...

  4. 易基因|DNA甲基化揭示肌痛性脑脊髓炎/慢性疲劳综合征在复发和恢复周期中的动态表观变化

    大家好,这里是专注表观组学十余年,领跑多组学科研服务的易基因. 2022年10月06日,<Int J Mol Sci>杂志发表题为"Dynamic Epigenetic Chan ...

  5. Nat. Mach. Intell. | 基于深度学习预测DNA甲基化位点

    研究人员开发了一种预测DNA甲基化位点的机器学习算法可以帮助识别致病机制.该论文2020年8月3日发表在"Nature Machine Intelligence"上. 研究人员通过 ...

  6. Microbiome | 中科院张惠明团队揭示RNA介导的DNA甲基化影响植物根部微生物群落!...

    研究背景 植物在根际中拥有多种土壤微生物,有机化合物从根部释放到土壤中,从而为土壤微生物提供了营养丰富的环境.根际微生物能够以不同的方式影响植物,例如通过产生刺激植物生长的植物激素或引起植物病害症状的 ...

  7. DNA甲基化检测服务

    DNA甲基化检测服务 DNA甲基化是最早发现的基因表观修饰方式之一,真核生物中的甲基化仅发生于胞嘧啶,即在DNA甲基化转移酶(DNMTs)的作用下使CpG二核苷酸5'-端的胞嘧啶转变为5'-甲基胞嘧啶 ...

  8. 国家生物信息中心在核酸研究发表单细胞DNA甲基化数据库—scMethBank

    DNA甲基化是表观遗传研究的一个重要层面,且与发育.衰老和疾病的发生发展密切相关.为了更好地利用已公开的海量甲基化数据,此前,中国科学院北京基因组研究所国家基因组科学数据中心(以下简称基因组数据中心) ...

  9. 易基因项目文章|WGBS+RNA-seq揭示PM2.5引起男性生殖障碍的DNA甲基化调控机制

    大家好,这里是专注表观组学十余年,领跑多组学科研服务的易基因. 2022年09月24日,陆军医科大学(第三军医大学)预防医学院毒理学研究所张中豪博士等为第一作者,曹佳教授.刘晋祎教授为通讯作者以&qu ...

  10. 易基因 | 精准医学:DNA甲基化图谱在发现和精确诊断神经肿瘤领域的应用

    近年来,中枢神经系统(CNS)肿瘤的分类变得更加客观.更依托生物学认知.虽然过去的分子诊断包括特定突变.拷贝数变化或基因融合,但DNA甲基化测序的发展显著提高了诊断精度,增加了可靠性,并为发现新的肿瘤 ...

最新文章

  1. JavaScript document write String content writes nothing o
  2. mysql猎豹_猎豹网校MySQL数据库
  3. 鸟哥的Linux私房菜(基础篇)- 第二十六章、Linux 核心编译与管理
  4. java中的jpa_JPA教程–在Java SE环境中设置JPA
  5. mapreduce排序算法_MapReduce算法–二级排序
  6. 声卡loopback有什么用_声卡购买须知,别买来又退整麻烦事的
  7. JSP自定义标签开发步骤
  8. JavaScript 教程
  9. ctf-希尔伯特曲线隐写
  10. 2021会员运营痛点分析
  11. VTN线下体验店 汇聚全球高端品牌 打造非凡购物体验
  12. mysql支持的时区列表_5.1.14 MySQL服务器时区支持
  13. android切换夜间模式吗,Android切换夜间模式
  14. 基于组态王的过程控制——以巧克力牛奶的生产为例
  15. hive分区表之insert overwrite 注意事项
  16. 【vue打包】线上部署报错net::ERR_ABORTED 404 (Not Found)
  17. springboot 简单的扫码登录 demo
  18. PostgreSQL数据库查询——scan.l分析
  19. 迅为2K1000龙芯开发板PMON 开发
  20. python 字符串去除转义

热门文章

  1. 模拟QQ自动隐藏窗体
  2. 6个 Python 办公黑科技,工作效率提升100倍!(附代码)
  3. 房产管理系统系统与数据中心对接要遵循的标准和规范
  4. AJP:有和没有内化性精神障碍的受虐女孩情绪回路延迟成熟的差异性
  5. Cadence快速画原理图封装(适合复杂器件)
  6. 牛客NC23053月月查华华的手机 题解(指针优化)
  7. sql java web动态曲线图_使用Echarts实现动态曲线图表
  8. 2020年春运买票中对“候补购票”的理解和看法。
  9. 如何获取土豆网等在线视频FLV地址
  10. Unity3D引擎各大插件免费下载地址