原标题:少妇体验2020年零代码4分苼信套路拆解+全文复现
解螺旋公众号·陪伴你科研的第2236天
Microenvironment不知道合不合你们胃口。会了代码的学员可能会不喜欢这种零代码的文章认为鼡上代码都发不了好文章,还零代码你这种想法从某种角度来说也没错,毕竟如果纯生信10+的文章那不上代码可能不行,但是找个好的切入点利用各种数据库零代码解决“温饱”,应该还是可以的这篇文章是2020年发表在Frontiers in
点击Expression DIY下方的Stage plot,进入下方页面后输入基因CXCL1并选择肿瘤为KIRC,然后设置相应颜色如下:
接着点击Plot,下方就会出来相应图片:
这个图片跟文章是一摸一样了按照上面的步骤依次对其他CXC家族成員进行分析,然后保存所有图片后使用AI进行拼图Figure 4又被不知不觉搞定了,好像速度有点快接下来继续往下。
Figure 5还是用GEPIA做的图表示的是肾細胞癌患者无病生存曲线(DFS)中不同表达CXC趋化因子的预后价值,那接着我们还没关闭的GEPIA继续往下:
按照上图红色框中的内容进行设置然后点擊Plot,这样CXCL1的图就又出来了(奇怪
一样的步骤再做其他CXC家族的成员,然后进行拼图即可接下来我们再往下:
OMG,又是生存曲线还只是把DFS換成了OS,你们别信啊生信不是这么简单的东西,只是刚好这个文章就用GEPIA弄了好几个图而已我们再来看看操作,还是刚才的页面把DFS换荿了OS,如下:
然后点击Plot这样CXCL1的OS生存曲线又出来了:
这样Figure 6也完成了,诶行吧,今天就到这里吧我们打到了Figure 6,下回从Figure 7开始
在左边选择腫瘤Kidney,然后选择相应的TCGA RCC数据集如下:
点击Query by gene,到达下面页面然后进行相应设置:
这个结果跟文章似乎稍微有点出入,可以多试试几个参數不过这个结果是我试了几次后最接近文章结果的结果了,考虑到cBioportal曾经改版过似乎这个有差异也算正常,既然说到cBioportal就不得不说一下,新版的网站取消了network的模块所以这7B和7E是没办法做了≡(▔﹏▔)≡,那我们接着继续到7C吧;
点击SEARCH进入下一个页面:
按上图输入CXC家族选择human后點击SEARCH,到达下一页面:
在Gene List的位置输入CXC家族然后Tissue Order选中dataset为KIRC肾透明细胞癌,其他选择默认参数然后点击Plot,得到下方图片:
这个图又跟文章的結果一致了有没有人狠话不多,接着7D;
这个图其实我在训练营就讲过了来过训练营的学员可能有印象,这个图是GeneMANIA (http://www.genemania.org)网站的图片我们打開网站:
在左上角输入CXC家族成员,点击SEARCH结果如下:
点击左下角的小圆圈可以给网络添加相应的颜色:
可以选择排名靠前的通路,也可以選择符合自己研究方向的通路当然还可以对网络排列进行修饰:
相应的操作都在图片坐标,多点点选择自己喜欢的style,又不花钱点点沒坏处嘛,这样我们7D也就解决了接下来的Figure 8是 肾细胞癌中不同表达的CXC趋化因子和50个最常改变的邻近基因的富集分析,其中50个邻近基因也是從cBioportal得到的结果也就是基于Figure 7E,网站改版后暂时无法复现接下来看到Table
点击中间的SEARCH,在出来的页面拉到下方然后如下图设置:
输入CXC家族成員之后,点击提交得到结果如下:
将表格整理就可以放入文章中,点击表格中红框的数字可以直接跳转到相应详情页面,这里就不给夶家展开啦最后到了Figure 9;
我们以CXCL1为例,选择TIMER首页的GENE模块点击进入下方页面进行设置:
设置完成后点击Submit提交,得到了下方的结果:
按照上媔的步骤依次分析剩下的CXC成员然后拼成Figure 9,这样就大功告成啦!等等等等好像还有一个表格:
这个表格其实也可以在TIMER进行分析,我们点擊TIMER模块中的Survival然后进行如下设置(这里我以CXCL1和CXCL2为例):
设置完成后,在右边就可以看到相应的分析结果了:
然后整理成表格就完成啦!这樣我们文章的复现操作就此结束!
这篇零代码的文章用的很多数据库都是在我们生信全书上篇中讲解过的网站,有了好工具关键的还昰怎么组成一个好的故事,当然思路也是很重要,当然文章复现过程中还有很多细节需要大家去进一步细化这些都是打磨的工作,相信难不倒大家一篇文章复现完了,也不是就此结束我们可以思考下,还有哪些网站的哪些分析可以加进来让数据更加丰富呢?
欢迎夶家关注解螺旋生信频道-挑圈联靠公号~
点下“在看”多根头发