SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
quantifying RNAseq reads for a genbank isoforme that is not in Reference Annotation ftmir Bioinformatics 0 02-11-2016 02:43 AM
plant annotation pipeline for assembled RNAseq? vebaev RNA Sequencing 27 08-19-2014 12:43 AM
Removing overlapping genes from annotation for RNAseq read count DRAT Bioinformatics 2 04-11-2014 03:53 AM
ERANGE-RNASeq with ensemble annotation Matthias Bioinformatics 0 02-23-2011 05:05 AM

Reply
 
Thread Tools
Old 01-23-2017, 10:46 AM   #1
Dani
Member
 
Location: Bay Area

Join Date: Dec 2016
Posts: 13
Default RNA class annotation RNASeq

Hi,

I am trying to to this pretty common "quantification" representation of different RNA classes in my RNA-Seq data set (meaning mRNA, tRNA, rRNA, ncRNA etc, like for example here http://journalofextracellularvesicle...#CIT0025_26760), but I can't seem to find any tools that do that - can anybody help with that?

Thanks,

Danielle
Dani is offline   Reply With Quote
Old 02-15-2017, 10:44 AM   #2
Dani
Member
 
Location: Bay Area

Join Date: Dec 2016
Posts: 13
Default

Btw. found that biotypes in ngi_visualization does the job if anybody has the same question=)
Dani is offline   Reply With Quote
Old 02-16-2017, 12:03 AM   #3
Michael.Ante
Senior Member
 
Location: Vienna

Join Date: Oct 2011
Posts: 120
Default

Hi Danielle,

I hope it's not to late for a suggestion. The Genecode and ENSEMBL annotations include in their gtf files per gene a field called "gene_biotype".
If you use htseq-count (for instance), you can switch from counting the reads per gene_id to counting the reads per gene_biotype:

Code:
 htseq-count -f bam -r pos -i gene_biotype my_alignment.bam my_genecode.gtf
Before doing so, check if the gtf and your bam file base on the same assembly (e.g. hg19) and have the same naming for the chromosomes' name.

Cheers,

Michael
Michael.Ante is offline   Reply With Quote
Old 02-16-2017, 08:02 AM   #4
Dani
Member
 
Location: Bay Area

Join Date: Dec 2016
Posts: 13
Default

Hi Michael,

that's a great suggestion, thanks - unfortunately, it complains:
Error occured when processing GFF file (line 8 of file /mnt/users/dga
Feature ENSG00000223972.5 does not contain a 'gene_biotype' attribu
[Exception type: ValueError, raised in count.py:53]
I am using: gencode.v25.primary_assembly.annotation.gtf which has the same assembly and naming as my bam.
Does that mean that the gtf files have to be "cleaned" up to not contain any empty biotypes before using?

Thanks,

Danielle
Dani is offline   Reply With Quote
Old 02-16-2017, 01:04 PM   #5
Dani
Member
 
Location: Bay Area

Join Date: Dec 2016
Posts: 13
Default

Sorry, realized that it should be -i gene_type =)
Dani is offline   Reply With Quote
Reply

Tags
rna class, rna-seq analysis

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 11:48 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO