Exploring the new world of the genome with DNA microarrays.

Thousands of genes are being discovered for the first time by sequencing the genomes of model organisms, an exhilarating reminder that much of the natural world remains to be explored at the molecular level. DNA microarrays provide a natural vehicle for this exploration. The model organisms are the first for which comprehensive genome-wide surveys of gene expression patterns or function are possible. The results can be viewed as maps that reflect the order and logic of the genetic program, rather than the physical order of genes on chromosomes. Exploration of the genome using DNA microarrays and other genome-scale technologies should narrow the gap in our knowledge of gene function and molecular biology between the currently-favoured model organisms and other species.

Large-scale identification of secreted and membrane-associated gene products using DNA microarrays.

Membrane-associated and secreted proteins are an important class of proteins and include receptors, transporters, adhesion molecules, hormones and cytokines. Although algorithms have been developed to recognize potential amino-terminal membrane-targeting signals or transmembrane domains in protein sequences, their accuracy is limited and they require knowledge of the entire coding sequence, including the N terminus, which is not currently available for most of the genes in most organisms, including human. Several experimental approaches for identifying secreted and membrane proteins have been described, but none have taken a comprehensive genomic approach. Furthermore, none of these methods allow easy classification of clones from arrayed cDNA libraries, for which large-scale gene-expression data are now becoming available through the use of DNA microarrays. We describe here a rapid and efficient method for identifying genes that encode secreted or membrane proteins. mRNA species bound to membrane-associated polysomes were separated from other mRNAs by sedimentation equilibrium or sedimentation velocity. The distribution of individual transcripts in the 'membrane-bound' and 'cytosolic' fractions was quantitated for thousands of genes by hybridization to DNA microarrays. Transcripts known to encode secreted or membrane proteins were enriched in the membrane-bound fractions, whereas those known to encode cytoplasmic proteins were enriched in the fractions containing mRNAs associated with free and cytoplasmic ribosomes. On this basis, we identified over 275 human genes and 285 yeast genes that are likely to encode previously unrecognized secreted or membrane proteins.

Promoter-specific binding of Rap1 revealed by genome-wide maps of protein-DNA association.

We determined the distribution of repressor-activator protein 1 (Rap1) and the accessory silencing proteins Sir2, Sir3 and Sir4 in vivo on the entire yeast genome, at a resolution of 2 kb. Rap1 is central to the cellular economy during rapid growth, targeting 294 loci, about 5% of yeast genes, and participating in the activation of 37% of all RNA polymerase II initiation events in exponentially growing cells. Although the DNA sequence recognized by Rap1 is found in both coding and intergenic sequences, the binding of Rap1 to the genome was highly specific to intergenic regions with the potential to act as promoters. This global phenomenon, which may be a general characteristic of sequence-specific transcriptional factors, indicates the existence of a genome-wide molecular mechanism for marking promoter regions.

Genome-wide analysis of DNA copy-number changes using cDNA microarrays.

Gene amplifications and deletions frequently contribute to tumorigenesis. Characterization of these DNA copy-number changes is important for both the basic understanding of cancer and its diagnosis. Comparative genomic hybridization (CGH) was developed to survey DNA copy-number variations across a whole genome. With CGH, differentially labelled test and reference genomic DNAs are co-hybridized to normal metaphase chromosomes, and fluorescence ratios along the length of chromosomes provide a cytogenetic representation of DNA copy-number variation. CGH, however, has a limited ( approximately 20 Mb) mapping resolution, and higher-resolution techniques, such as fluorescence in situ hybridization (FISH), are prohibitively labour-intensive on a genomic scale. Array-based CGH, in which fluorescence ratios at arrayed DNA elements provide a locus-by-locus measure of DNA copy-number variation, represents another means of achieving increased mapping resolution. Published array CGH methods have relied on large genomic clone (for example BAC) array targets and have covered only a small fraction of the human genome. cDNAs representing over 30,000 radiation-hybrid (RH)-mapped human genes provide an alternative and readily available genomic resource for mapping DNA copy-number changes. Although cDNA microarrays have been used extensively to characterize variation in human gene expression, human genomic DNA is a far more complex mixture than the mRNA representation of human cells. Therefore, analysis of DNA copy-number variation using cDNA microarrays would require a sensitivity of detection an order of magnitude greater than has been routinely reported. We describe here a cDNA microarray-based CGH method, and its application to DNA copy-number variation analysis in breast cancer cell lines and tumours. Using this assay, we were able to identify gene amplifications and deletions genome-wide and with high resolution, and compare alterations in DNA copy number and gene expression.

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.

