Microarray analysis techniques

Microarray analysis techniques
Example of an approximately 40,000 probe spotted oligo microarray with enlarged inset to show detail.

Microarray analysis techniques are used in interpreting the data generated from experiments on DNA, RNA, and protein microarrays, which allow researchers to investigate the expression state of a large number of genes - in many cases, an organism's entire genome - in a single experiment. Such experiments generate a very large volume of genetic data that can be difficult to analyze, especially in the absence of good gene annotation.

Contents

Introduction

Microarray data analysis involves several distinct steps, as outlined below. Changing any one of the steps has the potential to change the outcome of the analysis, so the MAQC Project[1] was created to identify a set of standard strategies. Companies exist that use the MAQC protocols to perform a complete analysis.[2]

Creating raw data

Most microarray manufacturers, such as Affymetrix and Agilent[3], provide commercial data analysis software with microarray equipment such as plate readers.

Background correction

Depending on the type of array, signal related to nonspecific binding of the fluorophore can be subtracted to achieve better results. One approach involves subtracting the average signal intensity of the area between spots. A variety of tools for background correction and further analysis are available from TIGR,[4] Agilent (GeneSpring),[5] and Ocimum Biosolutions (Genowiz).[6]

Quality control

Entire arrays may have obvious flaws detectable by visual inspection, pairwise comparisons to arrays in the same experimental group, or by analysis of RNA degradation[7]. Results may improve by removing these arrays from the analysis entirely.

Spot filtering

Visual identification of local artifacts, such as printing or washing defects, may likewise suggest the removal of individual spots. This can take a substantial amount of time depending on the quality of array manufacture. In addition, some procedures call for the elimination of all spots with an expression value below a certain intensity threshold.

Aggregation and normalization

Comparing two different arrays, or two different samples hybridized to the same array generally involves making adjustments for systematic errors introduced by differences in procedures and dye intensity effects. Dye normalization for two color arrays is often achieved by local regression. LIMMA provides a set of tools for background correction and scaling, as well an option to average on-slide duplicate spots.[8]. A common method for evaluating how well normalized an array is, is to plot an MA plot of the data.

Raw Affy data contains about twenty probes for the same RNA target. Half of these are "mismatch spots", which do not precisely match the target sequence. These can theoretically measure the amount of nonspecific binding for a given target. Robust Multichip Average (RMA) is a normalization approach that does not take advantage of these mismatch spots, but still must summarize the perfect matches through median polish.[9] The median polish algorithm, although robust, behaves differently depending on the number of samples analyzed.[10] Quantile normalization, also part of RMA, is one sensible approach to normalize a batch of arrays in order to make further comparisons meaningful.

The current Affymetrix MAS5 algorithm, which uses both perfect match and mismatch probes, continues to enjoy popularity and do well in head to head tests.[11]

Factor Analysis for Robust Microarray Summarization (FARMS)[12] is a model-based technique for summarizing array data at perfect match probe level. It is based on a factor analysis model for which a Bayesian maximum a posteriori method optimizes the model parameters under the assumption of Gaussian measurement noise. According to the Affycomp benchmark[13] FARMS outperformed all other summarizations methods with respect to sensitivity and specificity.

Identification of significant differential expression

Many strategies exist to identify which array probes show an unusual level of over expression or under expression. The simplest one is to call "significant" any probe that differs by an average of at least twofold between treatment groups. More sophisticated approaches are often related to t-tests or other mechanisms that take both effect size and variability into account. Curiously, the p-values associated with particular genes do not reproduce well between replicate experiments, and lists generated by straight fold change perform much better.[14][15] This represents an extremely important observation, since the point of performing experiments has to do with predicting general behavior. The MAQC group recommends using a fold change assessment plus a non-stringent p-value cutoff, further pointing out that changes in the background correction and scaling process have only a minimal impact on the rank order of fold change differences, but a substantial impact on p-values.

Pattern recognition

Commercial systems for gene network analysis such as Ingenuity[16] and Pathway studio[17] create visual representations of differentially expressed genes based on current scientific literature. Non-commercial tools such as GenMAPP also aid in organizing and visualizing gene network data procured from one or several microarray experiments. A wide variety of microarray analysis tools are available through Bioconductor written in the R programming language. The frequently cited SAM Excel module and other microarray tools[18] are available through Stanford University. Another set is available from Harvard and MIT[19].

