Microarray analysis techniques
Microarray analysis techniques are used in interpreting the data generated from experiments on DNA (Gene chip analysis), 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.[1] Such experiments can generate very large amounts of data, allowing researchers to assess the overall state of a cell or organism. Data in such large quantities is difficult – if not impossible – to analyze without the help of computer programs.
Introduction
Microarray data analysis is the final step in reading and processing data produced by a microarray chip. Samples undergo various processes including purification and scanning using the microchip, which then produces a large amount of data that requires processing via computer software. It involves several distinct steps, as outlined in the image below. Changing any one of the steps will change the outcome of the analysis, so the MAQC Project[2] was created to identify a set of standard strategies. Companies exist that use the MAQC protocols to perform a complete analysis.[3]
Techniques
Most microarray manufacturers, such as
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 as an option to average on-slide duplicate spots.[5] A common method for evaluating how well normalized an array is, is to plot an MA plot of the data. MA plots can be produced using programs and languages such as R and MATLAB.[6][7]
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 Multi-array Average (RMA)[8] 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 array probes that 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
Clustering
Clustering is a data mining technique used to group genes having similar expression patterns. Hierarchical clustering, and k-means clustering are widely used techniques in microarray analysis.
Hierarchical clustering
Hierarchical clustering is a statistical method for finding relatively
- Single linkage (minimum method, nearest neighbor)
- Average linkage (UPGMA)
- Complete linkage (maximum method, furthest neighbor)
Different studies have already shown empirically that the Single linkage clustering algorithm produces poor results when employed to gene expression microarray data and thus should be avoided.[18][19]
K-means clustering
K-means clustering is an algorithm for grouping genes or samples based on pattern into K groups. Grouping is done by minimizing the sum of the squares of distances between the data and the corresponding cluster
Pattern recognition
Commercial systems for gene network analysis such as Ingenuity
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
Significance analysis of microarrays (SAM)
Significance analysis of microarrays (SAM) is a statistical technique, established in 2001 by Virginia Tusher, Robert Tibshirani and Gilbert 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 of thousands of genes in a single hybridization experiment. The data generated is considerable, and a method for sorting out what is significant and what isn't is essential. SAM is distributed by Stanford University in an R-package.[31]
SAM identifies statistically significant genes by carrying out gene specific
Basic protocol
- Perform microarray experiments — DNA microarray with oligo and cDNA primers, SNP arrays, protein arrays, etc.
- Input Expression Analysis in Microsoft Excel — see below
- Run SAM as a Microsoft Excel Add-Ins
- Adjust the Delta tuning parameter to get a significant # of genes along with an acceptable false discovery rate (FDR)) and Assess Sample Size by calculating the mean difference in expression in the SAM Plot Controller
- List Differentially Expressed Genes (Positively and Negatively Expressed Genes)
Running SAM
- SAM is available for download online at http://www-stat.stanford.edu/~tibs/SAM/ for academic and non-academic users after completion of a registration step.
- SAM is run as an Excel Add-In, and the SAM Plot Controller allows Customization of the False Discovery Rate and Delta, while the SAM Plot and SAM Output functionality generate a List of Significant Genes, Delta Table, and Assessment of Sample Sizes
- Permutationsare calculated based on the number of samples
- Block Permutations
the number of permutations is set by the user when imputing correct values for the data set to run SAM
Response formats
Types:[32]
- Quantitative — real-valued (such as heart rate)
- One class — tests whether the mean gene expression differs from zero
- Two class — two sets of measurements
- Unpaired — measurement units are different in the two groups; e.g. control and treatment groups with samples from different patients
- Paired — same experimental units are measured in the two groups; e.g. samples before and after treatment from the same patients
- Multiclass — more than two groups with each containing different experimental units; generalization of two class unpaired type
- Survival — data of a time until an event (for example death or relapse)
- Time course — each experimental units is measured at more than one time point; experimental units fall into a one or two class design
- Pattern discovery — no explicit response parameter is specified; the user specifies eigengene (principal component) of the expression data and treats it as a quantitative response
Algorithm
SAM calculates a test statistic for relative difference in gene expression based on permutation analysis of expression data and calculates a false discovery rate. The principal calculations of the program are illustrated below.[32][33][34]
The so constant is chosen to minimize the coefficient of variation of di. ri is equal to the expression levels (x) for gene i under y experimental conditions.
Fold changes (t) are specified to guarantee genes called significant change at least a pre-specified amount. This means that the absolute value of the average expression levels of a gene under each of two conditions must be greater than the fold change (t) to be called positive and less than the inverse of the fold change (t) to be called negative.
The SAM algorithm can be stated as:
- Order test statistics according to magnitude [33][34]
- For each permutation compute the ordered null (unaffected) scores [33][34]
- Plot the ordered test statistic against the expected null scores [33][34]
- Call each gene significant if the absolute value of the test statistic for that gene minus the mean test statistic for that gene is greater than a stated threshold [34]
- Estimate the false discovery rate based on expected versus observed values [33][34]
Output
- Significant gene sets
- Positive gene set — higher expression of most genes in the gene set correlates with higher values of the phenotype y
- Negative gene set — lower expression of most genes in the gene set correlates with higher values of the phenotype y
SAM features
- Data from Oligo or cDNA arrays, SNP array, protein arrays, etc. can be utilized in SAM[33][34]
- Correlates expression data to clinical parameters[35]
- Correlates expression data with time[32]
- Uses data permutation to estimates False Discovery Rate for multiple testing[33][34][35][38]
- Reports local false discovery rate (the FDR for genes having a similar di as that gene)[32] and miss rates [32][33]
- Can work with blocked design for when treatments are applied within different batches of arrays[32]
- Can adjust threshold determining number of gene called significant[32]
Error correction and quality control
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.[39] Results may improve by removing these arrays from the analysis entirely.
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,[40] Agilent (GeneSpring),[41] and Ocimum Bio Solutions (Genowiz).[42]
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.
See also
- Microarray databases
- Significance analysis of microarrays
- Transcriptomics
- Proteomics
References
- ^ PMID 16199517.
- ^ Dr. Leming Shi, National Center for Toxicological Research. "MicroArray Quality Control (MAQC) Project". U.S. Food and Drug Administration. Retrieved 2007-12-26.
- ^ "GenUs BioSystems - Services - Data Analysis". Retrieved 2008-01-02.
- ^ "Agilent | DNA Microarrays". Archived from the original on December 22, 2007. Retrieved 2008-01-02.
- ^ "LIMMA Library: Linear Models for Microarray Data". Retrieved 2008-01-01.
- PMID 25690415.
- ^ "Create intensity versus ratio scatter plot of microarray data - MATLAB mairplot". MathWorks. Retrieved 2023-11-24.
- PMID 12925520.
- PMID 12538238.
- PMID 21070630.
- PMID 17646307.
- PMID 16473874.
- ^ "Affycomp III: A Benchmark for Affymetrix GeneChip Expression Measures".
- ^ PMID 16964229.
- S2CID 8192240.
- ISBN 978-0-387-29362-2.
- S2CID 760277.
- ^ PMID 24564555.
- ^ PMID 19038021.
- ^ "Home". biostat.ucsf.edu.
- ^ "Ingenuity Systems". Retrieved 2007-12-31.
- ^ "Ariadne Genomics: Pathway Studio". Archived from the original on 2007-12-30. Retrieved 2007-12-31.
- ^ "FunRich: Functional Enrichment Analysis". Retrieved 2014-09-09.
- ^ ["Significance Analysis of Microarrays". Retrieved 2007-12-31.]
- ^ "Software - Broad". Retrieved 2007-12-31.
- ^ "BioCarta - Charting Pathways of Life". Retrieved 2007-12-31.
- PMID 23443684.
- ^ "DBI Web". Archived from the original on 2007-07-05. Retrieved 2007-12-31.
- ^ "SCOPE". Archived from the original on 2011-08-17. Retrieved 2007-12-31.
- ^ "RssGsc". Retrieved 2008-10-15.
- ^ "SAM: Significance Analysis of Microarrays". tibshirani.su.domains. Retrieved 2023-11-24.
- ^ a b c d e f g h i Chu, G., Narasimhan, B, Tibshirani, R, Tusher, V. "SAM "Significance Analysis of Microarrays" Users Guide and technical document." [1]
- ^ PMID 17317331.
- ^ a b c d e f g h i <Zhang, S. (2007). "A comprehensive evaluation of SAM, the SAM R-package and a simple modification to improve its performance." BMC Bioinformatics 8: 230.
- ^ PMID 11309499.
- PMID 17612399.
- PMID 16872483.
- PMID 15921534.
- PMID 16076888.
- ^ "J. Craig Venter Institute -- Software". Retrieved 2008-01-01.
- ^ "Agilent | GeneSpring GX". Retrieved 2008-01-02.
- ^ "Ocimum Biosolutions | Genowiz". Archived from the original on 2009-11-24. Retrieved 2009-04-02.
External links
- ArrayExplorer - Compare microarray side by side to find the one that best suits your research needs
- FARMS - Factor Analysis for Robust Microarray Summarization, an R package —software
- StatsArray - Online Microarray Analysis Services —software
- ArrayMining.net - web-application for online analysis of microarray data —software
- FunRich - Perform gene set enrichment analysis —software
- Comparative Transcriptomics Analysis in Reference Module in Life Sciences
- SAM download instructions
- GeneChip® Expression Analysis-Data Analysis Fundamentals (by Affymetrix)
- Duke data_analysis_fundamentals_manual