How to check the allele strand of variants?
1
0
Entering edit mode
5.1 years ago
seta ★ 1.9k

Hi all,

I have a list of variants obtained from various GWAS studies, it is important for me that all variant alleles (alternate allele) located on the forward strand. Could you please kindly tell me how I can check this issue to make sure forward strand allele for all variants?

Thanks

variant • 3.5k views
ADD COMMENT
0
Entering edit mode

Do you mean forward relative to the gene in which the variants are, or relative to the chromosome?

Nevertheless, variants (in GWAS) are most likely on the "sense" direction ('+') of a chromosome. It is also unclear what you aim to achieve with these variants (why is the strand important?) and in which format your data is (plink? vcf?).

ADD REPLY
0
Entering edit mode

C: Inconsistencies between snp coordinates and alleles across references

Ensembl variants are always mapped to the positive strand. dbSNP maps to either but this is usually indicated in the record.

ADD REPLY
7
Entering edit mode
5.1 years ago
Emily 23k
  1. Do the paper authors tell you which strand they've quoted their alleles on?
  2. Look up the variants in public databases. Ensembl always quotes the forward strand alleles. If the paper talks about G/T and the database talks about C/A, then they've given the reverse strand alleles – flip them. Unfortunately if the paper talks about A/T and the database talks about T/A, then you can't always tell if they've used reverse strand, or just flipped reference and alternative.
  3. Look at the gene the paper discussed with reference to the variant. If they are talking about a reverse strand gene, they'll usually be talking about reverse strand variant alleles.
  4. If it's still not clear, email the paper authors to ask.

This is not an easy problem, because we're talking about human behaviour, rather than proper data conventions. This is why the GWAS catalog (and we at Ensembl when we import their data) have a policy of quoting exactly the alleles given in the paper, rather than trying to convert. If we claim to convert then it's expected that we will always convert correctly, whereas if we don't convert then we aren't responsible when we misinterpret a G/C or A/T.

ADD COMMENT
0
Entering edit mode

Thank you very much for your complete answer.

ADD REPLY
0
Entering edit mode

If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
Upvote|Bookmark|Accept

ADD REPLY
0
Entering edit mode

Thanks :)

ADD REPLY

Login before adding your answer.

Traffic: 1855 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6