Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Minimum amount of data needed for reliable results?

    We are trying to do analysis for whole genome metagenomic data taken from the surface of the eye. Each sample has millions of reads generated, but of those reads at most only 1-2% are bacterial reads. We are wondering if there is some information/resources about the amount of data available related to the reliability of the results eg. finding out the taxonomic information for bacteria down to species level that are greater than 1% relative abundance. Currently we are aligning the data to whole genome bacterial sequences, but there are many multi-mapping locations which many of which may be false positives. We have tried using Metaphlan2 to do alignment which uses a custom catalog of unique markers for different clades, but usually only several hundred reads will be mapped back -- many of the samples report very low/no species present. Specifically, we are wondering methods to do analysis for whole genome metagenomic sequences where the amount of data is very low. Any help is greatly appreciated.

    Daniel

  • #2
    You might try removing human sequence, then assembling the rest and BLASTing the contigs against nt/nr/RefSeq microbial. Assuming the contigs are longer than read length, they will give you more reliable hits. What kind of depth do you have for the bacteria? You can find that out with a kmer-frequency histogram, after human reads are removed.

    Comment


    • #3
      Thank you for the quick response. The idea of assembling the reads into contigs before alignment makes sense, I will let me supervisor know. Do you know of good software to do this? I have tried Velvet in the past but did not use it extensively.

      I forgot to mention that we have removed human sequences, although the revised reference genome that you created seems like it would be especially useful for us.

      Could you give more information on how to estimate the depth of the bacteria? There is generally less than 100,000 bacterial reads per sample out of 20-30 million initial reads (before any trimming/contaminant removal).

      Comment


      • #4
        I suggest Spades or Megahit for metagenome assembly. 100k is not many reads; you might not have sufficient depth for assembly. But in that case, you may get a better assembly by combining all bacterial reads from all samples and assembling together. Then you can quantify by mapping to the combined assembly.

        For human removal, the raw human genome is fine in your case (bacteria). The masked version is mainly to allow decontamination of eukaryotes, which have shared sequence with human; bacteria basically don't.

        Comment


        • #5
          You can create rarefaction curves to see if what you have is likely sufficient to describe the metagenomic profile.

          The basic process is to remove reads and see if your calculation of the species diversity is similar. A low complexity sample will plateau at a low coverage, while the diversity of a high complexity sample will just keep increasing substantially with more reads.

          Comment


          • #6
            I had some luck using IMSA in a similar situation to the one you describe, but only after human read removal

            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