生信宝典之傻瓜式 (五) - 文献挖掘查找指定基因调控网络

如何快速查找指定基因的调控网络介绍了使用在线查询数据库 (http://evexdb.org/)对PubMed和PubMed Central中发表文章的摘要和全文为依据进行文本挖掘探寻基因直接可能的相互作用的工具。反响很好,但现在网站似乎出了点问题,获得的相互作用细节信息不能展开了(推测可能是使用的JS库无法加载)。有朋友留言推荐 Cytoscape literature search,一个存在历史挺久的Cytoscape插件,通过给定关键字搜索文献,并且基于搜索结果构建互作网络,帮助研究者快速搜索和提取基因之间,蛋白之间可能的联系,兼容2.7和3.1版本,下载32,742次。

作为Cytoscape的插件,安装很方便,基本的Cytoscape使用见教程。

安装完之后,从菜单栏Apps-Agilent Literature Search启动,使用界面如下。

左侧一般输入一个或多个基因 (若输入多个则每一行输入一个),右边限制一个环境,可以是物种,也可以是某种疾病如lung cancer,或某个过程stem cell。下面的选项还可以选择是否使用别名 (选择后我们输入的pou5f1就被转成了oct4, otf4等),限定物种,限定相互作用的判断 (个人一般使用relaxed)。具体每个参数的含义详见后面解释。

前面输入的内容都会在Query Editor中转换为逻辑查询表达式的形式,方便查看搜索的内容是否符合自己的需要,也可以自行修改,比如我们把stem cell改为AND连接。

点击蓝色箭头就可以启动搜索。搜索到的文献展示在左下角,可点击跳转到PubMed,右键删除某一项。

右侧展示的是挖掘出的调控网络,可以根据属性进行一些修饰、美化和查询。

网络查看和美化

首先是调大字体 (Label Font Size),设置搜索出的基因和挖掘出的相互作用基因不同的颜色 (Fill Color-searchTerm-Discrete Mapping), 删掉与核心基因没有连线的点,Apply preferred layout重新调整布局 (工具栏上的刷新按钮)。

如果觉得线太过扭曲,可以Layout-Clear all edge bends,然后再点击Layout-Bundle edges使连接看上去圆润。初步修饰下,效果如下:

如果常用,每次调样式也比较麻烦,可以把样式保存起来,点击样式旁的三道杠,选择Copy style重新命令 (若不导出,关闭后就不见了),然后File-Export导出。下次查询好之后,再File-Import导入就好。也可后台回复 style,获取我们这个简单的样式,以此为基础修改。

每个节点,点击右键,按下图点选菜单,可以看到支持其相互作用的文章句子节选,方便快速阅读和理解潜在的调控关系。

选项解释

寻求帮助

在上述查询界面,按F1或点Help可以打开帮助页面如下。当前截图显示的是通过View菜单选择搜索的数据库,现在只有PubMedUSPTO (美国商标专利数据库)可选。

Max Engine Matches: 限制每个字符串在每个库最大查询到的结果数目。

Use Aliases: 选定后,将会根据Concept Lexicon中限定的物种寻找左侧输入框输入的内容的别名。查询时,有一个别名匹配上就可以。

Use context: 是否使用Context面板 (右侧输入框)限定查询。

Concept Lexicon Limits Search: 如果需要把搜索结果限制在某个物种,则勾选。

Concept Lexicon: 通常是物种相关的选项,对Use aliases的判断和搜索结果提取有效,但不用于限制查询结果。所以如果要在查询时就限制物种,则需要再右侧输入框输入物种的名字,会加快查询速度。

Interaction Lexicon: 限制判断相互作用的严格程度。对于每个包含搜索关键字的句子,都会来判断里面是都包含interaction lexicon收录的动词,如activate, enhance, cause等。这些关键词可以修改,有严格版和宽松版。

Load and Save: 搜索结果可以存储和再次导入。

选项进一步解释和自定义

如果Windows下,LiteratureSearch的配置文件在目录C:\Users\sxbd\CytoscapeConfiguration\app-data\com.agilent.labs.als.AgilentLiteratureSearch-3.1.1\data下。(把sxbd改为您的用户名)

Interaction Lexicon:前面提到的limit, relax, empty每一个的效果都记录在文件interaction-lexicon-map.txt中,文件内容如下

limited    data/strictVerbNames.txt
relaxed    data/verbNames.txt
empty    data/emptyVerbNames.txt

每个不同的参数表示使用的关键词列表不同,empty表示不进行限定,只要两个词出现在一个句子中就认为有作用。

strict表示严格限定,默认要求句子中必须含有收录的15个单词中的一个才认为存在相互作用 (在文件strictVerbNames.txt)。

relaxed默认要求句子中必须含有收录的75个单词中的一个才认为存在相互作用,涉及促进、抑制、结合、催化等对应的英文单词和变种,在使用过程中,我们也可以不断完善、添加更多词汇到verbNames.txt中,以获得更多关注的相互作用。

Concept Lexicon

这个由文件concept-lexicon-map.txt控制,默认收录了常见物种的KEGG注释信息、基因的别名信息。

Arabidopsis thaliana    data/.uc_Arabidopsis_thaliana    
Bos taurus    data/.uc_Bos_taurus
Caenorhabditis elegans    data/.uc_Caenorhabditis_elegans
Danio rerio    data/.uc_Danio_rerio
Drosophila melanogaster    data/.uc_Drosophila_melanogaster
Escherichia coli    data/.uc_Escherichia_coli
Homo sapiens    data/.uc_Homo_sapiens
Mus musculus    data/.uc_Mus_musculus
Rattus norvegicus    data/.uc_Rattus_norvegicus
Saccharomyces cerevisiae    data/.uc_Saccharomyces_cerevisiae

以人的数据为例,前面是KEGG编号、对应的描述,后面是基因的每个名字一行,方便使用别名搜索。如果不知道基因别名或想下载KEGG注释,直接用这提供的就好。

6.3.5.8    aminodeoxychorismate synthase    adc synthase    4-amino-4-deoxychorismate synthase    pabb
6.3.5.9    hydrogenobyrinic acid a, c-diamide synthase (glutamine-hydrolysing)    cobb
dynamin    dynamin-1    dynamin1
dynamin-2    dynamin2
epsin    epsin1    epsin-1
nf-kappaB    nfkappaB    nfkb1    nfkb    nf-kappa B    nfkappa B
frizzled    fz    fzd    fzd7
dsh    disheveled    dishevelled    dsh1    dvl1l1    dvl1
bcatenin    beta-catenin    beta catenin))

