What is the correct way to import transcript abundances from 3' Quantseq results that calculated by kallisto?
2
0
Entering edit mode
duhwa.lee • 0
@e491bd42
Last seen 12 weeks ago
Austria

Hi,

I am planning to import the transcript abundances from tsv files that were calculated by kallisto using Quantseq 3' mRNA seq data.

But I don't know what is a correct argument for tximport function for Quantseq results, especially for countFromAbundance.

I assume it will be "no", but I don't know well clearly. Could you explain it to me?

Best,

DH

Txi_gene <- tximport(path,
type = "kallisto",
tx2gene = Tx,
txOut = FALSE, #determines whether your data represented at transcript or gene level
countsFromAbundance = "no",
ignoreTxVersion = TRUE)

tximport tximportData Quantseq kallisto • 211 views
3
Entering edit mode
ATpoint ★ 2.1k
@atpoint-13662
Last seen 11 hours ago
Germany

Yes, option no is recommended for end-tagged data, see vignette https://bioconductor.org/packages/release/bioc/vignettes/tximport/inst/doc/tximport.html#3%E2%80%99_tagged_RNA-seq and this previous question tximport of 3 tagged RNA-seq counts from Salmon for downstream analysis with DESeq2

0
Entering edit mode

0
Entering edit mode
@james-w-macdonald-5106
Last seen 10 hours ago
United States

You probably don't want to use countsFromAbundance, because the transcript length is probably less relevant than for normal bulk RNA-Seq. But that is just a wild guess. You might consider aligning with STAR or whatever conventional aligner you like, and then look at the count distribution across the transcripts. If it's all just within say the first 100 nt of each transcript, then the transcript length is probably not relevant any longer, and using any transcript length information is likely to introduce bias.