Question: Abundance Counting for Non-Model Organism for edgeR
gravatar for sdbaney
3 months ago by
sdbaney0 wrote:

I am working with a non-model organism (no annotated genome or the option to use an annotation file in my abundance counting). I performed de novo transcriptome assembly of the Illumina HiSeq4000 reads.

I originally used StringTie as my abundance counter because I didn't need a gtf annotation file and it created one I could use for the differential expression analysis but I am now seeing that because of its feature to predict isoforms, I am having trouble matching up their IDs with my assembly IDs (there being more stringtie IDs than assembly because of its novel prediction).

Is there any way to disable this feature in StringTie? Is there another program I can use to create a gtf file so that I can use featureCounts or just count abundances that will not require a gtf?

edger abundance rna-seq • 191 views
ADD COMMENTlink modified 3 months ago by ATpoint22k • written 3 months ago by sdbaney0
gravatar for ATpoint
3 months ago by
ATpoint22k wrote:

You can use salmon to quantify reads against a transcriptome in fasta format. From there on you can summarize transcript count to the gene level with tximport.

ADD COMMENTlink written 3 months ago by ATpoint22k

Thank you, I think that’s what I’m going to end up doing. As stringtie seems to be causing more headache than anticipated.

ADD REPLYlink written 3 months ago by sdbaney0
Please log in to add an answer.


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