Equalizer reduces SNP bias in Affymetrix microarrays

BMC Bioinformatics

David Quigley¹ 

[1] Helen Diller Family Comprehensive Cancer Center, University of California at San Francisco, San Francisco 94158, CA, USA

关键词: Single nucleotide polymorphism;    eQTL;    Affymetrix;    Microarray;   

Others  :  1230600
DOI  :  10.1186/s12859-015-0669-y

 received in 2015-04-06, accepted in 2015-07-09,  发布年份 2015

【 摘 要 】

Background

Gene expression microarrays measure the levels of messenger ribonucleic acid (mRNA) in a sample using probe sequences that hybridize with transcribed regions. These probe sequences are designed using a reference genome for the relevant species. However, most model organisms and all humans have genomes that deviate from their reference. These variations, which include single nucleotide polymorphisms, insertions of additional nucleotides, and nucleotide deletions, can affect the microarray’s performance. Genetic experiments comparing individuals bearing different population-associated single nucleotide polymorphisms that intersect microarray probes are therefore subject to systemic bias, as the reduction in binding efficiency due to a technical artifact is confounded with genetic differences between parental strains. This problem has been recognized for some time, and earlier methods of compensation have attempted to identify probes affected by genome variants using statistical models. These methods may require replicate microarray measurement of gene expression in the relevant tissue in inbred parental samples, which are not always available in model organisms and are never available in humans.

Results

By using sequence information for the genomes of organisms under investigation, potentially problematic probes can now be identified a priori. However, there is no published software tool that makes it easy to eliminate these probes from an annotation. I present equalizer, a software package that uses genome variant data to modify annotation files for the commonly used Affymetrix IVT and Gene/Exon platforms. These files can be used by any microarray normalization method for subsequent analysis. I demonstrate how use of equalizer on experiments mapping germline influence on gene expression in a genetic cross between two divergent mouse species and in human samples significantly reduces probe hybridization-induced bias, reducing false positive and false negative findings.

Conclusions

The equalizer package reduces probe hybridization bias from experiments performed on the Affymetrix microarray platform, allowing accurate assessment of germline influence on gene expression.

【 授权许可】

   
2015 Quigley.

【 图 表 】

Fig. 1.

Fig. 2.

Fig. 3.

【 参考文献 】

• [1]Brem RB, Yvert G, Clinton R, Kruglyak L. Genetic dissection of transcriptional regulation in budding yeast. Science. 2002; 296(5568):752-5.
• [2]Chesler EJ, Lu L, Shou S, Qu Y, Gu J, Wang J, Hsu HC, Mountz JD, Baldwin NE, Langston MA et al.. Complex trait analysis of gene expression uncovers polygenic and pleiotropic networks that modulate nervous system function. Nat Genet. 2005; 37(3):233-42.
• [3]Chen Y, Zhu J, Lum PY, Yang X, Pinto S, MacNeil DJ, Zhang C, Lamb J, Edwards S, Sieberts SK et al.. Variations in DNA elucidate molecular networks that cause disease. Nature. 2008; 452(7186):429-35.
• [4]Quigley DA, To MD, Perez-Losada J, Pelorosso FG, Mao JH, Nagase H, Ginzinger DG, Balmain A. Genetic architecture of mouse skin inflammation and tumour susceptibility. Nature. 2009; 458(7237):505-8.
• [5]Doss S, Schadt EE, Drake TA, Lusis AJ. Cis-acting expression quantitative trait loci in mice. Genome Res. 2005; 15(5):681-91.
• [6]Alberts R, Terpstra P, Li Y, Breitling R, Nap JP, Jansen RC. Sequence polymorphisms cause many false cis eQTLs. PLoS One. 2007; 2(7):e622.
• [7]Ciobanu DC, Lu L, Mozhui K, Wang X, Jagalur M, Morris JA, Taylor WL, Dietz K, Simon P, Williams RW. Detection, validation, and downstream analysis of allelic variation in gene expression. Genetics. 2010; 184(1):119-28.
• [8]Ramasamy A, Trabzuni D, Gibbs JR, Dillman A, Hernandez DG, Arepalli S, Walker R, Smith C, Ilori GP, Shabalin AA et al.. Resolving the polymorphism-in-probe problem is critical for correct interpretation of expression QTL studies. Nucleic Acids Res. 2013; 41(7):e88.
• [9]Dannemann M, Lorenc A, Hellmann I, Khaitovich P, Lachmann M. The effects of probe binding affinity differences on gene expression measurements and how to deal with them. Bioinformatics. 2009; 25:2772.
• [10]Dannemann M, Lachmann M, Lorenc A. ‘maskBAD’--a package to detect and remove Affymetrix probes with binding affinity differences. BMC Bioinformatics. 2012; 13:56. BioMed Central Full Text
• [11]Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26(6):841-2.
• [12]Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010; 26(19):2363-7.
• [13]Quigley DA, To MD, Kim IJ, Lin KK, Albertson DG, Sjolund J, Perez-Losada J, Balmain A. Network analysis of skin tumor progression identifies a rewired genetic architecture affecting inflammation and tumor susceptibility. Genome Biol. 2011; 12(1):R5. BioMed Central Full Text
• [14]Keane TM, Goodstadt L, Danecek P, White MA, Wong K, Yalcin B, Heger A, Agam A, Slater G, Goodson M et al.. Mouse genomic variation and its effect on phenotypes and gene regulation. Nature. 2011; 477(7364):289-94.
• [15]Wong K, Bumpstead S, Van Der Weyden L, Reinholdt LG, Wilming LG, Adams DJ, Keane TM. Sequencing and characterization of the FVB/NJ mouse genome. Genome Biol. 2012; 13(8):R72. BioMed Central Full Text
• [16]Sjolund J, Pelorosso FG, Quigley DA, DelRosario R, Balmain A. Identification of Hipk2 as an essential regulator of white fat development. Proc Natl Acad Sci U S A. 2014; 111(20):7373-8.
• [17]Dejager L, Libert C, Montagutelli X. Thirty years of Mus spretus: a promising future. Trends Genet. 2009; 25(5):234-41.
• [18]1000 Genomes Project Consortium, Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56–65.
• [19]Li L, Fridley B, Kalari K, Jenkins G, Batzler A, Safgren S, Hildebrandt M, Ames M, Schaid D, Wang L. Gemcitabine and cytosine arabinoside cytotoxicity: association with lymphoblastoid cell expression. Cancer Res. 2008; 68(17):7050-8.
• [20]Bach FH, van Rood JJ. The major histocompatibility complex--genetics and biology. (First of three parts). N Engl J Med. 1976; 295(15):806-13.
• [21]R Core Team. R: A Language and Environment for Statistical Computing. Vienna; 2012.
• [22]Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007; 8(1):118-27.