The Stanford Microarray Database

The Stanford Microarray Database (SMD) stores raw and normalized data from microarray experiments, and provides web interfaces for researchers to retrieve, analyze and visualize their data. The two immediate goals for SMD are to serve as a storage site for microarray data from ongoing research at Stanford University, and to facilitate the public dissemination of that data once published, or released by the researcher. Of paramount importance is the connection of microarray data with the biological data that pertains to the DNA deposited on the microarray (genes, clones etc.). SMD makes use of many public resources to connect expression information to the relevant biology, including SGD [Ball,C.A., Dolinski,K., Dwight,S.S., Harris,M.A., Issel-Tarver,L., Kasarskis,A., Scafe,C.R., Sherlock,G., Binkley,G., Jin,H. et al. (2000) Nucleic Acids Res., 28, 77–80], YPD and WormPD [Costanzo,M.C., Hogan,J.D., Cusick,M.E., Davis,B.P., Fancher,A.M., Hodges,P.E., Kondu,P., Lengieza,C., Lew-Smith,J.E., Lingner,C. et al. (2000) Nucleic Acids Res., 28, 73–76], Unigene [Wheeler,D.L., Chappey,C., Lash,A.E., Leipe,D.D., Madden,T.L., Schuler,G.D., Tatusova,T.A. and Rapp,B.A. (2000) Nucleic Acids Res., 28, 10–14], dbEST [Boguski,M.S., Lowe,T.M. and Tolstoshev,C.M. (1993) Nature Genet., 4, 332–333] and SWISS-PROT [Bairoch,A. and Apweiler,R. (2000) Nucleic Acids Res., 28, 45–48] and can be accessed at http://genome-www.stanford.edu/microarray.

Delineating developmental and metabolic pathways in vivo by expression profiling using the RIKEN set of 18,816 full-length enriched mouse cDNA arrays

We have systematically characterized gene expression patterns in 49 adult and embryonic mouse tissues by using cDNA microarrays with 18,816 mouse cDNAs. Cluster analysis defined sets of genes that were expressed ubiquitously or in similar groups of tissues such as digestive organs and muscle. Clustering of expression profiles was observed in embryonic brain, postnatal cerebellum, and adult olfactory bulb, reflecting similarities in neurogenesis and remodeling. Finally, clustering genes coding for known enzymes into 78 metabolic pathways revealed a surprising coordination of expression within each pathway among different tissues. On the other hand, a more detailed examination of glycolysis revealed tissue-specific differences in profiles of key regulatory enzymes. Thus, by surveying global gene expression by using microarrays with a large number of elements, we provide insights into the commonality and diversity of pathways responsible for the development and maintenance of the mammalian body plan.

An orchestrated gene expression component of neuronal programmed cell death revealed by cDNA array analysis

Programmed cell death (PCD) during neuronal development and disease has been shown to require de novo RNA synthesis. However, the time course and regulation of target genes is poorly understood. By using a brain-biased array of over 7,500 cDNAs, we profiled this gene expression component of PCD in cerebellar granule neurons challenged separately by potassium withdrawal, combined potassium and serum withdrawal, and kainic acid administration. We found that hundreds of genes were significantly regulated in discreet waves including known genes whose protein products are involved in PCD. A restricted set of genes was regulated by all models, providing evidence that signals inducing PCD can regulate large assemblages of genes (of which a restricted subset may be shared in multiple pathways).

GenEST, a powerful bidirectional link between cDNA sequence data and gene expression profiles generated by cDNA-AFLP

The release of vast quantities of DNA sequence data by large-scale genome and expressed sequence tag (EST) projects underlines the necessity for the development of efficient and inexpensive ways to link sequence databases with temporal and spatial expression profiles. Here we demonstrate the power of linking cDNA sequence data (including EST sequences) with transcript profiles revealed by cDNA-AFLP, a highly reproducible differential display method based on restriction enzyme digests and selective amplification under high stringency conditions. We have developed a computer program (GenEST) that predicts the sizes of virtual transcript-derived fragments (TDFs) of in silico-digested cDNA sequences retrieved from databases. The vast majority of the resulting virtual TDFs could be traced back among the thousands of TDFs displayed on cDNA-AFLP gels. Sequencing of the corresponding bands excised from cDNA-AFLP gels revealed no inconsistencies. As a consequence, cDNA sequence databases can be screened very efficiently to identify genes with relevant expression profiles. The other way round, it is possible to switch from cDNA-AFLP gels to sequences in the databases. Using the restriction enzyme recognition sites, the primer extensions and the estimated TDF size as identifiers, the DNA sequence(s) corresponding to a TDF with an interesting expression pattern can be identified. In this paper we show examples in both directions by analyzing the plant parasitic nematode Globodera rostochiensis. Various novel pathogenicity factors were identified by combining ESTs from the infective stage juveniles with expression profiles of ∼4000 genes in five developmental stages produced by cDNA-AFLP.

