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
                  Best Practices for Single-Cell Sequencing Analysis
                  by seqadmin



                  While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
                  06-06-2024, 07:15 AM
                • seqadmin
                  Latest Developments in Precision Medicine
                  by seqadmin



                  Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                  Somatic Genomics
                  “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                  05-24-2024, 01:16 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, Yesterday, 08:58 AM
                0 responses
                9 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-12-2024, 02:20 PM
                0 responses
                15 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-07-2024, 06:58 AM
                0 responses
                182 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-06-2024, 08:18 AM
                0 responses
                231 views
                0 likes
                Last Post seqadmin  
                Working...
                X