Question: Abundance Counting for Non-Model Organism for edgeR
gravatar for sdbaney
12 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 • 350 views
ADD COMMENTlink modified 12 months ago by ATpoint35k • written 12 months ago by sdbaney0
gravatar for ATpoint
12 months ago by
ATpoint35k 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 12 months ago by ATpoint35k

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 12 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: 1927 users visited in the last hour