Characterization of Two Divergent Endo-β-1,4-Glucanase cDNA Clones Highly Expressed in the Nonclimacteric Strawberry Fruit

Two cDNAs clones (Cel1 and Cel2) encoding divergent endo-β-1,4-glucanases (EGases) have been isolated from a cDNA library obtained from ripe strawberry (Fragaria x ananassa Duch) fruit. The analysis of the amino acid sequence suggests that Cel1 and Cel2 EGases have different secondary and tertiary structures and that they differ in the presence of potential N-glycosylation sites. By in vitro translation we show that Cel1 and Cel2 bear a functional signal peptide, the cleavage of which yields mature proteins of 52 and 60 kD, respectively. Phylogenetic analysis revealed that the Cel2 EGase diverged early in evolution from other plant EGases. Northern analysis showed that both EGases are highly expressed in fruit and that they have different temporal patterns of accumulation. The Cel2 EGase was expressed in green fruit, accumulating as the fruit turned from green to white and remaining at an elevated, constant level throughout fruit ripening. In contrast, the Cel1 transcript was not detected in green fruit and only a low level of expression was observed in white fruit. The level of Cel1 mRNA increased gradually during ripening, reaching a maximum in fully ripe fruit. The high levels of Cel1 and Cel2 mRNA in ripe fruit and their overlapping patterns of expression suggest that these EGases play an important role in softening during ripening. In addition, the early expression of Cel2 in green fruit, well before significant softening begins, suggests that the product of this gene may also be involved in processes other than fruit softening, e.g. cell wall expansion.

Purification of the Trehalase GMTRE1 from Soybean Nodules and Cloning of Its cDNA. GMTRE1 Is Expressed at a Low Level in Multiple Tissues1

Trehalose (α-d-glucopyranosyl-1,1-α-d-glucopyranoside), a disaccharide widespread among microbes and lower invertebrates, is generally believed to be nonexistent in higher plants. However, the recent discovery of Arabidopsis genes whose products are involved in trehalose synthesis has renewed interest in the possibility of a function of trehalose in higher plants. We previously showed that trehalase, the enzyme that degrades trehalose, is present in nodules of soybean (Glycine max [L.] Merr.), and we characterized the enzyme as an apoplastic glycoprotein. Here we describe the purification of this trehalase to homogeneity and the cloning of a full-length cDNA encoding this enzyme, named GMTRE1 (G. max trehalase 1). The amino acid sequence derived from the open reading frame of GMTRE1 shows strong homology to known trehalases from bacteria, fungi, and animals. GMTRE1 is a single-copy gene and is expressed at a low but constant level in many tissues.

Purification and cDNA Cloning of Isochorismate Synthase from Elicited Cell Cultures of Catharanthus roseus

Isochorismate is an important metabolite formed at the end of the shikimate pathway, which is involved in the synthesis of both primary and secondary metabolites. It is synthesized from chorismate in a reaction catalyzed by the enzyme isochorismate synthase (ICS; EC 5.4.99.6). We have purified ICS to homogeneity from elicited Catharanthus roseus cell cultures. Two isoforms with an apparent molecular mass of 64 kD were purified and characterized. The Km values for chorismate were 558 and 319 μm for isoforms I and II, respectively. The isoforms were not inhibited by aromatic amino acids and required Mg2+ for enzyme activity. Polymerase chain reaction on a cDNA library from elicited C. roseus cells with a degenerated primer based on the sequence of an internal peptide from isoform II resulted in an amplification product that was used to screen the cDNA library. This led to the first isolation, to our knowledge, of a plant ICS cDNA. The cDNA encodes a protein of 64 kD with an N-terminal chloroplast-targeting signal. The deduced amino acid sequence shares homology with bacterial ICS and also with anthranilate synthases from plants. Southern analysis indicates the existence of only one ICS gene in C. roseus.

