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.