Fundamental patterns underlying gene expression profiles: Simplicity from complexity

Analysis of previously published sets of DNA microarray gene expression data by singular value decomposition has uncovered underlying patterns or “characteristic modes” in their temporal profiles. These patterns contribute unequally to the structure of the expression profiles. Moreover, the essential features of a given set of expression profiles are captured using just a small number of characteristic modes. This leads to the striking conclusion that the transcriptional response of a genome is orchestrated in a few fundamental patterns of gene expression change. These patterns are both simple and robust, dominating the alterations in expression of genes throughout the genome. Moreover, the characteristic modes of gene expression change in response to environmental perturbations are similar in such distant organisms as yeast and human cells. This analysis reveals simple regularities in the seemingly complex transcriptional transitions of diverse cells to new states, and these provide insights into the operation of the underlying genetic networks.

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.

Genome-wide gene expression profiles of the developing mouse hippocampus

We have analyzed the developmental molecular programs of the mouse hippocampus, a cortical structure critical for learning and memory, by means of large-scale DNA microarray techniques. Of 11,000 genes and expressed sequence tags examined, 1,926 showed dynamic changes during hippocampal development from embryonic day 16 to postnatal day 30. Gene-cluster analysis was used to group these genes into 16 distinct clusters with striking patterns that appear to correlate with major developmental hallmarks and cellular events. These include genes involved in neuronal proliferation, differentiation, and synapse formation. A complete list of the transcriptional changes has been compiled into a comprehensive gene profile database (http://BrainGenomics.Princeton.edu), which should prove valuable in advancing our understanding of the molecular and genetic programs underlying both the development and the functions of the mammalian brain.

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.

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.

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.

Temporal mapping of gene expression levels during the differentiation of individual primary hematopoietic cells

A hierarchical order of gene expression has been proposed to control developmental events in hematopoiesis, but direct demonstration of the temporal relationships between regulatory gene expression and differentiation has been difficult to achieve. We modified a single-cell PCR method to detect 2-fold changes in mRNA copies per cell (dynamic range, 250–250,000 copies/cell) and used it to sequentially quantitate gene expression levels as single primitive (CD34+,CD38−) progenitor cells underwent differentiation to become erythrocytes, granulocytes, or monocyte/macrophages. Markers of differentiation such as CD34 or cytokine receptor mRNAs and transcription factors associated with their regulation were assessed. All transcription factors tested were expressed in multipotent progenitors. During lineage-specific differentiation, however, distinct patterns of expression emerged. SCL, GATA-2, and GATA-1 expression sequentially extinguished during erythroid differentiation. PU.1, AML1B, and C/EBPα expression profiles and their relationship to cytokine receptor expression in maturing granulocytes could be distinguished from similar profiles in monocytic cells. These data characterize the dynamics of gene expression accompanying blood cell development and define a signature gene expression pattern for specific stages of hematopoietic differentiation.

Lipophosphoglycan from Leishmania suppresses agonist-induced interleukin 1β gene expression in human monocytes via a unique promoter sequence

Leishmania are parasites that survive within macrophages by mechanism(s) not entirely known. Depression of cellular immunity and diminished production of interleukin 1β (IL-1β) and tumor necrosis factor α are potential ways by which the parasite survives within macrophages. We examined the mechanism(s) by which lipophosphoglycan (LPG), a major glycolipid of Leishmania, perturbs cytokine gene expression. LPG treatment of THP-1 monocytes suppressed endotoxin induction of IL-1β steady-state mRNA by greater than 90%, while having no effect on the expression of a control gene. The addition of LPG 2 h before or 2 h after endotoxin challenge significantly suppressed steady-state IL-1β mRNA by 90% and 70%, respectively. LPG also inhibited tumor necrosis factor α and Staphylococcus induction of IL-1β gene expression. The inhibitory effect of LPG is agonist-specific because LPG did not suppress the induction of IL-1β mRNA by phorbol 12-myristate 13-acetate. A unique DNA sequence located within the −310 to −57 nucleotide region of the IL-1β promoter was found to mediate LPG’s inhibitory activity. The requirement for the −310 to −57 promoter gene sequence for LPG’s effect is demonstrated by the abrogation of LPG’s inhibitory activity by truncation or deletion of the −310 to −57 promoter gene sequence. Furthermore, the minimal IL-1β promoter (positions −310 to +15) mediated LPG’s inhibitory activity with dose and kinetic profiles that were similar to LPG’s suppression of steady-state IL-1β mRNA. These findings delineated a promoter gene sequence that responds to LPG to act as a “gene silencer,” a function, to our knowledge, not previously described. LPG’s inhibitory activity for several mediators of inflammation and the persistence of significant inhibitory activity 2 h after endotoxin challenge suggest that LPG has therapeutic potential and may be exploited for therapy of sepsis, acute respiratory distress syndrome, and autoimmune diseases.

Altered gene expression in the host brain caused by a trematode parasite: Neuropeptide genes are preferentially affected during parasitosis

Schistosome parasites adjust the physiology and behavior of their intermediate molluscan hosts to their own benefit. Previous studies demonstrated effects of the avian-schistosome Trichobilharzia ocellata on peptidergic centers in the brain of the intermediate snail host Lymnaea stagnalis. In particular, electrophysiological properties and peptide release of growth- and reproduction-controlling neuroendocrine neurons were affected. We now have examined the possibility that the expression of genes that control physiology and behavior of the host might be altered during parasitosis. A cDNA library of the brain of parasitized Lymnaea was constructed and differentially screened by using mRNA from the brain of both parasitized and nonparasitized snails. This screening yielded a number of clones, including previously identified cDNAs as well as novel neuronal transcripts, which appear to be differentially regulated. The majority of these transcripts encode neuropeptides. Reverse Northern blot analysis confirmed that neuropeptide gene expression is indeed affected in parasitized animals. Moreover, the expression profiles of 10 transcripts tested showed a differential, parasitic stage-specific regulation. Changes in expression could in many cases already be observed between 1.5 and 5 hr postinfection, suggesting that changes in gene expression are a direct effect of parasitosis. We suggest that direct regulation of neuropeptide gene expression is a strategy of parasites to induce physiological and behavioral changes in the host.

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.

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.

Regulation of T cell receptor (TCR) β gene expression by CD3 complex signaling in immature thymocytes: Implications for TCRβ allelic exclusion

During αβ thymocyte development, clonotype-independent CD3 complexes are expressed at the cell surface before the pre-T cell receptor (TCR). Signaling through clonotype-independent CD3 complexes is required for expression of rearranged TCRβ genes. On expression of a TCRβ polypeptide chain, the pre-TCR is assembled, and TCRβ locus allelic exclusion is established. We investigated the putative contribution of clonotype-independent CD3 complex signaling to TCRβ locus allelic exclusion in mice single-deficient or double-deficient for CD3ζ/η and/or p56lck. These mice display defects in the expression of endogenous TCRβ genes in immature thymocytes, proportional to the severity of CD3 complex malfunction. Exclusion of endogenous TCRβ VDJ (variable, diversity, joining) rearrangements by a functional TCRβ transgene was severely compromised in the single-deficient and double-deficient mutant mice. In contrast to wild-type mice, most of the CD25+ double-negative (DN) thymocytes of the mutant mice failed to express the TCRβ transgene, suggesting defective expression of the TCRβ transgene similar to endogenous TCRβ genes. In the mutant mice, a proportion of CD25+ DN thymocytes that failed to express the transgene expressed endogenous TCRβ polypeptide chains. Many double-positive cells of the mutant mice coexpressed endogenous and transgenic TCRβ chains or more than one endogenous TCRβ chain. The data suggest that signaling through clonotype-independent CD3 complexes may contribute to allelic exclusion of the TCRβ locus by inducing the expression of rearranged TCRβ genes in CD25+ DN thymocytes.

Induction of ARF tumor suppressor gene expression and cell cycle arrest by transcription factor DMP1

Expression of the DMP1 transcription factor, a cyclin D-binding Myb-like protein, induces growth arrest in mouse embryo fibroblast strains but is devoid of antiproliferative activity in primary diploid fibroblasts that lack the ARF tumor suppressor gene. DMP1 binds to a single canonical recognition site in the ARF promoter to activate gene expression, and in turn, p19ARF synthesis causes p53-dependent cell cycle arrest. Unlike genes such as Myc, adenovirus E1A, and E2F-1, which, when overexpressed, activate the ARF-p53 pathway and trigger apoptosis, DMP1, like ARF itself, does not induce programmed cell death. Therefore, apart from its recently recognized role in protecting cells from potentially oncogenic signals, ARF can be induced in response to antiproliferative stimuli that do not obligatorily lead to apoptosis.

Loss of the circadian clock-associated protein 1 in Arabidopsis results in altered clock-regulated gene expression

Little is known about plant circadian oscillators, in spite of how important they are to sessile plants, which require accurate timekeepers that enable the plants to respond to their environment. Previously, we identified a circadian clock-associated (CCA1) gene that encodes an Myb-related protein that is associated with phytochrome control and circadian regulation in plants. To understand the role CCA1 plays in phytochrome and circadian regulation, we have isolated an Arabidopsis line with a T DNA insertion that results in the loss of CCA1 RNA, of CCA1 protein, and of an Lhcb-promoter binding activity. This mutation affects the circadian expression of all four clock-controlled genes that we examined. The results show that, despite their similarity, CCA1 and LHY are only partially redundant. The lack of CCA1 also affects the phytochrome regulation of gene expression, suggesting that CCA1 has an additional role in a signal transduction pathway from light, possibly acting at the point of integration between phytochrome and the clock. Our results indicate that CCA1 is an important clock-associated protein involved in circadian regulation of gene expression.