A gene expression database for the molecular pharmacology of cancer.

We used cDNA microarrays to assess gene expression profiles in 60 human cancer cell lines used in a drug discovery screen by the National Cancer Institute. Using these data, we linked bioinformatics and chemoinformatics by correlating gene expression and drug activity patterns in the NCI60 lines. Clustering the cell lines on the basis of gene expression yielded relationships very different from those obtained by clustering the cell lines on the basis of their response to drugs. Gene-drug relationships for the clinical agents 5-fluorouracil and L-asparaginase exemplify how variations in the transcript levels of particular genes relate to mechanisms of drug sensitivity and resistance. This is the first study to integrate large databases on gene expression and molecular pharmacology.

Relation of gene expression phenotype to immunoglobulin mutation genotype in B cell chronic lymphocytic leukemia.

The most common human leukemia is B cell chronic lymphocytic leukemia (CLL), a malignancy of mature B cells with a characteristic clinical presentation but a variable clinical course. The rearranged immunoglobulin (Ig) genes of CLL cells may be either germ-line in sequence or somatically mutated. Lack of Ig mutations defined a distinctly worse prognostic group of CLL patients raising the possibility that CLL comprises two distinct diseases. Using genomic-scale gene expression profiling, we show that CLL is characterized by a common gene expression "signature," irrespective of Ig mutational status, suggesting that CLL cases share a common mechanism of transformation and/or cell of origin. Nonetheless, the expression of hundreds of other genes correlated with the Ig mutational status, including many genes that are modulated in expression during mitogenic B cell receptor signaling. These genes were used to build a CLL subtype predictor that may help in the clinical classification of patients with this disease.

Yeast actin-binding proteins: evidence for a role in morphogenesis.

Three yeast actin-binding proteins were identified using yeast actin filaments as an affinity matrix. One protein appears to be a yeast myosin heavy chain; it is dissociated from actin filaments by ATP, it is similar in size (200 kD) to other myosins, and antibodies directed against Dictyostelium myosin heavy chain bind to it. Immunofluorescence experiments show that a second actin-binding protein (67 kD) colocalizes in vivo with both cytoplasmic actin cables and cortical actin patches, the only identifiable actin structures in yeast. The cortical actin patches are concentrated at growing surfaces of the yeast cell where they might play a role in membrane and cell wall insertion, and the third actin-binding protein (85 kD) is only detected in association with these structures. This 85-kD protein is therefore a candidate for a determinant of growth sites. The in vivo role of this protein was tested by overproduction; this overproduction causes a reorganization of the actin cytoskeleton which in turn dramatically affects the budding pattern and spatial growth organization of the yeast cell.

Diverse effects of beta-tubulin mutations on microtubule formation and function.

We have used in vitro mutagenesis and gene replacement to construct five new cold-sensitive mutations in TUB2, the sole gene encoding beta-tubulin in the yeast Saccharomyces cerevisiae. These and one previously isolated tub2 mutant display diverse phenotypes that have allowed us to define the functions of yeast microtubules in vivo. At the restrictive temperature, all of the tub2 mutations inhibit chromosome segregation and block the mitotic cell cycle. However, different microtubule arrays are present in these arrested cells depending on the tub2 allele. One mutant (tub2-401) contains no detectable microtubules, two (tub2-403 and tub2-405) contain greatly diminished levels of both nuclear and cytoplasmic microtubules, one (tub2-104) contains predominantly nuclear microtubules, one (tub2-402) contains predominantly cytoplasmic microtubules, and one (tub2-404) contains prominent nuclear and cytoplasmic microtubule arrays. Using these mutants we demonstrate here that cytoplasmic microtubules are necessary for nuclear migration during the mitotic cell cycle and for nuclear migration and fusion during conjugation; only those mutants that possess cytoplasmic microtubules are able to perform these functions. We also show that microtubules are not required for secretory vesicle transport in yeast; bud growth and invertase secretion occur in cells which contain no microtubules.

