MetaPhlAn 4.2.2 work with humann v3.0.1

Microbial community profiling HUMAnN
1

I used MetaPhlAn version 4.2.2 to run the alignment successfully and generated an abundance file. Now, I hope to use this file as input and run humann v3.0.1 to obtain the correspondence between species, abundance, and functional profiles. However, I have been unable to achieve this. My codes were as follows:

/datapool/software/anaconda3/envs/phylophlan/bin/humann3 -i B38_03.merge.fq.gz -o B38_03.hum --threads 20 --metaphlan /datapool/software/anaconda3/bin/metaphlan --metaphlan-options “–index mpa_vJan25_CHOCOPhlAnSGB_202503 -t rel_ab --mapout B38_03.merge_metaphlan_bowtie2.txt --db_dir /datapool/software/anaconda3/lib/python3.7/site-packages/metaphlan/metaphlan_databases” --diamond /datapool/software/anaconda3/envs/phylophlan/bin/diamond

The above didn’t work. So I try to run the MetaPhlAn firstly, and put the output to humann3 as input.
/datapool/software/anaconda3/envs/phylophlan/bin/humann3 -i B38_03.merge.fq.gz -o B38_03.hum --threads 20 --taxonomic-profile B38_03.mp --metaphlan /datapool/software/anaconda3/bin/metaphlan --metaphlan-options “–index mpa_vJan25_CHOCOPhlAnSGB_202503 --db_dir /datapool/software/anaconda3/lib/python3.7/site-packages/metaphlan/metaphlan_databases”

But the codes also failed. The errors as follows:
ERROR: The MetaPhlAn2 taxonomic profile provided was not generated with the database version v30 . Please update your version of MetaPhlAn2 to v3.0.

Thank you.

2 Likes
2

I met same problem to be solved, my problem is that humann database could not use the metaphlan_bug_list.tsv file as input.

3

Sorry that that error message is not very informative. HUMAnN 4.0.0.alpha.1 is currently configured to work with the Oct22 SGB set. By default MetaPhlAn 4.2 is now profiling the Jan25 SGB set, so HUMAnN is not able to work with the corresponding profile. The post-alpha release of HUMAnN 4.0 will support the Jan25 SGBs.