Seqanswers Leaderboard Ad



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

  • Wolfgang Huber
    Dear Yohann

    thanks for the feedback. All of that behaviour is intended, and the rationale behind it is described in the DEXSeq Paper. Briefly, reads that touch multiple counting bins provide evidence for the presence of each of the bins, therefore they are counted for each of the bins. The evidence is not independent, but since the testing in DEXSeq is marginal (bin by bin), the dependence is not a problem (in the same way that the dependence in expression of different genes is not a problem for gene-by-gene testing methods). Therefore the sum of bins counts is typically larger than the gene count.

    Second, exons are split up by the preparation script into multiple parts (bins) if the GTF file has different boundaries for them. This is not always pretty, and could probably be coarse-grained in some cases (and you are welcome to do your own manual or automated curation of counting bins to this end!).

    Hope this.

    Kind regards

    Leave a comment:

  • Yohann
    started a topic DEXSeq - Counting with HT-seq at the exon level

    DEXSeq - Counting with HT-seq at the exon level


    I'm currently looking for differently expressed isoforms and got curious about the behaviour of HT-seq when counting exons. Basically, I wanted to check if those two methods would give me really close results when estimating genes expression :
    • counting at the gene level with HT-seq (HTg)
    • counting at the exon level, then summing all the exons per gene (HTe)

    The correlation is not that great (see attached pdf) and there is a global trend of higher counts from my HTe method. Some of the highlighted genes have crazy differences between the two methods :
    ensembl_gene_id value_HTg value_HTe ratio
    ENSG00000205336 21 6806 0.003231967
    ENSG00000165795 73 21996 0.003364095
    When looking in a genome browser for ENSG00000205336, I can count 21 mapping reads : it fits with HTg !
    I believe that if a read is mapping on a splicing junction, it will be counted 2 times when using HT-seq at the exon level and may explain some of the differences.

    In the first steps of the DEXSeq analysis, we have to process a GTF file (from Ensembl, for example) to obtain a GFF with "collapsed" exons from different transcripts of the same gene. For my example gene, the script "" generates really small "exonic_part", some have a length of 1bp !

    GTF for ENSG00000205336
    GFF for ENSG00000205336

    Is it expected ?
    I think that each of those exonic parts will be treated as an exon when doing the DEXSeq analysis, could it be a problem ?

    Thanks for your help!

    EDIT :
    If found that thread is pretty similar to my question :
    Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

    It seems that I can't sum the different exonic parts to estimate the gene value as a read can be counted multiple times.
    Attached Files
    Last edited by Yohann; 05-07-2014, 09:01 AM. Reason: found related thread

Latest Articles


  • seqadmin
    Recent Advances in Sequencing Analysis Tools
    by seqadmin

    The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
    05-06-2024, 07:48 AM





Topics Statistics Last Post
Started by seqadmin, Yesterday, 02:06 PM
0 responses
Last Post seqadmin  
Started by seqadmin, 05-14-2024, 07:03 AM
0 responses
Last Post seqadmin  
Started by seqadmin, 05-10-2024, 06:35 AM
0 responses
Last Post seqadmin  
Started by seqadmin, 05-09-2024, 02:46 PM
0 responses
Last Post seqadmin