Functional Characterization of the S. cerevisiae Genome by Gene Deletion and Parallel Analysis
10.1126/science.285.5429.901
Crossref journal-article
American Association for the Advancement of Science (AAAS)
Science (221)
Abstract

The functions of many open reading frames (ORFs) identified in genome-sequencing projects are unknown. New, whole-genome approaches are required to systematically determine their function. A total of 6925 Saccharomyces cerevisiae strains were constructed, by a high-throughput strategy, each with a precise deletion of one of 2026 ORFs (more than one-third of the ORFs in the genome). Of the deleted ORFs, 17 percent were essential for viability in rich medium. The phenotypes of more than 500 deletion strains were assayed in parallel. Of the deletion strains, 40 percent showed quantitative growth defects in either rich or minimal medium.

Bibliography

Winzeler, E. A., Shoemaker, D. D., Astromoff, A., Liang, H., Anderson, K., Andre, B., Bangham, R., Benito, R., Boeke, J. D., Bussey, H., Chu, A. M., Connelly, C., Davis, K., Dietrich, F., Dow, S. W., El Bakkoury, M., Foury, F., Friend, S. H., Gentalen, E., … Davis, R. W. (1999). Functional Characterization of the S. cerevisiae Genome by Gene Deletion and Parallel Analysis. Science, 285(5429), 901–906.

Authors 52
  1. Elizabeth A. Winzeler (first)
  2. Daniel D. Shoemaker (additional)
  3. Anna Astromoff (additional)
  4. Hong Liang (additional)
  5. Keith Anderson (additional)
  6. Bruno Andre (additional)
  7. Rhonda Bangham (additional)
  8. Rocio Benito (additional)
  9. Jef D. Boeke (additional)
  10. Howard Bussey (additional)
  11. Angela M. Chu (additional)
  12. Carla Connelly (additional)
  13. Karen Davis (additional)
  14. Fred Dietrich (additional)
  15. Sally Whelen Dow (additional)
  16. Mohamed El Bakkoury (additional)
  17. Françoise Foury (additional)
  18. Stephen H. Friend (additional)
  19. Erik Gentalen (additional)
  20. Guri Giaever (additional)
  21. Johannes H. Hegemann (additional)
  22. Ted Jones (additional)
  23. Michael Laub (additional)
  24. Hong Liao (additional)
  25. Nicole Liebundguth (additional)
  26. David J. Lockhart (additional)
  27. Anca Lucau-Danila (additional)
  28. Marc Lussier (additional)
  29. Nasiha M'Rabet (additional)
  30. Patrice Menard (additional)
  31. Michael Mittmann (additional)
  32. Chai Pai (additional)
  33. Corinne Rebischung (additional)
  34. Jose L. Revuelta (additional)
  35. Linda Riles (additional)
  36. Christopher J. Roberts (additional)
  37. Petra Ross-MacDonald (additional)
  38. Bart Scherens (additional)
  39. Michael Snyder (additional)
  40. Sharon Sookhai-Mahadeo (additional)
  41. Reginald K. Storms (additional)
  42. Steeve Véronneau (additional)
  43. Marleen Voet (additional)
  44. Guido Volckaert (additional)
  45. Teresa R. Ward (additional)
  46. Robert Wysocki (additional)
  47. Grace S. Yen (additional)
  48. Kexin Yu (additional)
  49. Katja Zimmermann (additional)
  50. Peter Philippsen (additional)
  51. Mark Johnston (additional)
  52. Ronald W. Davis (additional)
