稀有物种加测bcr/tcr vdj How to create a Cell Ranger compatible V(D)J reference? cellranger

10x官网给出了如下几个答案

How to create a Cell Ranger compatible V(D)J reference?

Question: I want to run Cell Ranger's V(D)J pipeline on an organism different from human or mouse. How can I create a compatible V(D)J reference?

Answer: To generate a Cell Ranger compatible reference for the V(D)J pipeline using the Ensembl database, you will need the reference genome FASTA and its corresponding gene annotations in GTF format for your organism. Then, you can run the Cell Ranger's cellranger mkvdjref utility to generate a compatible reference package.

To generate your reference using the IMGT database, use the fetch-imgt script provided with Cell Ranger v2.1+ to fetch the data from IMGT. Then run cellranger mkvdjrefto generate the reference.

Click here for more details on building a Cell Ranger compatible reference for the V(D)J pipeline. Also see a related KB article on building a custom reference for V(D)J using IMGT tool.


Note: Currently, in Cloud platform we do not support running Cell Ranger V(D)J in de novo mode without a reference. We are planning to add support for this option in Cloud in a future release. De novo mode without a reference is enabled in cellranger command line. 


Related Article: Usage of Cell Ranger VDJ denovo mode


Product: Single Cell Immune Profiling

Was this article helpful?

https://kb.10xgenomics.com/hc/en-us/articles/115002699972-How-to-create-a-Cell-Ranger-compatible-V-D-J-reference-

Recently viewed articles

  • Loupe Browser investigation for GEX + V(D)J data for differences in cell numbers
  • Web summary metrics for V(D)J Troubleshooting

Can I use Cell Ranger to analyze cell hashing data?

Question: I have single-cell gene expression samples prepared with a cell hashing protocol. Can I use Cell Ranger to analyze these data?

Answer: 10x Genomics now supports Feature Barcoding technology for Cell Multiplexing (CellPlex). This technology utilizes species-agnostic, lipid-based tagging compatible with cells or nuclei to load up to 12 samples onto one GEM well. Please see this page for more information. We also provide the cellranger multi pipeline in Cell Ranger (6.0 or later) to analyze 10x Genomics Cell Multiplexing data.

If you have used custom antibody-based hashtag oligos (HTOs) for sample multiplexing, you may use the cellranger multi pipeline in Cell Ranger v6+ by providing a reference for custom multiplexing oligos.  However, because the data generated from a custom workflow for multiplexing (e.g. CITE-seq) has not been validated with Cell Ranger’s algorithms, we cannot support custom workflows for multiplexing. Cell Ranger multi-web summaries will also provide a warning to customers to let them know a custom multiplexing reagent has been detected.  Users should exercise caution when interpreting results generated with custom reagents. 

If you encounter issues generating HTO libraries or any other custom reagents, we recommend contacting the vendor (e.g. New York Genome Center or BioLegend) directly.
 

Related articles:

How to analyze 5' multiplexed GEX+VDJ+FB(custom) data on multiple samples

I used antibody tags for cell surface protein capture and cell hashing with Single Cell 3' chemistry. How can I use Cell Ranger to analyze my data?

Was this article helpful?

 

How do I get the read counts for each barcode?

Question: Is there a way to get the read counts for each barcode in addition to UMIs?

Answer: Most customers only want the UMI counts because it corrects for amplification bias. If you are interested in the read counts, then you can extract them from the possorted_genome_bam.bam file with some custom coding. The Linux command shown below requires samtools, a copy of which can be found in your Cell Ranger installation.

Before starting, you should source the following file so Linux knows where to find samtools. Please be sure to change the bolded part below based on the location of the Cell Ranger version you have installed:

source /path/to/cellranger-x.y.z/sourceme.bash

Afterwards, the following command will compute reads per barcode. Try it out with the first 1000 reads to see if the output matches expectation. If the test run checks out, then you can remove the bolded | head -n 1000 to process the entire bam file.

samtools view possorted_genome_bam.bam | head -n 1000 | grep CB:Z: | sed 's/.*CB:Z:\([ACGT]*\).*/\1/' | sort | uniq -c > reads_per_barcode

