Quantitative phenotypic analysis of yeast deletion mutants using a highly parallel molecular bar-coding strategy.

A quantitative and highly parallel method for analysing deletion mutants has been developed to aid in determining the biological function of thousands of newly identified open reading frames (ORFs) in Saccharomyces cerevisiae. This approach uses a PCR targeting strategy to generate large numbers of deletion strains. Each deletion strain is labelled with a unique 20-base tag sequence that can be detected by hybridization to a high-density oligonucleotide array. The tags serve as unique identifiers (molecular bar codes) that allow analysis of large numbers of deletion strains simultaneously through selective growth conditions. Hybridization experiments show that the arrays are specific, sensitive and quantitative. A pilot study with 11 known yeast genes suggests that the method can be extended to include all of the ORFs in the yeast genome, allowing whole genome analysis with a single selective growth condition and a single hybridization.

Systematic variation in gene expression patterns in human cancer cell lines.

We used cDNA microarrays to explore the variation in expression of approximately 8,000 unique genes among the 60 cell lines used in the National Cancer Institute's screen for anti-cancer drugs. Classification of the cell lines based solely on the observed patterns of gene expression revealed a correspondence to the ostensible origins of the tumours from which the cell lines were derived. The consistent relationship between the gene expression patterns and the tissue of origin allowed us to recognize outliers whose previous classification appeared incorrect. Specific features of the gene expression patterns appeared to be related to physiological properties of the cell lines, such as their doubling time in culture, drug metabolism or the interferon response. Comparison of gene expression patterns in the cell lines to those observed in normal breast tissue or in breast tumour specimens revealed features of the expression patterns in the tumours that had recognizable counterparts in specific cell lines, reflecting the tumour, stromal and inflammatory components of the tumour tissue. These results provided a novel molecular characterization of this important group of human cell lines and their relationships to tumours in vivo.

Genome-wide mapping with biallelic markers in Arabidopsis thaliana.

Single-nucleotide polymorphisms, as well as small insertions and deletions (here referred to collectively as simple nucleotide polymorphisms, or SNPs), comprise the largest set of sequence variants in most organisms. Positional cloning based on SNPs may accelerate the identification of human disease traits and a range of biologically informative mutations. The recent application of high-density oligonucleotide arrays to allele identification has made it feasible to genotype thousands of biallelic SNPs in a single experiment. It has yet to be established, however, whether SNP detection using oligonucleotide arrays can be used to accelerate the mapping of traits in diploid genomes. The cruciferous weed Arabidopsis thaliana is an attractive model system for the construction and use of biallelic SNP maps. Although important biological processes ranging from fertilization and cell fate determination to disease resistance have been modelled in A. thaliana, identifying mutations in this organism has been impeded by the lack of a high-density genetic map consisting of easily genotyped DNA markers. We report here the construction of a biallelic genetic map in A. thaliana with a resolution of 3.5 cM and its use in mapping Eds16, a gene involved in the defence response to the fungal pathogen Erysiphe orontii. Mapping of this trait involved the high-throughput generation of meiotic maps of F2 individuals using high-density oligonucleotide probe array-based genotyping. We developed a software package called InterMap and used it to automatically delimit Eds16 to a 7-cM interval on chromosome 1. These results are the first demonstration of biallelic mapping in diploid genomes and establish means for generalizing SNP-based maps to virtually any genetic organism.

A radiation hybrid map of the zebrafish genome.

Recent large-scale mutagenesis screens have made the zebrafish the first vertebrate organism to allow a forward genetic approach to the discovery of developmental control genes. Mutations can be cloned positionally, or placed on a simple sequence length polymorphism (SSLP) map to match them with mapped candidate genes and expressed sequence tags (ESTs). To facilitate the mapping of candidate genes and to increase the density of markers available for positional cloning, we have created a radiation hybrid (RH) map of the zebrafish genome. This technique is based on somatic cell hybrid lines produced by fusion of lethally irradiated cells of the species of interest with a rodent cell line. Random fragments of the donor chromosomes are integrated into recipient chromosomes or retained as separate minichromosomes. The radiation-induced breakpoints can be used for mapping in a manner analogous to genetic mapping, but at higher resolution and without a need for polymorphism. Genome-wide maps exist for the human, based on three RH panels of different resolutions, as well as for the dog, rat and mouse. For our map of the zebrafish genome, we used an existing RH panel and 1,451 sequence tagged site (STS) markers, including SSLPs, cloned candidate genes and ESTs. Of these, 1,275 (87.9%) have significant linkage to at least one other marker. The fraction of ESTs with significant linkage, which can be used as an estimate of map coverage, is 81.9%. We found the average marker retention frequency to be 18.4%. One cR3000 is equivalent to 61 kb, resulting in a potential resolution of approximately 350 kb.

Functional analysis of the genes of yeast chromosome V by genetic footprinting.

