Tweet
I saw this in the latest issue of nature biotech: http://www.nature.com/nbt/journal/va.../nbt.3122.html. Here is the site for it: http://ccb.jhu.edu/software/stringtie/. Has anyone used this on their data yet? If so what did you think? I'm trying to estimate isoform abundance in my RNA-seq data.
-
-
It looks pretty good from the paper, and I'd definitely give StringTie a shot, but I've not tested the software yet. However, if you're looking to try out a new tool, I'd also take a look at the Bayesembler (http://genomebiology.com/2014/15/10/501). It seems to have a very nice and principled probabilistic model that it uses to jointly infer what transcripts are present and quantify their abundance.Last edited by robp; 02-26-2015, 05:46 PM. -
Bayesembler is nice, and also works with large quantities of data rapidly.
I haven't successfully tried Stringtie yet, think an older version aborted with errors.Comment
-
Hi folks,
Sorry for bumping an old thread but I am facing this weird output from StringTie and thought you could spot the issue.
I am running StringTie on STAR returned BAM files (I have used the option to insert XS tag in SAM/BAM fields)
I have now run StringTie on 6 samples and in each I have the same scenario -
FPKM value as "inf" or "-nan"
StringTie returns a GTF file and this an example of the 9th col for one transcript -
About 98% of the transcripts have FPKM as "inf" and rest fraction as "-nan"gene_id "STRG.67"; transcript_id "STRG.67.1"; reference_id "ENST00000447513"; ref_gene_id "ENSG00000157911"; ref_gene_name "PEX10"; cov "1.525000"; FPKM "inf";
This is the cmd I am using
I have checked the STAR BAM files and they seem ok: >85% reads with unique alignment./path/to/stringtie-1.0.4/stringtie \
Aligned.sortedByCoord.out.bam \
-o StringTie_out/StringTie.gtf \
-p 16 \
-G /path/to/Homo_sapiens.GRCh37.73_Only-REF.gtf \
-B \
-e
Any hunches? I am sure I am doing something stupid but can't figure out what.Comment
-
Also to add, for the same BAM (STAR returned) on which Cufflinks successfully returned FPKM values, StringTie is giving me "Inf" and "-nan"
:-(Comment
-
hi everyone,
I figured out the stupidity and that was -not reading the manual properly :-(
The Manual page says -
..If you run StringTie, and all FPKM values that were produced appear to have the value inf then it is very likely that the HI tag starts at 1. An easy command that you can run...
So STAR out BAM files have the HI tag starting from 1 and hence.Comment
-
Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
...07-09-2026, 11:10 AM -
Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.
There is no single reason why many patients don’t respond to treatment as expected. Cancer is...07-08-2026, 05:17 AM -
Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.07-01-2026, 11:43 AM
Comment