Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Salmon and RSEM TPM estimate divergence?

    Hello,

    I have been exploring which transcript quantification method to use for isoform level analysis. To date I have been performing gene level analysis with a STAR-htseqcounts-edgeR pipeline. However, I am now interested in quantifying isoform levels. To this end, I have mapped with STAR as before, but allowed the maximum number of reported mappings and reported out the transcriptome covering reads with the option available in the latest version of STAR. I have then used the default and appropriate read orientation settings for RSEM, salmon, and eXpress to see how much these tools agree with a single library of about 130million paired end human reads. When I compare the Salmon NumReads with either eXpress effective counts or RSEM estimated counts I find very good agreement with high Pearson correlation, (Salmon/RSEM read count figure attached) Although RSEM and Salmon are in much greater agreement linearly (Pearson correlation ~0.97 (RSEM/Salmon) vs ~0.70 (RSEM/eXpress, Salmon/eXpress), while Salmon reads are monotonically closer to eXpress values than RSEM ones (Spearman correlation ~0.92 (Salmon/eXpress) vs ~0.82 (Salmon/RSEM). However, it seems with Salmon and RSEM that both the FPKM, and more importantly, TPM values reported (another figure attached) disagree in a way that appears to be largely dependent on transcript length. Is there something in either model that significantly changes the way TPM is calculated to bias the transcript counts based on length when the read approximations produced by the two models are so similar? I wasn't able to identify such a bias from the documentation and literature.

    I should say I do notice that RSEM tends to estimate much higher TPM than either eXpress or Salmon for transcripts around 100-200bp (read length is 51bp); however, the agreement with eXpress TPM values for anything larger than this is profound - while the Salmon and RSEM values diverge as a function of decreasing transcript length as seen in the second attached figure. Any discussion would be much appreciated.

    Warm regards,
    David
    Attached Files

  • #2
    What version of Salmon did you use? Have you posted to the google users board for sailfish/salmon? Rob is usually very good about getting answers to questions.

    I agree with you somethings seems odd with the TPM values. Assuming this was done after the STAR alignments you are then just looking at the difference in the two programs. Perhaps there is a setting that is not performing as expected or not enabled by default. We noticed some odd behavior in salmon 2.7 when using stranded libraries but it tracked back to the command line implementation.

    Comment


    • #3
      Hi seversond12,

      I just came across this post. Did you ever track down the source of these differences? I'd be interested in knowing. Also, was this salmon result with / without bias correction?

      Thanks!
      Rob

      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, Yesterday, 06:37 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      9 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      49 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      67 views
      0 likes
      Last Post seqadmin  
      Working...
      X