This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Practical Haplotype Graph v2 not finding correct paths

Hello, I am using Practical Haplotype Graph v2.2.85.134 to build a pangenome graph using six diploid plant species (12 haplotypes). I was able to go through their Build and Load module. Then, I simulated 10 million WGS Illumina Novaseq paired end reads from two of my haplotypes (australasica_primary and fallglo_primary) and mapped it to the pangenome using their Imputation module. The idea was that the most haplotype paths in the imputed vcf file should originate from those two haplotypes however that's not the case. Below are the number of haplotype paths from each of my haplotypes.

  • 5,232 australasica_alternate
  • 325 australasica_primary
  • 11,390 australis_alternate
  • 14,105 australis_primary
  • 1,439 fallglo_alternate
  • 10,086 fallglo_primary
  • 12,617 fortune_alternate
  • 5 inodora_alternate
  • 17 inodora_primary

I understand there are a lot of moving parts to this pipeline and I would be happy to provide more details if requested. Any suggestion is highly appreciated. Thank you.

pangenome phg graph
• 1,757 views • link
written by beantkapoor16 • beantkapoor16

2 answers

There have been a couple of bugs found that affect imputation accuracy, one in build-kmer-index and one in find-paths. The find-paths bug only affects haploid path finding. Those should be fixed in the next few days. I will post here when the fix is available.

• 0 views • link
written by pjb39 • pjb39

Hello, just wanted to check in if the code has been fixed yet? Thank you.

• 0 views • link
written by beantkapoor16 • beantkapoor16

The timing for your question is excellent. The fix was released yesterday.

• 0 views • link
written by pjb39 • pjb39

That's great. I will test it and let you know how it goes. Thank you.

• 0 views • link
written by beantkapoor16 • beantkapoor16

So, I just tested it and it's performing much better now. Most of the haplotype paths are coming from the expected haplotypes (australasica primary and fallglo primary).

  • 302 australasica_alternate
  • 32,242 australasica_primary
  • 60 australis_alternate
  • 19 australis_primary
  • 1616 fallglo_alternate
  • 11,797 fallglo_primary
  • 7,662 fortune_alternate
  • 32 inodora_alternate
  • 33 inodora_primary
  • 6,805 wilking_alternate
  • 6,972 wilking_primary

My follow up question is how can we further improve imputation accuracy with our PHG database? Are there any parameters related to Anchorwave aligning step that we should try changing? Any suggestion is greatly appreciated. Thank you!

• 0 views • link
written by beantkapoor16 • beantkapoor16

Log in to answer this question.

More posts like this

No similar posts found.