×

信实 复现 是的 代码 文章

生信实操丨是的你没看错!0代码复现4分生信文章的7图1表

jnlyseo998998 jnlyseo998998 发表于2023-03-25 01:11:04 浏览27 评论0

抢沙发发表评论

撰写: S.Yi 来源:小张聊科研平台的“ i生信”公众号,微信公众号搜索“ i生信”即可关注/扫描关注见文末

#生信分析# #生信发文# #零代码# #数据库#

又是一篇0代码0实验的4+分生信文章,给还奋战在医学院校的科研僧们打了一只兴奋剂!话不多说,现在就和小编一起复现发表在Journal of Translational Medicine期刊上的这篇文章——The role of ARHGAP9: clinical implication and potential function in acute myeloid leukemia(图1)。

图1

这篇文章用8个免费在线的癌症数据库说明了ARHGAP9在急性髓细胞白血病的临床意义和潜在功能,对于忙碌的临床医生和求毕业的医学生来说简直是黑暗中的一丝光亮。

①在Figure1中(图2),作者用HPA和CCLE数据库中分析得到:AML细胞系中ARHGAP9过表达。

图2

展开全文

>> Figure 1a:打开HPA数据库()

图3

>> Figure 2c:打开CCLE数据库)。

图4

②在Figure2中(图5),作者用UALCAN、GEPIA和GEO数据库中分析得到:AML病人中ARHGAP9过表达。

图5

>> Figure 2a:打开UALCAN数据库,。

图6

>> Figure 2b:打开GEPIA数据库,。

图7

>> Figure 2c: 再在Expression DIY板块,选择BOX PLOTS,输入ARHGAP9,LAML数据库(图8),即可得到Figure 2c

图8

>> Figure 2d&e:在NCBI网站中,)

图9

>> 在Profile graph中输入226906_s_at(也就是ARHGAP9,可以通过旁边View data for GPL570查找基因对应编码)再点击Sample Values,得到ARHGAP9的表达量,复制数值与对应分组对应(图10),再作图为Figure 2d。同理可复制Figure 2e

图10

③在Table 1中(图11),说明了ARHGAP9与临床患者资料的相关性。作者在TCGA数据库中下载了151例患者ARHGAP9的表达量和临床资料,并统计为表格。

图11

>> 打开)

图12

>> 下载RNA-seq数据

TCGA-LAML数据库中一共有200例,有RNA-seq数据的只有151例(图13)

图13

点击进入后,选择HTSeq-FPKM这种RNA-seq结果呈现方式,再点击Add all files to carts,点击Download-Cart,即可获取每个病例的RNA-seq结果,搜索ENSG00000123329即可得到ARHGAP9表达量。注:下载后的文件名称需与Case UUID相对应(图14)

图14

>> 下载临床信息。

在此网站输入TCGA-LAML后(图12),点击RNA-seq进入Case页面(图13),随意点击其中一个Clinical,包含了两个文件,一个是此Case的临床信息,另一个是在TCGA-LAML中所有Case的临床资料,下载有200 cases的文件(图15)

图15

将下载的txt文件用Excel打开(图16),然后与之前下载的RNA-seq结果一一对应(RNAseq结果有Case ID,),即可通过统计得到Table 1。

图16

④在Figure 3&4中(图17),使用GEPIA和LinkedOmics数据库说明了ARHGAP9在AML中的预后价值。

图17

>> Figure 3a:打开GEPIA数据库,。

图18

>> Figure 3b:打开LinkedOmics数据库,)。

图19

>> Figure 3c:按照Table 1中的方式打开TCGA数据库,将临床资料与ARHGAP9的mRNA表达水平一一对应,列入Excel中(图20)。

然后通过Graphpad prism的X-Y功能得到生存分析图.

详细教程见官网链接:

图20

>> Figure 4:由于我们已获得临床资料和其对应ARHGAP9的mRNA表达,因此可整理不同临床分组的ARHGAP9表达值,然后再通过Graphpad prism的X-Y功能得到生存分析图。

⑤在Figure 5中(图21),作者通过火山图和韦恩图展示了ARHGAP9高表达组和ARHGAP9低表达组的转录组差异特征。

通常我们采取R语言进行火山图和韦恩图的绘制,在这里小编给大家推荐几款用于可视化做图的在线网址,助力大家轻松科研。

图21

>> ImageGP:打开网址)。

图22

>> LinkedOmics:首先选择肿瘤类型—LAML,接着选择数据类型—RNAseq结果,下一步填写基因名称ARHGAP9,选择靶数据—RNAseq类型,最终选择Spearman correlation,点击Submit Query,即可得到火山图(图23)。还可下载数据Download data进一步绘制韦恩图。

图23

⑥在Figure 6中(图24),作者进一步分析了ARHGAP9表达高、低两组之间具有特异性差异的转录基因的生物学功能通路。

图24

>> Figure 6a-c:打开Metascape网址:)。

详细内容可参考:

图25

⑦在Figure 7中(图26),作者进一步分析了ARHGAP9表达高、低两组之间具有特异性差异的转录基因的交互网络,寻找到与其密切相关并作用于AML的KIF20A基因。

图26

>> Figure 7a&b:打开String网址:。

图27

>> Figure 7c&d:在GEPIA数据库中,检索评分高的基因与AML的关系,鉴定得到KIF20A基因与之相关。

小编有话说,整篇文章通过8个开源的数据库,确定了ARHGAP9在AML疾病中的重要作用,并通过分析ARHGAP9高、低表达两组的转录组特异性差异,找到与之相关的KIF20A基因。那后续即可通过细胞实验&临床检测来验证推测,那又是一篇文章啦。临床医生和求毕业的医学生发3分文章也不是那么难呀!

原文链接: