8000 GitHub - snayfach/MIDAS at v1.2.0
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

snayfach/MIDAS

BranchesTags

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

234 Commits
bin
bin
 
 
docs
docs
 
 
images
images
 
 
midas
midas
 
 
scripts
scripts
 
 
test
test
 
 
.gitignore
.gitignore
 
 
LICENSE
LICENSE
 
 
README.md
README.md
 
 
setup.py
setup.py
 
 

Repository files navigation

Metagenomic Intra-Species Diversity Analysis System (MIDAS)

MIDAS is an integrated pipeline that leverages >30,000 reference genomes to estimate bacterial species abundance and strain-level genomic variation, including gene content and SNPs, from shotgun metagnomes.

Applications

  1. Profile bacterial species abundance: rapidly estimate the abundance of 5,952 bacterial species
  2. Strain-level pan-genome profiling: estimate the gene content of populations based on mapping to genes from reference genomes
  3. Single nucleotide polymorphism prediction: identify single-nucleotide polymorphisms (SNPs) of populations based on mapping to reference genomes
  4. Phylogenetic inference: reconstruct the phylogeny of strains from metagenomes and reference genomes
  5. Population genetic inference: quantify strain-level diversity, differentiation, and selection within and between metagenomes

Table of Contents

  1. [Step-by-step tutorial] (docs/tutorial.md)
  2. [Install or update MIDAS] (docs/install.md)
  3. [Download default MIDAS database] (docs/ref_db.md)
  4. [Build your own custom database] (docs/build_db.md)
  5. Scripts to run MIDAS on a single sample:
  • [Estimate species abundance] (docs/species.md)
  • [Predict pan-genome gene content] (docs/cnvs.md)
  • [Call single nucleotide polymorphisms] (docs/snvs.md)
  1. Scripts to merge MIDAS results across samples:
  • [Merge species abundance] (docs/merge_species.md)
  • [Merge gene content] (docs/merge_cnvs.md)
  • [Merge SNPs] (docs/merge_snvs.md)
  1. Example scripts for analyzing gene content and SNPs:
  • [Strain tracking] (docs/strain_tracking.md)
  • [Population diversity] (docs/snp_diversity.md)
  • [Core-genome phylogenetic trees] (docs/snp_trees.md)
  • [Gene content dynamics] (docs/compare_genes.md)

Citation

If you use this tool, please cite: S Nayfach, B Rodriguez-Mueller, N Garud, and KS Pollard. "An integrated metagenomics pipeline for strain profiling reveals novel patterns of transmission and global biogeography of bacteria". Genome Research 2016. doi: 10.1101/gr.201863.115

Pipeline

**An integrated pipeline to estimate bacterial species abundance and strain-level genomic variation from shotgun metagnomes** **A) Metagenome species profiling.** Reads from a metagenomic sample are aligned against a database of phylogenetic marker genes and are assigned to species groups. Mapped reads are used to estimate the genome-coverage and relative abundance of 5,952 genome-clusters. **B) Metagenome pan-genome profiling.** A pan-genome database is dynamically constructed based on the subset of species that are present at high coverage (e.g. >1x) in the metagenome. Reads are mapped to the gene database using Bowtie2. Mapped reads are used to infer gene copy number and gene presence/absence. **C) Single-nucleotide variant prediction.** A representative genome database is constructed, as described in (B). Reads are globally aligned to the genome database using Bowtie2. Mapped reads are used to identify variants, predict consensus alleles, and estimate allele frequencies. **D) Merge results.** For each species, results are merged across one or more samples to generate several outputs, including: a gene presence/absence matrix, an allele frequency matrix, an approximate maximum-likelihood phylogenetic tree.

Examples

**Comparative genomics of *Bacteroides ovatus* strains across host microbiomes** **A)** Presence or absence of genes in the *Bacteroides ovatus* pangenome across human faecal metagenomes. Column colors indicate whether a gene is core (blue; occurs in >95% of samples), auxiliary (red; occurs in 1-95% of samples ), or absent (green; occurs in < 1% of samples). **B)** Gene set enrichment analysis identifies functions overrepresented in the core genome, auxiliary genome, and genes that only occur in reference genomes.

About

An integrated pipeline for estimating strain-level genomic variation from metagenomic data

dx.doi.org/10.1101/gr.201863.115

Resources

Readme

License

GPL-3.0 license
Activity

Stars

127 stars

Watchers

12 watching

Forks

52 forks
Report repository

Releases 8

Version 1.3.2 Latest
Nov 26, 2017
+ 7 releases

Packages

No packages published

Contributors 3

Languages
0