Please keep in mind that this command can take a while to run since the computer needs to stream through the entire BAM file, and sort millions of barcodes.

Disclaimer: This article and code-snippet are provided for instructional purposes only. 10x Genomics does not support or guarantee the code.

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/640400.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

一种行之有效的防错策略:在支付系统中实施防呆设计的实践

聊个支付人都会碰到的问题:资损防控。做支付如果还没有碰到过资损,那就是做得时间还不够久。资损防控是一个很大的话题,需要开几篇文章才能讲完,今天只从一件小事入手聊一个简单而又行之有效的防错策略:防呆设计的实践…

论文阅读_tinyllama_轻量级大模型

英文名称: TinyLlama: An Open-Source Small Language Model中文名称: TinyLlama: 一个开源的小型语言模型链接: http://arxiv.org/abs/2401.02385v1代码: https://github.com/jzhang38/TinyLlama作者: Peiyuan Zhang, Guangtao Zeng, Tianduo Wang, Wei Lu机构: 新加坡科技与设…

OOM异常

OOM异常详解-CSDN博客 1.什么是OOM? OOM,java.lang.OutOfMemoryError 错误,也就是java内存溢出错误。一般当jvm虚拟机内存不够用的时候,就会抛出该错误。 2.OOM的各种情况

【C++】list容器迭代器的模拟实现

list容器内部基本都是链表形式实现,这里的迭代器实现的逻辑需要注意C语言中指针的转换。 list容器如同数据结构中的队列,通常用链式结构进行存储。在这个容器中,我们可以模仿系统的逻辑,在头结点后设置一个“ 哨兵 ”,…

14.2搭建ASP运行环境

14.2搭建ASP运行环境 可以使用WinXP来搭建ASP的运行环境,不过,WinXP不适宜作为商业用Web服务器。目前比较常用的Web服务器操作系统是Windows2003 Server,在Windows2003 Server上安装完IIS6.0并进行相关设置后,即完成了ASP运行环境…

【lettuce-排行榜】

背景&#xff1a; 这次游戏中台采用lettuce的zset完成游戏内的本服和跨服排行榜&#xff0c;因此写一下案例。 pom.xml <dependency><groupId>io.lettuce</groupId><artifactId>lettuce-core</artifactId><version>6.2.4.RELEASE</ve…

【Redis】Redis基础

Redis基础 初识Redis 认识NoSQL SQL&#xff1a;结构化查询语言 > 关系型数据库 NoSQL&#xff1a;非关系型数据库 SQL与NoSQL的差异&#xff1a; 数据结构 SQL结构化&#xff1a;表的信息依赖于表的结构NoSQL非结构化&#xff1a;存储的信息为KV形式 数据关联 SQL关联…

Kubernetes实战(十八)-Pod亲和与反亲和调度

1 调度Pod的主要方式 Pod调度到指定Node的方式主要有4种&#xff1a; nodeName调度&#xff1a;直接在Pod的yaml编排文件中指定nodeName&#xff0c;调度到指定name的节点上。nodeSelector调度&#xff1a;直接在Pod的yaml编排文件中指定nodeSelector&#xff0c;调度到带有指…

Java和Redis实现一个简单的热搜功能

1. 前言 我们有一个简单的需求&#xff1a; 搜索栏展示当前登陆的个人用户的搜索历史记录&#xff0c;删除个人历史记录。用户在搜索栏输入某字符&#xff0c;则将该字符记录下来 以zset格式存储的redis中&#xff0c;记录该字符被搜索的个数以及当前的时间戳 &#xff08;用…

Qt5.12.0 与 VS2017 在 .pro文件转.vcxproj文件

一、参考资料 stackoverflow qt - How to generate .sln/.vcproj using qmake - Stack Overflowhttps://stackoverflow.com/questions/2339832/how-to-generate-sln-vcproj-using-qmake?answertabtrending#tab-topqt - 如何使用 qmake 生成 .sln/.vcproj - IT工具网 (coder.wo…

Ubuntu安装向日葵远程控制, Team

