Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

bustools workflow for snRNA-seq detects many more genes/UMIs than cellranger #57

Open
jrknoedler opened this issue Apr 29, 2020 · 0 comments

Comments

@jrknoedler
Copy link

I'm a big fan of kallisto for bulk RNA-seq so decided to establish a kallisto | bustools workflow for my single nucleus RNA-seq data (admittedly partially because my university's cluster start throwing OOM-kill events when I was running cellranger, which had worked fine previously, for no obvious reason). I followed the tutorials for building a cDNA index by downloading the files directly from the Pachter lab's github repository and running the kallisto bus allignment, then bustools capture and count for the intron and cDNA alignments following the tutorial. So far, so good. However, there's an odd result - I'm getting something like 30% more unique genes/cell with Kallisto/bus than with cellranger-3.1 (something like 4,500 versus 3,000). Now, I don't want to complain, especially since this has already resulted in my being able to find markers that previously weren't showing up, but I also go by the credo that if something sounds too good to be true it usually is. Is this likely to be due to pseudoalignment to a cDNA/intron index performing better than cellranger? Has anyone seen anything like this before?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant