Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • No peak in BBNorm kmer-frequency histogram

    Hi,

    I'm working on de novo assembly of an insect genome. Our paired-end libraries were made from 10 ng of sheared DNA using a Rubicon ThruPLEX kit with 9 cycles of PCR. The bioanalyzer trace shows a mean insert size of 476 bp, and we sequenced around 120 million 125 base read-pairs from this library. We're expecting a genome size of around 500 Mbp (but that's a pretty rough guess).

    I'm having trouble getting contiguous assemblies. Among others, I've tried edena (L50 = 287 bp as reported by BBTools stats.sh) and velvet (L50 = 482 bp). I'm new to de novo assembly but 9 cycles of PCR sounds like a lot to me, and I'm wondering if the library complexity is too low. Also, fastqc reports a GC content of around 30 % which I know can exacerbate PCR bias.

    To troubleshoot, I'm looking at the before and after kmer-frequency histograms generated by BBNorm during normalisation (below). I can't see a peak in either histogram, but I'm not sure what that means. Can anyone help me interpret these plots or suggest further troubleshooting steps?

    In case it's relevant, the processing I did before assembly is: quality trimming (Q < 30 at 3´ end) and adaptor trimming (TruSeq indexed adaptor and TruSeq universal adaptor) with cutadapt; contaminant filtering using PhiX, sequencing_artifacts and adapters_no_transposase references (BBDuk); normalisation and error correction to target k-mer coverage of 57 (BBNorm).

    Please let me know if any more information would help.

    Thanks for reading,

    Tom
    Attached Files

  • #2
    It looks to me like low library complexity due to overamplification. It's hard to say, though. Note that there is a peak at ~16x. In order to see it, you have to rotate the image to the left by 45 degrees... the fact that it is such a weak peak indicates a very wide spread of coverage, indicative of overamplification, or an extremely high error rate. Or severe contamination, which can also be a problem with low-input amplified libraries. Do the reads generally BLAST to related insects?

    Comment


    • #3
      Hi Brian, thanks for the reply. Contamination sounds quite possible, I just BLASTed a random subset of the reads and got human, macaque, trees, zebrafish etc. as well as the occasional hit on other insects. Uh oh.

      Our server is going offline tonight but I'll do a more systematic investigation tomorrow and post the results.

      Comment


      • #4
        I blastn-ed 1000 R1 and 1000 R2 reads from this library against the 'nt' database. For R1, I got 544 hits with an evalue < 1. 493 of them had usable taxon identifiers. From that I got 89 plant hits (18%), 82 mammalian (17%), 65 insects (13%), 63 nematodes and 57 fish (and some other stuff). R2 numbers were similar.

        I don't know if evalue is the best way to look at BLAST results for NGS reads (i.e. short queries), but either way it looks like contamination to me.

        Thanks for the hint.
        Last edited by TomHarrop; 01-10-2017, 07:39 PM. Reason: more concise

        Comment


        • #5
          I'd suggest to do the contamination analysis is a more systematic way using biobloomtools with a couple of the different top hit plants, mammals, insect, ... genomes you got from Blast. It will probably take some time, but I'd be rather surprised if you really have so many different contaminations - as long as the person doing your library preps isn't also a dedicated gardener or fisherman

          Comment


          • #6
            Yes, to be honest, this does sound strange. Normally, contamination comes from 1 or 2 sources... a grab-bag of taxa is very unusual. Are you getting 100% identity to anything, or just weak hits?

            Comment


            • #7
              Originally posted by Brian Bushnell View Post
              Yes, to be honest, this does sound strange. Normally, contamination comes from 1 or 2 sources... a grab-bag of taxa is very unusual. Are you getting 100% identity to anything, or just weak hits?
              And in addition: Do you get your complete read seq aligned or are your hits rather the tiny 20-40bp local alignment crap Blast may output if there is nothing more suitable?

              EDIT: Just saw that you used Blastn. I'd suggest megablast here.
              Last edited by WhatsOEver; 12-22-2016, 05:54 AM.

              Comment


              • #8
                Thanks for the replies. Sorry about the slow response, I missed the email notification over the holidays.

                You're correct, the hits are mostly less than 60 bp, not the full read. I did try megablast but I don't get any hits (well, 11 out of 1000 reads had hits, about half to insects).

                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
                12 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
                68 views
                0 likes
                Last Post seqadmin  
                Working...
                X