Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Julien Roux
    Member
    • Dec 2011
    • 24

    DEXSeq gene level counts

    Dear DEXSeq authors,
    I have a question related to the following paragraph in the DEXSeq vignette:
    Note that a read that mapped to several exons of a gene is counted for each of these exons by the dexseq_count.py script. The table returned geneCountTable will hence count the read several time for the gene, which may skew downstream analyses in subtle ways. Hence, we recommend to use geneCountTable with care and use more sophisticated counting schemes where appropriate.
    I am wondering if this problem is also affecting the tests for differential expression of exons in DEXSeq? I have the feeling that the reads overlapping multiple exons are also counted several times when it comes to estimate gene expression levels (to control for overall differences in expression at the gene level when testing exons). Is this accounted for in the package functions?

    Thanks for your help
    Julien
  • Simon Anders
    Senior Member
    • Feb 2010
    • 995

    #2
    It's not really an issue because if a gene is expressed more strongle in one sample than in another, all exon counts, including those that result from double-counting due to overlaps, will be scaled up in a proportional manner, so that the expression ratio between the samples is still estimated correctly.

    Comment

    • Julien Roux
      Member
      • Dec 2011
      • 24

      #3
      Hi Simon,
      Thanks for your answer. However I am not fully convinced that it is not an issue...

      In my dataset at least 40% of the reads overlap an exon junction. If you count these read twice when estimating gene baseline levels, it can lead to major estimation errors for the expression ratio between conditions (especially when there is a large difference between conditions).

      An easy fix would be to provide DEXSeq the correct gene count table. Is there a way to do this and use it in the function "testGeneForDEU"?
      Thanks for your help
      Julien
      Last edited by Julien Roux; 11-27-2012, 03:03 PM.

      Comment

      • Simon Anders
        Senior Member
        • Feb 2010
        • 995

        #4
        As DEXSeq performs a hypothesis test, we only need to ensure that the "null" case is not affected. The null hypothesis in DEXSeq is that the relative proportions that the differnt isoforms contribute to overall gene expression does not depend on the condition, while the overall expression strength may change.

        If now a gene is expressed twice as strongly in one condition than in the other, all per-exon counts go up by a factor of two. The number of reads that fall onto two exons goes up by two, as well, so that the ratio of each per-exon count to the sum of all per-exon counts stays as constant. This is all that we need. There is no need for the sum of all per-exon counts to reflect gene expression strength faithfully.

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