Error in picrust2_pipeline.py command, possibly related to EPA-NG.

[Kirin_C]

Hello all!

I encountered the following problem when using the picrust2_pipeline.py command of PICRUSt2. Can anyone help me take a look：

picrust2_pipeline.py \
  -s dna-sequences.fasta \
  -i feature-table.biom \
  -o picrust2_out \
  -p 4

Error running this command:
place_seqs.py --study_fasta dna-sequences.fasta --ref_dir /home/ubuntu_user/miniconda3/envs/fresh-picrust2/lib/python3.8/site-packages/picrust2/default_files/prokaryotic/pro_ref --out_tree picrust2_out/out.tre --processes 1 --intermediate picrust2_out/intermediate/place_seqs --min_align 0.8 --chunk_size 5000 --placement_tool epa-ng

Standard error of the above failed command:

Error running this command:
epa-ng --tree /home/ubuntu_user/miniconda3/envs/fresh-picrust2/lib/python3.8/site-packages/picrust2/default_files/prokaryotic/pro_ref/pro_ref.tre --ref-msa picrust2_out/intermediate/place_seqs/ref_seqs_hmmalign.fasta --query picrust2_out/intermediate/place_seqs/study_seqs_hmmalign.fasta --chunk-size 5000 -T 1 -m /home/ubuntu_user/miniconda3/envs/fresh-picrust2/lib/python3.8/site-packages/picrust2/default_files/prokaryotic/pro_ref/pro_ref.model -w picrust2_out/intermediate/place_seqs/epa_out --filter-acc-lwr 0.99 --filter-max 100

Standard output of the above failed command:
INFO Selected: Output dir: picrust2_out/intermediate/place_seqs/epa_out/
INFO Selected: Query file: picrust2_out/intermediate/place_seqs/study_seqs_hmmalign.fasta
INFO Selected: Tree file: /home/ubuntu_user/miniconda3/envs/fresh-picrust2/lib/python3.8/site-packages/picrust2/default_files/prokaryotic/pro_ref/pro_ref.tre
INFO Selected: Reference MSA: picrust2_out/intermediate/place_seqs/ref_seqs_hmmalign.fasta
INFO Selected: Filtering by accumulated threshold: 0.99
INFO Selected: Maximum number of placements per query: 100
INFO Selected: Automatic switching of use of per rate scalers
INFO Selected: Preserving the root of the input tree
INFO Selected: Specified model file: /home/ubuntu_user/miniconda3/envs/fresh-picrust2/lib/python3.8/site-packages/picrust2/default_files/prokaryotic/pro_ref/pro_ref.model
INFO Selected: Reading queries in chunks of: 5000
INFO Selected: Using threads: 1
INFO     ______ ____   ___           _   __ ______
        / ____// __ \ /   |         / | / // ____/
       / __/  / /_/ // /| | ______ /  |/ // / __
      / /___ / ____// ___ |/_____// /|  // /_/ /
     /_____//_/    /_/  |_|      /_/ |_/ \____/ (v0.3.8)
INFO Using model parameters:
INFO    Rate heterogeneity: GAMMA (4 cats, mean),  alpha: 0.453141 (user),  weights&rates: (0.25,0.0250674) (0.25,0.220229) (0.25,0.782933) (0.25,2.97177)
        Base frequencies (user): 0.229585 0.22008 0.298596 0.251739
        Substitution rates (user): 1.00319 2.79077 1.5301 0.87441 3.83966 1
INFO Output file: picrust2_out/intermediate/place_seqs/epa_out/epa_result.jplace

The picrust2_out folder outputted also only contains 'intermediate'.

[Robyn Wright]

Hi there,

Was that the entire output? And do you have the output file from epa-ng there? Could you send some more information about the system that you're using, and also see if maybe running each of the steps separately works, or helps to narrow down what the issue might be? It looks like you're using the older verson of PICRUSt2, so that would be the instructions here.

Best wishes,

Robyn

[Simona Polizzi]

Hi! I had the same problem while running the entire pipeline and even with the place-seqs.py one. I'm sure I'm using the current release of picrust2. Did you find any solution? Thank you!

[Simona Polizzi]

Hi, I had the same issue but I'm sure I'm using the current version of picrust2 and the same error occurs while using test ASVs (so the problem is not about the sequences i'm using), have you found the solution??

Il giorno giovedì 15 maggio 2025 alle 16:24:50 UTC+2 roby...@gmail.com ha scritto:

[Robyn Wright]

Hi there,

As I mentioned above, could you try running the steps separately, as detailed here? And then please send the epa-ng output. 

Thanks,

Robyn

[Huy Nguyen]

Hi Robyn, I have the same issue when I ran Picrust2 on the example dataset. I tried to run steps separately as well but it also returned similar error. The log of epa-ng output from my run is (from epa_info.log)

INFO Selected: Output dir: placement_working_bac/epa_out/
INFO Selected: Query file: placement_working_bac/study_seqs_hmmalign.fasta
INFO Selected: Tree file: /mnt/c/Users/h.nguyen/Documents/Data/MiSeq35/picrust2_github/picrust2-2.6.2/picrust2/default_files/bacteria/bac_ref/bac_ref.tre
INFO Selected: Reference MSA: placement_working_bac/ref_seqs_hmmalign.fasta

INFO Selected: Filtering by accumulated threshold: 0.99
INFO Selected: Maximum number of placements per query: 100
INFO Selected: Automatic switching of use of per rate scalers
INFO Selected: Preserving the root of the input tree

INFO Selected: Specified model file: /mnt/c/Users/h.nguyen/Documents/Data/MiSeq35/picrust2_github/picrust2-2.6.2/picrust2/default_files/bacteria/bac_ref/bac_ref.model

INFO Selected: Reading queries in chunks of: 5000
INFO Selected: Using threads: 1
INFO     ______ ____   ___           _   __ ______
        / ____// __ \ /   |         / | / // ____/
       / __/  / /_/ // /| | ______ /  |/ // / __  
      / /___ / ____// ___ |/_____// /|  // /_/ /  
     /_____//_/    /_/  |_|      /_/ |_/ \____/ (v0.3.8)
INFO Using model parameters:

INFO    Rate heterogeneity: GAMMA (4 cats, mean),  alpha: 0.408013 (user),  weights&rates: (0.25,0.0178872) (0.25,0.187716) (0.25,0.739808) (0.25,3.05459)
        Base frequencies (user): 0.209289 0.232667 0.311531 0.246513
        Substitution rates (user): 1.04872 2.92815 1.63887 0.831209 3.68952 1
INFO Output file: placement_working_bac/epa_out/epa_result.jplace

I first tried on my own data set but it failed, so I ran on the example dataset, but it also failed with similar error at this step (running place_seqs.py, and error in epa)
Any suggestion on what I can do?

I am running on picrust2 on WSL on Windows 11.
Best,