Question: transcript count after string tie merge
0
gravatar for rthapa
18 months ago by
rthapa0
rthapa0 wrote:

Hi, I am trying to get transcript counts with stringtie merge command. I have set of gtf output, when I checked the transcript counts of individual gtf files, I got some values of transcript counts but when I used strigtie merge command to merge all the stringtie files, I got merged output file but the transcript counts are all 0 for all the transcripts. I used the following command lines:

stringtie --merge -p 8 -o stringtie_mergedA.gtf stringtieoutputgtf.txt

where stringtieoutputgtf.txt looks like:

/scratch/user/rthapa/RNA-seq/stringtieS1Routput.gtf
/scratch/user/rthapa/RNA-seq/stringtieS2Routput.gtf
/scratch/user/rthapa/RNA-seq/stringtieS3Routput.gtf
/scratch/user/rthapa/RNA-seq/stringtieS4Routput.gtf
/scratch/user/rthapa/RNA-seq/stringtieS5Routput.gtf
/scratch/user/rthapa/RNA-seq/stringtieoutputS6R.gtf
/scratch/user/rthapa/RNA-seq/stringtieoutputS7R.gtf
/scratch/user/rthapa/RNA-seq/stringtieoutputS8R.gtf
/scratch/user/rthapa/RNA-seq/stringtieoutputS9R.gtf
/scratch/user/rthapa/RNA-seq/stringtieoutputS10R.gtf

python prepDE.py -i sample1.txt

sample1.txt looks like stringtie_mergedA.gtf /scratch/user/rthapa/RNA-seq/stringtie_mergedA.gtf

Could anyone suggest me what might have gone wrong? Thanks

rna-seq stringtie • 980 views
ADD COMMENTlink modified 18 months ago by Kevin Blighe45k • written 18 months ago by rthapa0
1
gravatar for Kevin Blighe
18 months ago by
Kevin Blighe45k
Kevin Blighe45k wrote:

stringtie --merge is only supposed to be used to create a consensus GTF reference from multiple GTF files. The idea of doing this is that it helps with downstream analyses. For example: how could you accurately do differential expression analysis on multiple samples, each with different transcript co-ordinates in their respective GTF files? With a consensus GTF, this becomes a whole lot easier.

Once you generate the consensus GTF, you can then re-estimate abundances by re-running StringTie with the -e command-line parameter on all of your samples.

Take a look at the very helpful manual: StringTie Transcript assembly and quantification for RNA-Seq

Cross-posting: transcript count after string tie merge

Kevin

ADD COMMENTlink written 18 months ago by Kevin Blighe45k
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: 1473 users visited in the last hour