Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • MeDIP-seq versus ChIP-seq analyses?

    I've shortened my previous post in hopes of getting a response:

    - I am mapping MeDIP sequence reads (Illumina GAII) with Bowtie and using MACS for peak finding.
    - Question 1: Should only unique reads be mapped?
    - Question 2: How does the pattern/shape of peaks compare between MeDIP-seq and ChiP-seq?

    Thanks.

    jjw
    Last edited by jjw14; 08-20-2010, 11:44 AM. Reason: No responses. Shortened post.

  • #2
    Bumping shorter version of original post. Apologies if this is bad form. Just looking for a bit of advice. Much appreciated.
    jjw

    Comment


    • #3
      I would think that only unique reads should be mapped because if you map reads which match to multiple locations then you would not be sure about any resulting patterns.

      Comment


      • #4
        MeDIP-Seq is just another form of ChIP-Seq so in crude terms you can use the same techniques as you'd apply to ChIP-Seq to analyse MeDIP data. What I suspect you really want to know if if you should expect that the peaks you see in MeDIP should look like nice small independent events that many of the peak finders expect.

        In our experience we see two kinds of enrichment in MeDIP data, there are small, highly enriched events which are very easy to identify and then there are slightly more diffuse, lower enrichment events which can be harder to pick out. The majority of reads will sit in the highly enriched peaks, and looking at these is a good place to start, but there may be some value in looking at the lower levels of enrichment.

        Comment


        • #5
          Thanks very much. Your description of the "independent peaks" or peak "frequency" for MeDIP versus ChIP is exactly the type information I was interested in.

          I'm using MACS for peak detection. Do you feel there is a more appropriate algorithm for MeDIP peak detection, especially for the areas of low enrichment you describe?

          Much appreciated,
          jjw

          Comment


          • #6
            Thanks for the response. I appreciate it.
            jjw

            Originally posted by mattanswers View Post
            I would think that only unique reads should be mapped because if you map reads which match to multiple locations then you would not be sure about any resulting patterns.

            Comment


            • #7
              We do our analysis within SeqMonk (which we write, so we're somewhat biased!) rather than using an external peak caller.

              Normally we'd use the contig based probe generator to do the cluster detection. This allows you to set an enrichment threshold and build clusters where your reads exceed this threshold. There are then facilities for merging adjacent peaks or rejecting enriched regions which are too short. We find that this kind of approach lends itself well to more unusual ChIP experiments where the shape of the peaks might not be what you'd expect from a traditional 'transcription factor' type enrichment.

              Comment


              • #8
                Great information. I had actually downloaded the SeqMonk dmg to my Mac about two weeks ago, but wasn't sure if I had enough RAM (4 GB) to run it, so I defaulted to using the available tools on our campus Linux cluster. I see the SeqMonk readme file describes how to change the memory limits in OS X. I may also see if I can get SeqMonk installed on our Linux cluster.

                I'll try the sample data file provided on your web site with some of the cluster detection methods you describe. Thank you for providing this resource, and I'll gladly acknowledge your assistance in any publication/report that (hopefully!) is produced.

                jjw

                Comment


                • #9
                  The standard SeqMonk install assumes a 32-bit machine with 2GB RAM. You can configure it to use more than this, but for modestly sized experiments (up to 10 lanes or so) this shouldn't be necessary. The software is designed to run on desktop machines - I wouldn't see much advantage to installing it on a node of a cluster.

                  Comment

                  Latest Articles

                  Collapse

                  • 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
                  • seqadmin
                    The Impact of AI in Genomic Medicine
                    by seqadmin



                    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                    02-26-2024, 02:07 PM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, 03-14-2024, 06:13 AM
                  0 responses
                  32 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-08-2024, 08:03 AM
                  0 responses
                  71 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-07-2024, 08:13 AM
                  0 responses
                  80 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-06-2024, 09:51 AM
                  0 responses
                  68 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X