Specialized software tools for statistical analysis to determine the extent of over- or under-expression of a gene in a microarray experiment relative to a reference state have also been developed to aid in identifying genes or gene sets associated with particular phenotypes. One such method of analysis, known as Gene Set Enrichment Analysis (GSEA), uses a Kolmogorov-Smirnov-style statistic to identify groups of genes that are regulated together[20]. This third-party statistics package offers the user information on the genes or gene sets of interest, including links to entries in databases such as NCBI's GenBank and curated databases such as Biocarta[21] and Gene Ontology. Related system, PAINT[22] and SCOPE[23] performs a statistical analysis on gene promoter regions, identifying over and under representation of previously identified transcription factor response elements. Another statistical analysis tool is Rank Sum Statistics for Gene Set Collections (RssGsc), which uses rank sum probability distribution functions to find gene sets that explain experimental data[24]. A further approach is contextual meta-analysis, i.e. finding out how a gene cluster responds to a variety of experimental contexts. Genevestigator is a public tool to perform contextual meta-analysis across contexts such as anatomical parts, stages of development, and response to diseases, chemicals, stresses, and neoplasms.

References

  1. ^ Dr. Leming Shi, National Center for Toxicological Research. "MicroArray Quality Control (MAQC) Project". U.S. Food and Drug Administration. http://www.fda.gov/nctr/science/centers/toxicoinformatics/maqc/. Retrieved 2007-12-26. 
  2. ^ "GenUs BioSystems - Services - Data Analysis". http://www.genusbiosystems.com/services-data.shtml. Retrieved 2008-01-02. 
  3. ^ "Agilent | DNA Microarrays". http://www.chem.agilent.com/Scripts/PCol.asp?lPage=494. Retrieved 2008-01-02. [dead link]
  4. ^ "J. Craig Venter Institute -- Software". http://www.tigr.org/software/microarray.shtml. Retrieved 2008-01-01. 
  5. ^ "Agilent | GeneSpring GX". http://www.chem.agilent.com/scripts/pds.asp?lpage=27881. Retrieved 2008-01-02. 
  6. ^ "Ocimum Biosolutions | Genowiz". http://www.ocimumbio.com/web/bioinformatic/prod_details.asp?prod_id=20&prodType=1. Retrieved 2009-04-02. 
  7. ^ Wilson CL, Miller CJ (2005). "Simpleaffy: a BioConductor package for Affymetrix Quality Control and data analysis". Bioinformatics 21 (18): 3683–5. doi:10.1093/bioinformatics/bti605. PMID 16076888. 
  8. ^ "LIMMA Library: Linear Models for Microarray Data". http://bioinf.wehi.edu.au/limma/. Retrieved 2008-01-01. 
  9. ^ Bolstad BM, Irizarry RA, Astrand M, Speed TP (2003). "A comparison of normalization methods for high density oligonucleotide array data based on variance and bias". Bioinformatics 19 (2): 185–93. doi:10.1093/bioinformatics/19.2.185. PMID 12538238. 
  10. ^ Giorgi FM, Bolger AM, Lohse M, Usadel B (2010). "Algorithm-driven Artifacts in median polish summarization of Microarray data". BMC Bioinformatics 11: 553. doi:10.1186/1471-2105-11-553. PMC 2998528. PMID 21070630. http://www.pubmedcentral.nih.gov/articlerender.fcgi?tool=pmcentrez&artid=2998528. 
  11. ^ Lim WK, Wang K, Lefebvre C, Califano A (2007). "Comparative analysis of microarray normalization procedures: effects on reverse engineering gene networks". Bioinformatics 23 (13): i282–8. doi:10.1093/bioinformatics/btm201. PMID 17646307. 
  12. ^ Sepp Hochreiter, Djork-Arné Clevert, and Klaus Obermayer, 2006. A new summarization method for affymetrix probe level data. Bioinformatics, 22(8), 943-949. [1]
  13. ^ http://affycomp.jhsph.edu/
  14. ^ Shi L, Reid LH, Jones WD, et al. (2006). "The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements". Nat. Biotechnol. 24 (9): 1151–61. doi:10.1038/nbt1239. PMID 16964229. 
  15. ^ Guo L, Lobenhofer EK, Wang C, et al. (2006). "Rat toxicogenomic study reveals analytical consistency across microarray platforms". Nat. Biotechnol. 24 (9): 1162–9. doi:10.1038/nbt1238. PMID 17061323. 
  16. ^ "Ingenuity Systems". http://www.ingenuity.com/. Retrieved 2007-12-31. 
  17. ^ "Ariadne Genomics: Pathway Studio". http://www.ariadnegenomics.com/products/pathway-studio/. Retrieved 2007-12-31. 
  18. ^ ["Significance Analysis of Microarrays". http://www-stat.stanford.edu/~tibs/SAM/. Retrieved 2007-12-31. ]
  19. ^ "Software - Broad". http://www.broad.mit.edu/tools/software.html. Retrieved 2007-12-31. 
  20. ^ Subramanian A, Tamayo P, Mootha VK, et al. (2005). "Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles". Proc. Natl. Acad. Sci. U.S.A. 102 (43): 15545–50. doi:10.1073/pnas.0506580102. PMC 1239896. PMID 16199517. http://www.pubmedcentral.nih.gov/articlerender.fcgi?tool=pmcentrez&artid=1239896. 
  21. ^ "BioCarta - Charting Pathways of Life". http://www.biocarta.com/. Retrieved 2007-12-31. ]
  22. ^ "DBI Web". http://www.dbi.tju.edu/dbi/staticpages.php?page=tools&menu=37. Retrieved 2007-12-31. 
  23. ^ "SCOPE". http://genie.dartmouth.edu/scope/. Retrieved 2007-12-31. 
  24. ^ "RssGsc". http://rssgsc.sourceforge.net/. Retrieved 2008-10-15. 


External links


Wikimedia Foundation. 2010.

Игры ⚽ Поможем сделать НИР

Look at other dictionaries:

  • DNA microarray — A DNA microarray (also commonly known as gene chip, DNA chip, or biochip) is a collection of microscopic DNA spots attached to a solid surface. Scientists use DNA microarrays to measure the expression levels of large numbers of genes… …   Wikipedia

  • Factor analysis — is a statistical method used to describe variability among observed, correlated variables in terms of a potentially lower number of unobserved, uncorrelated variables called factors. In other words, it is possible, for example, that variations in …   Wikipedia

  • Significance analysis of microarrays — (SAM) is a statistical technique, established in 2001 by Tusher, Tibshirani and Chu, for determining whether changes in gene expression are statistically significant. With the advent of DNA microarrays it is now possible to measure the expression …   Wikipedia

  • Chemical compound microarray — A chemical compound microarray is a collection of organic chemical compounds spotted on a solid surface, such as glass and plastic. This microarray format is very similar to DNA microarray, protein microarray and antibody microarray. In chemical… …   Wikipedia

  • Principal component analysis — PCA of a multivariate Gaussian distribution centered at (1,3) with a standard deviation of 3 in roughly the (0.878, 0.478) direction and of 1 in the orthogonal direction. The vectors shown are the eigenvectors of the covariance matrix scaled by… …   Wikipedia

  • Principal components analysis — Principal component analysis (PCA) is a vector space transform often used to reduce multidimensional data sets to lower dimensions for analysis. Depending on the field of application, it is also named the discrete Karhunen Loève transform (KLT),… …   Wikipedia

  • Gene expression profiling — Heat maps of gene expression values show how experimental conditions influenced production (expression) of mRNA for a set of genes. Green indicates reduced expression. Cluster analysis has placed a group of down regulated genes in the upper left… …   Wikipedia

  • Bisulfite sequencing — Figure 1: Outline of bisulfite conversion of sample sequence of genomic DNA. Nucleotides in blue are unmethylated cytosines converted to uracils by bisulfite, while red nucleotides are 5 methylcytosines resistant to conversion …   Wikipedia

  • DecisionDx-UM — is a prognostic test that accurately determines the metastatic risk associated with ocular melanoma tumors of the eye. Ocular melanoma is a term commonly used to describe tumors of the uveal tract such as uveal melanoma, choroidal melanoma,… …   Wikipedia

  • Oncogenomics — is relatively new sub field of genomics, which applies high throughput technologies to characterize genes associated with cancer. Oncogenomics is synonymous with cancer genomics . Cancer is a genetic disease caused by accumulation of mutations to …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”