Hi
Does anybody know if htseq-count works on UTR? I tried it, only got gene list, there is no count in the output.
Does anybody know if htseq-count works on UTR? I tried it, only got gene list, there is no count in the output.
-
-
Read the instructions for htseq-count:
-t <feature type>, --type=<feature type>
feature type (3rd column in GFF file) to be used, all features of other type are ignored (default, suitable for RNA-Seq and Ensembl GTF files: exon)
By using the -t, --type= argument you can tell what features you want htseq-count to use. The default is exon, so you cant use the defaults. The question then is do you want reads mapping only to the UTR or to the entire transcript, including the UTR? That will determine what feature type you give it.
So say you want reads mapping only to UTRs. Whatever your UTRs are named in your GFF file, use that with the -t or --type= arguments when you run htseq-count. -
Hi chadn737,
I used -t feature in the command line and set feature to utr, please see the command i used:
htseq-count -s no -t utr -m intersection-nonempty file.sam utr.gff > urtcount.txt
the output looks like:
"06100073560ik"
"06100071803ik"
"0610008608Rik"
"0610008008Rik"
"0610008007Rik"
"0610008007Rik"
"0610009760Rik"
"06100094371ik"
"06100094166ik"
"06100094858ik"
"0610009020Rik"
"06100091744ik"
Only gene_id is there and no count. BTW, we used NCBI reference and modified it according to ensembl gff file format.
THanks,Comment
-
Look in your gff file to find out specifically how the UTRs are labeled. I'm not familiar with your specific species/gff file. For instance in the TAIR10 GFF they are named as three_prime_UTR or five_prime_UTR. Assuming your gff file even has the UTRs listed? You may need to modify the names using sed or alternatively count reads for 5' UTRs and 3' UTRs seperately.Comment
-
Hi Chand737,
The GFF used had been modified according to ensembl format as following
1 bed2gff utr5 8794025 8794051 . - . transcript_id "NM_27671"; utr5_number "0016"; gene_id "NM_27671"
1 bed2gff utr5 8872242 8872295 . - . transcript_id "NM_27671"; utr5_number "0017"; gene_id "NM_27671"
1 bed2gff utr5 8989267 8989330 . - . transcript_id "NM_27671"; utr5_number "0018"; gene_id "NM_27671"
1 bed2gff utr5 9193260 9193341 . - . transcript_id "NM_27671"; utr5_number "0019"; gene_id "NM_27671"
1 bed2gff utr5 9288214 9289811 . - . transcript_id "NM_27671"; utr5_number "0020"; gene_id "NM_27671"
total 9 columns. each columns separated by tab. Starting from transcript_id to gene_id is one column (the last column).
I am using mouse genome. 5'UTR
Thanks,
LiLast edited by emilyjia2000; 04-06-2012, 06:51 AM.Comment
-
Why does your output from your last post:Originally posted by emilyjia2000 View Post
"06100073560ik"
"06100071803ik"
"0610008608Rik"
"0610008008Rik"
"0610008007Rik"
"0610008007Rik"
"0610009760Rik"
"06100094371ik"
"06100094166ik"
"06100094858ik"
"0610009020Rik"
"06100091744ik"
Not match the gene_ids of your gff file?
That suggests to me that the problem has to do with naming of something in one of your files.Last edited by chadn737; 04-06-2012, 11:10 AM.Comment
-
Sorry I have two annotation files, one is the posted, another is the one with geneid as listed in the output. two files have the same format.
Sorry I didn't make it clear.Comment
-
It looks like you edited this post while or after I replied. In this gff file your feature type is listed as "utr5" not utr.Originally posted by emilyjia2000 View Post
So try using -t utr5 instead of -t utr or alternatively --type=utr5.Comment
-
Something else to check. At the very end of your output from htseq-count there should be the following categories:
no_feature
ambiguous
too_low_aQual
not_aligned
alignment_not_unique
How many reads fall into these categories? Assuming that everything in your sam and gff files is right, and the only problem is the feature type you give to htseq-count, then you should see read counts in these categories. Particularly since the majority of your reads are unlikely to be coming from UTR regions.
If these categories are empty, then something else is wrong.Comment
-
Thanks Chand737, I tried both features: -t utr or -t utr5, both of them delivered the same results.
As you pointed out, I checked the category as follows:
no_feature 122533035
ambiguous 1417743
too_low_aQual 0
not_aligned 0
alignment_not_unique 20439904
Could you please let me know what's wrong with it?
ThanksComment
-
Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?
Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating...Yesterday, 09:55 AM -
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment