ChengyangWang的个人博客分享 http://blog.sciencenet.cn/u/ChengyangWang

博文

画泡泡图,发CNS

已有 3687 次阅读 2018-3-28 14:49 |个人分类:生物paper|系统分类:科普集锦

 本文转载自嘉因微信公众号,已获得授权。查看最新文章,敬请关注嘉因,微信ID:rainbow-genome

作者:小丫   来源: 嘉因

小伙伴们看了《把Cell爆文的思路、方法用到自己的研究中》都想知道:找到关键转录因子的这个图是怎么画的?

Figure S5. Motif Dynamics, Related to Figure 3 (A) Enrichment scores of transcription factor binding motifs for the DHSs in human and mouse embryos.

出自Gao, L., Wu, K., Liu, Z., Yao, X., Yuan, S., & Tao, W., et al. (2018). Chromatin accessibility landscape in human early embryos and its association with evolution.Cell. 详见《把Cell爆文的思路、方法用到自己的研究中

 

彩色版

 

 

Figure 4 d. TF motif analysis of different regions in the human genome CGI in DMVs, non-CGI in DMVs, CGI in all promoters, and non-CGI in all promoters. The scale of the circles represents motif enrichment. Red represents homeobox TFs.

出自Li, Y., Zheng, H., Wang, Q., Zhou, C., Wei, L., Liu, X., Zhang, W., Zhang, Y., Du, Z., Wang, X. and Xie, W., 2018. Genome-wide analyses reveal a role of Polycomb in promoting hypomethylation of DNA methylation valleys. Genome Biology, 19(1), p.18. 详情见《任意两组RNA-seq变身,国自然得AAA

 

还有用颜色代表表达量的:

 

Figure 4 Identification of candidate regulators for early development. a, TF motifs identified from distal ATAC-seq peaks. Only TFs expressed at least at one stage (FPKM ≥ 5) and motif enrichment P value <1 × 10−10 at least at one stage were included.

出自Wu, J., Huang, B., Chen, H., Yin, Q., Liu, Y., Xiang, Y., ... & Li, W. (2016). The landscape of accessible chromatin in mammalian preimplantation embryos. Nature, 534(7609), 652.详见《表观遗传是怎样遗传的?| 颉伟讲怎样用ATAC-seq研究表观遗传机制的视频

 

 

 

到底怎么画呢?方法来啦!

 

上面三篇paper用的都是HOMER,又一款motif分析利器!

The script findMotifsGenome.pl in HOMER software was used to analyze transcription factors binding motifs enrichment in DHSs (Heinz et al., 2010).

 

http://homer.ucsd.edu/homer/ngs/peakMotifs.html

 

用HOMER的findMotifsGenome.pl算出P-value,然后跟RNA-seq的表达值FPKM一起作为输入,在R里面画图。

 

熟悉clusterProfiler的小伙伴儿,一定觉得这些图跟下面这个图很像,详见《最靠谱的富集分析,超炫的展示方式,TCGA也是他的粉丝》:

 

 

击鼓传花,3,2,1,开始:

  • 怎样用HOMER找motif?@生信技能树

  • 怎样用HOMER算出的P-value画出CNS级别的泡泡图?@Y叔

 


想用
ChIP-seqATAC-seq研究感兴趣的基因?想整合ChIP-seq、ATAC-seq、eCLIP/RIP-seq、RNA-seq数据寻找线索?找嘉因生物吧!从实验、测序,到多种数据整合分析,为您一站式解决。(点击文中蓝字了解详情)




嘉因公众号定位:客户共性问题解答,生信学习资源导航,高通量实验导购 | 为您提供高通量实验-测序-分析-验证一站式解决方案


 


电话:021-61539657

Email:marketing@rainbow-genome.com
地址:上海市杨浦区国伟路135号10号楼305室





https://wap.sciencenet.cn/blog-3372875-1106111.html

上一篇:我画的图只能发Plos One?Genome Research的图长啥样?
下一篇:神技能!转录调控必读
收藏 IP: 124.77.56.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-19 23:36

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部