Generation of longer cDNA fragments from serial analysis of gene expression tags for gene identification

We have developed a technique called the generation of longer cDNA fragments from serial analysis of gene expression (SAGE) tags for gene identification (GLGI), to convert SAGE tags of 10 bases into their corresponding 3′ cDNA fragments covering hundred bases. A primer containing the 10-base SAGE tag is used as the sense primer, and a single base anchored oligo(dT) primer is used as an antisense primer in PCR, together with Pfu DNA polymerase. By using this approach, a cDNA fragment extending from the SAGE tag toward the 3′ end of the corresponding sequence can be generated. Application of the GLGI technique can solve two critical issues in applying the SAGE technique: one is that a longer fragment corresponding to a SAGE tag, which has no match in databases, can be generated for further studies; the other is that the specific fragment corresponding to a SAGE tag can be identified from multiple sequences that match the same SAGE tag. The development of the GLGI method provides several potential applications. First, it provides a strategy for even wider application of the SAGE technique for quantitative analysis of global gene expression. Second, a combined application of SAGE/GLGI can be used to complete the catalogue of the expressed genes in human and in other eukaryotic species. Third, it can be used to identify the 3′ cDNA sequence from any exon within a gene. It can also be used to confirm the reality of exons predicted by bioinformatic tools in genomic sequences. Fourth, a combined application of SAGE/GLGI can be applied to define the 3′ boundary of expressed genes in the genomic sequences in human and in other eukaryotic genomes.

Toward functional genomics in bacteria: Analysis of gene expression in Escherichia coli from a bacterial artificial chromosome library of Bacillus cereus

As the study of microbes moves into the era of functional genomics, there is an increasing need for molecular tools for analysis of a wide diversity of microorganisms. Currently, biological study of many prokaryotes of agricultural, medical, and fundamental scientific interest is limited by the lack of adequate genetic tools. We report the application of the bacterial artificial chromosome (BAC) vector to prokaryotic biology as a powerful approach to address this need. We constructed a BAC library in Escherichia coli from genomic DNA of the Gram-positive bacterium Bacillus cereus. This library provides 5.75-fold coverage of the B. cereus genome, with an average insert size of 98 kb. To determine the extent of heterologous expression of B. cereus genes in the library, we screened it for expression of several B. cereus activities in the E. coli host. Clones expressing 6 of 10 activities tested were identified in the library, namely, ampicillin resistance, zwittermicin A resistance, esculin hydrolysis, hemolysis, orange pigment production, and lecithinase activity. We analyzed selected BAC clones genetically to identify rapidly specific B. cereus loci. These results suggest that BAC libraries will provide a powerful approach for studying gene expression from diverse prokaryotes.

Identification of diverse nerve growth factor-regulated genes by serial analysis of gene expression (SAGE) profiling

Neurotrophic factors such as nerve growth factor (NGF) promote a wide variety of responses in neurons, including differentiation, survival, plasticity, and repair. Such actions often require changes in gene expression. To identify the regulated genes and thereby to more fully understand the NGF mechanism, we carried out serial analysis of gene expression (SAGE) profiling of transcripts derived from rat PC12 cells before and after NGF-promoted neuronal differentiation. Multiple criteria supported the reliability of the profile. Approximately 157,000 SAGE tags were analyzed, representing at least 21,000 unique transcripts. Of these, nearly 800 were regulated by 6-fold or more in response to NGF. Approximately 150 of the regulated transcripts have been matched to named genes, the majority of which were not previously known to be NGF-responsive. Functional categorization of the regulated genes provides insight into the complex, integrated mechanism by which NGF promotes its multiple actions. It is anticipated that as genomic sequence information accrues the data derived here will continue to provide information about neurotrophic factor mechanisms.

Saccharomyces Genome Database provides tools to survey gene expression and functional analysis data

Upon the completion of the Saccharomyces cerevisiae genomic sequence in 1996 [Goffeau,A. et al. (1997) Nature, 387, 5], several creative and ambitious projects have been initiated to explore the functions of gene products or gene expression on a genome-wide scale. To help researchers take advantage of these projects, the Saccharomyces Genome Database (SGD) has created two new tools, Function Junction and Expression Connection. Together, the tools form a central resource for querying multiple large-scale analysis projects for data about individual genes. Function Junction provides information from diverse projects that shed light on the role a gene product plays in the cell, while Expression Connection delivers information produced by the ever-increasing number of microarray projects. WWW access to SGD is available at genome-www.stanford.edu/Saccharomyces/.

High-sensitivity array analysis of gene expression for the early detection of disseminated breast tumor cells in peripheral blood

