Dynamic modeling of gene expression data

We describe the time evolution of gene expression levels by using a time translational matrix to predict future expression levels of genes based on their expression levels at some initial time. We deduce the time translational matrix for previously published DNA microarray gene expression data sets by modeling them within a linear framework by using the characteristic modes obtained by singular value decomposition. The resulting time translation matrix provides a measure of the relationships among the modes and governs their time evolution. We show that a truncated matrix linking just a few modes is a good approximation of the full time translation matrix. This finding suggests that the number of essential connections among the genes is small.

Mechanism for fetal globin gene expression: Role of the soluble guanylate cyclase–cGMP-dependent protein kinase pathway

Despite considerable concerns with pharmacological stimulation of fetal hemoglobin (Hb F) as a therapeutic option for the β-globin disorders, the molecular basis of action of Hb F-inducing agents remains unclear. Here we show that an intracellular pathway including soluble guanylate cyclase (sGC) and cGMP-dependent protein kinase (PKG) plays a role in induced expression of the γ-globin gene. sGC, an obligate heterodimer of α- and β-subunits, participates in a variety of physiological processes by converting GTP to cGMP. Northern blot analyses with erythroid cell lines expressing different β-like globin genes showed that, whereas the β-subunit is expressed at similar levels, high-level expression of the α-subunit is preferentially observed in erythroid cells expressing γ-globin but not those expressing β-globin. Also, the levels of expression of the γ-globin gene correlate to those of the α-subunit. sGC activators or cGMP analogs increased expression of the γ-globin gene in erythroleukemic cells as well as in primary erythroblasts from normal subjects and patients with β-thalassemia. Nuclear run-off assays showed that the sGC activator protoporphyrin IX stimulates transcription of the γ-globin gene. Furthermore, increased expression of the γ-globin gene by well known Hb F-inducers such as hemin and butyrate was abolished by inhibiting sGC or PKG activity. Taken together, these results strongly suggest that the sGC–PKG pathway constitutes a mechanism that regulates expression of the γ-globin gene. Further characterization of this pathway should permit us to develop new therapeutics for the β-globin disorders.

The effects of aging on gene expression in the hypothalamus and cortex of mice

A better understanding of the molecular effects of aging in the brain may help to reveal important aspects of organismal aging, as well as processes that lead to age-related brain dysfunction. In this study, we have examined differences in gene expression in the hypothalamus and cortex of young and aged mice by using high-density oligonucleotide arrays. A number of key genes involved in neuronal structure and signaling are differentially expressed in both the aged hypothalamus and cortex, including synaptotagmin I, cAMP-dependent protein kinase C β, apolipoprotein E, protein phosphatase 2A, and prostaglandin D. Misregulation of these proteins may contribute to age-related memory deficits and neurodegenerative diseases. In addition, many proteases that play essential roles in regulating neuropeptide metabolism, amyloid precursor protein processing, and neuronal apoptosis are up-regulated in the aged brain and likely contribute significantly to brain aging. Finally, a subset of these genes whose expression is affected by aging are oppositely affected by exposure of mice to an enriched environment, suggesting that these genes may play important roles in learning and memory.

Hypoxia-induced gene expression profiling in the euryoxic fish Gillichthys mirabilis

Hypoxia is important in both biomedical and environmental contexts and necessitates rapid adaptive changes in metabolic organization. Mammals, as air breathers, have a limited capacity to withstand sustained exposure to hypoxia. By contrast, some aquatic animals, such as certain fishes, are routinely exposed and resistant to severe environmental hypoxia. Understanding the changes in gene expression in fishes exposed to hypoxic stress could reveal novel mechanisms of tolerance that may shed new light on hypoxia and ischemia in higher vertebrates. Using cDNA microarrays, we have studied gene expression in a hypoxia-tolerant burrow-dwelling goby fish, Gillichthys mirabilis. We show that a coherent picture of a complex transcriptional response can be generated for a nonmodel organism for which sequence data were unavailable. We demonstrate that: (i) although certain shifts in gene expression mirror changes in mammals, novel genes are differentially expressed in fish; and (ii) tissue-specific patterns of expression reflect the different metabolic roles of tissues during hypoxia.

Plastid redox state and sugars: Interactive regulators of nuclear-encoded photosynthetic gene expression

