Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • reps and analysis

    Hi all!

    I am doing a 16S amplicon analysis on soil samples with 5 different treatments.
    I have 9 replicates per sample.
    I got from illumina base space a file containing all the reads for all different taxonomic levels. Right now I am working on the class one.
    I normalized the reads to the total number of reads. now I would like to verify the variability within my treatment between the reps, and obtain a single data set to work with per traetment.

    Could you suggest me how is the best way to proceed with the analysis testing the variability, for istance, and then combine the sequences coming from the replicates into a unique sample sequence data set?
    Could you suggest me the best method to analyze the variability between my replicates?

    Do you have ideas?

    Thank you very much!

  • #2
    Generally you don't want to discard replicate information. All/most of the downstream analyses that you want to do can and should be done with all of your replicates intact.

    See the QIIME tutorial here:
    http://nbviewer.ipython.org/github/b...tutorial.ipynb

    A PCoA plot stratified by your five different treatments might be a quick and easy way to convince yourself of the reproducibility of your replicates.

    Comment


    • #3
      Hi,
      Thank you fanli for replying.
      I have a couple of doubts still about reps treatment for 16S sequnecing.
      I didn't mean to discard the replicates data. I would like just to ask how to deal with them.
      I see the PCoA plot as a tool to evaluate the reproducibility of the reps. What will you do then?
      would you combine the replicates reads ? would you sum the replicates reads? how will you proceed in order to obtain treatments sets?

      Thank you a lot!

      Cheers

      Comment


      • #4
        I suppose you could combine the replicate reads if you wanted a single list of OTUs for each treatment. What exactly are you trying to compare with your data? If you want to compare your treatments, it makes more sense to compute diversity/population statistics for each replicate and then compare those between treatments (e.g. a 9-vs-9 comparison).

        Comment


        • #5
          HI FAnli again,
          I would like to compare the treatments, taking into account the intrinsic variability within the samples (due to DNA extraction or due to field).
          I thought to have a single data set per treatment. maybe it is not the right way.
          Which tool will you use for the 9vs9 comparison?
          thnks

          Comment


          • #6
            Hi Vale16S,

            I've been using QIIME to compare alpha- and beta-diversity statistics between treatment groups. See http://qiime.org/scripts/compare_alpha_diversity.html.

            From the man page:
            For example, if your mapping file had a category called ‘Treatment’ that separated your samples into three groups (Treatment=’Control’, Treatment=’Drug’, Treatment=‘2xDose’), passing ‘Treatment’ to this script would cause it to compare (Control,Drug), (Control,2xDose), (2xDose, Drug) alpha diversity values. By default the two-sample t-test will be nonparametric (i.e. using Monte Carlo permutations to calculate the p-value), though the user has the option to make the test a parametric t-test
            So you'd be doing a 9vs9 t-test, essentially.

            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
            66 views
            0 likes
            Last Post seqadmin  
            Working...
            X