科学网

 找回密码
  注册

tag 标签: motif

相关帖子

版块 作者 回复/查看 最后发表

没有相关内容

相关日志

[转载]ggseqlogo:一个通用序列logos绘制R包
chinapubmed 2019-6-2 10:24
序列logos已经成为研究基因组中基本序列模式的重要可视化方法。尽管如此,仍然缺乏提供这些可视化通常需要的多功能性软件包。 ggseqlogo是一个以ggplot2包为基础的R包,旨在解决这个问题。ggseqlogo以一种高度可定制的方式提供了出版质量DNA、RNA和蛋白质序列logos,具有包括多logo绘制,定性和定量配色方案,logos注释和与其他图整合的特征。该包可以直观地使用,能够无缝整合到R分析流程中。 连续配色方案 logo注释 多图整合 下载: https://cran.r-project.org/web/packages/ggseqlogo 时间:20170720 参考: ggseqlogo: a versatile R package for drawing sequence logos
个人分类: 软件|3431 次阅读|0 个评论
找到了motif,怎样展示结果?
ChengyangWang 2018-1-31 15:55
本文转载自嘉因微信公众号,已获得授权。查看最新文章,敬请关注嘉因,微信ID:rainbow-genome 作者:小丫 来源: 嘉因 motif系列问题答疑: 去哪找motif? 史上最全物种转录因子、motif数据库footprintDB 这段DNA上有我关心的motif吗? 点鼠标就能找启动子区的motif | meme-FIMO motif scan结果怎样看? 互补链上的motif有意义吗? motif结果怎样展示到文章里? 今天回答这个问题 依稀记得Y叔介绍过一个工具,能用关键字搜到research paper里的图,我当初还收藏了,可是收藏的文章太多,找不到哇!翻biobabble历史消息,用了几个词搜也没搜到。想起有一次Y叔说在公众号后台回复“目录”,就能看所有文章的目录。其中有个“神器系列”看起来像是汇总工具的,其中第二个就是: 在这个看脸(哦不对,看图)的时代,做科研的你怎么能错过这个帮你找一张漂亮脸蛋的神器! 经常浏览优质公众号,平时看的帖子当时未必需要,需要的时候能想起来去哪找。嘉因的技术贴,最近我整理了一些,放在左下角阅读原文的目录里,或者加最熟悉嘉因技术贴的嘉因表观小助手,她的微信ID:epigenomics http://www.viziometrics.org/ 搜motif,没有我喜欢的 搜ChIP-seq motif 这样展示就很好 它的原理是给图加标签,用户还能帮忙校对。 2016年发表的paper,code都给你了。
个人分类: 转录调控|10671 次阅读|0 个评论
[转载]点鼠标就能找启动子区的motif | meme-FIMO
ChengyangWang 2018-1-23 10:05
本文转载自嘉因微信公众号,已获得授权。查看最新文章,敬请关注嘉因,微信ID:rainbow-genome 作者:小丫 来源: 嘉因 《 转录因子调控了谁?挖全天下所有高通量实验数据,只为您解决这一个问题 》给出了从基因组水平找靶基因的策略和方法: Plan A:RNA-seq。直接调控?间接调控?傻傻分不清 Plan B:ChIP-seq。最直接,最有效, 与RNA-seq整合分析更准确 Plan C:ATAC-seq + motif分析。没有ChIP级别抗体,也能准确找到靶基因 Plan D:基于motif预测,这是个垫底儿的 《 哪个蛋白质调控我感兴趣的基因?怎样筛选?基于分析或实验的可行方案V2.1 》给出了研究某个基因受谁调控的策略和方法: Plan A:基于大量ChIP-seq公共数据挖掘 Plan B:motif分析预测 Plan C:ATAC-seq结合motif分析 这些策略适用于所有物种的转录调控研究,但小伙伴们一结合自己的课题,疑问就多了起来。我们一个一个解决。 今天解决这个问题: 一段DNA序列上有没有某个转录因子的motif? 话说找motif的工具有很多哇!不好讲哪个好哪个差,谁能找出你心仪的motif谁就是好猫! 方法一 《 转录因子调控了谁? 》里的Plan D教大家用JASPAR + UCSC Table Browser查某个motif在全基因组范围内分布在哪些基因附近。好烦,我不要全基因组范围,我只要看这一个基因,见方法二和三。 方法二 《 互补链上的motif有意义吗? 》里的补充技能:用JASPAR2016在线工具,找一段DNA序列上是否存在指定的motif,输入序列长度受限。JASPAR2018提供R包TFBSTools做这件事,要会用R,不接地气。 方法三 meme是个强大的工具箱(此meme非 彼meme ),以后会以问题为导向,分多篇介绍里面的工具。今天用 FIMO ,'Find Individual Motif Occurrences', 判断一段DNA上是否存在某个motif 。还能同时上传多个motif和多条DNA序列,批量判断多个基因上游是否存在某几个motif。 首先,去JASPAR2018下载MEME格式的motif文件。 http://jaspar.genereg.net/matrix/MA0106.3/ 然后,打开FIMO,http://meme-suite.org/tools/fimo 1. Upload刚刚从JASPAR下载好的meme格式的motif文件; 2. 把感兴趣的DNA序列以fasta格式粘贴到Input the sequences区域,或者upload序列文件; 3. Start search,然后把地址加入收藏栏; 4. 需要运行一段时间,即使浏览器被关掉了,还能从收藏栏里打开结果。 运行结束,页面会变成这样: 两种接地气的查看方式,和一种高大上的展示方式: 一、以表格的形式查看哪个基因上有motif。 点开FIMO HTML.output,显示P0.0001的Top1000个位点。Alt ID就是你upload的fasta文件里“”右边的字,建议写基因名。 二、图形化展示输入序列的哪个位置有这个motif,展示效果类似于方法一。 1. 鼠标右键点击FIMO GFF output,复制gff文件地址; 2. 打开网址http://genome.ucsc.edu/cgi-bin/hgCustom,往UCSC genome browser里添加一个新的track,在URLs的框里粘贴gff文件地址,submit。 三、高大上的展示方式。 看这篇: 如何画类似MEME的注释序列 方法四、五、六、七。。。 看这篇: 研究单个基因的生物信息学分析工具(大全)
个人分类: 转录调控|8445 次阅读|0 个评论
[转载]互补链上的motif有意义吗?
ChengyangWang 2018-1-17 16:31
本文转载自嘉因微信公众号,已获得授权。查看最新文章,敬请关注嘉因,微信ID:rainbow-genome 作者:小丫 来源: 嘉因 有位小伙伴儿问:JASPAR里预测转录因子结合的strand-1是不是没有意义? 搜到好多论坛上都说strand-1无意义,只看strand+1就好。是这样吗? 下图清晰的描述了两种情形: a) 基因上游存在某一TF的motif,推测该TF可能结合到这里,调控该基因转录。 b) 基因上游互补链上存在TF的motif,这个TF还能调控这个基因的转录吗? 认为strand-1无意义,估计是考虑到如果motif的方向反了,TF的方向也跟着反了,调控的不就是相反方向的基因了? 事实上是这样的吗? 高中物理老师告诉我们:你想的都是错的,只有算过才是对的。 2016年发表在BMC Genomics上的文章已经算过了。 The orientation of transcription factor binding site motifs in gene promoter regions: does it matter? 开篇那个模式图就是这篇文章的Figure 1。感谢羊咩咩同学搜到这篇文章。 作者分析了拟南芥293个非回文序列(non-palindromic)转录因子结合位点、10个核心启动子motif的方向对基因表达调控的影响。 结论:方向不重要,结合才重要。 Conclusions : Our results suggest that for the motifs considered here, either no specific orientation rendering them functional across all their instances exists with orientational requirements instead depending on gene-locus specific additional factors, or that the binding orientation of transcription factors may generally not be relevant , but rather the event of binding itself. 回复20180117,直达文章链接。 更有趣的:核受体NR,以二聚体形式结合DNA,motif是各种方向的回文序列: 上图出自:Martin, Luc Tremblay, Jacques J.. (2010). Nuclear Receptors in Leydig Cell Gene Expression and Function. Biology of reproduction. 83. 3-14. 10.1095/biolreprod.110.083824. 补充技能: 在线预测转录因子结合位点,用JASPAR2016,入口在这里: http://jaspar2016.genereg.net/cgi-bin/jaspar_db.pl 从主页搜索感兴趣的转录因子,进入如下页面,在motif左侧的小方框内打钩,然后把DNA序列粘贴到右下角的大方框里,scan。 2016版本的这个在线工具只能上传20000bp以内的DNA序列。 JASPAR2018考虑到现在实验猫都会写代码了,就省掉了在线分析工具。如果想在更长的序列内查找更多的motif,就用 TFBSTools这个R包,包治百病。 总之,预测一个基因启动子区是否有某个转录因子的结合位点时,看到strand+1,意味着motif跟基因同向;strand-1,意味着motif在互补链上。无论motif在哪条链上,这个转录因子都有可能调控你的基因。 转录调控研究策略: 转录因子调控了谁? 哪个蛋白质调控我感兴趣的基因? 用ATAC-seq找关键转录因子、研究转录因子-靶基因 变异会影响转录?SNP影响转录因子结合?RegulomeDB non-coding区的SNP影响了转录?HaploReg
个人分类: 转录调控|2955 次阅读|0 个评论

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

GMT+8, 2024-6-16 13:10

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部