Feedback regulation of photosynthesis by carbon metabolites has long been recognized, but the underlying cellular mechanisms that control this process remain unclear. By using an Arabidopsis cell culture, we show that a block in photosynthetic electron flux prevents the increase in transcript levels of chlorophyll a/b-binding protein and the small subunit of Rubisco that typically occurs when intracellular sugar levels are depleted. In contrast, the expression of the nitrate reductase gene, which is induced by sugars, is not affected. These findings were confirmed in planta by using Arabidopsis carrying the firefly luciferase reporter gene fused to the plastocyanin and chlorophyll a/b-binding protein 2 gene promoters. Transcription from both promoters increases on carbohydrate depletion. Blocking photosynthetic electron transport with 3-(3′, 4′-dichlorophenyl)-1,1′-dimethylurea prevents this increase in transcription. We conclude that plastid-derived redox signaling can override the sugar-regulated expression of nuclear-encoded photosynthetic genes. In the sugar-response mutant, sucrose uncoupled 6 (sun6), plastocyanin-firefly luciferase transcription actually increases in response to exogenous sucrose rather than decreasing as in the wild type. Interestingly, plastid-derived redox signals do not influence this defective pattern of sugar-regulated gene expression in the sun6 mutant. A model, which invokes a positive inducer originating from the photosynthetic electron transport chain, is proposed to explain the nature of the plastid-derived signal.

Regulation of c-maf gene expression by Pax6 in cultured cells

c-Maf is a bZip transcription factor expressed in developmental and cellular differentiation processes. Recently, a c-maf knockout mouse model, showing abnormal lens development, has been reported. In order to study the regulation mechanisms of c-maf gene expression during the differentiation process we have cloned and functionally characterized the rat c-maf (maf-2) gene. The rat c-maf gene is an intronless gene, covering a length of 3.5 kb. Transient transfection analysis of the 5′-flanking region of the c-maf gene using luciferase as the reporter gene shows that Pax6, a master transcription factor for lens development, strongly activates the c-maf promoter construct. Endogenous c-maf is also activated by the Pax6 expression vector. Electrophoresis mobility shift assay and DNase I footprinting analysis show that at least three Pax6-binding sites are located in the 5′-flanking and 5′-non-coding regions of the rat c-maf gene. The c-maf gene was also markedly activated by its own product, c-Maf, through the MARE (Maf recognition element), suggesting that a positive autoregulatory mechanism controls this gene. In situ hybridization histochemical detection of Pax6 and c-Maf in the E14 lens showed that both mRNAs are expressed in the lens equator where lens epithelial cells are differentiating to lens fiber cells. These results suggest that a Pax6/c-Maf transcription factor cascade is working in lens development.

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

The TRANSFAC system on gene expression regulation

