Metagenomic data analysis discussion NEON Workshop
Download as PPTX, PDF2 likes1,874 views
This document discusses issues and considerations around analyzing metagenomic sequencing data. It notes that 50% or more of sequences from soil metagenomes are unknown and cannot be annotated. While assembly and binning methods exist, they have limitations like dependence on sequencing depth and need for unique bins. The document considers what types of signals or targets could be detected in metagenomic data and discusses tools like MG-RAST for analysis but notes they have limitations and oversimplifications. It suggests future tool development will be challenging and meeting all demands impossible, and that data management may be a bigger challenge than the analysis itself.
Metagenomic data analysis discussion NEON Workshop
- 1. Discussion on Metagenomic Data for NEON Workshop Adina Howe July 15, 2014 Boulder, CO
- 2. What I do with my sequencing data • Compare it to known references (BLAST, alignment mapping, HMMs) – Database gaps, 50%+ in soil metagenomes unknown – Annotation errors • Assembly (site specific reference) – Dependent on sequencing depth – Largest soil datasets, 1 billion reads, 10% assembled • Co-occurrence / binning – Abundance or nucleotide based – Sequencing coverage dependent, need “unique” bins – Particularly problematic for short reads
- 3. Sentinels in metagenomic data? • Presence/Absence? Quantification? • Phyla? Strain? Functional level? • Strong signals of similarity or differences in both “annotatable” and unknown sequences • Information about targets for which we need more references (by other methods) • Indexes (alpha diversity, copiotroph/oligotroph, % polymorphism, recA/16S gene count)
- 4. What I do with my sequencing data • Compare it to known references (BLAST, alignment mapping, HMMs) – Database gaps, 50%+ in soil metagenomes unknown – Annotation errors • Assembly (site specific reference) – Well developed pipelines and tutorials – Almost too many choices • Co-occurrence / binning – Some tools available, still in development largely – Eventually need to link to knowns for validation (or experimental validation)
- 5. What has been done – MG-RAST • Free • Easy • Growing in usage by many • Not used by the “leading wedge” (at least as more than a blast engine) • Free data / annotation storage and maintenance • NEON – do you want to build your own? Does MG-RAST suffice?
- 6. MG-RAST
- 10. Fine tuning is not possible • What happens when your sequence matches two known genes equally? • What happens if there are multiple taxa related to one function? • What if you want more than ordination (e.g., primer evaluation?) • What if you want to change your distance matrix for ordination? • How do (not) you normalize / rarify your data? • What if you want more sensitivity than a BLAT analysis (e.g., amplicons) • Computationally-optimized, not biologically – Denitrification, Nitrate reductase, nitrate ammonification all distinct categories in Level 3 – If you want to merge, not trivial – especially in dealing with abundance estimations
- 11. Looking forward • Development of all-pleasing computational, statistical, visualization tool will be impossible • NEON community demands will more than likely be easy to use and aimed at reducing (perhaps overly simplified) complex data into comparable metrics (TBD) • NEON data is actually not a big dataset and current tools should (imho) suffice • Data management and associated “metadata” may be the larger challenge (in the face of rapidly changing datasets)