This is a test version of Biostars. For the public version, visit https://www.biostars.org.
scRNAseq - infer type of single-cell chemistry from fastq

Hello all,

it's a known thing there's lots of variants of single cell experiments - e.g. kb --list from kallisto-bustools lists at least 10 different variants. If you're processing a reasonably big selection of scRNA-seq experiments, it would make sense to try and infer it automatically.

Is there by any chance a tool that can guess the chemistry from the reads?

Thank you in advance, as always.

scrna-seq single cell rna-seq
• 1,881 views • link
updated by bosun1988 • bosun1988 • written by predeus predeus

Perhaps the only thing possible with 10x data (assuming sequencing recommendations were followed strictly) is to look at the length of reads and check.

• 0 views • link
written by GenoMax GenoMax

Thank you. 10x is easier, yes - although 5' is still tricky.

• 0 views • link
written by predeus predeus

1 answer

Hi,

I don't know if you have solved this yourself already, but I wrote this package:

https://github.com/Nusob888/fasterqParseR/

Hope it helps, at the moment it will only take "SRR" run names as input filenames (it isn't written to handle complex underscore naming).

A workaround is to remove all the underscores and append SRR at the front. It will rename those formats to 10X cell ranger names so will correct the names back if you need it cell ranger (which I assume you don't since that has an auto detect feature).

The one bug I have noticed is that citeseq assays break it. It will not return a chemistry for non-whitelist barcodes.

Additionally it will correct R1, R2 and I1 files, since I noticed for many uploads to SRA, the order is incorrect which also breaks kb-python

• 0 views • link
written by bosun1988 • bosun1988

Log in to answer this question.

More posts like this