Normalization of Illumina Infinium whole-genome SNP data improves copy number estimates and allelic intensity ratios

BackgroundIllumina Infinium whole genome genotyping (WGG) arrays are increasingly being applied in cancer genomics to study gene copy number alterations and allele-specific aberrations such as loss-of-heterozygosity (LOH). Methods developed for normalization of WGG arrays have mostly focused on diploid, normal samples. However, for cancer samples genomic aberrations may confound normalization and data interpretation. Therefore, we examined the effects of the conventionally used normalization method for Illumina Infinium arrays when applied to cancer samples.ResultsWe demonstrate an asymmetry in the detection of the two alleles for each SNP, which deleteriously influences both allelic proportions and copy number estimates. The asymmetry is caused by a remaining bias between the two dyes used in the Infinium II assay after using the normalization method in Illumina's proprietary software (BeadStudio). We propose a quantile normalization strategy for correction of this dye bias. We tested the normalization strategy using 535 individual hybridizations from 10 data sets from the analysis of cancer genomes and normal blood samples generated on Illumina Infinium II 300 k version 1 and 2, 370 k and 550 k BeadChips. We show that the proposed normalization strategy successfully removes asymmetry in estimates of both allelic proportions and copy numbers. Additionally, the normalization strategy reduces the technical variation for copy number estimates while retaining the response to copy number alterations.ConclusionThe proposed normalization strategy represents a valuable tool that improves the quality of data obtained from Illumina Infinium arrays, in particular when used for LOH and copy number variation studies.

[1]  Hanna Göransson,et al.  Screening for copy‐number alterations and loss of heterozygosity in chronic lymphocytic leukemia—A comparative study of four differently designed, high resolution microarray platforms , 2008, Genes, chromosomes & cancer.

[2]  Emmanuel Barillot,et al.  Spatial normalization of array-CGH data , 2006, BMC Bioinformatics.

[3]  Cisca Wijmenga,et al.  High-resolution copy number analysis of paraffin-embedded archival tissue using SNP BeadArrays. , 2007, Genome research.

[4]  Rafael A Irizarry,et al.  Exploration, normalization, and genotype calls of high-density oligonucleotide SNP array data. , 2006, Biostatistics.

[5]  Michael G. Barnes,et al.  Experimental comparison and cross-validation of the Affymetrix and Illumina gene expression analysis platforms , 2005, Nucleic acids research.

[6]  Terence P. Speed,et al.  A comparison of normalization methods for high density oligonucleotide array data based on variance and bias , 2003, Bioinform..

[7]  S. P. Fodor,et al.  Genotyping over 100,000 SNPs on a pair of oligonucleotide arrays , 2004, Nature Methods.

[8]  S. Dudoit,et al.  Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. , 2002, Nucleic acids research.

[9]  John Quackenbush.  Microarray data normalization and transformation , 2002, Nature Genetics.

[10]  E. S. Venkatraman,et al.  A faster circular binary segmentation algorithm for the analysis of array CGH data , 2007, Bioinform..

[11]  K. Gunderson,et al.  High-resolution genomic profiling of chromosomal aberrations using Infinium whole-genome genotyping. , 2006, Genome research.

[12]  S. Sen,et al.  Aneuploidy and cancer , 2000, Current opinion in oncology.

[13]  Joseph T. Glessner,et al.  PennCNV: an integrated hidden Markov model designed for high-resolution copy number variation detection in whole-genome SNP genotyping data. , 2007, Genome research.

[14]  Weihua Chang,et al.  Whole-genome genotyping with the single-base extension assay , 2005, Nature Methods.

[15]  Cheng Li,et al.  dChipSNP: significance curve and clustering of SNP-array-based loss-of-heterozygosity data , 2004, Bioinform..

[16]  L. Chin,et al.  A comparison of DNA copy number profiling platforms. , 2008, Cancer research.

[17]  T. LaFramboise,et al.  SNP arrays in heterogeneous tissue: highly accurate collection of both germline and somatic genetic information from unpaired single tumor samples. , 2008, American journal of human genetics.

[18]  Rabab Kreidieh Ward,et al.  BMC Bioinformatics Methodology article A stepwise framework for the normalization of array CGH data , 2005 .

[19]  Simon Tavaré,et al.  Statistical issues in the analysis of Illumina data , 2008, BMC Bioinformatics.

[20]  K. Gunderson,et al.  A genome-wide scalable SNP genotyping assay using microarray technology , 2005, Nature Genetics.

[21]  Terry Speed,et al.  Normalization of cDNA microarray data. , 2003, Methods.

[22]  M. Ringnér,et al.  Segmentation-based detection of allelic imbalance and loss-of-heterozygosity in cancer cells using whole genome SNP arrays , 2008, Genome Biology.

[23]  B J Williams,et al.  Comparative genomic hybridization. , 1996, Methods in molecular medicine.

[24]  Johan Staaf,et al.  Normalization of array-CGH data: influence of copy number imbalances , 2007, BMC Genomics.