Amygdala-enriched genes identified by microarray technology are restricted to specific amygdaloid subnuclei

Microarray technology represents a potentially powerful method for identifying cell type- and regionally restricted genes expressed in the brain. Here we have combined a microarray analysis of differential gene expression among five selected brain regions, including the amygdala, cerebellum, hippocampus, olfactory bulb, and periaqueductal gray, with in situ hybridization. On average, 0.3% of the 34,000 genes interrogated were highly enriched in each of the five regions, relative to the others. In situ hybridization performed on a subset of amygdala-enriched genes confirmed in most cases the overall region-specificity predicted by the microarray data and identified additional sites of brain expression not examined on the microarrays. Strikingly, the majority of these genes exhibited boundaries of expression within the amygdala corresponding to cytoarchitectonically defined subnuclei. These results define a unique set of molecular markers for amygdaloid subnuclei and provide tools to genetically dissect their functional roles in different emotional behaviors.

Statistical modeling of large microarray data sets to identify stimulus-response profiles

A statistical modeling approach is proposed for use in searching large microarray data sets for genes that have a transcriptional response to a stimulus. The approach is unrestricted with respect to the timing, magnitude or duration of the response, or the overall abundance of the transcript. The statistical model makes an accommodation for systematic heterogeneity in expression levels. Corresponding data analyses provide gene-specific information, and the approach provides a means for evaluating the statistical significance of such information. To illustrate this strategy we have derived a model to depict the profile expected for a periodically transcribed gene and used it to look for budding yeast transcripts that adhere to this profile. Using objective criteria, this method identifies 81% of the known periodic transcripts and 1,088 genes, which show significant periodicity in at least one of the three data sets analyzed. However, only one-quarter of these genes show significant oscillations in at least two data sets and can be classified as periodic with high confidence. The method provides estimates of the mean activation and deactivation times, induced and basal expression levels, and statistical measures of the precision of these estimates for each periodic transcript.

Recursive partitioning for tumor classification with gene expression microarray data

Precise classification of tumors is critically important for cancer diagnosis and treatment. It is also a scientifically challenging task. Recently, efforts have been made to use gene expression profiles to improve the precision of classification, with limited success. Using a published data set for purposes of comparison, we introduce a methodology based on classification trees and demonstrate that it is significantly more accurate for discriminating among distinct colon cancer tissues than other statistical approaches used heretofore. In addition, competing classification trees are displayed, which suggest that different genes may coregulate colon cancers.

Screening insertion libraries for mutations in many genes simultaneously using DNA microarrays

We describe a method to screen pools of DNA from multiple transposon lines for insertions in many genes simultaneously. We use thermal asymmetric interlaced–PCR, a hemispecific PCR amplification protocol that combines nested, insertion-specific primers with degenerate primers, to amplify DNA flanking the transposons. In reconstruction experiments with previously characterized Arabidopsis lines carrying insertions of the maize Dissociation (Ds) transposon, we show that fluorescently labeled, transposon-flanking fragments overlapping ORFs hybridize to cognate expressed sequence tags (ESTs) on a DNA microarray. We further show that insertions can be detected in DNA pools from as many as 100 plants representing different transposon lines and that all of the tested, transposon-disrupted genes whose flanking fragments can be amplified individually also can be detected when amplified from the pool. The ability of a transposon-flanking fragment to hybridize declines rapidly with decreasing homology to the spotted DNA fragment, so that only ESTs with >90% homology to the transposon-disrupted gene exhibit significant cross-hybridization. Because thermal asymmetric interlaced–PCR fragments tend to be short, use of the present method favors recovery of insertions in and near genes. We apply the technique to screening pools of new Ds lines using cDNA microarrays containing ESTs for ≈1,000 stress-induced and -repressed Arabidopsis genes.

Purification and cDNA Cloning of Maize Poly(ADP)-Ribose Polymerase

