Question: How to obtain transcript level TPM using stringtie?
0
gravatar for kousi31
6 days ago by
kousi310
kousi310 wrote:

Hi all,

Stringtie reports TPM at gene level, but I need TPM at transcript level. Is there a way to obtain transcript level TPM using stringtie?

rna-seq stringtie • 127 views
ADD COMMENTlink modified 3 days ago by ruisergioluis20 • written 6 days ago by kousi310
2
gravatar for manuel.belmadani
6 days ago by
Canada
manuel.belmadani840 wrote:

It does provide transcript level TPM since 2015:

10/19/2015 - v1.1.0 release This StringTie release includes the following updates: major memory usage improvements due to: changes of internal data structures, collapsing reads aligned in the same place and filtering of spurious spliced alignments within large bundles -- most RNA-seq data samples use much less than 1Gb of memory now. TPM is now also reported for transcripts and genes (besides FPKM)

Also from https://ccb.jhu.edu/software/stringtie/index.shtml?t=manual:

for each RNA-Seq sample, run StringTie using the -B/-b and -e options in order to estimate transcript abundances and generate read coverage tables for Ballgown.

ADD COMMENTlink modified 6 days ago • written 6 days ago by manuel.belmadani840
1
gravatar for ruisergioluis
3 days ago by
ruisergioluis20 wrote:

The TPMs values are located in the GTF file, returned by StringTie. If you want a simple two-column file with Transcript ID and TPM, you can simply run the following command:

grep -P "\ttranscript\t"  StringTie_sampleName.gtf  | cut -f9 | awk '{gsub("\"","",$0);gsub(";","",$0);print $4,$12}'

I didn't test it but should be fine.

ADD COMMENTlink modified 3 days ago • written 3 days ago by ruisergioluis20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 560 users visited in the last hour