Thank you. I will try your suggestion and let you know.
• 0 views
•
link
written
by
Naresh D J •
How to separate htseq-count table to coding RNAs and non-coding RNAs
Hi,
I have generated the raw read counts for genes from RNA-seq data using htseq-count. Now I want to separate the this table into coding RNAs and non-coding RNAs.
I am new to the NGS data analysis.
Can anyone help me or suggest me ideas how to do it?
Thank you,
Naresh
• 6,229 views
•
link
updated
by
Ram
•
written
by
Naresh D J •
•
written
by
Naresh D J •
2 answers
What is your organism model?
If you are using some genome from ensembl, and used the gtf file with the set of anotations in the HTSeq-count, you can import all the tables with counts in txt files inside a data.frame in R.
With bioconductor do:
biocLite("biomaRt")
With this package you can get from the ensembl, a dataset based in several filters, for example, the biotype (if it is coding or noncoding).
Then you can simple merge the two tables based in the ensembl ID's and separate them based in your criteria. If you do not want to use R, the ensemble has a graphic web interface in http://www.ensembl.org/biomart, although I recommend R, because will be more easy later to create better graphics and statistics.
Links:
P.S.: biomaRt also can handle Uniprot and HapMap databases
• 0 views
•
link
updated
by
Ram
•
written
by
tiago211287
•
written
by
tiago211287
To do so, you need a file with a relation (range of bases) of the sequences that are coding and not coding. Mapping reads to the reference genome or transcriptome is not aware of this information
• 0 views
•
link
written
by
Antonio R. Franco
Log in to answer this question.
More posts like this
-
Can I convert HTSeq count into RPKM or TPM value or standard unit of RNA-Seq
written by wathunya.h •
-
Which RNA type I should considered as long non-coding RNAs
written by zizigolu
-
HT-Seq count data coding gene
written by Rob •
-
Generating read counts for my microRNA reads for a non model invertebrate
written by karthic
-
Inserting Noise into RNA-seq raw counts
written by saamar.rajput
-
differentially expression analysis using raw read counts
written by ashkan
-
Where to download processed time course RNA-seq data in human?
written by Naresh D J •
-
Help with DEseq2 design matrix and DE expression analysis
written by Naresh D J •
-
How to run the BLAST program with the un-mapped reads in bam format from tophat
written by Naresh D J •
-
Why cufflinks does not report raw read count anymore? Which version of it do report?
written by bibilong1111 •
What you mean by coding and non-coding RNAs? Do you mean separating counts for coding and non-coding transcripts ? Or do you mean separating counts for coding (exonic) and non-coding (intronic, UTRs) regions for a given transcript?
@Ashutosh Pandey, yes I want to separate the counts for coding and non-coding transcripts.
For separation of coding and non-coding regions there is a tool RSeQC.
Well RSeQC will give you the numbers or fractions of reads aligned to different genic features but it won't separate them. Anyways, what you need is the annotation of transcripts (genes) based on their biotypes. If these are ENSEMBL genes or gene IDs then you can use Biomart (http://www.ensembl.org/biomart) to download the "Biotype" for each gene and then annotate ENSEMBL genes in the count file as protein-coding, rRNA, tRNA, snoRNA, miRNA etc.
Thank you. I will try your suggestion and let you know.