Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Genome assembly and Rna-seq mapping

    I every one, is my first time here but I follow for many months. My question is if someone have experience with genome assembly and post mapping of ran-seq data.
    I use velvet, celera and clc to make de novo assembly of relative difficult bacteria, and i have relative goods statistics with velvet and celera (revised with published draft genomes). Then I check the assemblies mapping rna-seq data of the same organism and I obtain lows mapping percentages (near 50%) but when I check clc assembly this up to 90%. I try to optimise the clc assembly but if I drecrease the data or make more stringent trimming, the assembly decrease in quality.

    I use 250pb PE in miseq plataform.

    Best regards !
    Cristian.

  • #2
    Hi,

    for bacteria, I was using Edena http://http://www.genomic.ch/edena.php with good results

    Maybe you could give a try...
    Last edited by SylvainL; 10-31-2014, 05:38 AM.

    Comment


    • #3
      Originally posted by SylvainL View Post
      Hi,

      for bacteria, I was using Edena http://http://www.genomic.ch/edena.php with good results

      Maybe you could give a try...
      Thxs for reply. I will try. I obtain high rna-seq reads mapping with a5 pipeline, but i can't increase genome assembly statistics.

      Comment


      • #4
        Originally posted by SylvainL View Post
        Hi,

        for bacteria, I was using Edena http://http://www.genomic.ch/edena.php with good results

        Maybe you could give a try...
        I try Edena but i can't get the same length of my reads, If you can help me I appreciated.

        Comment


        • #5
          Hi,

          what do you mean by not getting the same length of your reads? In my case, when the Qscores were bad at the end of the reads, I had to trimm them... You can easily do this step with fastx toolkit (fastx_trimmer in this case) or seqtk (trimfq)...

          Let me know at which step you are stuck...

          Comment


          • #6
            Originally posted by SylvainL View Post
            Hi,

            what do you mean by not getting the same length of your reads? In my case, when the Qscores were bad at the end of the reads, I had to trimm them... You can easily do this step with fastx toolkit (fastx_trimmer in this case) or seqtk (trimfq)...

            Let me know at which step you are stuck...
            The error is when start edena:

            Rapid file(s) examination... 158 220
            [err] All reads within a file must be the same length.

            I make pre-processing with bbmap. Maybe the problem is paired end data?

            Comment


            • #7
              I am also having the same error "[err] All reads within a file must be the same length". I am not sure I understand the rationale behind this. Once the raw reads e.g. 100 or 150nt length are trimmed for quality, adapters etc., read length becomes variable.

              Comment


              • #8
                For bacterial assemblies SPAdes (http://bioinf.spbau.ru/spades) should be in your list of programs to try.

                Comment


                • #9
                  I try spades, and really can't obtain good assembly in my specific data (I try with other bacteria data and I had got good results). This is the principal reason why I try others assemblers.

                  Comment


                  • #10
                    Originally posted by fahmida View Post
                    I am also having the same error "[err] All reads within a file must be the same length". I am not sure I understand the rationale behind this. Once the raw reads e.g. 100 or 150nt length are trimmed for quality, adapters etc., read length becomes variable.
                    In this case, I would advice you to look at the minimum length of your reads and trim all of them to have this minimum length, or depending of your FastQC report, I would go without adapter trimming... Really it is worthwhile trying Edena. As example for a total de novo assembly of Staphylococcus aureus, I got 12 contigs (which stopped because of the rRNA operons).

                    s.
                    Last edited by SylvainL; 11-05-2014, 11:41 PM.

                    Comment


                    • #11
                      Thxs for your help. Finally I run edena, but I can't improve rna-seq mapping (~50%) comparing with clc, a5 and spades (~90%). I think that the next step is choice assembly with relative good statistics and at least ~90% rna-seq map.

                      Regards.

                      Comment


                      • #12
                        Originally posted by freestile View Post
                        Thxs for your help. Finally I run edena, but I can't improve rna-seq mapping (~50%) comparing with clc, a5 and spades (~90%). I think that the next step is choice assembly with relative good statistics and at least ~90% rna-seq map.

                        Regards.
                        When you say your rna-seq mapping is low do you mean only 50% of your rna-seq reads are mapped to your assembly? If yes, did you try a blast on some unmapped reads to see if they really come from your bacterium (or closed to)?

                        Another question: after rna-seq mapping to your asembly, do you have some regions which are not covered at all?

                        These RNAseq were performed with rRNA depletion? Do you have rRNA operons on your assembly?

                        Comment


                        • #13
                          Originally posted by SylvainL View Post
                          When you say your rna-seq mapping is low do you mean only 50% of your rna-seq reads are mapped to your assembly? If yes, did you try a blast on some unmapped reads to see if they really come from your bacterium (or closed to)?

                          Another question: after rna-seq mapping to your asembly, do you have some regions which are not covered at all?

                          These RNAseq were performed with rRNA depletion? Do you have rRNA operons on your assembly?
                          Yes, only 50% and I blast some unmapped reads and correspond to bacteria genes.

                          I have to check te second question

                          And yes I make rRNA depletion in library preparation.

                          Regards !

                          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
                          7 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, Yesterday, 06:07 PM
                          0 responses
                          7 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 03-22-2024, 10:03 AM
                          0 responses
                          49 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 03-21-2024, 07:32 AM
                          0 responses
                          66 views
                          0 likes
                          Last Post seqadmin  
                          Working...
                          X