上游:cellranger

cellranger更新到5啦(全新使用教程) - 云+社区 - 腾讯云 (tencent.com)

cellranger更新到6.0啦(全新使用教程) - 简书 (jianshu.com)


下载安装软件

官网: https://support.10xgenomics.com/single-cell-gene-expression/software/downloads/latest

#官网简单注册获得cellranger下载地址

wget -O cellranger-6.1.2.tar.gz "https://cf.10xgenomics.com/releases/cell-exp/cellranger-6.1.2.tar.gz?Expires=1644969388&Policy=eyJTdGF0ZW1lbnQiOlt7IlJlc291cmNlIjoiaHR0cHM6Ly9jZi4xMHhnZW5vbWljcy5jb20vcmVsZWFzZXMvY2VsbC1leHAvY2VsbHJhbmdlci02LjEuMi50YXIuZ3oiLCJDb25kaXRpb24iOnsiRGF0ZUxlc3NUaGFuIjp7IkFXUzpFcG9jaFRpbWUiOjE2NDQ5NjkzODh9fX1dfQ__&Signature=h4FL5tsIzSePQ92iKlyuD1s4ybXSl68zgxD-5MBkpQBIS8liLWaKPcOPDi-6yURDXPbEsujPlvif2CiWzPdIP5ZlxSZ-s2n9Ki73Wo2mYq-MmVelLTWhO0cEznzDNZAo-aZJKro008WJQ~zzhKVSDYy7kKC4O64pmcpl1OTj5dBrL64zZrv0bLc1J~W26dbx1IU983P2VTVA2l-gVLSKwpoIGooWYPTdXOCIYzp77iEnPfkdMuZ~M7nfm6JycmSKh3XPaEkkWlNFEX9stBCYKxrxV0Z9gfNXNg10mZ5vITB1DLe2YcLg-V9Kwbh-~dZgTgVn81oZsIJt9Zp~3FFQZA__&Key-Pair-Id=APKAI7S6A5RYOXBWRPDA"

或者

curl -o cellranger-6.1.2.tar.gz "https://cf.10xgenomics.com/releases/cell-exp/cellranger-6.1.2.tar.gz?Expires=1645022614&Policy=eyJTdGF0ZW1lbnQiOlt7IlJlc291cmNlIjoiaHR0cHM6Ly9jZi4xMHhnZW5vbWljcy5jb20vcmVsZWFzZXMvY2VsbC1leHAvY2VsbHJhbmdlci02LjEuMi50YXIuZ3oiLCJDb25kaXRpb24iOnsiRGF0ZUxlc3NUaGFuIjp7IkFXUzpFcG9jaFRpbWUiOjE2NDUwMjI2MTR9fX1dfQ__&Signature=LnkvS5UP2QM4FR7DPlL0lgkBQROtRWmv35P715ML5AZ1~2mcAY7x9OHlB76hjnYmpHbIcAq0Mtvi8mtaRebtE6zTdkaTrIdmiZSFOY7C4Bt3aQzWmNWLA7kyFCbjtktpaCRaSAbjTmOmRC3gEq9fr4EZM5m2nvYXbN0CnRtV1Rk0iyGTBTXZuyjgRKcddT1iYDmUJ1znFWpxE9Dpx-goxIfSOTLOTW3IGfD2Q-7TE-6iZbOrknmFQ5nYBaSqAMGZy~H8Xop-yyNU3PBiClseIlE6iJ1pKHJy5S~Rxfjyoiom30tHp9dVLCvkX4-UzBdGbzFYx-neWa6b8u~ddNvajA__&Key-Pair-Id=APKAI7S6A5RYOXBWRPDA"

#解压软件(82base)

tar -xzvf cellranger-6.1.2.tar.gz

#添加到环境变量

vim ~/.bashrc

export PATH=/opt/cellranger-6.0.2:$PATH

source ~/.bashrc

下载ref

#下载小鼠参考基因组

wget https://cf.10xgenomics.com/supp/cell-exp/refdata-gex-mm10-2020-A.tar.gz

#下载人参考基因组

wget https://cf.10xgenomics.com/supp/cell-exp/refdata-gex-GRCh38-2020-A.tar.gz

#解压参考基因组

cd /data2/ff/ref/cellranger.mm10

tar -xzvf refdata-gex-mm10-2020-A.tar.gz

###

如果你是直接拿到的R1/R2的fastq文件,那么就直接上cellranger count

如果你是I1/R1/R2的数据,那就麻烦还要跑个cellranger mkfastq

cellranger count定量(82 base)

#图中的四对样本来自于同一群细胞,我想将这四个文件一起一次性跑完

#cellranger count当中的--sample参数只识别的时候只识别S1之前的字段

cat CC5F0-1_S1_L003_R1_001.fastq.gz CC5F0-2_S1_L003_R1_001.fastq.gz CC5F0-3_S1_L003_R1_001.fastq.gz CC5F0-4_S1_L003_R1_001.fastq.gz > CC5F0_S1_L003_R1_001.fastq.gz

#运行

#单样品版

#换目录

cd /data2/ff/project/LJ.22.02.sc

cellranger count -- id=KOAC \  # 给你这次的任务取一个响亮的名字,随意,id指定输出文件存放目录名

--transcriptome=/data2/ff/ref/cellranger.mm10/refdata-gex-mm10-2020-A \  #参考基因组

--fastqs=/data2/ff/project/LJ.22.02.sc \   #你的fastq.gz文件保存的路径,完整

--sample=KOAC \  # 你的fastq.gz文件名当中S1之前的字段,作为软件识别的标志

--nosecondary

####循环版本

cd /data2/ff/project/LJ.22.02.sc

for i in KOAC KONC WTAC WTNC

do

#定量
cellranger count \

--id=${i} \  # 给你这次的任务取一个响亮的名字,随意

--transcriptome=/data2/ff/ref/cellranger.mm10/refdata-gex-mm10-2020-A \  #参考基因组

--fastqs=/data2/ff/project/LJ.22.02.sc/${i} \   #你的fastq.gz文件保存的路径,完整

--sample=${i} \  # 你的fastq.gz文件名当中S1之前的字段,作为软件识别的标志

done


#--nosecondary  ## nosecondary 只获得表达矩阵,不进行后续的降维、聚类和可视化分析(因为后期会自行用R包去做)

输出结果


filtered_gene_bc_matrices:
是重要的一个目录,下面又包含了 barcodes.tsv、features.tsv、matrix.mtx,是下游Seurat、Scater、Monocle等分析的输入文件。

web_summary.html:质控比对报告(一般认为外显子的比对率要在60%以上)。

index标记样本,

barcode标记细胞,其次,barcodes数量时要大于细胞数量的(以保证每个细胞都会有barcode来进行区分)。

UMI标记转录本,Unique Molecular Identifier,由4-10个随机核苷酸组成,在mRNA反转录后,进入到文库中,每一个mRNA随机连上一个UMI,根据PCR结果可以计数不同的UMI,最终统计mRNA的数量。它的主要作用是,处理PCR 扩增偏差,因为起始文库很小时需要的PCR扩增次数就越多,因为越容易引入扩增误差。

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容