The TRANSFAC database on transcription factors and their DNA-binding sites and profiles (http://www.gene-regulation.de/) has been quantitatively extended and supplemented by a number of modules. These modules give information about pathologically relevant mutations in regulatory regions and transcription factor genes (PathoDB), scaffold/matrix attached regions (S/MARt DB), signal transduction (TRANSPATH) and gene expression sources (CYTOMER). Altogether, these distinct database modules constitute the TRANSFAC system. They are accompanied by a number of program routines for identifying potential transcription factor binding sites or for localizing individual components in the regulatory network of a cell.

The Mouse Gene Expression Database (GXD)

The Gene Expression Database (GXD) is a community resource of gene expression information for the laboratory mouse. By combining the different types of expression data, GXD aims to provide increasingly complete information about the expression profiles of genes in different mouse strains and mutants, thus enabling valuable insights into the molecular networks that underlie normal development and disease. GXD is integrated with the Mouse Genome Database (MGD). Extensive interconnections with sequence databases and with databases from other species, and the development and use of shared controlled vocabularies extend GXD’s utility for the analysis of gene expression information. GXD is accessible through the Mouse Genome Informatics web site at http://www.informatic s.jax.org/ or directly at http://www.informatics.jax.org/me nus/expression_menu.shtml.

RECODE: a database of frameshifting, bypassing and codon redefinition utilized for gene expression

The RECODE database is a compilation of ‘programmed’ translational recoding events taken from the scientific literature and personal communications. The database deals with programmed ribosomal frameshifting, codon redefinition and translational bypass occurring in a variety of organisms. The entries for each event include the sequences of the corresponding genes, their encoded proteins for both the normal and alternate decoding, the types of the recoding events involved, trans-factors and cis-elements that influence recoding. The database is freely available at http://recode.genetics.utah.edu/.

Fatty acids and hypolipidemic drugs regulate peroxisome proliferator-activated receptors α- and γ-mediated gene expression via liver fatty acid binding protein: A signaling path to the nucleus

Peroxisome proliferator-activated receptor α (PPARα) is a key regulator of lipid homeostasis in hepatocytes and target for fatty acids and hypolipidemic drugs. How these signaling molecules reach the nuclear receptor is not known; however, similarities in ligand specificity suggest the liver fatty acid binding protein (L-FABP) as a possible candidate. In localization studies using laser-scanning microscopy, we show that L-FABP and PPARα colocalize in the nucleus of mouse primary hepatocytes. Furthermore, we demonstrate by pull-down assay and immunocoprecipitation that L-FABP interacts directly with PPARα. In a cell biological approach with the aid of a mammalian two-hybrid system, we provide evidence that L-FABP interacts with PPARα and PPARγ but not with PPARβ and retinoid X receptor-α by protein–protein contacts. In addition, we demonstrate that the observed interaction of both proteins is independent of ligand binding. Final and quantitative proof for L-FABP mediation was obtained in transactivation assays upon incubation of transiently and stably transfected HepG2 cells with saturated, monounsaturated, and polyunsaturated fatty acids as well as with hypolipidemic drugs. With all ligands applied, we observed strict correlation of PPARα and PPARγ transactivation with intracellular concentrations of L-FABP. This correlation constitutes a nucleus-directed signaling by fatty acids and hypolipidemic drugs where L-FABP acts as a cytosolic gateway for these PPARα and PPARγ agonists. Thus, L-FABP and the respective PPARs could serve as targets for nutrients and drugs to affect expression of PPAR-sensitive genes.

Identification and dissection of an enhancer controlling epithelial gene expression in skin

Keratins 14 and 5 are the structural hallmarks of the basal keratinocytes of the epidermis and outer root sheath (ORS) of the hair follicle. Their genes are controlled in a tissue-specific manner and thus serve as useful tools to elucidate the regulatory mechanisms involved in keratinocyte-specific transcription. Previously we identified several keratinocyte-specific DNase I hypersensitive sites (HSs) in the 5′ regulatory sequences of the K14 gene and showed that a 700-bp regulatory domain encompassing HSs II and III can confer epidermal and ORS-specific gene expression in transgenic mice in vivo. Although HS II harbored much of the transactivation activity in vitro, it was not sufficient to restrict expression to keratinocytes in vivo. We now explore the HS III regulatory element. Surprisingly, this element on its own confers gene expression to the keratinocytes of the inner root sheath (IRS) of the hair follicle, whereas a 275-bp DNA fragment containing both HSs II and III shifts the expression from the IRS to the basal keratinocytes and ORS in vivo. Electrophoretic mobility-shift assays and mutational studies of HSs III reveal a role for CACCC-box binding proteins, Sp1 family members, and other factors adding to the list of previously described factors that are involved in keratinocyte-specific gene expression. These studies highlight a cooperative interaction of the two HSs domains and strengthen the importance of combinatorial play of transcription factors that govern keratinocyte-specific gene regulation.

An efficient system for conditional gene expression in embryonic stem cells and in their in vitro and in vivo differentiated derivatives

We have developed a universally applicable system for conditional gene expression in embryonic stem (ES) cells that relies on tamoxifen-dependent Cre recombinase-loxP site-mediated recombination and bicistronic gene-trap expression vectors that allow transgene expression from endogenous cellular promoters. Two vectors were introduced into the genome of recipient ES cells, successively: (i) a bicistronic gene-trap vector encoding the β-galactosidase/neoR fusion protein and the Cre-ERT2 (Cre recombinase fused to a mutated ligand-binding domain of the human estrogen receptor) and (ii) a bicistronic gene-trap vector encoding the hygroR protein and the human alkaline phosphatase (hAP), the expression of which is prevented by tandemly repeated stop-of-transcription sequences flanked by loxP sites. In selected clones, hAP expression was shown to be regulated accurately by 4′hydroxy-tamoxifen. Strict hormone-dependent expression of hAP was achieved (i) in vitro in undifferentiated ES cells and embryoid bodies, (ii) in vivo in virtually all the tissues of the 10-day-old chimeric fetus (after injection of 4′hydroxy-tamoxifen to foster mothers), and (iii) ex vivo in primary embryonic fibroblasts isolated from chimeric fetuses. Therefore, this approach can be applied to drive conditional expression of virtually any transgene in a large variety of cell types, both in vitro and in vivo.

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.

Effect of rabies virus infection on gene expression in mouse brain

A variety of molecular genetic approaches were used to study the effect of rabies virus (RV) infection on host gene expression in mouse brain. The down-regulation of gene expression was found to be a major effect of RV infection by using subtraction hybridization. However, a combination of techniques identified approximately 39 genes activated by infection. These included genes involved in regulation of cell metabolism, protein synthesis, synaptic activity, and cell growth and differentiation. Northern blot analysis to monitor temporal activation of several of these genes following infection revealed essentially two patterns of activation: (i) an early response with up-regulation beginning within 3 days after infection and correlating with transcription of RV nuclear protein; and (ii) a late response with enhanced expression occurring at days 6–7 after infection and associated with peak RV replication. The gene activation patterns and the known functions of their products suggest that a number of host genes may be involved in the replication and spread of RV in the brain.