Thursday, June 22, 2017

SNP allele discrepancies due to dbSNP updates and Orientation - SNPedia

https://www.snpedia.com/index.php/Orientation

By standard convention, the locations of DNA variations (SNPs) are based on their chromosomal position. This number changes every time a new reference human genome assembly is released; the current "build" is GRCh38, and the previous one, which is still used by many sources for a variety of reasons, is GRCh37 (also known as hg19).

Genes are read (transcribed) in either the forward or reverse direction, as numbered along the chromosome. If a new build comes along that flips a large segment of a chromosome, the gene direction will change. As a result, at different times, as well as in different publications or different databases, the same SNP can be defined as being on the forward (plus) or as being on the reverse (minus) strand. In terms of the nucleotides for that SNP, the pairing of A with T,and C with G, in the DNA double helix means that an A on the plus strand by definition is a T on the minus strand, and vice versa, and a C on the plus strand means a G on the minus strand (and vice versa). The way a SNP is defined is based on its flanking sequences, so barring processing errors at dbSNP, the major and minor alleles of a SNP should not change between builds, however, the position number along the chromosome will change, and whether it is on the plus or minus strand may change.

SNPedia currently deals with these complexities using two fields, Orientation, and, StabilizedOrientation, which are usually seen on each SNP (rs#)page at the top of the sidebar on the right of the page. Orientation indicates the orientation reported in the most current build (currently GRCh38), which is reported below the genotypes in the Reference field. The chromosome number and nucleotide position in that reference build is shown next. StabilizedOrientation is the orientation that is relevant to the genotypes that have been defined in SNPedia for each SNP, allowing Promethease to minimize confusion despite how reference builds may change over time or how companies may report genotypes in a person's raw data. [Both fields are carefully controlled and you should not edit them; if you believe either is incorrect for a given entry, report it to us and we will investigate.]

Companies follow their own protocols, which are often different. 23andMe currently reports all genotype data based on the plus strand of GRCh37, whether or not dbSNP defined the SNP as being on the plus strand (in that build or any other build); they explain this for their customers here, here and here.

This often leads to confusion since a 23andMe customer may see a genotype in their raw data that will not match the genotypes defined by dbSNP, in SNPedia, or in their Promethease report. This will most often happen when the StabilizedOrientation is minus. In these cases, the alleles need to be "flipped" to match:

  • A->T
  • T->A
  • C->G
  • G->C

This can still lead to some ambiguous flips.

In a Promethease report, this is done automatically as Promethease correctly flips genotypes from each company's orientation into the orientation used by SNPedia.


discussion at reddit[edit]

Changes to move away from Orientation are coming soon https://www.reddit.com/r/promethease/comments/5vhh14/orientation_gut_check/

https://www.reddit.com/r/promethease/comments/3ayg64/orientation_confusion/

Excerpt from a real email exchange[edit]

I noticed that there are errors in my report. SNP rs1056836 was flagged by Genetic Genie as abnormal (GG) came back as normal (CC) on your report. I confirmed that the GG is correct in my raw data on 23andMe. Also, rs651852 is TT in my raw data and AA on my Promethease report. This is just in searching 4 SNPs that came back abnormally. Do you have any insight on this?

This is due to orientation issues. SNPedia and Promethease are more correct than any other source that you have seen.
I'm really new to all of this, and I'm not quite understanding. Are you saying that the alleles reported don't matter in these cases and that the interpretation of "normal" by Promethease is more correct than Genetic Genie's interpretation as "abnormal"?
The alleles matter, but you need to take orientation into account. Promethease goes to great lengths to do this; other sites don't.


No comments:

Post a Comment

20–33% PwP have mild cognitive impairment (MCI) at time of diagnosis, but 83 percent of those still living with the disease after 20 years have dementia. | Parkinson’s Dementia: Stages, Behavior, Symptoms, and More

Dementia is a condition that causes a decline in thinking, reasoning, and problem-solving.  About 20–33% PwP have  mild cognitive impairment...