单细胞RNA测序(scRNA-seq)cellranger count的细胞定量和aggr整合

本文详细介绍了单细胞RNA测序(scRNA-seq)中cellranger count的细胞定量步骤,包括参考基因组和注释文件的准备、测序数据的处理、细胞定量结果的分析以及多个结果的整合。内容涵盖shell脚本的使用,以实现自动化处理和数据分析。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

单细胞RNA测序(scRNA-seq)基础知识可查看以下文章:

单细胞RNA测序(scRNA-seq)工作流程入门

单细胞RNA测序(scRNA-seq)细胞分离与扩增

单细胞RNA测序(scRNA-seq)SRA数据下载及fastq-dumq数据拆分

单细胞RNA测序(scRNA-seq)Cellranger流程入门和数据质控

细胞定量是scRNA-seq重要的分析步骤,主要是进行细胞与基因的定量, cell ranger将比对、质控、定量都封装了起来,使用起来也相当便捷。

1. 参考基因组和注释文件准备

1.1 参考基因组、注释下载和参考基因组索引构建

参考基因组、注释下载注释和参考基因组索引参考以下文章:

单细胞RNA测序(scRNA-seq)构建人类参考基因组注释

构建完成后,GRCh38( cellranger mkref --genome 参数内容为指定输出文件夹)目录会出现以下内容。

# genes中gtf文件为.gz格式则解压
gzip -d genes.gtf.gz

在这里插入图片描述

# 查看genes/genes.gtf信息
grep -v '^#' genes.gtf |less -S

genes.gtf

1.2 注释基因信息更新

  1. 在fasta/genome.fa的fasta文件基础上增加序列信息;
  2. 在genes/genes.gtf的gtf文件基础上增加注释信息(根据GTF文件格式);
  3. 重新执行运行cellranger mkref流程

2 cellranger count细胞定量

2.1 测序数据fastq文件命名说明

# 一般命名方式
[Sample Name]_S1_L00[Lane Number]_[Read Type]_001.fastq.gz

# 其中Read Type有以下三种:
# I1: Sample index read (optional)
# R1: Read 1
# R2: Read 2

# 查看rawdata目录下fastq.gz文件
ls -lh rawdata 

rawdata目录文件列表:

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

生信与基因组学

每一份鼓励是我坚持下去动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值