Question: Illumina exome beadchip duplicated positions
gravatar for miaowzai
4 months ago by
United States
miaowzai60 wrote:

I have a plink format data (bim/bed/fam) from data provider. As far as I know, the sequencing was done by Illumina exome SNP beadchip. I didn't call the genotypes.

There are some variant IDs startwith "exm". They do have rs numbers when I look'em up. What are these SNPs? Why use "exm"?

Also, I found some duplicated monomorphic positions, like:

1 exm1771279 2435759 0 G

1 exm2277033 2435759 0 C

I checked the genotypes for both calls and the file says all people have GG for exm1771279 and all people have CC for exm2277033. However, they are at the same site. How is this possible?


beadchip snp illumina exome • 189 views
ADD COMMENTlink modified 4 months ago • written 4 months ago by miaowzai60

It is caused by the allele is from different strand. GG actually == CC they are from Watson and Crick strand, respectively.

Now you know the above situation, the follow situation should also be paid attention:

When merging two datasets, it is clearly very important that the two sets of SNPs are concordant in terms of positive or negative strand. Whereas some mismatches will be easy to spot as more than two alleles will be observed in the merged dataset, other instances will not be so easy to spot, i.e. for A/T and C/G SNPs.

Be sure to Flip DNA strand for SNPs when it is need in your analysis.

Plink only use ID to be the key. Here exm1771279 and exm2277033 are allowed in the analysis. but usually, for me, I will transfer exm id to rs id, and therefore, in my anlaysis, both of these two allele will be mapped to one rs SNP and therefore one of them will be removed in my further analysis. It is very important to do that.

ADD REPLYlink modified 4 months ago • written 4 months ago by Shicheng Guo6.1k

This makes a lot of sense! Thank you so much.

ADD REPLYlink written 4 months ago by miaowzai60
Please log in to add an answer.


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