Microarray analysis exercises 1

Starting Data     Processed Data

You'll be using Excel to do most of the mathematical analyses, since this will show the exact formulas used to perform every step of the analysis pipeline. As a result, you'll need to use Excel functions and be familiar with some Excel conventions. See the Excel help for the details.

Preliminary information: Image analysis and calculation of expression value

1. As described in Su et al., 2002, human tissue samples were hybridized on Affymetrix (one-color) arrays and chips were scanned. For each tissue, at least two independendent samples were hybridized to separate chips.
2. Scanned images were quantified (including measurement of background) using standard software.
3. Data from a probeset (a series of oligos designed to a specific gene target) were used to calculate an expression values for that probeset using standard Affymetrix algorithms.
4. See the manuals from Affymetrix for more information about these processes, and the Statistical Algorithms Description Document for the actual equations used.
5. Note that these analysis protocols are generally specific to the chip type and its manufacturer.

1. Why normalize? Chips may have been hybridized to different amounts of RNA, for different amounts of time, with different batches of solutions, etc. Normalization should remove systematic biases and make any comparisons between chips more meaningful.
2. Download the starting data for the exercises.
3. Look at the expression data for a selected series of experiments
• Down the first column are the Affymetrix probeset IDs, each corresponding to a target gene.
• Each other column shows a tissue name, followed by expression values in some arbitrary unit.
• There should be two chips each for adult and fetal brain and liver (a total of 8 chips).
• Note that since this is a one-color array, the expression measures are absolute values and not ratios. If you're more used to 2-color (cDNA) arrays, this would be similar to (but not exactly the same as) 4 chips, each with one channel/dye for fetal tissue and another for adult tissue.
4. Calculate the median of all expression values from each chip.
• Use the Excel "MEDIAN" function to calculate the median at the bottom of each column. Ex: =MEDIAN(B2:B12627)
5. To perform global median normalization, on a new sheet ("norm") of the same file, scale each column of data so that each median is 100.
• Divide each expression signal in a chip by the median for the chip and multiply by 100. Ex: =(raw!B2/raw!B$12628)*100 [if B12628 contains the median for the chip]
• To confirm that this was correctly done, calculate the median of each normalized chip.
• Global mean normalization is also possible but more influenced by outliers. If you want to compare mean and median yourself [after the class], try both methods and then compare the results with scatterplots.