Feb 08, 2016
Metabolic profiling with HUMAnN
- Curtis Huttenhower
- VERVE Net
- Huttenhower Lab
Protocol Citation: Curtis Huttenhower 2016. Metabolic profiling with HUMAnN. protocols.io https://dx.doi.org/10.17504/protocols.io.d9799m
License: This is an open access protocol distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Protocol status: Working
Created: December 07, 2015
Last Modified: November 20, 2017
Protocol Integer ID: 2079
Keywords: metabolic profiling of metagenomic dataset, metabolic profiling with humann, metagenomic dataset, metabolic profiling, metabolic profile, taxonomic biomarker discovery with lefse, taxonomic biomarker discovery, multiple samples taxonomic biomarker discovery, using humann, graphlan visualization, metaphlan, humann, graphlan, metaphlan output merge
Abstract
In this protocol, we perform a metabolic profiling of metagenomic datasets using HUMAnN applied to the 20 samples taxonomically profiled in the other protocols contained in this collection.
Once the metabolic profile has been generated, the following protocols can be performed:
- MetaPhlAn output merge and visualizations
- GraPhlAn visualization of single and multiple samples
- Taxonomic biomarker discovery with LEfSe
Guidelines
Metabolic profiles of the 20 HMP samples with HUMAnN can be generated with the following steps:
- Place the translated BLAST results using KEGG gene identifiers in the
input
directory (optionally can be gzipped or bzipped, several other formats can be enabled editing the settings in the
SConstruct
file).
- Run the
scons
command, optionally parallelizing multiple analyses using the
-j
flag. Results will be placed in the "output" directory.
We provide the HUMAnN output for users who want to perform the downstream analysis pipeline but avoid the computational intensive steps above. The 20 samples profiled with HUMAnN are available here.
Using HUMAnN output, you can perform the metabolic counterpart of the taxonomic pipeline presented in this tutorial. With a table of metabolomic abundances, the previous protocols can be performed with little to no modification. We report below the first command to obtain the merged table of metabolic abundances.
$ mkdir output $ utils/merge_metaphlan_tables.py humann_profiling/*.txt > output/merged_humann_abundance_table.txt
Troubleshooting
Place the translated BLAST results using KEGG gene identifiers in the
input
directory.
Note
Optionally can be gzipped or bzipped, several other formats can be enabled editing the settings in the
SConstruct
Run the
scons
command, optionally parallelizing multiple analyses using the
-j
flag. Results will be placed in the "output" directory.
We provide the HUMAnN output for users who want to perform the downstream analysis pipeline but avoid the computational intensive steps above. The 20 samples profiled with HUMAnN are available here.
Using HUMAnN output, you can perform the metabolic counterpart of the taxonomic pipeline presented in this tutorial. We report below the first command to obtain the merged table of metabolic abundances.
Command