296
edits
Changes
no edit summary
==How to obtain PGx allele definitions from literature==PGx alleles are collected and distributed through various channels* PGx alleles on JSON-LD format from the [https://api.pharmgkb.org PharmGKB API]* PGx alleles on Excel-style formats from the PGRN Translational Pharmacogenetics Project (TPP), published through the [https://api.pharmgkb.org PharmGKB API]* PGx alleles as VCF files from [https://www.pharmvar.org PharmVar]PGx allele definitions are given in either GRCh37 or GRCh38 reference coordinates. The process of changing from GRCh37 to GRCh38 for PharmGKB is only partially completed.We suspect that there are serious discrepancies between PharmVar and PharmGKB sources (which is strange as PharmGKB gets data form PharmVar). ==How to define PGx allelesfor next generation sequencing==
PGx alleles are defined as collections of one or more SNPs, INDELs or structural variants. When a patient is sequenced by next generation sequencing ([[NGS]]) we may typically observe more variants than those which are included in any individual PGx allele definitions.
[[File:Variant tree outline.png|thumb|The 16 possible haplotypes for a four loci, decomposed variant calling]]
We illustrate some of the problems that we encountered when trying to match patient haplotypes to the PGx allele definitions, by a four loci PGx gene
===The SNP array method===
[[File:Variant tree allele snp definition.png|thumb|PGx alleles defined as collections of variants, with no requirement on loci that are not part of the definition, will assign the same PGx allele to several different haplotypes]]
This definition only requires matches for variants explicitly included in PGx allele definitions.
*But the presence of additional variants will have no effect on reported PGx alleles
===The PharmCAT method===
This definition requires matches also for variants not explicitly included in PGx allele definitions.
[[File:Variant tree allele pharmcat definition.png|thumb|PGx alleles defined as complete haplotypes classifies the patient uniquely]]
(Note that in practice PharmCAT lets the user decide which allele definitions to use in their [https://github.com/PharmGKB/PharmCAT/wiki/NamedAlleleMatcher-101 NamedAlleleMatcher])
===Which definition should we stick to?===
{| class="wikitable"
|-