【原】玩转RNA-seq数据也可以不需要linux ？

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./geo/query/acc.cgi?acc=GSE109480 ： 

IDNAME	NAME
GSM2944116	Control 1
GSM2944117	Control 2
GSM2944118	Control 3
GSM2944119	Control 4
GSM2944120	Control 5
GSM2944121	KPC 1
GSM2944122	KPC 2
GSM2944123	KPC 3
GSM2944124	KPC 4
GSM2944125

提供表达矩阵下载：

下载文件，进入R处理它们：

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./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，你的数据有无限的可能！

另外，我这周六日在广州和深圳有演讲：