Genetic footprinting was used to assess the phenotypic effects of Ty1 transposon insertions in 268 predicted genes of chromosome V of Saccharomyces cerevisiae. When seven selection protocols were used, Ty1 insertions in more than half the genes tested (157 of 268) were found to result in a detectable reduction in fitness. Results could not be obtained for fewer than 3 percent of the genes tested (7 of 268). Previously known mutant phenotypes were confirmed, and, for about 30 percent of the genes, new mutant phenotypes were identified.

The transcriptional program in the response of human fibroblasts to serum.

The temporal program of gene expression during a model physiological response of human cells, the response of fibroblasts to serum, was explored with a complementary DNA microarray representing about 8600 different human genes. Genes could be clustered into groups on the basis of their temporal patterns of expression in this program. Many features of the transcriptional program appeared to be related to the physiology of wound repair, suggesting that fibroblasts play a larger and richer role in this complex multicellular response than had previously been appreciated.

Whole genome analysis: experimental access to all genome sequenced segments through larger-scale efficient oligonucleotide synthesis and PCR.

The recent ability to sequence whole genomes allows ready access to all genetic material. The approaches outlined here allow automated analysis of sequence for the synthesis of optimal primers in an automated multiplex oligonucleotide synthesizer (AMOS). The efficiency is such that all ORFs for an organism can be amplified by PCR. The resulting amplicons can be used directly in the construction of DNA arrays or can be cloned for a large variety of functional analyses. These tools allow a replacement of single-gene analysis with a highly efficient whole-genome analysis.

Allelic diversification at the class II DQB locus of the mammalian major histocompatibility complex.

The allelic diversity at HLA class II loci either arose after the divergence of hominoid lineages or, alternatively, the polymorphism was present before speciation and has been maintained by selection. Here, we report the use of oligonucleotide primers to amplify, by the polymerase chain reaction, and sequence the polymorphic second exon of the DQB locus from 11 species, spanning more than 40 million years of mammalian evolution. Phylogenetic analysis reveals that of the four human DQB allelic types (DQB1-B4), three (DQB1-3) were found in chimpanzee and gorilla and two (DQB3 and -4) were identified in the rhesus monkey, suggesting that some of these types are 5-20 million years old. The ratio of replacement to silent substitutions was calculated between members of the same allelic type from different species. These results suggest that the evolution of the DQB3 allelic type is more constrained than that of the DQB1 allelic type; both evolve more slowly than the DXB locus, a linked but presumably nonexpressed locus. Further, the clustering of allelic subtypes by species in the phylogenetic tree indicates that allelic diversification has occurred subsequent to the divergence of hominoids. Finally, some haplotype combinations of DQA and DQB alleles are common to several hominoid species and may have been maintained for at least 5 million years.

An automated multiplex oligonucleotide synthesizer: development of high-throughput, low-cost DNA synthesis.

An automated oligonucleotide synthesizer has been developed that can simultaneously and rapidly synthesize up to 96 different oligonucleotides in a 96-well microtiter format using phosphoramidite synthesis chemistry. A modified 96-well plate is positioned under reagent valve banks, and appropriate reagents are delivered into individual wells containing the growing oligonucleotide chain, which is bound to a solid support. Each well has a filter bottom that enables the removal of spent reagents while retaining the solid support matrix. A seal design is employed to control synthesis environment and the entire instrument is automated via computer control. Synthesis cycle times for 96 couplings are < 11 min, allowing a plate of 96 20-mers to be synthesized in < 5 hr. Oligonucleotide synthesis quality is comparable to commercial machines, with average coupling efficiencies routinely > 98% across the entire 96-well plate. No significant well-to-well variations in synthesis quality have been observed in > 6000 oligonucleotides synthesized to date. The reduced reagent usage and increased capacity allow the overall synthesis cost to drop by at least a factor of 10. With the development of this instrument, it is now practical and cost-effective to synthesize thousands to tens of thousands of oligonucleotides.

Mapping expressed sequence tag sites on yeast artificial chromosome clones of Arabidopsis thaliana DNA.

We describe a method for efficient parallel mapping of expressed sequence tag (EST) sites onto yeast artificial chromosome (YAC) clones. The strategy involves an initial YAC clone pooling scheme that minimizes the number of required PCR amplifications. This is followed by parallel analysis of PCR amplicons of EST sequences. Using this method, we have screened 600 EST sites in combinatorial pools of 3449 YAC clones that contain Arabidopsis thaliana DNA inserts. The presence of these genes on YACs was detected by amplifying EST sequences with PCR and analyzing the reaction products by agarose gel electrophoresis. Of the 600 ESTs, 271 were found to map to individual YACs. Software tools are presented that allow for the automated analysis of this electrophoresis data. Suggestions for the scale-up of this method to map large genomes are discussed.

Use of a fluorescent-PCR reaction to detect genomic sequence copy number and transcriptional abundance.

