Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • ChIP-Seq Integration & Interpretation

    Hi. I have three ChIP-Seq data for the same sample. Peak calling was done separately. How to combine the peaks of these three experiments together?

    One more question, How to associate genes with these binding sites?

    Thanks in advance!

  • #2
    these are technical replicates, right? So you could simply do the peak calling on the combined data. If you want to do the peak calling first the easiest - maybe not the best - way is to define an overlapping set of peaks.

    About the association, there are some tools around that can help you (i guess Galaxy). On the other hand that is not a too complicated job, isn't it? Calculating the distances to some gene features and stuff.

    Comment


    • #3
      I think this is what you are looking for:


      There are some other ideas on my blog at:
      Everything Ethan knows about biology - The Ethan-ome


      On the first question, I agree combining the data is a good way to go. You can also get the common peaks between the samples by using the BedTools intersectBed program. It is also available on Galaxy if you are not comfortable with the command line.
      --------------
      Ethan

      Comment


      • #4
        Are the replicates purely technical (mutliple lanes of the same library) or was the experiment repeated three times (three different libraries)? This makes a big difference. If it is just three lanes of the same library, then I agree the data should be merged and peaks called once. However without replicating the experiment, you will never be able to get reliable results -- you have no clue as to the variation inherent in the experiment, and whether your library represents an outlier.

        If the experiment was repeated three times, then you should not combine the reads, but rather work with the three peaksets. Start by checking their overlaps (3-way venn diagram). To minimize false positives, you could take only the peaks that all three replicates agree on, or be less stringent and take those that overlap in at least two of the three libraries. Another approach is to take all the peaks, then go back to the aligned reads and re-score each peak for each replicate (whether or not it was identified as a peak in that replicate).

        See also this paper for ideas on how to combine peaksets: http://projecteuclid.org/DPubS?servi...oas/1318514284

        If you are using R and Bioconductor, there are a number of packages that can help with this type of analysis. For example DiffBind helps with overlap and re-scoring, and ChIPpeakAnno will map peaks to genes.
        Last edited by rory; 12-02-2011, 07:51 AM. Reason: typo

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Advancing Precision Medicine for Rare Diseases in Children
          by seqadmin




          Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
          12-16-2024, 07:57 AM
        • seqadmin
          Recent Advances in Sequencing Technologies
          by seqadmin



          Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

          Long-Read Sequencing
          Long-read sequencing has seen remarkable advancements,...
          12-02-2024, 01:49 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 12-17-2024, 10:28 AM
        0 responses
        26 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 12-13-2024, 08:24 AM
        0 responses
        42 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 12-12-2024, 07:41 AM
        0 responses
        28 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 12-11-2024, 07:45 AM
        0 responses
        42 views
        0 likes
        Last Post seqadmin  
        Working...
        X