References 48 Referenced 3,365
  1. 10.1038/ng0497-339
  2. 10.1016/S0378-1119(97)00140-6
  3. 10.1073/pnas.92.14.6479
  4. 10.1101/gad.8.9.1087
  5. 10.1073/pnas.94.1.190
  6. 10.1126/science.274.5295.2069
  7. 10.1093/nar/21.14.3329
  8. 10.1002/yea.320101310
  9. 10.1016/0378-1119(95)00144-U
  10. 10.1002/(SICI)1097-0061(19970915)13:11<1065::AID-YEA159>3.0.CO;2-K
  11. 10.1038/ng1296-450
  12. 10.1126/science.7618105
  13. To construct deletion strains two long oligonucleotide primers are synthesized each containing (3′ to 5′) 18 or 19 bases of homology to the antibiotic resistance cassette KanMX4 (U1 D1) a unique 20-bp tag sequence an 18-bp tag priming site (U2 or D2) and 18 bases of sequence complementary to the region upstream or downstream of the yeast ORF being targeted (including the start codon or stop codon; see ). These 74-mers are used to amplify the heterologous KanMX4 module which contains a constitutive efficient promoter from a related yeast strain Ashbya gosspii fused to the kanamycin resistance gene nptI (5). Because oligonucleotide synthesis is 3′ to 5′ and the fraction of full-size molecules decreases with increasing length improved targeting is achieved by performing a second round of PCR using primers bearing 45 bases of homology to the region upstream and downstream of a particular ORF. Transformation with the PCR product results in replacement of the targeted gene upon selection for G418 resistance. The unique 20-mer tag sequences are covalently linked to the sequence that targets them to the yeast genome creating a permanent association and genetic linkage between a particular deletion strain and the tag sequence. The use of two tags increases confidence in the analyses and the redundancy is useful in case one of the tags carries a mutation or performs poorly in hybridization assays. To verify correct integration of the deletion cassette genomic DNA was prepared from the resistant strains and used as template in PCR reactions using two primers common to the KanMX4 module (KanB (5′-CTGCAGCGAGGAGCCGTAAT-3′) and KanC (5′-TGATTTTGATGACGAGCGTAAT-3′) and four ORF-specific primers (A B C and D). A and D are from regions 200 to 400 bases upstream or downstream of the start codon whereas B and C are from within the ORF (see ). For verification both the A-KanB and the D-KanC PCR reactions were required to give the correct size product when analyzed by gel electrophoresis. If one of either the A-KanB or D-KanC reactions failed to yield a product the identification of the correctly-sized AD product could suffice. In addition haploid deletion strains were tested for the disappearance of the wild-type AB and CD products. All ORFs encoding proteins greater than 100 amino acids in size were initially selected for deletion. The deletion cassettes were designed to remove the entire coding sequence for a given ORF but to leave the start and stop codon intact. Although ∼10% of ORFs in S. cerevisiae overlap one another the positions of the deletions were not adjusted nor was any attempt made to avoid essential genes genes in which a previous deletion had been constructed or genes with a well-defined function. Genes represented multiple times in the genome (telomeric ORF Ty-elements) were usually not deleted as their targeted disruption would pose a challenge due to the conservation of upstream and downstream regions. Some smaller nonannotated ORFs (NORFs) will be deleted in the future. Transcripts from many of these NORFs have been detected in SAGE analysis warranting their inclusion in the study [
  14. 10.1016/S0092-8674(00)81845-0
  15. ]. All oligonucleotide primers (5 nmol scale) were synthesized on an automated multiplex oligonucleotide synthesizer [
  16. 10.1073/pnas.92.17.7912
  17. ] in batches of 96 primers. Scripts were written to automate the selection of primers. Primer sequences and ORF locations were chosen from the Stanford Genome Database (http: //genome-www.stanford.edu/Saccharomyces/) at various times over a 2-year period. The KanMX4 cassette was PCR-amplified and the resulting PCR products were sent to participating laboratories where 1 μg of PCR product was used to transform yeast by means of a variation on the standard lithium acetate procedure (. edu/group/yeast_deletion_project/protocols.html) in a 96-well format. Electronic records accessible over the World-Wide Web were kept for every strain constructed. MAT a haploid strains were given record numbers of less than 10 000 MATα haploid strains were given record numbers between 10 000 and 20 000 the heterozygous diploid between 20 000 and 30 000 and the homozygous diploid greater than 30 000. Each record consists of primer sequence information the results of the different diploid tests that were performed and notes about the phenotype. Data for completed strains are accessible from www-sequence.stanford.edu/group/deletion/index. html. These strains frozen in 15% glycerol can be obtained from Research Genetics (Huntsville AL) or EUROSCARF (Frankfurt Germany).
  18. Four different types of strains containing several genetic markers (Table 1) were generated for each ORF—two haploid strains one of each mating type and two diploid strains heterozygous and homozygous for the deletion loci and kanamycin marker. The homozygous diploid was constructed by mating the two haploid strains obtained from independent transformations. Typically the heterozygous diploid and one of the haploid strains were obtained by direct transformation while the other haploid strain and the homozygous diploid were obtained by sporulation and mating respectively. Essential genes were identified by 2:2 segregation of viability in tetrads derived from the heterozygous diploid. Strains were sporulated by patching them on a fresh GNA plate (5% D -glucose 3% Difco nutrient broth 1% Difco yeast extract 2% Difco Bacto agar) for 1 day at 30°C before transfer to liquid sporulation medium (1% potassium acetate and 0.005% zinc acetate supplemented with 0.1 mM uracil 0.15 mM histidine-HCL or 1.0 mM leucine as necessary). Sporulation cultures were incubated on a rollerwheel for 4 to 5 days at 25°C. If 2:2 segregation of viability was consistently observed in two independently-transformed heterozygous deletion strains the gene was designated essential. All confirmed diploid strains obtained through mating were required to pass two of three tests: a diploid budding pattern the ability to sporulate and the inability to mate. In a few cases essential genes overlapped other essential genes (20 pairs) or a gene whose viability status was unknown (four pairs) making it difficult to determine the cause of lethality.
  19. 10.1093/nar/25.1.28
  20. 10.1016/S1097-2765(00)80114-8
  21. To construct the pools each deletion strain was patched on YPD plates in the presence of 150 mg/liter G418. Approximately equal numbers of cells were harvested from the plate for each strain and combined. Aliquots of the pools were stored in the presence of 15% glycerol at–80°C.
  22. The UPTAG and DOWNTAG sequences were separately amplified from genomic DNA by means of PCR using primers B-U2-comp (5′ biotin-GTCGACCTGCAGCGTACG-3′) and U1 (5′-GATGTCCACGAGGTCTCT-3′) or primers B-D2-comp (5′-biotin-CGAGCTCGAATTCATCG-3′) and D1 (5′-CGGTGTCGGTCTCGTAG-3′). In both cases a twofold molar excess of biotinylated primers was used in the reactions. The amplified UPTAG and DOWNTAG sequences were combined and hybridized to high-density oligonucleotide arrays in 200 μl of 6×SSPE [1 M NaCl 66 mM NaH 2 PO –4 6.6 mM EDTA (pH 7.4)] containing 0.005% Triton X-100 (SSPE-T) 200 pmol U1 200 pmol U2 (5′-CGTACGCTGCAGGTCGAC-3′) 200 pmol D1 and 200 pmol primer D2 (5′-CGATGAATTCGAGCTCG-3′). The addition of complementary primers to the hybridization mix was shown to improve the signal-to-noise ratio [D. D. Shoemaker thesis Stanford University Stanford CA (1998)]. Samples were heated to 100°C for two min and then cooled on ice before being applied to the array. Samples were hybridized for 1 hour at 42°C. The arrays were washed two times with six changes of 6×SSPE-T. The arrays were then stained at 42°C for 10 min with 6×SSPE-T containing 2 μg/ml phycoerythrin-streptavidin (Molecular Probes) and 1 μg/ml acetylated bovine serum albumin washed two times with five changes 6×SSPE-T and scanned at an emission wavelength of 560 nm using an Affymetrix GeneChip Scanner. Of the strains analyzed 157 contained only a single tag sequence (UPTAG). Six of these strains were not detected in the hybridization mix. Of the strains represented with two-tag sequences 98.5% exhibited either an UPTAG or DOWNTAG bar code hybridization signal that was greater than threefold over background. Sequencing of 186 deletion regions tags showed that 25% of mutations in the tags or tag priming sites resulted in a nonfunctional tag that could not be amplified or detected by hybridization or both. In only 1.1% of cases was a complete lack of hybridization signal not associated with a mutation in the tag or tag priming site. Mutations were most often found in the tags or tag priming sites (0.85% per base) and were less frequent in the regions of yeast homology (0.25% per base) most likely due to selection against the mutated PCR products during the recombination event or to the two-step PCR strategy.
  23. Grids were aligned to the scanned images using the known feature dimensions of the array. The hybridization intensities for each of the elements in the grid were determined using the 75th percentile method in the Affymetrix GeneChip software package. Subsequent analysis of the hybridization intensities consisted of two steps: adjustment of data to achieve approximate equality of background and maximal signals on each array and analysis of the decrease (or increase) of the UPTAG and DOWNTAG signal strength over time. Equalization of signal strength relied on the fact that for most array sites the amount of tag DNA present did not vary over time. A consensus score for these sites was obtained from the first principal component of the logarithms of the signals. The logarithms of the hybridization intensity for each element on the array were linearly transformed to make each array's overall signal approximately equal to the consensus. The growth rate for each strain was determined by using the model log 2 (signal) = max( a + b t 0) + e where a and b are model parameters t represents the number of population doublings and e is a random error term. The growth rate is calculated as 1 + b. For strains that have dropped out of the pool ( a + bt < 0) the model describes the statistical distribution of background signals. When e is normally distributed background signals have a lognormal distribution. This appeared to be roughly true in our data. However we found that a small fraction of signals on each array are liable to have very high values much larger than can be accounted for by a purely lognormal background model. To obtain a degree of robustness against these occasional outliers we assumed that the error term e had a scaled t distribution with one degree of freedom. The use of this heavy-tailed distribution reduced the likelihood of false positive identifications of deficient strains due to occasional high signal levels at t = 0 at the expense of a possible reduction in the ability to detect marginally deficient strains.
  24. The analysis was not as accurate for strains with growth rates of less than 0.5 that of the wild type because generally only three data points were above background for these.
  25. 10.1074/jbc.273.6.3253
  26. A comprehensive study of chromosome V genes using genetic footprinting (3) provided an opportunity to validate the data: the results generally agreed with a few exceptions. For example of the 52 genes whose disruption had no effect on strain fitness under all conditions tested by genomic footprinting we detected a growth defect in deletants yel033w (0.68 R; 0.83 M) yel050c (0.73 R; 0.68 M) and yer028c (0.79 R; 0.69 M). The observed phenotype for yel033w probably results from interference with a neighboring gene ( HYP2 encoding translation initiation factor eIF-5A). In addition the hem14 deletant showed a strong growth defect in rich medium while genomic footprinting revealed a salt-specific defect but no defect in rich medium. Of the 11 genes that had been shown by genomic footprinting to have a severe growth defect in three cases our deletants appeared to have no discernible phenotype ( nrf1 gda1 pcl6 ). These differences could be due to our use of diploid versus haploid strains to the auxotrophies carried by our deletants or to our using 30°C versus 25°C as the growth temperature. The effect of temperature is probably the cause of discrepancies for the nrf1 strain which grows slowly at 25°C but grows faster than normal at 36.5°C. Among other disparities two ( YER082C and YEL026W ) of the 22 genes on chromosome V determined to be essential by tetrad analysis were wild-type by footprinting (perhaps because of cross-feeding in the Ty pool or because the products are required for germination) and the minimal medium–specific growth defects we detected for hom3 and ilv1 mutants were not detected by genomic footprinting.
  27. Cultures of BY4743 ho::KanMX4 were grown at 30°C in YPD or minimal media supplemented with histidine uracil lysine methionine and adenine in the presence of G418 and were harvested at mid-log phase. cDNA was prepared from 20 μg of polyadenylated RNA from each sample using a dT 21 primer and Superscript II reverse transcriptase (GibcoBRL) according to the manufacturer's recommendation. cDNA was fragmented using DNaseI (GibcoBRL) biotinylated using ddATP (NEN) and Terminal Transferase (Boehringer) and hybridized to yeast full-genome arrays (Affymetrix) as described in L. Wodicka et al. [ Nature Biotechnol. 15 1359 (1997)]. After scanning the average signal from each array was normalized to the average signal strength of all eight chips.
  28. 10.1002/(SICI)1097-0061(199811)14:15<1417::AID-YEA334>3.0.CO;2-N
  29. Contributing groups include all authors G. Valle S. Kelley J. Strathern and D. Garfinkel.
  30. EUROFAN projects B0 and B9 at www.mips.biochem.mpg.de/proj/eurofan/index.html; R. Niedenthal et al. Yeast in press.
  31. 10.1038/6791
  32. 10.1038/42711
  33. 10.1093/emboj/17.8.2235
  34. 10.1016/0014-5793(92)81325-G
  35. 10.1074/jbc.270.1.437
  36. 10.1242/jcs.101.3.709
  37. 10.1002/(SICI)1097-0061(19970315)13:3<281::AID-YEA74>3.0.CO;2-E
  38. 10.1002/j.1460-2075.1994.tb06759.x
  39. 10.1074/jbc.272.27.17045
  40. 10.1083/jcb.133.1.85
  41. 10.1074/jbc.271.44.27829
  42. 10.1139/g97-022
  43. 10.1101/gad.4.5.740
  44. 10.1128/jb.176.5.1488-1499.1994
  45. 10.1126/science.1990438
  46. 10.1073/pnas.91.11.5022
  47. 10.1002/(SICI)1097-0061(19980130)14:2<115::AID-YEA204>3.0.CO;2-2
  48. We thank D. Lashkari for establishing the oligonucleotide synthesis facility T. Nguyen M. Sigrist and K. Tanner for help in tetrad analysis S. Voegeli for DNA sequence analyses P. Koetter for distribution of deletion strains J. Rine for helpful advice and M.Cherry and K. Wolfe for files. E.A.W. is supported by the John Wasmuth fellowship in Genomic Analysis (HG00185-02). Supported by NIH grants HG01633 HG01627 HG00198 by an operating grant from the Medical Research Council of Canada by grants from the European Commission (BIO4-CT97-2294) by the Swiss Federal Office for Education and Science and by the region de Bruxelles-Capital Belgium.
Dates
Type When
Created 23 years ago (July 27, 2002, 5:37 a.m.)
Deposited 1 year, 7 months ago (Jan. 13, 2024, 3:55 a.m.)
Indexed 1 day, 3 hours ago (Aug. 22, 2025, 12:45 a.m.)
Issued 26 years ago (Aug. 6, 1999)
Published 26 years ago (Aug. 6, 1999)
Published Print 26 years ago (Aug. 6, 1999)
Funders 0

None

@article{Winzeler_1999, title={Functional Characterization of the S. cerevisiae Genome by Gene Deletion and Parallel Analysis}, volume={285}, ISSN={1095-9203}, url={http://dx.doi.org/10.1126/science.285.5429.901}, DOI={10.1126/science.285.5429.901}, number={5429}, journal={Science}, publisher={American Association for the Advancement of Science (AAAS)}, author={Winzeler, Elizabeth A. and Shoemaker, Daniel D. and Astromoff, Anna and Liang, Hong and Anderson, Keith and Andre, Bruno and Bangham, Rhonda and Benito, Rocio and Boeke, Jef D. and Bussey, Howard and Chu, Angela M. and Connelly, Carla and Davis, Karen and Dietrich, Fred and Dow, Sally Whelen and El Bakkoury, Mohamed and Foury, Françoise and Friend, Stephen H. and Gentalen, Erik and Giaever, Guri and Hegemann, Johannes H. and Jones, Ted and Laub, Michael and Liao, Hong and Liebundguth, Nicole and Lockhart, David J. and Lucau-Danila, Anca and Lussier, Marc and M’Rabet, Nasiha and Menard, Patrice and Mittmann, Michael and Pai, Chai and Rebischung, Corinne and Revuelta, Jose L. and Riles, Linda and Roberts, Christopher J. and Ross-MacDonald, Petra and Scherens, Bart and Snyder, Michael and Sookhai-Mahadeo, Sharon and Storms, Reginald K. and Véronneau, Steeve and Voet, Marleen and Volckaert, Guido and Ward, Teresa R. and Wysocki, Robert and Yen, Grace S. and Yu, Kexin and Zimmermann, Katja and Philippsen, Peter and Johnston, Mark and Davis, Ronald W.}, year={1999}, month=aug, pages={901–906} }