Question: RNA-Seq DEGenes Without GTF
0
gravatar for isu2017
20 months ago by
isu20170
isu20170 wrote:

I am working a specific species where we have a reference genome, and I was able to align our RNA-seq data to it. However, I am looking to create a count table to funnel our dataset to edgeR for DE analysis. I don’t have a GTF to use with something like HTSeq. Anybody have an idea of how to proceed?

rna-seq • 523 views
ADD COMMENTlink modified 20 months ago by RamRS27k • written 20 months ago by isu20170
2

It's a species that has a reference genome but that is not well annotated in terms of gene expression? If that is the case, I would have done:

  1. HISAT2 / StringTie to produce a reference transcriptome GTF
  2. Use HTseq with the new GTF to derive raw counts
  3. Feed raw counts into EdgeR
ADD REPLYlink written 20 months ago by Kevin Blighe59k
1
  1. ...
  2. ...
  3. ...
  4. extract the unannotated gtf features and blast them to get an annotation.
ADD REPLYlink modified 20 months ago • written 20 months ago by h.mon29k

Valeu, cara!

ADD REPLYlink written 20 months ago by Kevin Blighe59k
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: 1189 users visited in the last hour