Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • spabinger
    Member
    • Jun 2011
    • 13

    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
  • areyes
    Senior Member
    • Aug 2010
    • 165

    #2
    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



    Alejandro

    Comment

    • spabinger
      Member
      • Jun 2011
      • 13

      #3
      Hi Alejandro,

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

      Best
      Stephan

      Comment

      • townway
        Member
        • May 2009
        • 41

        #4
        how about several treatments and also samples without replicates, does DEXSeq handle as well? Thank you
        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



        Alejandro

        Comment

        • areyes
          Senior Member
          • Aug 2010
          • 165

          #5
          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, 01:15 AM.

          Comment

          • emilyjia2000
            Member
            • May 2011
            • 59

            #6
            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

            Comment

            • Simon Anders
              Senior Member
              • Feb 2010
              • 995

              #7
              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.

              Comment

              • emilyjia2000
                Member
                • May 2011
                • 59

                #8
                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?

                Comment

                • Simon Anders
                  Senior Member
                  • Feb 2010
                  • 995

                  #9
                  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.

                  Comment

                  • Lisa
                    Member
                    • Jan 2010
                    • 10

                    #10
                    DESeq in Galaxy

                    I am wondering if it is possible to put DESeq in Galaxy? Thanks a lot! Lisa

                    Comment

                    • emilyjia2000
                      Member
                      • May 2011
                      • 59

                      #11
                      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,

                      Comment

                      • jbrwn
                        Member
                        • Mar 2011
                        • 37

                        #12
                        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".

                        Comment

                        • emilyjia2000
                          Member
                          • May 2011
                          • 59

                          #13
                          Thanks Jbrwn,
                          after I erased "i ID" from the command, same error message.

                          Comment

                          • jbrwn
                            Member
                            • Mar 2011
                            • 37

                            #14
                            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.

                            Comment

                            Latest Articles

                            Collapse

                            • seqadmin
                              Pathogen Surveillance with Advanced Genomic Tools
                              by seqadmin




                              The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
                              03-24-2025, 11:48 AM
                            • seqadmin
                              New Genomics Tools and Methods Shared at AGBT 2025
                              by seqadmin


                              This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

                              The Headliner
                              The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
                              03-03-2025, 01:39 PM

                            ad_right_rmr

                            Collapse

                            News

                            Collapse

                            Topics Statistics Last Post
                            Started by seqadmin, 03-20-2025, 05:03 AM
                            0 responses
                            49 views
                            0 reactions
                            Last Post seqadmin  
                            Started by seqadmin, 03-19-2025, 07:27 AM
                            0 responses
                            57 views
                            0 reactions
                            Last Post seqadmin  
                            Started by seqadmin, 03-18-2025, 12:50 PM
                            0 responses
                            50 views
                            0 reactions
                            Last Post seqadmin  
                            Started by seqadmin, 03-03-2025, 01:15 PM
                            0 responses
                            201 views
                            0 reactions
                            Last Post seqadmin  
                            Working...