This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Single-nucleus RNA-seq processing (10x) with Pseudoaligners (Kallisto)

I currently use CellRanger on single nucleus RNA-seq data using the pre-mRNA reference that they advise. https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/advanced/references

I want to try using Kallisto Bustools for increased processing speed. However, I haven't seen it used yet with single nucleus data, where there are a high proportion of introns, which won't be presented in the standard cDNA reference that one uses to build an index.

Is there a reference transcriptome for human that has unspliced mRNA?

Or does anyone have any references for using Kallisto with nuclear RNA?

scrna-seq rna-seq kallisto
• 5,016 views • link
updated by Lior Pachter Lior Pachter • written by achamess achamess

1 answer

There is a tutorial for how to do that here: https://github.com/BUStools/getting_started/blob/master/kallisto_bus_mouse_nuclei_tutorial.ipynb

This is part of a series of tutorials for how to use kallisto | bustools for a variety of single-cell RNA-seq pre-processing tasks: https://www.kallistobus.tools/tutorials

• 0 views • link
updated by Ram Ram • written by Lior Pachter Lior Pachter

Thanks! Can the same cDNA/intron index be used with Kallisto for processing of bulk nuclear RNA-seq too?

• 0 views • link
written by achamess achamess

Unfortunately it's not that simple- bulk RNA-seq has reads from across transcripts and so quantification requires normalization by length. Having said that, if you are only interested in producing transcript compatibility counts then yes, this index can be used.

• 0 views • link
written by Lior Pachter Lior Pachter

Log in to answer this question.

More posts like this