Question: How to determine parent of origin of phased variants when all members are het
0
gravatar for QVINTVS_FABIVS_MAXIMVS
11 months ago by
USA SoCal
QVINTVS_FABIVS_MAXIMVS2.2k wrote:

I'm a little confused on the format of phased variants in a VCF file. I was under the assumption that A|B is equal to Paternal Haplotype | Maternal Haplotype but I believe it's Haplotype Block A | Haplotype Block B.

So in the case of

KID    DAD    MOM

0|1     0|1      0|1

How does one determine the parent of origin for the derived allele?

phase vcf • 531 views
ADD COMMENTlink modified 11 months ago by WouterDeCoster35k • written 11 months ago by QVINTVS_FABIVS_MAXIMVS2.2k
2
gravatar for WouterDeCoster
11 months ago by
Belgium
WouterDeCoster35k wrote:

I believe it's Haplotype Block A | Haplotype Block B.

Indeed. Wherever possible using read-based evidence variants are phased in blocks, but e.g. across the repetitive elements, this information is lost. So this does not correspond to maternal and paternal haplotypes and needs long distance phasing (statistical, or using different sequencing technologies) to get longer 'blocks'

Based on your example it's not possible to determine the parent of origin, but perhaps SNPs in the same block can help nailing this down.

ADD COMMENTlink modified 11 months ago • written 11 months ago by WouterDeCoster35k

Lame. I was afraid this was the case. The phasing leveraged nanopore long reads, so the blocks are really long. You would think that phasing software that uses inheritance information would encode the parent of origin. Thanks for the insight. To add, if anyone knows of a software that can assign parent of origin to ambiguous cases like this, please comment.

ADD REPLYlink modified 11 months ago • written 11 months ago by QVINTVS_FABIVS_MAXIMVS2.2k

phasing software that uses inheritance information

That's new information, you didn't specify this before.

It should be feasible to identify the parent for each block if those are sufficiently big. You just have to find variant(s) specific to one parent to assign the entire block.

ADD REPLYlink written 11 months ago by WouterDeCoster35k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1646 users visited in the last hour