Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Originally posted by F_KVH View Post
    Dear Simon,
    Thanks for your help. I installed again & now it works. Do you think I need to install xcode app?
    Best regards
    No, you'd only need something like xcode if you wanted to edit the source code. It wouldn't be needed just to run it.

    Comment


    • Dealing with Random Hexamer Bias?

      Originally posted by lletourn View Post
      The illumina RNA protocol uses random hexamers to amplify the RNA. The thing is they are not 100% random so the beginning looks skewed for base composition, but that's because of the amplification.

      For mapping it's no problem. For assembly it might confuse some assemblers. (When assembling I would trim the 5' of RNA, not for mapping)
      Our group recently used the HiSeq 2000 platorm to generate transcriptome data (single-end, 50 bp reads). I have noticed that Illumina transcriptome sequencing yields typically yield these errors ("Per base sequence content" and "Per base GC content") during FastQC analysis. You suggest that these can safely be ignored when mapping to a genome? I wasn't sure if there was a "best practices" approach to dealing with these biases. Is there a mapping algorithm that is preferred among those who are dealing with Illumina transcriptome sequencing data?

      (I am sorry if this question has been answered thoroughly elsewhere in the forum... I have only just joined, and despite trying to navigate the posts with the "Search" tool, I have not yet come across an answer).

      Comment


      • Originally posted by boilermaker View Post
        Our group recently used the HiSeq 2000 platorm to generate transcriptome data (single-end, 50 bp reads). I have noticed that Illumina transcriptome sequencing yields typically yield these errors ("Per base sequence content" and "Per base GC content") during FastQC analysis. You suggest that these can safely be ignored when mapping to a genome? I wasn't sure if there was a "best practices" approach to dealing with these biases.
        It would probably be better to ignore them than try to correct them, though if you posted the fastqc graphs it would be easier to say.
        Is there a mapping algorithm that is preferred among those who are dealing with Illumina transcriptome sequencing data?
        I prefer BBMap when mapping Illumina RNA-seq data. It's more robust to errors than other RNA-seq aligners, and doesn't require an annotation file. Oh, and I wrote it, but that's not why.

        Comment


        • Originally posted by Brian Bushnell View Post
          It would probably be better to ignore them than try to correct them, though if you posted the fastqc graphs it would be easier to say.
          Thank you Brian. I will certainly give BBMap a try (and thank you very much for scripting it!)

          I have attached a "typical" fastqc graphs (per base gc content, per base sequence content) from one of my datasets (most have profiles like this example).
          Attached Files

          Comment


          • Doesn't look ideal, but I can't think of a good way to improve it, assuming you've already trimmed adapters (which can alter the base composition).

            Comment


            • The base composition bias you are seeing is very typical for RNA-Seq, and has indeed been mentioned on numerous posts here. The most relevant publication about this can be found here: http://nar.oxfordjournals.org/content/38/12/e131.full.

              In essence, the bias is normally introduced by the random priming step in the RNA-Seq library preparation which is not quite as random as you would hope it was. Trimming the first positions of every read wouldn't make any difference since the sequence would still align to the very same position.

              Comment


              • Originally posted by fkrueger View Post
                The base composition bias you are seeing is very typical for RNA-Seq, and has indeed been mentioned on numerous posts here. The most relevant publication about this can be found here: http://nar.oxfordjournals.org/content/38/12/e131.full.

                In essence, the bias is normally introduced by the random priming step in the RNA-Seq library preparation which is not quite as random as you would hope it was. Trimming the first positions of every read wouldn't make any difference since the sequence would still align to the very same position.
                If I had a penny for every time this question has been asked (and answered)

                You (or Simon) really should make a sticky post with a few example plots about this (and k-mers). Will save many, a bunch of time and worry. Another suggestion would be to put a note on the FastQC page itself, so it would be visible to anyone downloading the software.
                Last edited by GenoMax; 04-30-2014, 04:10 PM.

                Comment


                • Originally posted by GenoMax View Post
                  If I had a penny for every time this question has been asked (and answered)

                  You (or Simon) really should make a sticky post with a few example plots about this (and k-mers). Will save many, a bunch of time and worry. Another suggestion would be to put a note on the FastQC page itself, so it would be visible to anyone downloading the software.
                  I don't think we, as normal users, can create sticky posts can we?

                  I've added some information about this topic to the FastQC help so at least there will be a bit more guidance in the next release. I'll also have to make a new video for that so I'll make sure to mention it there. This is something we talk about at some length in the RNA-Seq analysis courses I run since it is a true technical bias but just one we tend to ignore (mostly because of not having any other option).

                  It's maybe also worth noting that a similar bias now seems to be appearing in transposase fragmented libraries, so it's not just RNA-Seq libraries which see this.

                  Comment


                  • Originally posted by simonandrews View Post
                    I don't think we, as normal users, can create sticky posts can we?
                    I was thinking that you can create a post and then PM ECO (or one of the other moderators) to see if they can make it sticky.

                    Originally posted by simonandrews View Post
                    I've added some information about this topic to the FastQC help so at least there will be a bit more guidance in the next release. I'll also have to make a new video for that so I'll make sure to mention it there.
                    Would it be useful to make this information available in the "Documentation" section of the main FastQC page where you have the Good/bad data set examples? Something along the lines of "this type of nucleotide distribution is normal if you have RNAseq data". Creating a video is a good idea but having the information available on a page that people can glance through (while their FastQC download is ongoing?) can make a better impact.

                    Comment


                    • Originally posted by GenoMax View Post
                      Would it be useful to make this information available in the "Documentation" section of the main FastQC page where you have the Good/bad data set examples?
                      That's exactly what I did. When the package updates to the next version it will be there.

                      Comment


                      • Hi All,

                        has by chance anybody figured out a pipeline that would automatically produce/convert the FASTQC report as a single file -perhaps as a single PDF file?
                        In my eyes it would be great if one could run FASTQC from the command line and not have to deal with an HTML folder and multiple files.
                        Thanks in advance.

                        Comment


                        • Originally posted by luc View Post
                          Hi All,

                          has by chance anybody figured out a pipeline that would automatically produce/convert the FASTQC report as a single file -perhaps as a single PDF file?
                          In my eyes it would be great if one could run FASTQC from the command line and not have to deal with an HTML folder and multiple files.
                          Thanks in advance.
                          The next release of the program (which will be out by the end of next week if it kills me!) creates a single HTML file with embedded graphics so you will then be able to distribute just that file instead of having to keep the existing folder structure.

                          Comment


                          • Hi Simon,

                            thanks a lot. That sounds very practical.

                            Comment


                            • After a very long gestation I've finally released a new version of FastQC (v0.11.1). This is now available from the project web site.

                              This is a major release of the software which introduces a load of new features. Some of the big changes are:
                              • Added configurable warn/fail thresholds for all modules
                              • Allow modules to be selectively turned off
                              • Added a per-tile quality plot for Illumina libraries
                              • Added an adapter content plot
                              • Improved the duplication plot
                              • Improved the Kmer module
                              • Used embedded graphics in the HTML output so you can distribute a single file
                              • Added the ability to read data from stdin
                              • Changed how base grouping works to better accommodate long reads
                              • Dropped support for Solexa64 format (NB not Phred 64 which is still supported) to avoid mis-detection errors


                              We've done a fair bit of testing on the new version but I'm aware that there's a lot of new code in there so please report any problems either directly into our bug tracking system or via email to [email protected]

                              Comment


                              • 2-6-13: Version 0.11.1 released
                                Shouldn't that be

                                2-6-14: Version 0.11.1 released
                                Thank you for adding the example reports on the FastQC page!

                                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
                                11 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