Yeast cell cycle protein CDC48p shows full-length homology to the mammalian protein VCP and is a member of a protein family involved in secretion, peroxisome formation, and gene expression.

Yeast mutants of cell cycle gene cdc48-1 arrest as large budded cells with microtubules spreading aberrantly throughout the cytoplasm from a single spindle plaque. The gene was cloned and disruption proved it to be essential. The CDC48 sequence encodes a protein of 92 kD that has an internal duplication of 200 amino acids and includes a nucleotide binding consensus sequence. Vertebrate VCP has a 70% identity over the entire length of the protein. Yeast Sec18p and mammalian N-ethylmaleimide-sensitive fusion protein, which are involved in intracellular transport, yeast Pas1p, which is essential for peroxisome assembly, and mammalian TBP-1, which influences HIV gene expression, are 40% identical in the duplicated region. Antibodies against CDC48 recognize a yeast protein of apparently 115 kD and a mammalian protein of 100 kD. Both proteins are bound loosely to components of the microsomal fraction as described for Sec18p and N-ethylmaleimide-sensitive fusion protein. This similarity suggests that CDC48p participates in a cell cycle function related to that of N-ethylmaleimide-sensitive fusion protein/Sec18p in Golgi transport.

NDC1: a nuclear periphery component required for yeast spindle pole body duplication.

The spindle pole body (SPB) of Saccharomyces cerevisiae serves as the centrosome in this organism, undergoing duplication early in the cell cycle to generate the two poles of the mitotic spindle. The conditional lethal mutation ndc1-1 has previously been shown to cause asymmetric segregation, wherein all the chromosomes go to one pole of the mitotic spindle (Thomas, J. H., and D. Botstein. 1986. Cell. 44:65-76). Examination by electron microscopy of mutant cells subjected to the nonpermissive temperature reveals a defect in SPB duplication. Although duplication is seen to occur, the nascent SPB fails to undergo insertion into the nuclear envelope. The parental SPB remains functional, organizing a monopolar spindle to which all the chromosomes are presumably attached. Order-of-function experiments reveal that the NDC1 function is required in G1 after alpha-factor arrest but before the arrest caused by cdc34. Molecular analysis shows that the NDC1 gene is essential and that it encodes a 656 amino acid protein (74 kD) with six or seven putative transmembrane domains. This evidence for membrane association is further supported by immunofluorescent localization of the NDC1 product to the vicinity of the nuclear envelope. These findings suggest that the NDC1 protein acts within the nuclear envelope to mediate insertion of the nascent SPB.

Ultrastructure of the yeast actin cytoskeleton and its association with the plasma membrane.

We characterized the yeast actin cytoskeleton at the ultrastructural level using immunoelectron microscopy. Anti-actin antibodies primarily labeled dense, patchlike cortical structures and cytoplasmic cables. This localization recapitulates results obtained with immunofluorescence light microscopy, but at much higher resolution. Immuno-EM double-labeling experiments were conducted with antibodies to actin together with antibodies to the actin binding proteins Abp1p and cofilin. As expected from immunofluorescence experiments, Abp1p, cofilin, and actin colocalized in immuno-EM to the dense patchlike structures but not to the cables. In this way, we can unambiguously identify the patches as the cortical actin cytoskeleton. The cortical actin patches were observed to be associated with the cell surface via an invagination of plasma membrane. This novel cortical cytoskeleton-plasma membrane interface appears to consist of a fingerlike invagination of plasma membrane around which actin filaments and actin binding proteins are organized. We propose a possible role for this unique cortical structure in wall growth and osmotic regulation.

