HUMAnN 3.8 - How to get taxonomic results at the order level instead of genus?

Hello,

I’ve been using HUMAnN 3.8.
I got many unclassified results because I think the taxonomic resolution of genus is too precise.
So, I’m interested in obtaining taxonomic results at the order level rather than at the genus level.
Unfortunately it seems that is not possible? Can you please help me ?

Best regards,
Jérémy Tournayre

Dear Jeremy,

I’m not sure whether you can specify the taxonomic level in the HUMAnN directly. Nevertheless, you can ran MetaPhlAn separately and specify the level there:

--tax_lev 'o' (orders only)

Then you can ran HUMAnN and refer to the MetaPhlAn output:

--taxonomic-profile metaphlan_output

I hope it solves your problem!

1 Like

Hello,

I found solution here:

with humann_infer_taxonomy

humann_infer_taxonomy -l Genus -d uniref90-tol-lca -i humann_genefamilies.tsv -o Genus/humann_genefamilies.tsv"

It works but I don’t know why there are some stange things like this below.
Why do I have “unclassified” instead of “g__Prevotella” and “g__Prevotellaceae_unclassified” in the Genus file?

Original “humann_genefamilies.tsv” :

UniRef90_A0A010YBU9 43.4636486692
UniRef90_A0A010YBU9|g__Prevotella.s__Prevotella_sp_ne3005 15.4012377598
UniRef90_A0A010YBU9|g__Prevotella.s__Prevotella_ruminicola 12.1212121212
UniRef90_A0A010YBU9|g__Prevotellaceae_unclassified.s__Prevotellaceae_bacterium_MN60 6.0506050605
UniRef90_A0A010YBU9|g__Prevotella.s__Prevotella_sp_tc2_28 4.4004400440
UniRef90_A0A010YBU9|g__Prevotella.s__Prevotella_sp_tf2_5 3.2789872934
UniRef90_A0A010YBU9|g__Prevotella.s__Prevotella_aff_ruminicola_Tc2_24 2.2111663903

In Phylum:

UniRef90_A0A010YBU9 43.4636
UniRef90_A0A010YBU9|p__Bacteroidetes 43.4636

In Genus:

UniRef90_A0A010YBU9 43.4636
UniRef90_A0A010YBU9 unclassified  43.4636

Best regards,
Jérémy Tournayre