This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to identify which genome the variations in the pan-genome graph originate from

We want to use vg to construct a pan-genome. In order to detect variations in the subsequent pan-genome graph and determine which genomes they exist in, we intend to retain the genome source information for the paths within the pan-genome graph. We attempted to add genome information to the ID and info columns in the VCF input provided to vg. However, after completing the construction process, the ID and info columns in the resulting pan-genome graph were replaced with randomly generated identifiers, and the corresponding genome information couldn't be found. How should we proceed?

vg
• 3,124 views • link
written by Wenke • Wenke

2 answers

vg has no mechanism for preserving the VCF IDs, unfortunately. You will need to use position lookups (vg find) to compare your VCF and graph.

If your VCF is phased, vg will store the haplotypes in your VCF in the graph by way of the GBWT/GBZ which will allow you to do some times of queries on your samples.

• 0 views • link
written by glenn.hickey glenn.hickey

Thank you! What does "VCF is phased" mean ? Does it mean to add the allele types at mutation positions of different genomes?

• 0 views • link
written by Wenke • Wenke

A phased VCF specifies which variants co-occur on the same haplotype, whereas an unphased VCF only lists genotypes at each site, with no assertion of which combination is on each haplotype. In the VCF format, phasing is usually indicated by using a | symbol between the alleles instead of a /.

• 0 views • link
written by Jordan M Eizenga • Jordan M Eizenga

Is there any way to know which specific path in the resequencing alignment has been aligned

• 0 views • link
written by Wenke • Wenke

If you are referring to the sequence of node IDs that the read is aligned to, then yes. The GAM and GAF formats both specify the path. GAF is a text-based format and is pretty easy to inspect manually to see the path. GAM is based on Protobuf, but if you want to inspect some reads, you can convert it to JSON with vg view -a -j. However, the JSON option is not very efficient at genome scale.

• 0 views • link
written by Jordan M Eizenga • Jordan M Eizenga

Is there a way to construct a pan-genome graph using vg to create a GFA graph with "W"

• 0 views • link
written by Wenke • Wenke

If you use PGGB to build a pangenome from FASTAs ( so not from fasta + VCF) you get a resultant VCF out with the SNPs in each pangenome. Of course, theres' no guarantee that it will find the SNPs you're expecting, but there must be some overlap if they are present in the genomes.

• 0 views • link
written by colindaven colindaven

Log in to answer this question.

More posts like this