Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • ChIP-seq duplicate reads/ Poisson distribution

    Hi everybody

    I started analyzing my first ChIP-seq data set, it contains one ChIP-sample and one input sample. After mapping the reads to a reference genome by using Bowtie and additionally MAQs, around 70% of all reads were uniquely mapped to the reference, which should be a quiet good rate (I guess).
    For the input sample ~20 mio reads were left, and for the ChIP-sample ~17 mio reads.
    However, I found ~2 mio duplicated reads (matching the same chromosomal location) in the input sample and ~16 mio duplicated reads in the ChIP-sample, which might be due to amplification errors or library preparation.
    In literature I read that the expected number of reads matching the same position and strand can be modeled by the Poisson distribution. Is this assumption also true for ChIP-samples, where we enrich specific chromosomal locations and get rid of those, where the TF doesn't bind? Wouldn't we expect to find more more duplicated reads in ChIP-samples than in the Input samples

    To identify Peaks I used MACS that removes duplicated reads before calling the peaks. Does anyone know a good peak caller that parameters concerning duplicated reads can be adjusted by the user? I want to try to set a customized threshold for the number of duplicated reads depending on my duplicated read distribution and check the sequences of my peak regions. Fortunately the TFBS motif of my TF is already known, so I can verify my results.

    It would be great to get some comments or ideas, as I am an absolute beginner in NGS analyses...

    Besides that, thank you for the great forum, it's a great help

    Thanks a lot in advance
    Kathrin

  • #2
    Hello Kathrin,

    Originally posted by kathrin View Post
    After mapping the reads to a reference genome by using Bowtie and additionally MAQs, around 70% of all reads were uniquely mapped to the reference, which should be a quiet good rate (I guess).
    You may also try other mappers (but this is another thread). I usually get ~80-90% of "good quality alignments" in ChIP-seq experiments.

    Originally posted by kathrin View Post
    In literature I read that the expected number of reads matching the same position and strand can be modeled by the Poisson distribution. Is this assumption also true for ChIP-samples, where we enrich specific chromosomal locations and get rid of those, where the TF doesn't bind? Wouldn't we expect to find more more duplicated reads in ChIP-samples than in the Input samples
    Indeed ChIP-seq samples (but also Input samples) show reads distributed with a power law, which essentially means "rich gets richer". In other words, there are sequences that have higher chances to be enriched, and the more reads you have the more will match on those regions.
    I see your point but consider also that genomic DNA (the Input) may contain some genomic features (open chromatin) that may be constitutively "enriched", plus a number of repeated sequences you may not find in the IP. Also, TF binding sites + background noise cover a slice of your genome which is probably "wide" enough to sparse your duplicates.
    You may try do zap duplicates from your alignment using picard (assuming you have BAM files)

    Originally posted by kathrin View Post
    Does anyone know a good peak caller that parameters concerning duplicated reads can be adjusted by the user?
    You may try FindPeaks4. It has a -duplicatefilter flag that can be omitted. Also, it offers a huge set of options to model your experiment.

    HTH

    d

    Comment


    • #3
      I am used to getting only 30% unique reads when doing ChIP-Seq on Arabidopsis genome. The number goes up to around 80% when all alignments are counted.

      I would be very interested if dawe could elaborate on the read distribution, because our Input DNA does not span the genome evenly and particularly shows enrichment in exons. We have been wondering why this would be so ?

      Comment


      • #4
        Hi, check this paper: Zhang et al. Modeling ChIP sequencing in silico with applications. PLoS Comput Biol (2008) vol. 4 (8) pp. e1000158
        Indeed the Input is not random nor is flat. There are regions (open chromatin, fragile sites) that may be preferentially enriched. One should try to use naked DNA to isolate sequence bias (I believe there's a paper on this... I can't get the reference right now).

        d

        Comment


        • #5
          dawe, thank you for the reference and alerting us to this.

          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