Significance analysis of microarrays

Significance analysis of microarrays

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 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.

SAM identifies statistically significant genes by carrying out gene specific t-tests and computes a statistic "dj" for each gene "j", which measures the strength of the relationship between gene expression and a response variable 1,7,8. This analysis uses non-parametric statistics, since the data may not follow a normal distribution. The response variable describes and groups the data based on experimental conditions. In this method, repeated permutations of the data are used to determine if the expression of any gene is significant related to the response. The use of permutation-based analysis accounts for correlations in genes and avoids parametric assumptions about the distribution of individual genes. This is an advantage over other techniques (for example ANOVA and Bonferroni), which assume equal variance and/or independence of genes6.

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.

*AM 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

*Permutations are calculated based on the number of samples

*Block Permutations
**Blocks are batches of microarrays; for example for eight samples split into two groups (control and affected) there are 4!=24 permutations for each block and the total number of permutations is (24)(24)= 196. A minimum of 1000 permutations are recommended; 1,2 ,3 the number of permutations is set by the user when imputing correct values for the data set to run SAM

Response formats1

Types
**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; eg. 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 (principle component) of the expression data and treats it as a quantitative response

AM calculations 1, 7, 8

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 principle calculations of the program are illustrated below.

The "s"o constant is chosen to minimize the coefficient of variation of "di". r"i" 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 7, 8
#For each permutation compute the ordered null (unaffected) scores 7, 8
#Plot the ordered test statistic against the expected null scores 7, 8
#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 8
#Estimate the false discovery rate based on expected versus observed values 7, 8

AM 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

AM features

*Data from Oligo or cDNA arrays, SNP array, protein arrays,etc can be utilized in SAM 7, 8
*Correlates expression data to clinical parameters6
*Correlates expression data with time1
*Uses data permutation to estimates False Discovery Rate for multiple testing 5, 6, 7, 8
*Reports local false discovery rate (the FDR for genes having a similar di as that gene)1 and miss rates 1, 7
*Can work with blocked design for when treatments are applied within different batches of arrays1
*Can adjust threshold determining number of gene called significant1

External links

* [http://www-stat.stanford.edu/~tibs/SAM/|Where to download SAM]

References

#Chu, G., Narasimhan, B, Tibshirani, R, Tusher, V. "SAM "Significance Analysis of Microarrays" Users Guide and technical document." [http://www-stat.stanford.edu/~tibs/SAM/sam.pdf]
#Dinu, I. P., JD; Mueller, T; Liu, Q; Adewale, AJ; Jhangri, GS; Einecke, G; Famulski, KS; Halloran, P; Yasui, Y. (2007). "Improving gene set analysis of microarray data by SAM-GS." BMC Bioinformatics 8: 242.
#Jeffery, I. H., DG; Culhane, AC. (2006). "Comparison and evaluation of methods for generating differentially expressed gene lists from microarray data." BMC Bioinformatics 7: 359.
#Kooperberg, C., S. Sipione, et al. (2002). "Evaluating test statistics to select interesting genes in microarray experiments." Hum. Mol. Genet. 11(19): 2223–2232.
#Larsson, O. W., C; Timmons, JA. (2005). "Considerations when using the significance analysis of microarrays (SAM) algorithm." BMC Bioinformatics 6: 129.
#Tusher, V. G., R. Tibshirani, et al. (2001). "Significance analysis of microarrays applied to the ionizing radiation response." Proceedings of the National Academy of Sciences 98(9): 5116–5121.
#Zang, S., R. Guo, et al. (2007). "Integration of statistical inference methods and a novel control measure to improve sensitivity and specificity of data analysis in expression profiling studies." Journal of Biomedical Informatics 40(5): 552–560
#Zhang, S. (2007). "A comprehensive evaluation of SAM, the SAM R-package and a simple modification to improve its performance." BMC Bioinformatics 8: 230.


Wikimedia Foundation. 2010.

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

Look at other dictionaries:

  • 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… …   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

  • 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

  • Genexpressionsanalyse — Heatmaps von Genexpressionsdaten zeigen wie experimentelle Bedingungen die Produktion (Expression) von mRNA für eine Gruppe von Genen beeinflussen. Grün zeigt reduzierte Expression an, während rot für eine verstärkte Expression steht.… …   Deutsch Wikipedia

  • List of statistics topics — Please add any Wikipedia articles related to statistics that are not already on this list.The Related changes link in the margin of this page (below search) leads to a list of the most recent changes to the articles listed below. To see the most… …   Wikipedia

  • Sam — may be:;Used as a name: * Male given names (U.S. Census Bureau statistics 1990): [U.S. Census Bureau [http://www.census.gov/genealogy/names/dist.male.first U.S. List of male names] , 1990] **Samuel (0.306%) **Sam (0.092%) **Sammy (0.025%)… …   Wikipedia

  • Gilbert Chu — (Chinese: ) is an American Professor of Biochemistry and Medicine at Stanford University.He received a B.A. in Physics from Princeton University in 1967, a Ph.D in Physics from M.I.T. in 1973. a M.D from Harvard Medical School in 1980. He joined… …   Wikipedia

  • Methylated DNA immunoprecipitation — (MeDIP or mDIP) is a large scale (chromosome or genome wide) technique that is used to enrich for methylated DNA sequences. It consists of isolating methylated DNA fragments via an antibody raised against 5 methylcytosine (5mC). This technique… …   Wikipedia

  • Multiple comparisons — In statistics, the multiple comparisons or multiple testing problem occurs when one considers a set of statistical inferences simultaneously.[1] Errors in inference, including confidence intervals that fail to include their corresponding… …   Wikipedia

  • Health and Disease — ▪ 2009 Introduction Food and Drug Safety.       In 2008 the contamination of infant formula and related dairy products with melamine in China led to widespread health problems in children, including urinary problems and possible renal tube… …   Universalium

Share the article and excerpts

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