之前如何快速查找指定基因的调控网络文章下有朋友留言,非模式生物怎么查找,一个是利用生信宝典之傻瓜式(四)蛋白蛋白互作网络在线搜索中提到的在线工具STRING收录了2031物种。另外一个就是在这自定义需要的文件,使用此插件搜索。

原文发布于微信公众号 - 生信宝典(Bio_data)

原文发表时间:2018-01-19

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏腾讯NEXT学位

那些让编码效率起飞(前端)的工具了解一下

? | 导语 想晚上吃鸡?前端编码效率提升工具了解一下? 一、Bash篇(Mac) iTerm2 iTerm 2 is a terminal emulato...

1733
来自专栏Angular&服务

关于做angular4引入echarts图表

NgxEchartsService 是全局 echarts 对象的包装器。您可以直接获取本机echarts对象或使用包装器方法。例如:

5733
来自专栏SDNLAB

VXLAN与EVPN的结合使用

EVPN是近几年最热门的网络技术之一。如果你还没听过EVPN,你的网络技能可能已经落伍了,赶紧看看之前的《EVPN简介》吧!EVPN全称是Ethernet VP...

4243
来自专栏Golang语言社区

从零开始创建一个基于Go语言的web service

20个小时的时间能干什么?也许浑浑噩噩就过去了,也许能看一些书、做一些工作、读几篇博客、再写个一两篇博客,等等。而黑客马拉松(HackAthon),其实是一种自...

4439
来自专栏张戈的专栏

教你如何去掉友荐和无觅的隐藏外链和版权链接

最近,为了做无觅的 APP,装上了无觅的相关推荐,结果果断不给力,打包了 2 个星期还在打包,还能再坑点么? ? 蜗牛般的效率,暂且就不吐槽了。偶尔用站长工具检...

4358
来自专栏FreeBuf

新型勒索病毒软件GruxEr来袭:深度分析如何传播、加密及如何删除

本文旨在帮助您清除Gruxer ransomware感染,并通过此ransomware病毒恢复AES加密文件。 据报道,GruxEr名称的ransomwar...

2726
来自专栏杨建荣的学习笔记

假期前的数据库检查脚本之主备关系(r11笔记第46天)

快过年了,很多系统都要进入最后的检查和复验阶段,一方面在节假日前,提前发现问题总比过节的时候发现要好。另一方面如果出现故障的时候能及时进行处理,这个时候我们就需...

28810
来自专栏SAP最佳业务实践

SAP最佳业务实践:MM–不交货与库存调拨(135)-2业务处理

3 流程概览表 选项一:不含库存调拨订单的库存调拨 – 一次性库存调拨 流程步骤业务条件业务角色事务代码预期结果一次性库存调拨需调拨的库存仓库文员MIGO在系统...

3413
来自专栏Timhbw博客

Hexo-完全免费全平台搭建个人博客(2)-域名主题设置

2017-03-1011:01:58 发表评论 913℃热度 Hexo-完全免费全平台搭建个人博客(1)-整体搭建 上一篇文章把 Hexo 博客整体搭建一遍了...

37712
来自专栏张戈的专栏

另类SEO分享:利用JS封装iframe躲过搜索引擎的抓取

前言:很多博友不仔细看完内容就直接认为用 iframe 不好之类的云云,而实际上本文就是教你在必须使用 iframe 的时候,该如何躲过搜索引擎的抓取,避免不利...

5066

扫码关注云+社区