宏基因组基因定量分析
一、salmon简介
由于样品中包含的微生物丰度不同，相应的基因丰度也不同，类似于 RNAseq 中基因表达量的差异。宏基因组中同样可以对基因进行定量。利用 Salmon 软件可以对宏基因组基因丰度进行定量。salmon 是一款新的、极快的计数软件。它与 Kallisto 和 Sailfish 类似，可以不通过 mapping 而获得基因的 counts 值。Salmon 的结果可由 edgeR 或 DESeq2 等进行 counts值的下游分析。
网址：https://combine-lab.github.io/salmon/
案例：https://2017-cicese-metagenomics.readthedocs.io/en/latest/salmon_tutorial.html
二、软件安装
#安装salmon
conda install -y salmon=1.4.0 三、使用 salmon 定量
#创建索引目录
mkdir salmon
cd salmon
ln -s ../prokka/mg.filter.ffn .
#建索引
salmon index -t mg.filter.ffn -i mg
salmon quant -i mg -l A -1 /share/home/xiehs/18.mags/2/illumina/ERR4007992_1.fastq.gz -2 /share/home/xiehs/18.mags/2/illumina/ERR4007992_2.fastq.gz -o mg.quant
#多个样本分析后，文件夹*.quant就可以使用以下命令合并，后续做差异分析
#合并TPM
#salmon quantmerge --quants ../salmon/*.quant -o result/salmon/gene.TPM
#合并原始reads count值
#salmon quantmerge --quants ../salmon/*.quant --column NumReads -o result/salmon/gene.count
#sed -i '1 s/.quant//g' result/salmon/gene.*写在最后：有时间我们会努力更新的。大家互动交流可以前去论坛，地址在下面，复制去浏览器即可访问，弥补下公众号没有留言功能的缺憾。
bioinfoer.com有些板块也可以预设为大家日常趣事的分享等，欢迎大家来提建议。
本文参与 腾讯云自媒体同步曝光计划，分享自微信公众号。
原始发表：2023-01-14，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
推荐阅读