Question: Abundance Counting for Non-Model Organism for edgeR
gravatar for sdbaney
25 days 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 • 104 views
ADD COMMENTlink modified 25 days ago by ATpoint17k • written 25 days ago by sdbaney0
gravatar for ATpoint
25 days ago by
ATpoint17k 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 25 days ago by ATpoint17k

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 25 days 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: 1119 users visited in the last hour