Dynamics of Gene Expression Revealed by Comparison of Serial Analysis of Gene Expression Transcript Profiles from Yeast Grown on Two Different Carbon SourcesD⃞

We describe a genome-wide characterization of mRNA transcript levels in yeast grown on the fatty acid oleate, determined using Serial Analysis of Gene Expression (SAGE). Comparison of this SAGE library with that reported for glucose grown cells revealed the dramatic adaptive response of yeast to a change in carbon source. A major fraction (>20%) of the 15,000 mRNA molecules in a yeast cell comprised differentially expressed transcripts, which were derived from only 2% of the total number of ∼6300 yeast genes. Most of the mRNAs that were differentially expressed code for enzymes or for other proteins participating in metabolism (e.g., metabolite transporters). In oleate-grown cells, this was exemplified by the huge increase of mRNAs encoding the peroxisomal β-oxidation enzymes required for degradation of fatty acids. The data provide evidence for the existence of redox shuttles across organellar membranes that involve peroxisomal, cytoplasmic, and mitochondrial enzymes. We also analyzed the mRNA profile of a mutant strain with deletions of the PIP2 and OAF1 genes, encoding transcription factors required for induction of genes encoding peroxisomal proteins. Induction of genes under the immediate control of these factors was abolished; other genes were up-regulated, indicating an adaptive response to the changed metabolism imposed by the genetic impairment. We describe a statistical method for analysis of data obtained by SAGE.

Global analysis of gene expression in pulmonary fibrosis reveals distinct programs regulating lung inflammation and fibrosis

The molecular mechanisms of pulmonary fibrosis are poorly understood. We have used oligonucleotide arrays to analyze the gene expression programs that underlie pulmonary fibrosis in response to bleomycin, a drug that causes lung inflammation and fibrosis, in two strains of susceptible mice (129 and C57BL/6). We then compared the gene expression patterns in these mice with 129 mice carrying a null mutation in the epithelial-restricted integrin β6 subunit (β6−/−), which develop inflammation but are protected from pulmonary fibrosis. Cluster analysis identified two distinct groups of genes involved in the inflammatory and fibrotic responses. Analysis of gene expression at multiple time points after bleomycin administration revealed sequential induction of subsets of genes that characterize each response. The availability of this comprehensive data set should accelerate the development of more effective strategies for intervention at the various stages in the development of fibrotic diseases of the lungs and other organs.

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.

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.

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.

The pattern of gene expression in human CD15+ myeloid progenitor cells

We performed a genome-wide analysis of gene expression in primary human CD15+ myeloid progenitor cells. By using the serial analysis of gene expression (SAGE) technique, we obtained quantitative information for the expression of 37,519 unique SAGE-tag sequences. Of these unique tags, (i) 25% were detected at high and intermediate levels, whereas 75% were present as single copies, (ii) 53% of the tags matched known expressed sequences, 34% of which were matched to more than one known expressed sequence, and (iii) 47% of the tags had no matches and represent potentially novel genes. The correct genes were confirmed by application of the generation of longer cDNA fragments from SAGE tags for gene identification (GLGI) technique for high-copy tags with multiple matches. A set of genes known to be important in myeloid differentiation were expressed at various levels and used different spliced forms. This study provides a normal baseline for comparison of gene expression in myeloid diseases. The strategy of using SAGE and GLGI techniques in this study has broad applications to the genome-wide identification of expressed genes.

Importance of replication in microarray gene expression studies: Statistical methods and evidence from repetitive cDNA hybridizations

We present statistical methods for analyzing replicated cDNA microarray expression data and report the results of a controlled experiment. The study was conducted to investigate inherent variability in gene expression data and the extent to which replication in an experiment produces more consistent and reliable findings. We introduce a statistical model to describe the probability that mRNA is contained in the target sample tissue, converted to probe, and ultimately detected on the slide. We also introduce a method to analyze the combined data from all replicates. Of the 288 genes considered in this controlled experiment, 32 would be expected to produce strong hybridization signals because of the known presence of repetitive sequences within them. Results based on individual replicates, however, show that there are 55, 36, and 58 highly expressed genes in replicates 1, 2, and 3, respectively. On the other hand, an analysis by using the combined data from all 3 replicates reveals that only 2 of the 288 genes are incorrectly classified as expressed. Our experiment shows that any single microarray output is subject to substantial variability. By pooling data from replicates, we can provide a more reliable analysis of gene expression data. Therefore, we conclude that designing experiments with replications will greatly reduce misclassification rates. We recommend that at least three replicates be used in designing experiments by using cDNA microarrays, particularly when gene expression data from single specimens are being analyzed.

