Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • De-novo assembly of bacteria genomes - which tools?

    Hi all,

    I have been given the task of assembling three bacterial genomes.

    Each have been sequenced with 454 with a read length around 250nt. But also with Solexa mate-pair (2.5kb space - 35nt reads) because the initial assembly of the 454 data gave too many contigs (using the Newbler assembler resulting in a couple of thousand contigs per genome - these being GC rich).

    Which assembly software would be appropriate? And are there any recommendations on how to proceed with assembly having two types of data?

    Cheers,

    Martin

  • #2
    Originally posted by maasha View Post
    Which assembly software would be appropriate? And are there any recommendations on how to proceed with assembly having two types of data?
    You could try out release candidate 1 of MIRA 3 (or wait for rc2 which will be put online this weekend): http://chevreux.org/projects_mira.html (comes with extensive docs)

    On top of my head: Velvet and Euler-SR will also work with 454 and Solexa. I certainly missed some. Have also a look at: http://seqanswers.com/forums/showthread.php?t=43

    Regards,
    B.

    PS: I'm a bit biased as I wrote MIRA

    Comment


    • #3
      Thanks BaCh,

      So where to start? The Solexa Data or the 454? It occurs to me that the 454 - with the longer reads - are best for initial assembly - and then gaps can be closed with the mate-pair Solexa data.

      And can you feed the resulting contigs of one of these assemblies into the software as reads?


      Martin

      Comment


      • #4
        Originally posted by maasha View Post
        So where to start? The Solexa Data or the 454? It occurs to me that the 454 - with the longer reads - are best for initial assembly - and then gaps can be closed with the mate-pair Solexa data.
        If you can't feed both read types at the same time, 454 will indeed give you a better initial assembly. Though things may change slowly with 76mers and soon the double.

        Originally posted by maasha View Post
        And can you feed the resulting contigs of one of these assemblies into the software as reads?
        Most of the time not: contigs are normally pretty long sequences where assemblers have problems with (they're currently built to expect reads up to 1-2k bases).

        Try feeding the assembler with all reads at first. Compare what looks best to you. If that does not work at all, the resort to assembling Solexas, shred resulting contigs and feed those shreds alongside 454 reads to "normal" assemblers (Newbler, CABOG, MIRA, and others that grok 454).

        At least that's what I'd do.

        HOWEVER ...

        First, I am skeptical that the Solexas will greatly reduce contigs in your project. First, having "thousands" of contigs for a small bacterium with 454 data is *extremely* unusual (unless your average coverage is a low single digit number). Something doesn't feel right.

        Second, you wrote "GC" rich. Which might get problematic as Solexa has problem sometimes with GGC.G motifs. What I've often seen is this:

        refseq .....GGCGGCGGCxxxxxxxxGCCGCCGGC.......

        Now, the GGC motif in forward and reverse direction pretty much leads to a complete depletion of correct Solexa reads (they all have errors in the bases marked with x).

        B.

        Comment


        • #5
          Try Velvet http://www.ebi.ac.uk/~zerbino/velvet/

          has support for using both short (illumina) and long (454) reads in the assembly.

          Comment


          • #6
            @BaCh,

            Thanks,

            First, I am skeptical that the Solexas will greatly reduce contigs in your project. First, having "thousands" of contigs for a small bacterium with 454 data is *extremely* unusual (unless your average coverage is a low single digit number). Something doesn't feel right.
            Remember that the Solexa data is mate pair. So the question is also - is there any software that will take both single reads and mate pair reads as input?

            @ewilbanks,

            Thanks, velvet is on top of my list of tools to try (along with MIRA of cause).


            Martin

            Comment


            • #7
              Originally posted by maasha View Post
              Remember that the Solexa data is mate pair. So the question is also - is there any software that will take both single reads and mate pair reads as input?
              Any software able to treat paired-end should be very well able to handle single reads. It needs to anyway for cases where a mate is not present or unusable.

              B.

              Comment


              • #8
                Now I have battled extensively to assemble one bacterial genome sequenced with 454 and Solexa (mate pair).

                454 data: 478840 reads covering 117689186nt with a mean length of 246 and 33% GC content (the other two genomes are the ones with high GC content).

                Solexa data: Mate pair reads (distance ~2500nt) of 35nt and 44nt length in total 3308974 reads covering 130704473nt and GC content of 38% (can anyone explain this bias?).

                Expected genome size: 2.8Mb

                Using the Newbler assembler on the 454 data alone results in ~50 contigs.

                Using Velvet on both data type I could at the best get down to 459 contigs using the below:

                velveth test_both 31 -long M1_454.fna -shortPaired M1_solexa_paired.fna
                velvetg test_both/ -ins_length 2500 -cov_cutoff 10 -exp_cov 41 -max_coverage 300

                This is not too good (


                (Installing MIRA is not going well because of the libboost requirement - can't libboost be included in MIRA?)



                Martin

                Comment


                • #9
                  Originally posted by maasha View Post
                  Using the Newbler assembler on the 454 data alone results in ~50 contigs.
                  ...
                  This is not too good (
                  I beg to differ: this is already not too bad for a 40x+ coverage. Now, having a 3MB data set that assembles into 700 contigs ... that's "not too good"

                  Originally posted by maasha View Post
                  (Installing MIRA is not going well because of the libboost requirement - can't libboost be included in MIRA?)
                  Won't happen: including an almost standard library with a program is not a good idea, especially when all major distributions have it included. You could give the binary packages a try though, or don't they run on your machine?

                  Regards,
                  B.

                  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
                  25 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 07-16-2024, 05:49 AM
                  0 responses
                  39 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 07-15-2024, 06:53 AM
                  0 responses
                  45 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 07-10-2024, 07:30 AM
                  0 responses
                  42 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X