前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >玩转RNA-seq数据也可以不需要linux ?

玩转RNA-seq数据也可以不需要linux ?

作者头像
生信技能树
发布2019-09-17 17:15:02
6860
发布2019-09-17 17:15:02
举报
文章被收录于专栏:生信技能树

Hepatocytes direct the formation of a pro-metastatic niche in the liver. Nature 2019 Mar;567(7747):249-252.

看到这篇实验文章 [ PMID: 30842658 ],它里面用到了mRNA-seq技术:

We next performed mRNA sequencing on RNA isolated from the livers of control and KPC mice. We identified 275 differentially expressed genes, and found that genes upregulated in KPC mice were associated with immune-related processes.

但是我看了看文章的方法部分,并没有介绍linux环境或者R数据处理,只有纯粹的软件工具。

文章用到的数据公布在 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE109480

提供表达矩阵下载:

下载文件,其实我们会R的话,就可以进入R处理它们:

代码语言:javascript
复制
rm(list = ls())
options(stringsAsFactors = F)
setwd('files/')
el=lapply(list.files(pattern = '*.gz'),function(x){
  read.table(x,header = T,row.names = 1)[,9,drop=FALSE]
})
lapply(el, head)
exprSet=do.call(cbind,el)
colnames(exprSet)=paste0(c('control','KPC'),c(1:5,1:5))
# The reads per 1k bases of exon model per million mapped reads.
# (http://www.nature.com/nmeth/journal/v5/n7/abs/nmeth.1226.html)

差异分析后的热图

功能富集及注释

chemoattractant相关基因集的表达量分布

GSEA分析挑选感兴趣通路

Enrichment of IL-6–JAK–STAT3 signalling genes in the liver (n = 5 for control mice and NTB KPC mice).

写在后面

真正的粉丝看到这样的标题肯定会以为我们公众号被盗了,居然不推荐一个纯正的生信工程师学linux,开什么国际玩笑!

是的,本文并没有劝退大家离开linux的意思,相反,我会极力推广:

视频都在B站:

其实吧,如果只是为了获得表达矩阵,那当然,linux不重要,可是,表达矩阵只不过是转录组数据分析的冰山一角:

会linux,你的数据有无限的可能!

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2019-09-14,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 生信技能树 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档