SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
DEXSeq - more than two different conditions yylilly RNA Sequencing 22 04-25-2019 08:00 PM
DEXseq mathew Bioinformatics 15 05-23-2014 04:11 AM
Bug? duplicated genes in cufflinks output genes.expr silin284 Bioinformatics 3 05-17-2014 11:19 PM
DEXSeq error:Failed to set up model frames for genes zhuya5607 Bioinformatics 6 08-22-2012 01:34 AM
aggregate read cycle info from bam files osenan Bioinformatics 0 06-28-2012 06:43 AM

Reply
 
Thread Tools
Old 06-17-2013, 04:38 PM   #1
leda
Junior Member
 
Location: California

Join Date: Feb 2013
Posts: 7
Default dexseq aggregate genes

Hi,

I'm trying to understand the "-r no" option in the latest version of dexseq_prepare_annotation.py. It seems that the script only ignores an exonic part if it exactly overlaps an exonic part from another gene - is this correct? Couldn't gene-level differential expression also disrupt differential exon usage calls if an exonic part of one gene is a subset of an exonic part of another?

thanks,
Leda

----------

E.g. here are some lines from running the script on gencode 14 annotation. The first exonic part contains the subsequent ones and more (from a different gene), so I'm wary when DEXSeq tells me it's used differentially.

chr19 dexseq_prepare_annotation.py aggregate_gene 2252252 2269758 . - . gene_id "ENSG00000167476.5"
chr19 dexseq_prepare_annotation.py exonic_part 2269408 2269758 . - . transcripts "ENST00000593238.1"; exonic_part_number "012"; gene_id "ENSG00000167476.5"
...
chr19 dexseq_prepare_annotation.py aggregate_gene 2269519 2273487 . + . gene_id "ENSG00000104904.6"
chr19 dexseq_prepare_annotation.py exonic_part 2269519 2269519 . + . transcripts "ENST00000583542.2+ENST00000582888.2"; exonic_part_number "001"; gene_id "ENSG00000104904.6"
chr19 dexseq_prepare_annotation.py exonic_part 2269520 2269528 . + . transcripts "ENST00000583542.2+ENST00000322297.4+ENST00000582888.2"; exonic_part_number "002"; gene_id "ENSG00000104904.6"
chr19 dexseq_prepare_annotation.py exonic_part 2269529 2269529 . + . transcripts "ENST00000583542.2+ENST00000581150.1+ENST00000582888.2+ENST00000322297.4"; exonic_part_number "003"; gene_id "ENSG00000104904.6"
...
leda is offline   Reply With Quote
Reply

Tags
dexseq

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 12:21 AM.


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