No_name after humann_rename_table

Microbial community profiling HUMAnN
1

Hi,
I am sorry if I missed a solution to my problem, but I cannot find it in the previous forum.

I got gene families.tsv file from running Humann (v3.7) and ‘uniref90_201901b_full.dmnd’ was used as reference database. Then, I tried to rename the UniRef90 ID to EC or MetaCyc ID with the following code line. humann_rename_table --input merged_knead_2102420018_genefamilies.tsv --output ./rename_test.tsv --names metacyc-pwy, and I have got a message: Renamed 0 of 178911 entries (0.00%).

I even confirmed that UniRef90_A0A1F3ZUU1, which is the L-lacate dehydrogenase (EC 1.1.1.27), in my gene families.tsv file, whereas this was also renamed as No_name.

Is this phenomenon due to the extension of UniRef90 ID such as UniRef90_A0A1F3ZUU1 (extension in italics)? or I missed something?

I appreciate your help in advance.
Gihyeon.

2

The rename script is used to attach human-readable names (e.g. “Cyclin-dependent Kinase 1”) to systematic IDs (e.g. “UniRef90_ABC”). To group genes into broader functional categories you would use the regroup script instead. Note that you cannot directly regroup genes to pathways, but you can regroup them to ECs or MetaCyc reactions (or KOs, Pfams, GO terms, etc.).