Poly(ADP)-ribose polymerase (PADPRP) has been purified to apparent homogeneity from suspension cultures of the maize (Zea mays) callus line. The purified enzyme is a single polypeptide of approximately 115 kD, which appears to dimerize through an S-S linkage. The catalytic properties of the maize enzyme are very similar to those of its animal counterpart. The amino acid sequences of three tryptic peptides were obtained by microsequencing. Antibodies raised against peptides from maize PADPRP cross-reacted specifically with the maize enzyme but not with the enzyme from human cells, and vice versa. We have also characterized a 3.45-kb expressed-sequence-tag clone that contains a full-length cDNA for maize PADPRP. An open reading frame of 2943 bp within this clone encodes a protein of 980 amino acids. The deduced amino acid sequence of the maize PADPRP shows 40% to 42% identity and about 50% similarity to the known vertebrate PADPRP sequences. All important features of the modular structure of the PADPRP molecule, such as two zinc fingers, a putative nuclear localization signal, the automodification domain, and the NAD+-binding domain, are conserved in the maize enzyme. Northern-blot analysis indicated that the cDNA probe hybridizes to a message of about 4 kb.

Biosynthesis of Lipoic Acid in Arabidopsis: Cloning and Characterization of the cDNA for Lipoic Acid Synthase1

Lipoic acid is a coenzyme that is essential for the activity of enzyme complexes such as those of pyruvate dehydrogenase and glycine decarboxylase. We report here the isolation and characterization of LIP1 cDNA for lipoic acid synthase of Arabidopsis. The Arabidopsis LIP1 cDNA was isolated using an expressed sequence tag homologous to the lipoic acid synthase of Escherichia coli. This cDNA was shown to code for Arabidopsis lipoic acid synthase by its ability to complement a lipA mutant of E. coli defective in lipoic acid synthase. DNA-sequence analysis of the LIP1 cDNA revealed an open reading frame predicting a protein of 374 amino acids. Comparisons of the deduced amino acid sequence with those of E. coli and yeast lipoic acid synthase homologs showed a high degree of sequence similarity and the presence of a leader sequence presumably required for import into the mitochondria. Southern-hybridization analysis suggested that LIP1 is a single-copy gene in Arabidopsis. Western analysis with an antibody against lipoic acid synthase demonstrated that this enzyme is located in the mitochondrial compartment in Arabidopsis cells as a 43-kD polypeptide.

Overexpression of an Arabidopsis cDNA Encoding a Sterol-C241-Methyltransferase in Tobacco Modifies the Ratio of 24-Methyl Cholesterol to Sitosterol and Is Associated with Growth Reduction

Higher plants synthesize 24-methyl sterols and 24-ethyl sterols in defined proportions. As a first step in investigating the physiological function of this balance, an Arabidopsis cDNA encoding an S-adenosyl-l-methionine 24-methylene lophenol-C241-methyltransferase, the typical plant enzyme responsible for the production of 24-ethyl sterols, was expressed in tobacco (Nicotiana tabacum L.) under the control of a constitutive promoter. Transgenic plants displayed a novel 24-alkyl-Δ5-sterol profile: the ratio of 24-methyl cholesterol to sitosterol, which is close to 1 in the wild type, decreased dramatically to values ranging from 0.01 to 0.31. In succeeding generations of transgenic tobacco, a high S-adenosyl-l-methionine 24-methylene lophenol-C241-methyltransferase enzyme activity and, consequently, a low ratio of 24-methyl cholesterol to sitosterol, was associated with reduced growth compared with the wild type. However, this new morphological phenotype appeared only below the threshold ratio of 24-methyl cholesterol to sitosterol of approximately 0.1. Because the size of cells was unchanged in small, transgenic plants, we hypothesize that a radical decrease of 24-methyl cholesterol and/or a concomitant increase of sitosterol would be responsible for a change in cell division through as-yet unknown mechanisms.

Expressed Sequence Tags from a Root-Hair-Enriched Medicago truncatula cDNA Library1

The root hair is a specialized cell type involved in water and nutrient uptake in plants. In legumes the root hair is also the primary site of recognition and infection by symbiotic nitrogen-fixing Rhizobium bacteria. We have studied the root hairs of Medicago truncatula, which is emerging as an increasingly important model legume for studies of symbiotic nodulation. However, only 27 genes from M. truncatula were represented in GenBank/EMBL as of October, 1997. We report here the construction of a root-hair-enriched cDNA library and single-pass sequencing of randomly selected clones. Expressed sequence tags (899 total, 603 of which have homology to known genes) were generated and made available on the Internet. We believe that the database and the associated DNA materials will provide a useful resource to the community of scientists studying the biology of roots, root tips, root hairs, and nodulation.