Question: Once I Got The Snparray Genotyped , What'S The Next Step
0
gravatar for edufoju
4.9 years ago by
edufoju10
Spain
edufoju10 wrote:

I just genotyped an SNP array from .CEL file using the R package "crlmm", now I have a variable called crlmmResult, which is the .CEL file genotyped ( a big file with samples and the SNPs genotype (1, 2 or 3).

Can someone help me in what is the next step for processing the SNP array (which R packages to use)? If I want to know what SNPs are related to disease for example.

Also I don't understand exactly what the genotype mean in SNP, AA means no change and BB means there's a SNP?

Thanks.

ADD COMMENTlink modified 4.9 years ago by Neilfws48k • written 4.9 years ago by edufoju10
2
gravatar for Devon Ryan
4.9 years ago by
Devon Ryan88k
Freiburg, Germany
Devon Ryan88k wrote:

Since you already have everything in R, you might try the R interface to the Ensembl Variant Effect Predictor, which will give you much of the information you're after. Regarding the genotyping, AA is usually homozygous for the reference, whereas BB would be homozygous for the variant (there can, of course, be more than one variant).

ADD COMMENTlink written 4.9 years ago by Devon Ryan88k

I think this is the package I am searching for, thank you!. But how I convert the .CEL file to .VCF, which is the file I need to input to the package.

ADD REPLYlink written 4.9 years ago by edufoju10

Presumably there's a convenient way to convert the genotype calls to a VCF. From googling around a bit, it looks like you can at least make .ped and .map files (for plink), which can be converted to VCF. I suspect there's a better way to do that, but I'm not familiar with the crlmm package.

Having said that, since you're using an array (and therefore not looking for new mutations), you might consider the GWAS approach outlined in the link in neilfws' reply.

ADD REPLYlink written 4.9 years ago by Devon Ryan88k
0
gravatar for Neilfws
4.9 years ago by
Neilfws48k
Sydney, Australia
Neilfws48k wrote:

The next step is to read all of the crlmm documentation. It seems that you have read the genotyping vignette. Next, try the downstream analysis vignette, which explains the basics of conversion to a SNPset for GWAS.

ADD COMMENTlink written 4.9 years ago by Neilfws48k
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: 1400 users visited in the last hour