Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • ChIP seq-% alignment

    Hi all:

    Just wondering after running the ChIP-seq data with Pipeline which by default uses ELAND for alignment, what kind of % aligment is acceptable, indicating the experiment was working fine? We got 40-60% aligment, is it reasonable? Thanks!

  • #2
    @SPGSling

    What length reads are you aligning and against what reference?

    Cheers
    Roald

    ** Disclaimer: I work at CLC bio **

    Comment


    • #3
      Our read length is 36 bp and the reference genome is human. Thanks!

      Comment


      • #4
        @SPGSling

        Are you allowing the algorithm to align reads with multiple matches?
        If this dramatically increases the number of aligned reads it would indicate problems with overly amplified repeat regions.
        Have you looked at the contigs with the aligned reads and do the data seem to make sense?

        Cheers

        Roald

        ** Disclaimer: I work at CLC bio **

        Comment


        • #5
          We just run it with pipeline by default using ELAND to align. We are currently checking if the NM are repeats, primer dimers or other source of contamination. And we are also mapping the seq to the genome to look for peaks as well. Just want to get an idea what kind of %alignment people usually get as an indication of a good ChIP-seq experiment. Thanks a lot!

          Comment


          • #6
            40-60% using the default ELAND settings is a little on the low side for 36bp (paired end?). If you want to try alternate settings, bowtie (see the forum here) is very amenable to messing around with settings and also extremely fast!

            Have you looked at your primer dimer occurance rate? i.e. look through your sequences for the illumina adapter sequence reverse complement (the read you would get if for an adapter/primer dimer)

            Comment


            • #7
              I thought ChIP-seq was supposed to have a lower alignment rate than e.g. genomic DNA? We've just done some ChIP experiments and whilst we are very much still optimising the protocols our alignment was ~54% to human genome. The amount of reads containing adapter-dimer sequences (or slightly truncated ones) was <<<1%.

              Anyone else willing to share their alignment percentages for human ChIP alignments?

              Comment


              • #8
                Our Illumina data are from 1G, with 30-60 % aligned reads but this number does not at all indicate how succesful the ChIP is.

                Comment


                • #9
                  We saw chip seq on human data with about 40-50% mapped reads using eland as well..
                  Is there a quick way to check for adapter-dimer? Do you do it with eland itself?
                  --
                  bioinfosm

                  Comment


                  • #10
                    We were able to sort out all the reads generated by pipeline 1.4 (default setting) for our ChIP seq experiments. This is the breakdown of what we got:

                    % alignment 60-70%
                    % repeats 20-30%
                    % no match 3-6%

                    So now we know since ELAND exclude all the reads aligned to more than 2 sites. Most of the non-aligned reads in our experiment was repearts. We did get some no-match, we are in the process to verify if it's primer dimer or other source of contaminants.

                    Chipper, you are absolutely right. % aligned reads does not at all indicate how succesful the ChIP is! We use it as a quality assessment of our ChIP library preparation. How good was our ChIP, we have to wait till we can validate our peaks.

                    Comment


                    • #11
                      Originally posted by bioinfosm View Post
                      We saw chip seq on human data with about 40-50% mapped reads using eland as well..
                      Is there a quick way to check for adapter-dimer? Do you do it with eland itself?
                      hi bioinfosm,

                      To find out the % of adapter dimers in our reads i count the number of reads containing the sequences GATCGGAAGAGCTCGTATGCCGTCTTCTGCTTG or
                      GATCGGAAGAGCTCGTATGCCGTCTTCT in the s_N_export.txt file using the grep command: [grep -c "GATCGGAAGAGCTCGTATGCCGTCTTCT" s_N_export.txt] and then count the total number of reads [grep -c "" s_N_export.txt].

                      very quick and dirty..

                      Comment


                      • #12
                        Originally posted by elaney_k View Post
                        I thought ChIP-seq was supposed to have a lower alignment rate than e.g. genomic DNA? We've just done some ChIP experiments and whilst we are very much still optimising the protocols our alignment was ~54% to human genome. The amount of reads containing adapter-dimer sequences (or slightly truncated ones) was <<<1%.

                        Anyone else willing to share their alignment percentages for human ChIP alignments?
                        Hi!
                        In a ChIP-Seq we made months ago we had 64/65% IP aligned to hg18 as from Eland Summary.htm, or 50/51% from the ShortRead report.
                        The latter is the percentage related to the total reads, while the first is the percentage of the reads passing the filters.
                        In that experiment we had a 70% of reads passing the filtering (eland default, pipeline 1.3) out of ~ 14M of clusters per lane.
                        We hadn't extrapolated the quota of non-aligned reads that are actually repeats, but we will...hopefully.

                        Cheers
                        gabriele bucci

                        Comment


                        • #13
                          In case someone else reads this, this could be of use - FastQC (http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/) analyzes your fastq dataset and will report overrepresented sequences (among other useful info). A few times I've had datasets have an adapter contamination and FastQC reported it. It even identified which adapter - in the case I'm looking at now, TruSeq Adapter 3. It was <1.1%. I'm not sure what threshold makes a sequence overrepresented though.

                          Comment


                          • #14
                            Originally posted by captainentropy View Post
                            In case someone else reads this, this could be of use - FastQC (http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/) analyzes your fastq dataset and will report overrepresented sequences (among other useful info). A few times I've had datasets have an adapter contamination and FastQC reported it. It even identified which adapter - in the case I'm looking at now, TruSeq Adapter 3. It was <1.1%. I'm not sure what threshold makes a sequence overrepresented though.
                            FastQC is a good tool. We use it routinely as a command line pipe.
                            With the option --casava it also recognizes the several fastq.gz splitted by the new sequencer and output one single full sample QC.


                            /usr/bin/fastqc --casava -t 10 ${name}_R?_???.fastq.gz

                            HTH
                            gabriele bucci

                            Comment

                            Latest Articles

                            Collapse

                            • seqadmin
                              Exploring the Dynamics of the Tumor Microenvironment
                              by seqadmin




                              The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
                              07-08-2024, 03:19 PM
                            • seqadmin
                              Exploring Human Diversity Through Large-Scale Omics
                              by seqadmin


                              In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
                              06-25-2024, 06:43 AM

                            ad_right_rmr

                            Collapse

                            News

                            Collapse

                            Topics Statistics Last Post
                            Started by seqadmin, 07-19-2024, 07:20 AM
                            0 responses
                            29 views
                            0 likes
                            Last Post seqadmin  
                            Started by seqadmin, 07-16-2024, 05:49 AM
                            0 responses
                            42 views
                            0 likes
                            Last Post seqadmin  
                            Started by seqadmin, 07-15-2024, 06:53 AM
                            0 responses
                            50 views
                            0 likes
                            Last Post seqadmin  
                            Started by seqadmin, 07-10-2024, 07:30 AM
                            0 responses
                            43 views
                            0 likes
                            Last Post seqadmin  
                            Working...
                            X