Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Is it OK to use paired-end 150bp to sequence ATAC-seq libraries?

    Hi!

    I am planning to do ATAC-seq assays and I am wondering if it is fine to use paired-end 150bp to sequence ATAC-seq libraries?
    Just because the sequencing platform I am working with is only doing 150bp run with HiSeq3000...
    I have the impression that people are usually using PE sequencing with 50 to 75bp reads. Would the only problem be that I would sequence a lot of adapters (exept money lost, adapter trimming should be enough to solve this)? Or would I lost many "small" fragments?

    Thank you for your help!
    Vero

  • #2
    Don't use large insert libraries for ATAC-seq

    Hi,

    This is probably too late to be of any help to you, but for anyone else doing ATAC-seq, I highly recommend doing paired-end sequencing with 50x50bp or lower.

    Even with adapter trimming, I have had a very difficult time mapping a 76x76bp PE Illumina NextSeq ATAC-seq library using bowtie2. There is always a dip in mapped reads around the 76bp insert size (see attached example plot), and this has subsequently caused NucleoATAC to have trouble generating its models for nucleosome positioning analysis.

    And by have trouble, I mean the program crashes when running "nucleoatac run ..." making that analysis method unusable. There was an issue created in the NucleoATAC GitHub repo that discusses this problem a little bit (https://github.com/GreenleafLab/NucleoATAC/issues/62).

    Chris
    Attached Files

    Comment


    • #3
      You could always hard-clip your reads to 50 bases.

      --
      Phillip

      Comment


      • #4
        Originally posted by pmiguel View Post
        You could always hard-clip your reads to 50 bases.
        Hmm... that is something I hadn't thought of actually. I am going to give that a shot to see how the insert size distribution looks after that. I'll hard clip to 50 bases, then use trim_adapters (https://github.com/ParkerLab/atactk/.../trim_adapters) to trim, but throw out reads smaller than 30bp after trimming.

        I'd be interested if any other people had similar issues with an ATAC-seq library like this, and what solutions they would recommend.

        I did try various adapter trimmers on the 76x76bp data and compared the mapping results. For example, I tried trim_adapters, trimmomatic and cutadapt both with and without the trimming done by bcl2fastq (so 6 different possibilities total). Using bcl2fastq trimming, and then trim_adapters seemed to work the best for me, but there wasn't much of a difference.

        Most recently I tried remapping with very slow bowtie2 settings to see if that would help. For a Mus musculus sample with about 60 million paired-end reads, it took about 5 days to finish mapping (12 threads, ~2.6GHz per), but it did actually help a bit by getting me 8-10% more uniquely mapping concordant pairs, and more 70-80bp insert size reads to map.

        For reference, these are the bowtie2 settings I used: "-k 4 --dovetail --no-mixed -D 20 -R 3 -N 1 -L 8 -i S,1,0.50". Mostly the same as "--very-sensitive-local", just lowering the seed size (-L), and allowing for a mis-match (-N 1).

        --very-sensitive-local == -D 20 -R 3 -N 0 -L 20 -i S,1,0.50

        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
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        49 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