A common goal in many microbiome studies is to identify features (i.e., species, OTUs, gene families, etc.) that differ according to some study condition of interest. While often done, this is a difficult task, and in the Introduction to the Statistical Analysis of Microbiome Data in R post I touch on some of the reasons for this.
Below I provide scripts to implement the current default workflow for taxonomic and functional profiling using the Huttenhower Lab’s Biobakery Tool Suite used by the Microbial Metagenomics Analysis Center (MMAC) at CCHMC for paired-end data.
Below I provide scripts to implement several workflows for denoising 16s rRNA gene sequences used by the Microbial Metagenomics Analysis Center (MMAC) at CCHMC for paired-end data. These scripts are written to run on the CCHMC high-performance computing (HPC) cluster.
Below I provide scripts to implement the current default workflow for taxonomic profiling using Kraken2 and Bracken and functional profiling using HUMAnN2 used by the Microbial Metagenomics Analysis Center (MMAC) at CCHMC for paired-end data.
I recently read through Calgaro et. al. “Assessment of statistical methods from single cell, bulk RNA-seq, and metagenomics applied to microbiome data” where they examined the performance of statistical models developed for bulk RNA (RNA-seq), single-cell RNA-seq (scRNA-seq), and microbial metagenomics to:
During the Introduction to Metagenomics Summer Workshop we discussed denoising amplicon sequence variants and worked through Ben Callahan’s DADA2 tutorial. During that session, I mentioned several other approaches and algorithms for denoising or clustering amplicon sequence data including UNOISE3, DeBlur and Mothur.
A collaborator recently asked if I could help pull down a few thousand sequence files from the NCBI Sequence Read Archive (SRA) for a secondary analysis. This is a short post primarily to help me (and hopefully others) remember how to do this once you have a set of SRR IDs of interest.
This post is from a tutorial demonstrating the processing of amplicon short read data in R taught as part of the Introduction to Metagenomics Summer Workshop. It provides a quick introduction some of the functionality provided by phyloseq and follows some of Paul McMurdie’s excellent tutorials.
This post is also from the Introduction to Metagenomics Summer Workshop and provides a quick introduction to some common analytic methods used to analyze microbiome data. I thought it might be of interest to a broader audience so decided to post it here.