SUN: 参考&#xff1a; Ubuntu安装向日葵远程控制——Ubuntu18.04/20.04/22.04_ubuntu20 anzhuang xiangrikui-CSDN博客 xubuntu16.04下安装向日葵并设置开机自启_./sunloginclient_11.0.1.44968_amd64.deb: 2: syntax-CSDN博客xubuntu16.04下安装向日葵并设置开机自启_./sun…

《WebKit 技术内幕》学习之七(2): 渲染基础

2 网页层次和RenderLayer树 2.1 层次和RenderLayer对象 前面章节介绍了网页的层次结构&#xff0c;也就是说网页是可以分层的&#xff0c;这有两点原因&#xff0c;一是为了方便网页开发者开发网页并设置网页的层次&#xff0c;二是为了WebKit处理上的便利&#xff0c;也就是…

数据库复试—关系数据库标准语言SQL

数据库复试—关系数据库标准语言SQL SQL&#xff1a;结构化查询语言 以教材中的学生-课程数据库为例进行SQL基础语法的复习 数据库实验环境选择SQLServer 11 关系模式 学生表Student(Sno,Sname,Ssex,Sage,Sdept) 课程表Course(Cno,Cname,Cpno,Ccredit) 学生选课表SC&#xf…

【Java】学习一门开发语言,从TA的Hello World开始

欢迎来到《小5讲堂》 大家好&#xff0c;我是全栈小5。 这是《Java》序列文章&#xff0c;每篇文章将以博主理解的角度展开讲解&#xff0c; 特别是针对知识点的概念进行叙说&#xff0c;大部分文章将会对这些概念进行实际例子验证&#xff0c;以此达到加深对知识点的理解和掌握…

【Android12】Android Framework系列---Adb和PMS安装apk源码流程

Adb和PMS安装apk源码流程 adb install命令 通过adb install命令可以将apk安装到Android系统&#xff08;注意&#xff1a;特定类型的apk&#xff0c;比如persist类型是无法通过adb安装的&#xff09; 下述命令中adb解析install命令&#xff0c;并调用Android PackageManagerS…

eNSP学习——部分VLAN间互通、部分VLAN间隔离、VLAN内用户隔离(MUX-VLAN)

MUX VLAN&#xff08;Multiplex VLAN&#xff09;提供了一种通过VLAN进行网络资源控制 的机制。通过MUX VLAN提供的二层流量隔离的机制可以实现企业内部员 工之间互相通信&#xff0c;而企业外来访客之间的互访是隔离的。 特点&#xff1a; 一、主VLAN端口可以和所有VLAN通信 二…

java基础08-Java中的null和空字符串的区别是什么?

Java中的null和空字符串的区别是什么&#xff1f; null表示一个变量没有引用任何对象&#xff0c;而空字符串表示一个字符串变量引用了一个长度为0的字符串。 null&#xff1a;表示一个变量没有引用任何对象。也就是说&#xff0c;这个变量还没有被分配内存。如&#xff1a; …

基本控件(一)

目录 控件的定义&#xff1a;使用控件的好处&#xff1a; setupUiwidget.hwidget.cppmain.cpp 控件的定义&#xff1a; QT控件是在QT框架下实现的一组用户界面元素&#xff0c;包括按钮、标签、输入框、列表框、滚动条、菜单等等&#xff0c;它们用于快速构建各种应用程序的图…

【Linux】gcc与g++的认识

gcc与g是linux下的编译器&#xff0c; 我们都知道写好了源文件gcc一下就可以得到.exe文件 那么其中会经过几个过程呢&#xff1f;也就是源文件翻译过程&#xff0c; 预处理&#xff08;预编译&#xff09;编译汇编链接 本篇文章的侧重点是围绕链接进行理解C语言标准库 目录 …

Java项目:ssm框架基于spring+springmvc+mybatis框架的民宿预订管理系统设计与实现(ssm+B/S架构+源码+数据库+毕业论文)

一、项目简介 本项目是一套ssm827基于SSM框架的民宿预订管理系统设计与实现&#xff0c;主要针对计算机相关专业的正在做毕设的学生与需要项目实战练习的Java学习者。 包含&#xff1a;项目源码、数据库脚本等&#xff0c;该项目附带全部源码可作为毕设使用。 项目都经过严格调…