We present a fluorescent-PCR-based technique to assay genomic sequence copy number and transcriptional abundance. This technique relies on the ability to follow fluorescent PCR progressively in real time during the exponential phase of the reaction so that quantitative PCR is accomplished. We demonstrated the ability of this technique to quantitate both known deletions and amplifications of loci that have been measured previously by other methods, and to measure transcriptional abundance. Using an efficient variant of the fluorescent-PCR technology, we can monitor transcription semiquantitatively. The ability to detect all amplifications and deletions at any single copy locus by PCR makes this the technique of choice to assay genomic sequence copy number anomalies in birth defects and cancers. The ability to detect variations in transcript abundance enables this technique to fashion a time and tissue analysis of transcription.

Yeast microarrays for genome wide parallel genetic and gene expression analysis.

We have developed high-density DNA microarrays of yeast ORFs. These microarrays can monitor hybridization to ORFs for applications such as quantitative differential gene expression analysis and screening for sequence polymorphisms. Automated scripts retrieved sequence information from public databases to locate predicted ORFs and select appropriate primers for amplification. The primers were used to amplify yeast ORFs in 96-well plates, and the resulting products were arrayed using an automated micro arraying device. Arrays containing up to 2,479 yeast ORFs were printed on a single slide. The hybridization of fluorescently labeled samples to the array were detected and quantitated with a laser confocal scanning microscope. Applications of the microarrays are shown for genetic and gene expression analysis at the whole genome level.

High-throughput genomic and proteomic analysis using microarray technology.

BACKGROUND: High-density microarrays are ideally suited for analyzing thousands of genes against a small number of samples. The next step in the discovery process is to take the resulting genes of interest and rapidly screen them against thousands of patient samples, tissues, or cell lines to further investigate their involvement in disease risk or the response to medication. METHODS: We used a microarray technology platform for both single-nucleotide polymorphisms (SNPs) and protein expression. Each microarray contains up to 250 elements that can be customized for each application. Slides contained either a 16- or 96-microarray format (4000-24,000 elements per slide), allowing the corresponding number of samples to be rapidly processed in parallel. RESULTS: Results for SNP genotyping and protein profiling agreed with results of restriction fragment length polymorphism (RFLP) analysis or ELISA, respectively. Genotyping analyses, using the microarray technology, on large sample sets over multiple polymorphisms in the NAT2 gene were in full agreement with traditional methodologies, such as sequencing and RFLP analysis. The multiplexed protein microarray had correlation coefficients of 0.82-0.99 (depending on analyte) compared with ELISAs. CONCLUSIONS: The integrated microarray technology platform is adaptable and versatile, while offering the high-throughput capabilities needed for drug development and discovery applications.

A genetic linkage map for zebrafish: comparative analysis and localization of genes and expressed sequences.

Genetic screens in zebrafish (Danio rerio) have isolated mutations in hundreds of genes with essential functions. To facilitate the identification of candidate genes for these mutations, we have genetically mapped 104 genes and expressed sequence tags by scoring single-strand conformational polymorphisms in a panel of haploid siblings. To integrate this map with existing genetic maps, we also scored 275 previously mapped genes, microsatellites, and sequence-tagged sites in the same haploid panel. Systematic phylogenetic analysis defined likely mammalian orthologs of mapped zebrafish genes, and comparison of map positions in zebrafish and mammals identified significant conservation of synteny. This comparative analysis also identified pairs of zebrafish genes that appear to be orthologous to single mammalian genes, suggesting that these genes arose in a genome duplication that occurred in the teleost lineage after the divergence of fish and mammal ancestors. This comparative map analysis will be useful in predicting the locations of zebrafish genes from mammalian gene maps and in understanding the evolution of the vertebrate genome.

The nucleotide sequence of Saccharomyces cerevisiae chromosome V.

Here we report the sequence of 569,202 base pairs of Saccharomyces cerevisiae chromosome V. Analysis of the sequence revealed a centromere, two telomeres and 271 open reading frames (ORFs) plus 13 tRNAs and four small nuclear RNAs. There are two Tyl transposable elements, each of which contains an ORF (included in the count of 271). Of the ORFs, 78 (29%) are new, 81 (30%) have potential homologues in the public databases, and 112 (41%) are previously characterized yeast genes.

Distinctive gene expression patterns in human mammary epithelial cells and breast cancers.

cDNA microarrays and a clustering algorithm were used to identify patterns of gene expression in human mammary epithelial cells growing in culture and in primary human breast tumors. Clusters of coexpressed genes identified through manipulations of mammary epithelial cells in vitro also showed consistent patterns of variation in expression among breast tumor samples. By using immunohistochemistry with antibodies against proteins encoded by a particular gene in a cluster, the identity of the cell type within the tumor specimen that contributed the observed gene expression pattern could be determined. Clusters of genes with coherent expression patterns in cultured cells and in the breast tumors samples could be related to specific features of biological variation among the samples. Two such clusters were found to have patterns that correlated with variation in cell proliferation rates and with activation of the IFN-regulated signal transduction pathway, respectively. Clusters of genes expressed by stromal cells and lymphocytes in the breast tumors also were identified in this analysis. These results support the feasibility and usefulness of this systematic approach to studying variation in gene expression patterns in human cancers as a means to dissect and classify solid tumors.