In vivo detection of gene expression in liver by 31P nuclear magnetic resonance spectroscopy employing creatine kinase as a marker gene

In vivo assessment of gene expression is desirable to obtain information on the extent and duration of transduction of tissue after gene delivery. We have developed an in vivo, potentially noninvasive, method for detecting virally mediated gene transfer to the liver. The method employs an adenoviral vector carrying the gene for the brain isozyme of murine creatine kinase (CK-B), an ATP-buffering enzyme expressed mainly in muscle and brain but absent from liver, kidney, and pancreas. Gene expression was monitored by 31P magnetic resonance spectroscopy (MRS) using the product of the CK enzymatic reaction, phosphocreatine, as an indicator of transfection. The vector was administered into nude mice by tail vein injection, and exogenous creatine was administered in the drinking water and by i.p. injection of 2% creatine solution before 31P MRS examination, which was performed on surgically exposed livers. A phosphocreatine resonance was detected in livers of mice injected with the vector and was absent from livers of control animals. CK expression was confirmed in the injected animals by Western blot analysis, enzymatic assays, and immunofluorescence measurements. We conclude that the syngeneic enzyme CK can be used as a marker gene for in vivo monitoring of gene expression after virally mediated gene transfer to the liver.

Estrogen and thyroid hormone interaction on regulation of gene expression.

Estrogen receptor (ER) and thyroid hormone receptors (TRs) are ligand-dependent nuclear transcription factors that can bind to an identical half-site, AGGTCA, of their cognate hormone response elements. By in vitro transfection analysis in CV-1 cells, we show that estrogen induction of chloramphenicol acetyltransferase (CAT) activity in a construct containing a CAT reporter gene under the control of a minimal thymidine kinase (tk) promoter and a copy of the consensus ER response element was attenuated by cotransfection of TR alpha 1 plus triiodothyronine treatment. This inhibitory effect of TR was ligand-dependent and isoform-specific. Neither TR beta 1 nor TR beta 2 cotransfection inhibited estrogen-induced CAT activity, although both TR alpha and TR beta can bind to a consensus ER response element. Furthermore, cotransfection of a mutated TR alpha 1 that lacks binding to the AGGTCA sequence also inhibited the estrogen effect. Thus, the repression of estrogen action by liganded TR alpha 1 may involve protein-protein interactions although competition of ER and TR at the DNA level cannot be excluded. A similar inhibitory effect of liganded TR alpha 1 on estrogen induction of CAT activity was observed in a construct containing the preproenkephalin (PPE) promoter. A study in hypophysectomized female rats demonstrated that the estrogen-induced increase in PPE mRNA levels in the ventromedial hypothalamus was diminished by coadministration of triiodothyronine. These results suggest that ER and TR may interact to modulate estrogen-sensitive gene expression, such as for PPE, in the hypothalamus.

Artificial regulation of gene expression in Escherichia coli by RNase P.

Plasmids encoding various external guide sequences (EGSs) were constructed and inserted into Escherichia coli. In strains harboring the appropriate plasmids, the expression of fully induced beta-galactosidase and alkaline phosphatase activity was reduced by more than 50%, while no reduction in such activity was observed in strains with non-specific EGSs. The inhibition of gene expression was virtually abolished at restrictive temperatures in strains that were temperature-sensitive for RNase P (EC 3.1.26.5). Northern blot analysis showed that the steady-state copy number of EGS RNAs was several hundred per cell in vivo. A plasmid that contained a gene for M1 RNA covalently linked to a specific EGS reduced the level of expression of a suppressor tRNA that was encoded by a separate plasmid. Similar methods can be used to regulate gene expression in E. coli and to mimic the properties of cold-sensitive mutants.

