Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • MACS: ChIP and control with different tag size

    Hi,
    I am using MACS2 for ChIP-seq peak calling.

    I have a treatment file with original reads of 51bp. Before alignment, the reads were dynamically trimmed on 3'end for removing low quality bases (with Trimmomatic software). At the end, aligned reads ranges from 42 to 51 bp.

    Then, I have a control (input) file with original reads of 36 bp. Before alignment, the reads were dynamically trimmed on 3'end for removing low quality bases (with Trimmomatic software). At the end, aligned reads ranges from 30 to 36 bp.

    If I let MACS to automatically determine the tags size, MACS comes out with a wrong number. Moreover, since I think MACS scan only the first n reads, MACS does not give the average of the tags length.

    So, I have to specify the tag size with the option --tsize. However, it is impossible to set a tsize for the treatment and a different tsize for the input.

    What do you think is the best solution for this situation?

    Thanks very much!
    Albert

  • #2
    Why didn't you sequence the IP and input at the same time? I'm guessing that you also didn't construct the libraries at the same time and therefore have an uncontrolled batch effect completely confounding your results.

    Comment


    • #3
      In the real world, sometimes you have to use old data with new data, or data that comes from different labs.. any idea?

      Comment


      • #4
        Even if I fully agree with dpryan, it happened to me that researchers really wanted to use their old datasets with a new one. Except that in my case, both samples were done with the same chromatin then fragments sizes were similar. I just trimmed the longest reads to fit the shortest...
        Hope both samples (treatment and control) were obtained with the same chromatin (I mean same sonicated chromatin, not only from same celltype...)

        If not, maybe the best would be to use MACS without control for both samples and then you crossed the results yourself...

        Comment


        • #5
          Yes SylvainL, I am exacly in that situation: my wets wanted to use their old data.

          Is there any paper that I can give to my wets, where is clearly stated that treatment and input should comes from the same chromatin?

          Thanks guys for your help!

          Comment


          • #6
            They need a paper for that? They would have the exact same problem with their wet lab experiments. It'd be like taking an image of a Western blot from a paper and comparing the raw band intensities to one of your own Westerns. In fact, just use that example, since they should understand it.

            Comment


            • #7
              They told me that they think that they can make only one input from the same cell type, and use it as a reference input for all the future experiments from the same cell type..

              Comment


              • #8
                Funny, so they must think that they only ever need to run 1 Western as a baseline for comparisons then too...

                Comment


                • #9
                  So typical when wet-biologists do not understand a technique... I had to deal with some of these quite recently. Then, we finished the discussion by me telling them: if you want to perform only one Input, then prepare a huge batch of chromatin... They didn't like this comment when they realized they would have to prepare 1000 plates together

                  If they don't understand the example of dpryan (they may argue that it depends of the exposure of the blot, etc...), simply ask them to summarize all the Ct they obtained for their control genes in qPCR (all same cell type of course) let's say the last 4 years....

                  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
                  33 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-08-2024, 08:03 AM
                  0 responses
                  72 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-07-2024, 08:13 AM
                  0 responses
                  81 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