Early detection is an effective means of reducing cancer mortality. Here, we describe a highly sensitive high-throughput screen that can identify panels of markers for the early detection of solid tumor cells disseminated in peripheral blood. The method is a two-step combination of differential display and high-sensitivity cDNA arrays. In a primary screen, differential display identified 170 candidate marker genes differentially expressed between breast tumor cells and normal breast epithelial cells. In a secondary screen, high-sensitivity arrays assessed expression levels of these genes in 48 blood samples, 22 from healthy volunteers and 26 from breast cancer patients. Cluster analysis identified a group of 12 genes that were elevated in the blood of cancer patients. Permutation analysis of individual genes defined five core genes (P ≤ 0.05, permax test). As a group, the 12 genes generally distinguished accurately between healthy volunteers and patients with breast cancer. Mean expression levels of the 12 genes were elevated in 77% (10 of 13) untreated invasive cancer patients, whereas cluster analysis correctly classified volunteers and patients (P = 0.0022, Fisher's exact test). Quantitative real-time PCR confirmed array results and indicated that the sensitivity of the assay (1:2 × 108 transcripts) was sufficient to detect disseminated solid tumor cells in blood. Expression-based blood assays developed with the screening approach described here have the potential to detect and classify solid tumor cells originating from virtually any primary site in the body.

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).

Analysis of gene expression during myc oncogene-induced lymphomagenesis in the bursa of Fabricius

The transcriptional effects of deregulated myc gene overexpression are implicated in tumorigenesis in a spectrum of experimental and naturally occurring neoplasms. In follicles of the chicken bursa of Fabricius, myc induction of B-cell neoplasia requires a target cell population present during early bursal development and progresses through preneoplastic transformed follicles to metastatic lymphomas. We developed a chicken immune system cDNA microarray to analyze broad changes in gene expression that occur during normal embryonic B-cell development and during myc-induced neoplastic transformation in the bursa. The number of mRNAs showing at least 3-fold change was greater during myc-induced lymphomagenesis than during normal development, and hierarchical cluster analysis of expression patterns revealed that levels of several hundred mRNAs varied in concert with levels of myc overexpression. A set of 41 mRNAs were most consistently elevated in myc-overexpressing preneoplastic and neoplastic cells, most involved in processes thought to be subject to regulation by Myc. The mRNAs for another cluster of genes were overexpressed in neoplasia independent of myc expression level, including a small subset with the expression signature of embryonic bursal lymphocytes. Overexpression of myc, and some of the genes overexpressed with myc, may be important for generation of preneoplastic transformed follicles. However, expression profiles of late metastatic tumors showed a large variation in concert with myc expression levels, and some showed minimal myc overexpression. Therefore, high-level myc overexpression may be more important in the early induction of these lymphomas than in maintenance of late-stage metastases.

Comparative Analysis of the Regulation of Expression and Structures of Two Evolutionarily Divergent Genes for Δ1-Pyrroline-5-Carboxylate Synthetase from Tomato1

We isolated two tomato (Lycopersicon esculentum) cDNA clones, tomPRO1 and tomPRO2, specifying Δ1-pyrroline-5-carboxylate synthetase (P5CS), the first enzyme of proline (Pro) biosynthesis. tomPRO1 is unusual because it resembles prokaryotic polycistronic operons (M.G. García-Ríos, T. Fujita, P.C. LaRosa, R.D. Locy, J.M. Clithero, R.A. Bressan, L.N. Csonka [1997] Proc Natl Acad Sci USA 94: 8249–8254), whereas tomPRO2 encodes a full-length P5CS. We analyzed the accumulation of Pro and the tomPRO1 and tomPRO2 messages in response to NaCl stress and developmental signals. Treatment with 200 mm NaCl resulted in a >60-fold increase in Pro levels in roots and leaves. However, there was a <3-fold increase in the accumulation of the tomPRO2 message and no detectable induction in the level of the tomPRO1 message in response to NaCl stress. Although pollen contained approximately 100-fold higher levels of Pro than other plant tissues, there was no detectable increase in the level of either message in pollen. We conclude that transcriptional regulation of these genes for P5CS is probably not important for the osmotic or pollen-specific regulation of Pro synthesis in tomato. Using restriction fragment-length polymorphism mapping, we determined the locations of tomPRO1 and tomPRO2 loci in the tomato nuclear genome. Sequence comparison suggested that tomPRO1 is similar to prokaryotic P5CS loci, whereas tomPRO2 is closely related to other eukaryotic P5CS genes.

Parallel human genome analysis: microarray-based expression monitoring of 1000 genes.

Microarrays containing 1046 human cDNAs of unknown sequence were printed on glass with high-speed robotics. These 1.0-cm2 DNA "chips" were used to quantitatively monitor differential expression of the cognate human genes using a highly sensitive two-color hybridization assay. Array elements that displayed differential expression patterns under given experimental conditions were characterized by sequencing. The identification of known and novel heat shock and phorbol ester-regulated genes in human T cells demonstrates the sensitivity of the assay. Parallel gene analysis with microarrays provides a rapid and efficient method for large-scale human gene discovery.

Expression and analysis of presenilin 1 in a human neuronal system: localization in cell bodies and dendrites.

