SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to analyze RNA-seq data ? bioR RNA Sequencing 3 09-12-2013 08:08 AM
what is the best software to use to analyze illumina NGS data? bbsinfo Illumina/Solexa 3 08-21-2013 05:48 AM
Help to analyze Illumina HiSeq 2000 Human data kiradi Bioinformatics 4 12-09-2011 05:30 AM
New Tools to analyze Next Gen Seq data JMFA Bioinformatics 0 12-02-2010 09:26 AM
any one use grape to analyze solexa data zhidong Bioinformatics 11 05-22-2009 12:07 AM

Reply
 
Thread Tools
Old 08-08-2011, 06:11 AM   #1
spabinger
Member
 
Location: Europe

Join Date: Jun 2011
Posts: 13
Default DEXSeq - able to analyze paired (matched) data?

Hi everyone,

does anybody know if DEXSeq is able to analyze paired (matched) data - e.g.: same patient before and after treatment?

If not, is any software capable of handling the analysis of differential exone usage in RNA-Seq paired data (don't confuse it with paired end)?

Thanks in advance,
Stephan
spabinger is offline   Reply With Quote
Old 08-24-2011, 10:25 AM   #2
areyes
Senior Member
 
Location: Heidelberg

Join Date: Aug 2010
Posts: 165
Default

Hi Stephan,

Yes it is possible by changing the formula at the time of estimating dispersions and testing for differential exon usage. For more information check the vignette

http://bioconductor.org/packages/2.9...doc/DEXSeq.pdf

Alejandro
areyes is offline   Reply With Quote
Old 08-25-2011, 01:37 AM   #3
spabinger
Member
 
Location: Europe

Join Date: Jun 2011
Posts: 13
Default

Hi Alejandro,

thanks a lot for your reply. I will have a look at it.

Best
Stephan
spabinger is offline   Reply With Quote
Old 09-06-2011, 01:32 PM   #4
townway
Member
 
Location: Rockville

Join Date: May 2009
Posts: 40
Default

how about several treatments and also samples without replicates, does DEXSeq handle as well? Thank you
Quote:
Originally Posted by areyes View Post
Hi Stephan,

Yes it is possible by changing the formula at the time of estimating dispersions and testing for differential exon usage. For more information check the vignette

http://bioconductor.org/packages/2.9...doc/DEXSeq.pdf

Alejandro
townway is offline   Reply With Quote
Old 09-07-2011, 02:12 AM   #5
areyes
Senior Member
 
Location: Heidelberg

Join Date: Aug 2010
Posts: 165
Default

About several treatments, you can add them to the design data frame at the time of creating the ExonCountSet object. At the time of testing the analysis would tell you which exons are being affected due to changes on the treatment, nevertheless you would not know immediately which one the treatment that is causing the differential exon usage, then probably make pairwise comparisons. You could also adapt the formulas to test for the changes being introduced due to a certain treatment.

About the samples with no replicates, DEXSeq shares the same "motivation" as DESeq or edgeR. Meaning, use biological replicates to estimate biological variation using a negative binomial distribution. Though, you could use DEXSeq: the dispersion estimates will tend to 0 and you will have a lot of false positive results, but it does NOT really makes sense to make these kinds of analysis without biological replicates and a proper estimation of variation. Simon Anders and others have discussed lots about this in this forum!

Last edited by areyes; 09-07-2011 at 02:15 AM.
areyes is offline   Reply With Quote
Old 02-06-2012, 11:18 AM   #6
emilyjia2000
Member
 
Location: usa

Join Date: May 2011
Posts: 59
Default

I am new to DEXseq, I just wonder if I don't have replicates and have only one condition, is it possible to use DEXseq to count exon level reads?
Thanks,
emily
emilyjia2000 is offline   Reply With Quote
Old 02-06-2012, 12:17 PM   #7
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

Quote:
Originally Posted by emilyjia2000 View Post
I am new to DEXseq, I just wonder if I don't have replicates and have only one condition, is it possible to use DEXseq to count exon level reads?
DEXSeq is a tool to test for differential exon usgae, so using it with just a single sample is rather pointless. DEXSeq does come with a Python script to count reads for exons. However, I wonder what use the script's output might be for you.
Simon Anders is offline   Reply With Quote
Old 02-07-2012, 08:05 AM   #8
emilyjia2000
Member
 
Location: usa

Join Date: May 2011
Posts: 59
Default

Hi Simon,

THanks for the reply. I searched the tools that would help to generate exon expression level reads. for example, HTseq, Rsamtools, DEseq and DEXseq, and I am so confused and don't know which one I should choose to use in the analysis.
I have 20 RNA-seq samples (patient data), no replicates on only one condition, but we would like to know the exon level reads on each sample. Do you have any recommendation for the tools?
emilyjia2000 is offline   Reply With Quote
Old 02-07-2012, 08:27 AM   #9
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

There are many tools which take a list of genomic intervals (e.g., a list of exon positions) and count how many reads fall into each of the intervals. DEXSeq comes with such a script, GenomicRanges contains function for it, BEDtools might work etc.

I only really cannot fathom why you would want such an analysis. You will end up with a huge table with 20 columns and a few hundred thousand rows. Clearly, you cannot inspect all these rows by simply looking through them, as, I hope, you realize.

Hence you must have some idea what your next step will be, or what you hope to find in the data. The people in this forum might be much more able and willing to help you if you explained what your next planned step will be.
Simon Anders is offline   Reply With Quote
Old 02-08-2012, 02:20 AM   #10
Lisa
Member
 
Location: Ohio

Join Date: Jan 2010
Posts: 10
Default DESeq in Galaxy

I am wondering if it is possible to put DESeq in Galaxy? Thanks a lot! Lisa
Lisa is offline   Reply With Quote
Old 02-08-2012, 09:21 AM   #11
emilyjia2000
Member
 
Location: usa

Join Date: May 2011
Posts: 59
Default

Hello Simon,

Thanks for the suggestion, we would like to see transcript level, junction reads as well, not just exon level reads. But we are sort of stuck at the exon level reads. We used tophat to do the mapping.
Now I used HTseq-count to do the exon reads. it generated the error message when running it:
Message:
/usr/local/python-2.7/bin/htseq-count: Error: Please provide two arguments.
Call with '-h' to get usage information.

I am not very clear what the 2 arguments mentioned in the Error message?
command line to use:
htseq-count –t exonic_part –s no i ID -m intersection-nonempty myfile.sam Homo_sapiens.GRCh37.65_chr.gff > exoncount.txt

Do I have to sort/reorder sam files before working on htseq-count?

THanks,
emilyjia2000 is offline   Reply With Quote
Old 02-08-2012, 10:20 AM   #12
jbrwn
Member
 
Location: Denver, CO

Join Date: Mar 2011
Posts: 37
Default

Quote:
Originally Posted by emilyjia2000 View Post
Hello Simon,

Thanks for the suggestion, we would like to see transcript level, junction reads as well, not just exon level reads. But we are sort of stuck at the exon level reads. We used tophat to do the mapping.
Now I used HTseq-count to do the exon reads. it generated the error message when running it:
Message:
/usr/local/python-2.7/bin/htseq-count: Error: Please provide two arguments.
Call with '-h' to get usage information.

I am not very clear what the 2 arguments mentioned in the Error message?
command line to use:
htseq-count t exonic_part s no i ID -m intersection-nonempty myfile.sam Homo_sapiens.GRCh37.65_chr.gff > exoncount.txt

Do I have to sort/reorder sam files before working on htseq-count?

THanks,
the script doesn't like that you don't have a dash (-) before the "i".
jbrwn is offline   Reply With Quote
Old 02-08-2012, 11:04 AM   #13
emilyjia2000
Member
 
Location: usa

Join Date: May 2011
Posts: 59
Default

Thanks Jbrwn,
after I erased "i ID" from the command, same error message.
emilyjia2000 is offline   Reply With Quote
Old 02-08-2012, 11:40 AM   #14
jbrwn
Member
 
Location: Denver, CO

Join Date: Mar 2011
Posts: 37
Default

Quote:
Originally Posted by emilyjia2000 View Post
Thanks Jbrwn,
after I erased "i ID" from the command, same error message.
bummer.

for exon counting i have reverted to using there supplied script within the dexseq package (dexseq_count.py) which depends on htseq.

i assume you already have the package and used one script to write the flattened gff file. the scripts are located in:
DEXSeq/inst/python_scripts

the options are outlined in the script and account for paired-end reads, strandedness, and quality.
jbrwn is offline   Reply With Quote
Reply

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 06:40 AM.


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