Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 454 runAssembly contig lengths

    Hi

    I've made assemblies of several 454 reads (in sff and fasta format) by using two different programs: runAssembly (reads in sff format) and phrap (reads in fasta format).

    The contig length of the runAssembly program is much smaller than the one of phrap.

    Does anybody how to set the parameters/options of the runAssembly program in order to modify the contig length?
    I have tried the "notrim" option but it didn't have any effect on the contig length.

    Thanks for your advice!

  • #2
    I would be really careful in using phrap for assembly of 454 reads. The program was designed for assembly of Sanger reads, which have very different kind of errors than 454 reads. Even though the contigs may be longer in the phrap assembly, I would do some checks to see if they are correct (do you have a reference genome to compare the contigs to? PCR?).

    runAssembly (newbler) is so far (on of) the best assembly program for 454 reads (some people might disagree?). Be sure to feed it sff files, as they contain more information (flowgrams) than just fasta and quality files.

    Comment


    • #3
      We have good experience with the performance of 454 Newbler assembler. If you're doing denovo assembly however it might be a good idea to use an iterative approach. Dividing your reads into bins, assembling these separately and then assembling the contigs that you got. Especially useful if you have a bit too much coverage.

      The notrim option turns off the additional read trimming in the assembly (See the 454TrimStatus.txt that is generated with the assembly).

      Comment


      • #4
        Thank you for your kind replies!

        I appreciate any advice and suggestion because I am very new to bioinformatics and just have run my first assembly.

        >flxlex:
        I have compared the contigs to a reference genome. In some cases the phrap and in others the 454 contig is in agreement with the reference.
        Of course, I have used quality scores also for phrap (extracted from the sff file). What information does the "flowgrams" contain?

        >Tuxido:
        In another study, I am doing a de novo assembly. The sequencer made 800'000 reads. I have executed runAssembly on the whole set of reads. The output was contigs ranging upto 16'000bp. Would you still bin the reads, and if so, how many into each bin?

        Comment


        • #5
          That would depend on your genome size and your average read length (i.e. it depends on your coverage).
          So starting from I guess 50x coverage it might be worthwhile to give it a try. I would just try a few different bin settings and see which one gives you the best contigs. (I only did this once myself and used 18 bins on 500Mb of sequence data for a 4.5Mb genome, which improved assembly considerably).

          There's also a paper about this from Bas Dutilh
          http://bioinformatics.oxfordjournals...bstract/btp377 but he uses it when working with metagenomes.

          Comment


          • #6
            estimate genome size

            >Tuxido:
            Thanks for the link and advice!

            It's a de novo assemly, and I don't have any reference or knowledge of the genome size. Is there a way to estimate the genome size after a 454 run with 800k reads and an average read length of 354bp?

            Comment


            • #7
              Originally posted by DNAjunk View Post
              >flxlex:
              I have compared the contigs to a reference genome. In some cases the phrap and in others the 454 contig is in agreement with the reference.
              Of course, I have used quality scores also for phrap (extracted from the sff file). What information does the "flowgrams" contain?
              The flowgram contains the signal intensity for each flow obtained during sequencing. If you run the command

              sffinfo yourfile.sff

              you will see among the ouput something like this:

              Flowgram: 1.01 0.04 1.01 0.06 0.06 0.98 0.05 1.04
              2.33 1.16 1.10 0.06 0.21 0.89 0.07 2.00 0.12 0.08
              1.89 0.10 0.45 1.02 1.84 0.17 0.92 0.34 0.09 0.99

              With flow order TACG, this means that at the first flow, T, the signal intensity was just over 1, meaning most likely 1 T. The next flow, A, gave no signal. etc

              newbler uses this information for the assembly/quality scoring

              Comment


              • #8
                Originally posted by Tuxido View Post
                That would depend on your genome size and your average read length (i.e. it depends on your coverage).
                So starting from I guess 50x coverage it might be worthwhile to give it a try. I would just try a few different bin settings and see which one gives you the best contigs. (I only did this once myself and used 18 bins on 500Mb of sequence data for a 4.5Mb genome, which improved assembly considerably).
                Tuxido, what happens to repeated regions from the genome when you do the binning approach? How do you place them in your assembly? Newbler collapses these regions into single contigs.

                Comment


                • #9
                  I understood that with the old Newbler repeat reads where simply not used in the assembly. While with the upgrade they are now used once. I have no idea how these end up in the final contigs. We only did such an experiment once, and purely looking at contig number and contig length, binning seemed to work very well.

                  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, Yesterday, 06:53 AM
                  0 responses
                  12 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 07-10-2024, 07:30 AM
                  0 responses
                  34 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 07-03-2024, 09:45 AM
                  0 responses
                  204 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 07-03-2024, 08:54 AM
                  0 responses
                  213 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X