Mutations in the recently identified presenilin 1 gene on chromosome 14 cause early onset familial Alzheimer disease (FAD). Herein we describe the expression and analysis of the protein coded by presenilin 1 (PS1) in NT2N neurons, a human neuronal model system. PS1 was expressed using recombinant Semliki Forest virions and detected by introduced antigenic tags or antisera to PS1-derived peptides. Immunoprecipitation revealed two major PS1 bands of approximately 43 and 50 kDa, neither of which were N-glycosylated or O-glycosylated. Immunoreactive PS1 was detected in cell bodies and dendrites of NT2N neurons but not in axons or on the cell surface. PS1 was also detected in BHK cells, where it was also intracellular and colocalized with calnexin, a marker for the rough endoplasmic reticulum. A mutant form of PS1 linked to FAD did not differ from the wild-type protein at the light microscopic level. The model system described here will enable studies of the function of PS1 in human neurons and the role of mutant PS1 in FAD.

Analysis of lipopolysaccharide-response genes in B-lineage cells demonstrates that they can have differentiation stage-restricted expression and contain SH2 domains.

Bacterial lipopolysaccharide (LPS) is a potent stimulator of B-cell activation, proliferation, and differentiation. We examined the genetic response of B-lineage cells to LPS via trapping of expressed genes with a gene-trap retrovirus. This analysis showed that expression of only a small fraction of genes is altered during LPS stimulation of B-lineage cells. Isolation of the cellular portion of the trapped LPS-response genes via 5' RACE (rapid amplification of cDNA ends) cloning identified novel genes for all the cloned loci. These novel LPS-response genes were also found to have differentiation stage-restricted expression within the B-lymphoid lineage. That LPS-response genes in B cells also have differentiation stage-restricted expression suggests that these genes may be involved in the control of B-cell function and differentiation, since the known members of this class of genes have frequently been found to play a role in the function and differentiation of B-lineage cells. The isolation of novel members of this class of genes, including a gene that contains a putative SH2 domain, will further increase our understanding of the molecular events involved in the control of B-cell differentiation and function.

T-cell epitope analysis using subtracted expression libraries (TEASEL): application to a 38-kDA autoantigen recognized by T cells from an insulin-dependent diabetic patient.

Studies on circulating T cells and antibodies in newly diagnosed type 1 diabetic patients and rodent models of autoimmune diabetes suggest that beta-cell membrane proteins of 38 kDa may be important molecular targets of autoimmune attack. Biochemical approaches to the isolation and identification of the 38-kDa autoantigen have been hampered by the restricted availability of islet tissue and the low abundance of the protein. A procedure of epitope analysis for CD4+ T cells using subtracted expression libraries (TEASEL) was developed and used to clone a 70-amino acid pancreatic beta-cell peptide incorporating an epitope recognized by a 38-kDa-reactive CD4+ T-cell clone (1C6) isolated from a human diabetic patient. The minimal epitope was mapped to a 10-amino acid synthetic peptide containing a DR1 consensus binding motif. Data base searches did not reveal the identity of the protein, though a weak homology to the bacterial superantigens SEA (Streptococcus pyogenes exotoxin A) and SEB (Staphylococcus aureus enterotoxin B) (23% identity) was evident. The TEASEL procedure might be used to identify epitopes of other autoantigens recognized by CD4+ T cells in diabetes as well as be more generally applicable to the study low-abundance autoantigens in other tissue-specific autoimmune diseases.

Analysis of differential gene expression by display of 3' end restriction fragments of cDNAs.

We have developed an approach to study changes in gene expression by selective PCR amplification and display of 3' end restriction fragments of double-stranded cDNAs. This method produces highly consistent and reproducible patterns, can detect almost all mRNAs in a sample, and can resolve hidden differences such as bands that differ in their sequence but comigrate on a gel. Bands corresponding to known cDNAs move to predictable positions on the gel, making this a powerful approach to correlate gel patterns with cDNA data bases. Applying this method, we have examined differences in gene expression patterns during T-cell activation. Of a total of 700 bands that were evaluated in this study, as many as 3-4% represented mRNAs that are upregulated, while approximately 2% were down-regulated within 4 hr of activation of Jurkat T cells. These and other results suggest that this approach is suitable for the systematic, expeditious, and nearly exhaustive elucidation of subtle changes in the patterns of gene expression in cells with altered physiologic states.

Comparative expressed-sequence-tag analysis of differential gene expression profiles in PC-12 cells before and after nerve growth factor treatment.

Nerve growth factor-induced differentiation of adrenal chromaffin PC-12 cells to a neuronal phenotype involves alterations in gene expression and represents a model system to study neuronal differentiation. We have used the expressed-sequence-tag approach to identify approximately 600 differentially expressed mRNAs in untreated and nerve growth factor-treated PC-12 cells that encode proteins with diverse structural and biochemical functions. Many of these mRNAs encode proteins belonging to cellular pathways not previously known to be regulated by nerve growth factor. Comparative expressed-sequence-tag analysis provides a basis for surveying global changes in gene-expression patterns in response to biological signals at an unprecedented scale, is a powerful tool for identifying potential interactions between different cellular pathways, and allows the gene-expression profiles of individual genes belonging to a particular pathway to be followed.