Primary structure of hepatocyte nuclear factor/forkhead homologue 4 and characterization of gene expression in the developing respiratory and reproductive epithelium.

Members of the winged helix/forkhead family of transcription factors are believed to play a role in cell-specific gene expression. A cDNA encoding a member of this family of proteins, termed hepatocyte nuclear factor/forkhead homologue 4 (HFH-4), has been isolated from rat lung and rat testis cDNA libraries. This cDNA contains an open reading frame of 421 amino acids with a conserved DNA binding domain and several potential transactivating regions. During murine lung development, a single species of HFH-4-specific transcript (2.4 kb long) is first detected precisely at the start of the late pseudoglandular stage (embryonic day 14.5) and, by in situ hybridization, is specifically localized to the proximal pulmonary epithelium. The unique temporal and spatial pattern of HFH-4 gene expression in the developing lung defines this protein as a marker for the initiation of bronchial epithelial cell differentiation and suggests that it may play an important role in cell fate determination during lung development. In addition to expression in the pulmonary epithelium, RNA blot analysis reveals 2.4-kb HFH-4 transcripts in the testis and oviduct. By using mice with genetic defects in spermatogenesis, HFH-4 expression in the testis is found to be associated with the appearance of haploid germ cells and in situ hybridization studies indicate that HFH-4 expression is confined to stages I-VII of spermatogenesis. This pattern of HFH-4 gene expression during the early stages of differentiation of haploid germ cells suggests that HFH-4 may play a role in regulating stage-specific gene expression and cell-fate determination during lung development and in spermatogenesis.

Inhibition of gene expression in human cells through small molecule-RNA interactions

Small molecules that bind their biological receptors with high affinity and selectivity can be isolated from randomized pools of combinatorial libraries. RNA-protein interactions are important in many cellular functions, including transcription, RNA splicing, and translation. One example of such interactions is the mechanism of trans-activation of HIV-1 gene expression that requires the interaction of Tat protein with the trans-activation responsive region (TAR) RNA, a 59-base stem-loop structure located at the 5′ end of all nascent HIV-1 transcripts. Here we demonstrate the isolation of small TAR RNA-binding molecules from an encoded combinatorial library. We have made an encoded combinatorial tripeptide library of 24,389 possible members from d-and l-alpha amino acids on TentaGel resin. Using on-bead screening we have identified a small family of mostly heterochiral tripeptides capable of structure-specific binding to the bulge loop of TAR RNA. In vitro binding studies reveal stereospecific discrimination when the best tripeptide ligand is compared with diastereomeric peptide sequences. In addition, the most strongly binding tripeptide was shown to suppress transcriptional activation by Tat protein in human cells with an IC50 of ≈50 nM. Our results indicate that tripeptide RNA ligands are cell permeable, nontoxic to cells, and capable of inhibiting expression of specific genes by interfering with RNA-protein interactions.

Long-range disruption of gene expression by a selectable marker cassette

Recent studies have suggested that the retention of selectable marker cassettes (like PGK–Neo, in which a hybrid gene consisting of the phosphoglycerate kinase I promoter drives the neomycin phosphotransferase gene) in targeted loci can cause unexpected phenotypes in “knockout” mice due to disruption of expression of neighboring genes within a locus. We have studied targeted mutations in two multigene clusters, the granzyme B locus and the β-like globin gene cluster. The insertion of PGK–Neo into the granzyme B gene, the most 5′ gene in the granzyme B gene cluster, severely reduced the normal expression of multiple genes within the locus, even at distances greater than 100 kb from the mutation. Similarly, the insertion of a PGK–Neo cassette into the β-globin locus control region (LCR) abrogates the expression of multiple globin genes downstream from the cassette. In contrast, a targeted mutation of the promyelocyte-specific cathepsin G gene (which lies just 3′ to the granzyme genes in the same cluster) had minimal effects on upstream granzyme gene expression. Although the mechanism of these long distance effects are unknown, the expression of PGK–Neo can be “captured” by the regulatory domain into which it is inserted. These results suggest that the PGK–Neo cassette can interact productively with locus control regions and thereby disrupt normal interactions between local and long-distance regulatory regions within a tissue-specific domain.