Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • MNase-seq data analysis

    Hi,

    Has anyone ever played with the MNase-seq datasets?
    I got some mouse datasets running but there seems several points should take into consideration:
    1. repetitive region: as about 50% of genome is repetitive sequence and these can arise those multiple alignments during the mapping stage. Generally, only uniquely mapped reads were filtered out for the downstream analysis, which means we lost about 50% region's information of the genome. That's critical if we want to see how about the repeat regions happened in the cell. I don't know if anyone did similar analysis and added the multi-reads and how to process them.
    2. input signal: in general ChIP-seq data, the control might generates some peak signal and one of the reason arise this is due to the open chromatin structure. But in this case, say, if I notice the input sharp signal, still, that might be due to the open chromatin in that region, yet might also be generated due to other unknown reason. Has anyone also crush similar problem and how to process them?
    3. peak callor choice: there are several peak callers for the mononucleosome calling, I'm just wondering whether MACS can also accomplish this if I'm not only interested in the mononucleosome peaks.

    Welcome any kinds of communications and suggestions.
    Thanks!

    lix

  • #2
    Hello there

    Originally posted by lix View Post
    2. input signal: in general ChIP-seq data, the control might generates some peak signal and one of the reason arise this is due to the open chromatin structure. But in this case, say, if I notice the input sharp signal, still, that might be due to the open chromatin in that region, yet might also be generated due to other unknown reason. Has anyone also crush similar problem and how to process them?
    If you have sharp signals in your input you should find them in IP too. BTW, MNase input should be naked DNA, which is not the same input as in general ChIP-seq data

    Originally posted by lix View Post
    3. peak callor choice: there are several peak callers for the mononucleosome calling, I'm just wondering whether MACS can also accomplish this if I'm not only interested in the mononucleosome peaks.
    If you google for that you may find something. BTW, I'm having good results with DANPOS which is still unpublished. An alternative is NPS. Problem is that you'll need a lot of custom scripts to extract some general information (i.e. phase, distograms, TSS depletion...)

    HTH

    d

    Comment


    • #3
      Originally posted by dawe View Post
      Hello there



      If you have sharp signals in your input you should find them in IP too. BTW, MNase input should be naked DNA, which is not the same input as in general ChIP-seq data



      If you google for that you may find something. BTW, I'm having good results with DANPOS which is still unpublished. An alternative is NPS. Problem is that you'll need a lot of custom scripts to extract some general information (i.e. phase, distograms, TSS depletion...)

      HTH

      d
      Hi, in an attempt to use Danpos using 2 bed tag alignment files as input, I am unable to get it to work correctly, we get the following error:

      python2.7 danpos.py tagAlign1-tagAlign2 -k 1
      danpos version 2.0.1
      command:
      python danpos.py tagAlign1-tagAlign2 -k 1

      Namespace(bg=None, clonalcut=1e-10, count=None, distance=100, edge=0, extend=80, fs=None, gapfill=0, height=5, keep=1, lmd=300, mafrsz=250, mifrsz=50, name='result', nor='F', paired=0, path='chr22.Gm12878.nuc.tagAlign-chr22.K562.nuc.tagAlign', pcfer=0, smooth_width=20, span=10, statis='P', testcut=1e-05, width=40)
      time elasped: 0.36073589325 seconds

      normalizing wigs ...
      less than 2 datasets, no normalization to be done
      saving normalized wigs ...
      time elasped: 0.360930919647 seconds

      pooling each group ...
      Traceback (most recent call last):
      File "danpos.py", line 125, in <module>
      runDANPOS()
      File "danpos.py", line 122, in runDANPOS
      smooth_width=args.smooth_width,paired=args.paired)
      File "functions.py", line 262, in danpos
      pooledgroups[groupname]=groups[groupname].pop(filenames[0])
      IndexError: list index out of range

      Am I missing something in terms of being able to tell DANPOS that these are the two files I want to compare nucleosomal positioning data on?

      Please help, urgent

      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, Yesterday, 06:37 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      51 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      67 views
      0 likes
      Last Post seqadmin  
      Working...
      X