MOTIVATION
Preliminary results on the data produced using the Affymetrix large-scale genotyping platforms show that it is necessary to construct improved genotype calling algorithms. There is evidence that some of the existing algorithms lead to an increased error rate in heterozygous genotypes, and a disproportionately large rate of heterozygotes with missing genotypes. Non-random errors and missing data can lead to an increase in the number of false discoveries in genetic association studies. Therefore, the factors that need to be evaluated in assessing the performance of an algorithm are the missing data (call) and error rates, but also the heterozygous proportions in missing data and errors.
RESULTS
We introduce a novel genotype calling algorithm (GEL) for the Affymetrix GeneChip arrays. The algorithm uses likelihood calculations that are based on distributions inferred from the observed data. A key ingredient in accurate genotype calling is weighting the information that comes from each probe quartet according to the quality/reliability of the data in the quartet, and prior information on the performance of the quartet.
AVAILABILITY
The GEL software is implemented in R and is available by request from the corresponding author at nicolae@galton.uchicago.edu.
[1]
Felix Naef,et al.
Absolute mRNA concentrations from sequence-specific calibration of oligonucleotide arrays.
,
2003,
Nucleic acids research.
[2]
David Kulp,et al.
Model-P: a basecalling method for resequencing microarrays of diploid samples
,
2005,
ECCB/JBI.
[3]
Lee J. Bain,et al.
A Two-Sample Test of Equal Gamma Distribution Scale Parameters With Unknown Common Shape Parameter
,
1983
.
[4]
Cheng Li,et al.
Allele-Specific Amplification in Cancer Revealed by SNP Array Analysis
,
2005,
PLoS Comput. Biol..
[5]
Terence P. Speed,et al.
Genome analysis A genotype calling algorithm for affymetrix SNP arrays
,
2005
.
[6]
Jing Huang,et al.
Dynamic model based algorithms for screening and genotyping over 100K SNPs on oligonucleotide microarrays
,
2005,
Bioinform..
[7]
S. P. Fodor,et al.
Genotyping over 100,000 SNPs on a pair of oligonucleotide arrays
,
2004,
Nature Methods.