bivalent domains 如何计算 - 做专业的生信问答网站

登录

该问题已被锁定！

2: 关注

1124: 浏览

bivalent domains 如何计算

Chip_Seq

还有下面这种不同状态的分析，感觉像ChromHMM 分析的又好像不是。

图片是来自同济大学高绍荣团队发的文章,https://www.sciencedirect.com/science/article/pii/S1934590922002508?via%3Dihub

试着发邮件要了下画图数据分析过程和代码也没有回复

文章说计算来自一篇cell 文章https://www.sciencedirect.com/science/article/pii/S0092867406003801?via%3Dihub

挺想做上面这些分析的，哪位大佬能搞个详细教程呀。

好问题 0 评论收藏举报

2 回答

孟浩巍超级管理员用户来自于: 北京市
2023-06-02 17:36

这个所谓的bivalent domain指的就是一段染色质区域，处于即可能开放又可能抑制的状态。

比如你上面提到的，H3K4me3这个就是典型的active的 histone marker；

H3K27me3就是repressive的 histone marker；

然后如果这两个同时存在，就可以定义为bivalent domain。

看你上面的截图，就应该是先把基因组分区域，分区域的方法通常是分bin，然后你就可以计算这一个bin区间里的H3K4me3和H3K27me3的信号强度，两个都存在的时候就可以定义。

我看你上面还提到了chromHMM，用这个也可以，就是用chromHMM分出来H3K4me3和H3K27me3同时存在的那个类型就定义为bivalent domain即可。这个和直接分bin自己算coverage类似，没太大差别。

赞同 1 0评论

chatGPT机器人机器人用户来自于: IANA
2023-06-01 16:27

Bivalent domains是指同时存在H3K4me3和H3K27me3修饰的染色质区域，它们在基因调控中发挥重要的作用。下面简单介绍一些常用的计算方法。一、ChromHMM ChromHMM是一种常用的计算bivalent domains的工具。其基本思路是将染色质区域划分为不同的状态（如active promoter, enhancer, bivalent promoter等），再根据H3K4me3和H3K27me3的分布情况确定bivalent domains的位置。具体的计算流程可以参考ChromHMM官方文档。二、MACS2 MACS2是一种常用的ChIP-seq分析工具，它可以用来计算H3K4me3和H3K27me3的峰值。根据峰值的分布情况，可以确定bivalent domains的位置。具体的计算流程可以参考MACS2官方文档。三、HOMER HOMER是一种常用的染色质分析工具集，其中包含了多种计算bivalent domains的工具。比如，可以使用annotatePeaks.pl命令将H3K4me3和H3K27me3的信号合并起来，然后根据信号的分布情况确定bivalent domains的位置。具体的计算流程可以参考HOMER官方文档。以上是一些常用的计算bivalent domains的方法，具体的选择可以根据实验设计和数据分析的需求来确定。

赞同 0 0评论

关于作者

: tidyboy 注册会员
这家伙很懒，还没有设置简介

0: 回答

0: 文章

1: 问题

问题动态

发布时间: 2023-06-01 16:17

更新时间: 2023-06-02 17:36

关注人数: 2 人关注

相关问题

x <- x[keep.exprs,, keep.lib.sizes=FALSE] 请问一下这条命令该如何解读呢？: 604 浏览 1 关注 0 回答 2 评论

知道一批基因的具体位置，如何批量的从基因组中取出基因序列？: 978 浏览 3 关注 2 回答 0 评论

FPKM计算问题: 646 浏览 2 关注 1 回答 0 评论

GAPIT包导出的GWAS结果如何添加新的阈值线？以及GAPIT的结果文件中的nobs、H&B.P.Value、Effect分别是什么意思？: 1202 浏览 2 关注 1 回答 0 评论

大咖们好！如何像群主一样优秀自带关注: 948 浏览 4 关注 2 回答 0 评论

进行转录组数据分析时，进行cuffdiff后的输出文件gene_exp.diff中，一个基因出现了两个不同的表达量数据，应该如何处理？: 636 浏览 2 关注 1 回答 0 评论

ATAC-seq数据多样本call peak如何合并？: 694 浏览 2 关注 1 回答 0 评论

如何从NCBI上分别下载所有的RNA病毒和DNA病毒的序列: 1072 浏览 2 关注 2 回答 0 评论

如何下载指定文献的原始数据？？: 838 浏览 2 关注 1 回答 0 评论

R软件里Bray-Curtis指数计算: 491 浏览 1 关注 0 回答 0 评论

推荐内容

问老铁们，我想问下，如果homer找到了motif1，那我有可能反过来知道哪些peak含有motif1么。: 650 浏览 0 关注 1 回答 0 评论

问 ChIP-Seq deeptools热图信号比较为啥一般不做统计检验比较？: 1297 浏览 2 关注 2 回答 0 评论

问 CUT&TAG 的GC含量异常: 609 浏览 2 关注 1 回答 0 评论

问 ChIP-seq的IP的比对率低的可能原因: 1232 浏览 2 关注 1 回答 0 评论

问 annotatePeak peak基因组注释: 1050 浏览 2 关注 1 回答 0 评论

问 chromosome名称转换: 795 浏览 2 关注 1 回答 0 评论

问 aspera SRA 数据批量下载: 1158 浏览 2 关注 1 回答 0 评论

问 MEME motif analysis 软件的安装: 1056 浏览 2 关注 2 回答 0 评论

问是不是单细胞的chip-seq，cut run这些call peak是不需要对照的，我看代码里没有-c.所以Rna-seq的峰图也是不需要control来call peak的吗？: 780 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024