CIBERSORTx
CIBERSORT初探input输入 output 代码详解 原理解析免疫浸润 反卷积_YoungLeelight的博客-CSDN博客
https://blog.csdn.net/qq_52813185/article/details/127469758?csdn_share_tail=%7B%22type%22%3A%22blog%22%2C%22rType%22%3A%22article%22%2C%22rId%22%3A%22127469758%22%2C%22source%22%3A%22qq_52813185%22%7D
持续安利免疫浸润分析工具的小云又来啦~前边几期给大家分享了很多的工具,有同学留言问小云cibersort的教程呢?呐~这不就来了,重要的总要在后边登场。各位同学请放心,cibersort的教程,虽晚必到。
cibersort是目前引用次数最多的免疫细胞浸润估计分析工具,2015年首次发表于Nature method。CIBERSORT是基于线性支持向量回归(linear support vector regression)原理对人类免疫细胞亚型的表达矩阵进行去卷积的工具。用于芯片表达矩阵和测序表达矩阵,对未知混合物和含有相近的细胞类型的表达矩阵的去卷积分析优于其他方法 (LLSR,LLSR,PERT,RLR,MMAD,DSA) 。该方法基于已知参考数据集,默认提供22种免疫细胞亚型的基因表达特征集:LM22。
官网介绍:原文介绍 Bulk tissue deconvolution with single-cell reference profiles.Determining cell type abundance and expression from bulk tissues with digital cytometry | Nature BiotechnologyWe previously developed an approach for digital cytometry, called CIBERSORT, that enables estimation of cell type abundances from bulk tissue transcriptomes. We now introduce CIBERSORTx, a machine learning method that extends this framework to infer cell-type-specific gene expression profiles without physical cell isolation. By minimizing platform-specific variation, CIBERSORTx also allows the use of single-cell RNA-sequencing data for large-scale tissue dissection. We evaluated the utility of CIBERSORTx in multiple tumor types, including melanoma, where single-cell reference profiles were used to dissect bulk clinical specimens, revealing cell-type-specific phenotypic states linked to distinct driver mutations and response to immune checkpoint blockade.
LM22.txt获取方法:
在Cibersort论文中下载Supplementry table 1
https://www.nature.com/articles/nmeth.3337#MOESM207
CIBERSORT同时提供R包和在线网页工具,可以对用户的表达矩阵进行非负矩阵分解后计算不同细胞类型的比例。
CIBERSORT网页版链接为:https://cibersortx.stanford.edu/
需要说明的是,使用网页版需要先进行账号注册。
CIBERSORT需要两个准备文件,芯片或者测序表达矩阵以及参考数据集。
表达矩阵格式第一列是基因名,第一行是样品名,不能有重复基因名,第一列列名不能空白。矩阵中不能存在空白或NA值,不要对表达量取Log2。
网页版有LM22的参考数据集,只需要提供表达矩阵文件就可以了。
具体操作如下:
点击运行,等待结果即可。
Cibersort R 代码操作:
在R中新建R Script,复制以下网址中代码,保存为“Cibersort.R”
https://rdrr.io/github/singha53/amritr/src/R/supportFunc_cibersort.R
(4)以上三个文件需保存在同一文件夹,运行Cibersort的代码如下:
setwd("三个文件的文件夹")source('Cibersort.R')result <- CIBERSORT('LM22.txt','DATA.txt', perm = 1000, QN = T) #perm置换次数=1000,QN分位数归一化=TRUE
在同一文件夹下可以得到运算结果("CIBERSORT-Results.txt")
注意Cibersort结果的默认文件名为CIBERSORT-Results.txt,在同一文件夹下进行第二次运算会覆盖第一次得到的文件,建议在每一次运算之后对文件重命名。
运行结果示例:
基于结果表格,可以进行多种可视化展示,比如箱线图、柱状图等。
cibersort作为目前被引用次数最多的免疫浸润工具是不是也不难,如果大家感兴趣的话,可以做起来啦~~