Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Differential expression analysis on RPKMs - Contrasts and contrasts of contrasts

    Dear all,

    I'm looking for a tool (better would be a Bioconductor package) able to perform differential expression analysis on a table of RPKMs. I've seen so far only solutions implementing raw read counts, like DEGSeq, edgeR and BaySeq. And yeah, in this particular case I cannot transform RPKMs back to read counts.

    Do you have any idea? Furthermore, a great advantage would be the capability of calculating significance of "contrasts of contrasts" analysis, i.e. of tetrafactorial designs.

    E.g. having 4 conditions: WT control, WT treated, mutant control, mutant treated.

    Code:
    DE1 = WT treated vs. WT control
    DE2 = mutant treated vs. mutant control
    
    DE1DE2 = (mutant trated vs. mutant control) vs. (WT treated vs. WT control)
    Something like the great limma does for microarrays. But on RNASeq RPKMs.
    Thanks a lot for any hint!

  • #2
    Why can't you extract the raw read counts?

    Comment


    • #3
      It's unfortunately an issue sometimes, since in the beginning RPKM (Mortazavi et al., 2008) were considered a standard and some softwares (e.g. Cufflinks) by default output only FPKM and pairwise condition comparisons. There's an open (more like a closed now) debate on whether one should use RPKM or raw read counts for differential analysis, the latter being the prime choice for statistical reason (see for example the manuals of edgeR and the wiki of this website).
      As for your question, given only a table of RPKM, you can try the NOISeq R script (http://bioinfo.cipf.es/noiseq/doku.php?id=tutorial), if you can manage to convince it that your data is already rpkm-normalized. Good luck!

      Comment


      • #4
        Any valid statistical analysis which only get RPKM values as input will have very little power, because lacking information on whether a certain count value stems from a long gene with many counts or a short gene with few counts, it is bound to assume that all numbers stem from low-count and hence high-variance genes.

        I'd be careful with NOISeq, by the way. In the paper on NOISeq, the method is evaluated only on dataset that contain no replicates or only technical replicates, i.e., it is not shown that this method is able to correctly assess biological variability and hence control type-I error rate in a biologically meaningful way.

        Comment


        • #5
          NOIseq fails

          in my hands NOIseq on RAW read counts failed to get significant DE genes called (i.e. only VERY few called), similar to tophat/cuffdiff (only 2 DE genes called) if comparing 12 normal vs 12 colon cancer samples. DESeq, edgeR, BaySeq, SAMseq worked much better...

          with edgeR and limma (design matrix) you can analyse your design (wt-mock, wt-treat, mut-mock, mut-treat), but you need raw reads...

          Comment


          • #6
            Thanks everyone for your super-informative answers. I guess I will have to convince my collaborators to move towards raw read counts then.

            Comment


            • #7
              Check out a recent comparison of methods:

              Nookaew et. al. (2012), Nucleic Acids Research, Vol. 40, 10084-10097.

              NOISeq does not lack power and compares well with other methods as well as microarrays.

              “Although all the methods were in agreement in identifying significant GO terms related to growth (a consequence of the increased specific growth rate during batch cultivations), GO terms known to be relevant during fully respiratory growth are not all in agreement with the different methods. Specifically, edgeR showed some inconsistencies in capturing GO terms associated with fatty acid beta-oxidation terms (as well as DESeq), fatty acid metabolic process and TCA cycle, whereas baySeq weakly identify increased expression of ATP-coupled proton transport and ion transport. Interestingly, the results derived from NOISeq seem to give stronger signals that explain the known differences between batch and chemostat growth better than the results derived from the other methods.”

              The authors speculate that edgeR suffers from increased Type I error.

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Strategies for Sequencing Challenging Samples
                by seqadmin


                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                03-22-2024, 06:39 AM
              • seqadmin
                Techniques and Challenges in Conservation Genomics
                by seqadmin



                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                Avian Conservation
                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                03-08-2024, 10:41 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 03-27-2024, 06:37 PM
              0 responses
              12 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-27-2024, 06:07 PM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              53 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              68 views
              0 likes
              Last Post seqadmin  
              Working...
              X