a microarray-based screening procedure for detecting differentially represented yeast mutants rafael...
TRANSCRIPT
![Page 1: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/1.jpg)
A Microarray-Based Screening Procedure for Detecting Differentially
Represented Yeast Mutants
Rafael A. Irizarry
Department of Biostatistics, JHU
http://biostat.jhsph.edu/~ririzarr
![Page 2: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/2.jpg)
![Page 3: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/3.jpg)
kanRA
Transformation into deletion pool
Select for Ura+ transformantsGenomic DNA preparation
Circular pRS416
PCRCy5 labeled PCR products Cy3 labeled PCR products
Oligonucleotide array hybridization
B
EcoRI linearized PRS416
NHEJ Defective
MCS
CEN/ARS
URA3 ttaaaatt
CEN/ARS
URA3
UPTAG DOWNTAG
![Page 4: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/4.jpg)
Which mutants are NHEJ defective?
• Find mutants defective for transformation with linear DNA
• Dead in linear transformation (green)
• Alive in circular transformation (red)
• Look for spots with large log(R/G)
![Page 5: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/5.jpg)
• .
Y K U 7 0 N E J 1 Y K U 8 0
Y K U 7 0 N E J 1 Y K U 8 0
![Page 6: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/6.jpg)
Data
• 5718 mutants
• 3 replicates on each slide
• 5 Haploid slides, 4 Diploid slides
• Arrays are divided into 2 downtags, 3 uptag (2 of which replicate uptags)
![Page 7: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/7.jpg)
Average Red and Green Scatter Plot
![Page 8: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/8.jpg)
Average Red and Green MVA plot
![Page 9: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/9.jpg)
![Page 10: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/10.jpg)
Improvement to usual approach
• Take into account that some mutants are dead and some alive
• Use a statistical model to represent this• Mixture model?• With ratio’s we lose information about R
and G separately • Look at them separately (absolute
analysis)
![Page 11: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/11.jpg)
Histograms
![Page 12: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/12.jpg)
Using model we can attach uncertainty to tests
For example posterior z-test,
weighted average of z-tests with weights obtained using the posterior probability (obtained from EM)
Is Normal(0,1)
![Page 13: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/13.jpg)
QQ-Plot
![Page 14: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/14.jpg)
Uptag/Downtag Z-Scores
![Page 15: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/15.jpg)
Average Red and Green MVA Plot
![Page 16: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/16.jpg)
Average Red and Green Scatter Plot
![Page 17: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/17.jpg)
![Page 18: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/18.jpg)
![Page 19: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/19.jpg)
![Page 20: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/20.jpg)
![Page 21: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/21.jpg)
![Page 22: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/22.jpg)
ResultsTable
1 YMR106C 9.5 47 69.2 a a 1002 YOR005C 19.7 35 44.9 a d 1003 YLR265C 6.1 32 35.8 a m 1004 YDL041W 10.4 32 35.6 a m 1005 YIL012W 12.2 31 21.7 a a 1006 YIL093C 4.8 29 30.8 a a 1007 YIL009W 5.6 29 -23.5 a a 1008 YDL042C 12.9 29 32.1 a d 1009 YIL154C 1.8 28 91.3 m m 8210 YNL149C 1.7 27 93.4 m d 7111 YBR085W 2.5 26 -15.8 a a 8412 YBR234C 1.7 26 87.5 m d 7513 YLR442C 6.1 26 -100.0 a a 100
![Page 23: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/23.jpg)
Acknowledgements
• Siew Loon Ooi
• Jef Boeke
• Forrest Spencer
• Jean Yang
![Page 24: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/24.jpg)
END
![Page 25: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/25.jpg)
Summary
• Simple data exploration useful tool for quality assessment
• Statistical thinking helpful for interpretation
• Statistical models may help find signals in noise
![Page 26: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/26.jpg)
Acknowledgements
UC Berkeley StatBen BolstadSandrine DudoitTerry SpeedJean Yang
MBG (SOM)Jef BoekeSiew-Loon OoiMarina LeeForrest Spencer
BiostatisticsKarl BromanLeslie CopeCarlo CoulantoniGiovanni ParmigianiScott Zeger
Gene LogicFrancois Colin Uwe Scherf’s Group
PGATom Cappola Skip GarciaJoshua Hare
WEHIBridget HobbsNatalie Thorne
![Page 27: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/27.jpg)
![Page 28: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/28.jpg)
![Page 29: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/29.jpg)
![Page 30: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/30.jpg)
![Page 31: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/31.jpg)
Warning
• Absolute analyses can be dangerous for competitive hybridization slides
• We must be careful about “spot effect”
• Big R or G may only mean the spot they where on had large amounts of cDNA
• Look at some facts that make us feel safer
![Page 32: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/32.jpg)
Correlation between replicates
R1 R2 R3 G1 G2 G3
R1 1.00 0.95 0.95 0.94 0.90 0.90
R2 0.95 1.00 0.96 0.90 0.95 0.91
R3 0.95 0.96 1.00 0.91 0.92 0.95
G1 0.94 0.90 0.91 1.00 0.96 0.96
G2 0.90 0.95 0.92 0.96 1.00 0.97
G3 0.90 0.91 0.95 0.96 0.97 1.00
![Page 33: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/33.jpg)
Correlation between red, green, haploid, diplod, uptag, downtag
RHD RHU RDD RDU GHD GHU GDD GDU
RHD 1.00 0.59 0.56 0.32 0.95 0.58 0.54 0.37RHU 0.59 1.00 0.38 0.56 0.58 0.95 0.40 0.58RDD 0.56 0.38 1.00 0.58 0.54 0.39 0.92 0.64RDU 0.32 0.56 0.58 1.00 0.33 0.53 0.58 0.89GHD 0.95 0.58 0.54 0.33 1.00 0.62 0.56 0.39GHU 0.58 0.95 0.39 0.53 0.62 1.00 0.41 0.58GDD 0.54 0.40 0.92 0.58 0.56 0.41 1.00 0.73GDU 0.37 0.58 0.64 0.89 0.39 0.58 0.73 1.00
![Page 34: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/34.jpg)
BTW
The mean squared error across slides is about 3 times bigger than the mean squared error within slides
![Page 35: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/35.jpg)
Mixture Model
We use a mixture model that assumes:
• There are three classes:– Dead– Marginal– Alive
• Normally distributed with same correlation structure from gene to gene
![Page 36: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/36.jpg)
Random effect justification
Each x = (r1,…,r5,g1,…,g5) will have the following effects:
• Individual effect: same mutant same expression (replicates are alike)
• Genetic effect: same genetics same expression
• PCR effect : expect difference in uptag, downtag
![Page 37: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/37.jpg)
Does it fit?
![Page 38: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/38.jpg)
Does it fit?
![Page 39: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/39.jpg)
What can we do now that we couldn’t do before?
• Define a t-test that takes into account if mutants are dead or not when computing variance
• For each gene compute likelihood ratios comparing two hypothesis:
alive/dead vs.dead/dead or alive/alive
![Page 40: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/40.jpg)
QQ-plot for new t-test
![Page 41: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/41.jpg)
Better looking than others
![Page 42: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/42.jpg)
![Page 43: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/43.jpg)
![Page 44: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/44.jpg)
![Page 45: A Microarray-Based Screening Procedure for Detecting Differentially Represented Yeast Mutants Rafael A. Irizarry Department of Biostatistics, JHU rafa@jhu.edu](https://reader036.vdocument.in/reader036/viewer/2022062409/56649f3d5503460f94c5cfc3/html5/thumbnails/45.jpg)
1 YMR106C 9.5 47 69.2 a a 1002 YOR005C 19.7 35 44.9 a d 1003 YLR265C 6.1 32 35.8 a m 1004 YDL041W 10.4 32 35.6 a m 1005 YIL012W 12.2 31 21.7 a a 1006 YIL093C 4.8 29 30.8 a a 1007 YIL009W 5.6 29 -23.5 a a 1008 YDL042C 12.9 29 32.1 a d 1009 YIL154C 1.8 28 91.3 m m 8210 YNL149C 1.7 27 93.4 m d 7111 YBR085W 2.5 26 -15.8 a a 8412 YBR234C 1.7 26 87.5 m d 7513 YLR442C 6.1 26 -100.0 a a 100