Yeast: an experimental organism for modern biology.

The yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe have become popular and successful model systems for understanding eukaryotic biology at the cellular and molecular levels. The reasons for this success are experimental tractability, especially in applying classical and molecular genetic methods to associate genes with proteins and functions within the cell.

Homozygosity mapping: a way to map human recessive traits with the DNA of inbred children.

An efficient strategy for mapping human genes that cause recessive traits has been devised that uses mapped restriction fragment length polymorphisms (RFLPs) and the DNA of affected children from consanguineous marriages. The method involves detection of the disease locus by virtue of the fact that the adjacent region will preferentially be homozygous by descent in such inbred children. A single affected child of a first-cousin marriage is shown to contain the same total information about linkage as a nuclear family with three affected children. Calculations show that it should be practical to map a recessive disease gene by studying DNA from fewer than a dozen unrelated, affected inbred children, given a complete RFLP linkage map. The method should make it possible to map many recessive diseases for which it is impractical or impossible to collect adequate numbers of families with multiple affected offspring.

Diverse essential functions revealed by complementing yeast calmodulin mutants.

Calmodulin, a cytoplasmic calcium-binding protein, is indispensable for eukaryotic cell growth. Examination of 14 temperature-sensitive yeast mutants bearing one or more phenylalanine to alanine substitutions in the single essential calmodulin gene of yeast (CMD1) revealed diverse essential functions. Mutations could be classified into four intragenic complementation groups. Each group showed different characteristic functional defects in actin organization, calmodulin localization, nuclear division, or bud emergence. Phenylalanine residues implicated in calmodulin localization and nuclear division are located in the amino-terminal half of the protein, whereas those implicated in actin organization and bud emergence are located in the carboxyl-terminal half.

A yeast actin-binding protein is encoded by SAC6, a gene found by suppression of an actin mutation.

The protein encoded by SAC6, a gene that can mutate to suppress a temperature-sensitive defect in the yeast actin gene, has been identified as a 67-kilodalton actin-binding protein (ABP 67) that associates with all identifiable actin structures. This finding demonstrates the in vivo functional importance of the actin-ABP 67 interaction previously established in vitro and illustrates the use of suppressor analysis to identify physically interacting proteins.

Lethal disruption of the yeast actin gene by integrative DNA transformation.

A mutant allele of the chromosomal locus corresponding to the cloned actin gene of the yeast Saccharomyces cerevisiae has been constructed by DNA transformation with a hybrid plasmid which integrates into, and thereby disrupts, the protein-encoding sequences of the gene. In a diploid strain of yeast, disruption of the actin gene on one chromosome results in a mutation that segregates as a recessive lethal tightly linked to a selectable genetic marker on the integrated plasmid. The actin gene, therefore, must encode an essential function for yeast cell growth.

Many random sequences functionally replace the secretion signal sequence of yeast invertase.

In the process of protein secretion, amino-terminal signal sequences are key recognition elements; however, the relation between the primary sequence of an amino-terminal peptide and its ability to function as an export signal remains obscure. The limits of variation permitted for functional signal sequences were determined by replacement of the normal signal sequence of Saccharomyces cerevisiae invertase with essentially random peptide sequences. Since about one-fifth of these sequences can function as an export signal the specificity with which signal sequences are recognized must be very low.

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 of sporulation in budding yeast.

Diploid cells of budding yeast produce haploid cells through the developmental program of sporulation, which consists of meiosis and spore morphogenesis. DNA microarrays containing nearly every yeast gene were used to assay changes in gene expression during sporulation. At least seven distinct temporal patterns of induction were observed. The transcription factor Ndt80 appeared to be important for induction of a large group of genes at the end of meiotic prophase. Consensus sequences known or proposed to be responsible for temporal regulation could be identified solely from analysis of sequences of coordinately expressed genes. The temporal expression pattern provided clues to potential functions of hundreds of previously uncharacterized genes, some of which have vertebrate homologs that may function during gametogenesis.