Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Metagenome analysis - first steps

    Hi,
    I have a few general questions how to analyze metagenomes.

    Platform: Illuminia Highseq, 2x 150 bp PE reads with ~80 bp overlap
    Aim: Get taxonomic representation and functional profile (e.g. KEGG)
    Preferred tool: MEGAN 5

    My questions to get started:

    1. Should I first merge the paired reads by their overlap or should both be analyzed seperately? In MEGAN I could check a box for PE reads, so maybe not merging them serves a purpose. If I should merge, which tool do you recommend?

    2. Should I collapse my sequences before blasting them. I could do that to group identical reads to OTUs. They get a new header than, stating how often that OTU was present. That would greatly reduce the blast time. Question is, can downstream tools like MEGAN deal with that?

    3. Any recommended BLAST settings to get a good balance between accuracy and computational time?


    Sorry for these stupid questions. I am not a bioinformatician. I am normally a try and error learner. But if every computational step takes several days, I better go the right way from the start...

    Thanks,
    Sören

  • #2
    1. Combine the reads together, especially with that much overlap. Use FLASH, though I think Pandaseq does the same.

    2. Depends on how much of your reads are unique. If a lot of your read content is not unique, it may not give you much.

    I'm not sure what you mean by OTU's. OTU is Operational Taxonomic Unit and that comes up in the context of 16S data. fastx_collapser assigns a unique number to a unique read, but it's not an OTU.

    3. An option you may want to think about is MG-RAST, or if you feel lucky, you can attempt a metagenomic de novo assembly, extract the open reading frames with prodigal and blast those.

    Comment


    • #3
      Originally posted by sfranzenburg View Post
      1. Should I first merge the paired reads by their overlap or should both be analyzed seperately? In MEGAN I could check a box for PE reads, so maybe not merging them serves a purpose. If I should merge, which tool do you recommend?
      I won't comment on whether it is better to merge or not in your case, but we typically assemble the metagenome, then map to the assembly. The assembly gets annotated, and only reads that did not map to the assembly get merged, then annotated.

      If you do merge, I suggest BBMerge; it's faster and more accurate than FLASH - orders of magnitude fewer false positives with similar merging rates, in my tests (which were done on 2x150bp libraries with average 30bp overlap).

      Usage:

      bbmerge.sh in1=read1.fq in2=read2.fq out=merged.fq outu=unmerged.fq hist=histogram.txt

      Comment


      • #4
        Thanks to both of you, very helpful insight.

        @Brian: Do you do de novo assembly with Trinity? How does it handle closely related bacterial species. My samples are gut tissues and I expect a variety of Acetobacter species in there.

        Comment


        • #5
          I have no experience with Trinity. We generally use Soap for metagenome assemblies. As for handling closely-related strains... well, the assembly will be worse in that case Still, you should end up with a lot of contigs much longer than reads. That means you have much less data to process, and the annotation/blasting/whatever will be more informative.

          I've found normalization of reads helpful prior to metagenome assembly with Soap, especially when you have a lot of data - it yields a smaller assembly with less redundancy, better continuity, and more of the source reads map back to it. The link I gave also has a tool called BBNorm (bbnorm.sh) that can be used for normalization. But, that will tend to make it harder to resolve closely-related strains, so it depends on what you're trying to study.

          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
          8 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, Yesterday, 06:07 PM
          0 responses
          8 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
          67 views
          0 likes
          Last Post seqadmin  
          Working...
          X