Changes in global gene expression patterns during development and maturation of the rat kidney

We set out to define patterns of gene expression during kidney organogenesis by using high-density DNA array technology. Expression analysis of 8,740 rat genes revealed five discrete patterns or groups of gene expression during nephrogenesis. Group 1 consisted of genes with very high expression in the early embryonic kidney, many with roles in protein translation and DNA replication. Group 2 consisted of genes that peaked in midembryogenesis and contained many transcripts specifying proteins of the extracellular matrix. Many additional transcripts allied with groups 1 and 2 had known or proposed roles in kidney development and included LIM1, POD1, GFRA1, WT1, BCL2, Homeobox protein A11, timeless, pleiotrophin, HGF, HNF3, BMP4, TGF-α, TGF-β2, IGF-II, met, FGF7, BMP4, and ganglioside-GD3. Group 3 consisted of transcripts that peaked in the neonatal period and contained a number of retrotransposon RNAs. Group 4 contained genes that steadily increased in relative expression levels throughout development, including many genes involved in energy metabolism and transport. Group 5 consisted of genes with relatively low levels of expression throughout embryogenesis but with markedly higher levels in the adult kidney; this group included a heterogeneous mix of transporters, detoxification enzymes, and oxidative stress genes. The data suggest that the embryonic kidney is committed to cellular proliferation and morphogenesis early on, followed sequentially by extracellular matrix deposition and acquisition of markers of terminal differentiation. The neonatal burst of retrotransposon mRNA was unexpected and may play a role in a stress response associated with birth. Custom analytical tools were developed including “The Equalizer” and “eBlot,” which contain improved methods for data normalization, significance testing, and data mining.

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.

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.

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.

Patterns of Starchy Endosperm Acidification and Protease Gene Expression in Wheat Grains following Germination1

Cereal aleurone responses to gibberellic acid (GA3) include activation of synthesis of hydrolytic enzymes and acidification of the external medium. We have studied the effect of the pH of the incubation medium on the response of wheat (Triticum aestivum) aleurone cells to GA3. De-embryonated half grains show the capacity for GA3-activated medium acidification when incubation is carried out at pH 6.0 to 7.0 but not at lower pHs. In addition, the activating effect of GA3 on the expression of carboxypeptidase III and thiol protease genes is more efficient when the hormone treatment is carried out at neutral pH. In situ pH staining showed that starchy endosperm acidification takes place upon imbibition and advances from the embryo to the distal part of the grain. In situ hybridization experiments showed a similar pattern of expression of a carboxypeptidase III gene, which is up-regulated by GA3 in aleurone cells. However, aleurone gene expression precedes starchy endosperm acidification. These findings imply that in vivo GA perception by the aleurone layer takes place at neutral pH and suggest that the acidification of the starchy endosperm is regulated by GA3 in germinated wheat grains.

Embryonic expression patterns of the neural cell adhesion molecule gene are regulated by homeodomain binding sites.

During development of the vertebrate nervous system, the neural cell adhesion molecule (N-CAM) is expressed in a defined spatiotemporal pattern. We have proposed that the expression of N-CAM is controlled, in part, by proteins encoded by homeobox genes. This hypothesis has been supported by previous in vitro experiments showing that products of homeobox genes can both bind to and transactivate the N-CAM promoter via two homeodomain binding sites, HBS-I and HBS-II. We have now tested the hypothesis that the N-CAM gene is a target of homeodomain proteins in vivo by using transgenic mice containing native and mutated N-CAM promoter constructs linked to a beta-galactosidase reporter gene. Segments of the 5' flanking region of the mouse N-CAM gene were sufficient to direct expression of the reporter gene in the central nervous system in a pattern consistent with that of the endogenous N-CAM gene. For example, at embryonic day (E) 11, beta-galactosidase staining was found in postmitotic neurons in dorsolateral and ventrolateral regions of the spinal cord; at E14.5, staining was seen in these neurons throughout the spinal cord. In contrast, mice carrying an N-CAM promoter-reporter construct with mutations in both homeodomain binding sites (HBS-I and HBS-II) showed altered expression patterns in the spinal cord. At E11, beta-galactosidase expression was seen in the ventrolateral spinal cord, but was absent in the dorsolateral areas, and at E 14.5, beta-galactosidase expression was no longer detected in any cells of the cord. Homeodomain binding sites found in the N-CAM promoter thus appear to be important in determining specific expression patterns of N-CAM along the dorsoventral axis in the developing spinal cord. These experiments suggest that the N-CAM gene is an in vivo target of homeobox gene products in vertebrates.

Feedback-inducible nuclear-receptor-driven reporter gene expression in transgenic mice

Understanding nuclear receptor signaling in vivo would be facilitated by an efficient methodology to determine where a nuclear receptor is active. Herein, we present a feedback-inducible expression system in transgenic mice to detect activated nuclear receptor effector proteins by using an inducible reporter gene. With this approach, reporter gene induction is not limited to a particular tissue, and, thus, this approach provides the opportunity for whole-animal screens. Furthermore, the effector and reporter genes are combined to generate a single strain of transgenic mice, which enables direct and rapid analysis of the offspring. The system was applied to localize sites where the retinoic acid receptor ligand-binding domain is activated in vivo. The results identify previously discovered sources of retinoids in the embryo and indicate the existence of previously undiscovered regions of retinoic acid receptor signaling in vivo. Notably, the feedback-inducible nuclear-receptor-driven assay, combined with an independent in vitro assay, provides evidence for a site of retinoid synthesis in the isthmic mesenchyme. These data illustrate the potential of feedback-inducible nuclear-receptor-driven analyses for assessing in vivo activation patterns of nuclear receptors and for analyzing pharmacological properties of natural and synthetic ligands of potential therapeutic value.

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.

Cluster analysis and display of genome-wide expression patterns

A system of cluster analysis for genome-wide expression data from DNA microarray hybridization is described that uses standard statistical algorithms to arrange genes according to similarity in pattern of gene expression. The output is displayed graphically, conveying the clustering and the underlying expression data simultaneously in a form intuitive for biologists. We have found in the budding yeast Saccharomyces cerevisiae that clustering gene expression data groups together efficiently genes of known similar function, and we find a similar tendency in human data. Thus patterns seen in genome-wide expression experiments can be interpreted as indications of the status of cellular processes. Also, coexpression of genes of known function with poorly characterized or novel genes may provide a simple means of gaining leads to the functions of many genes for which information is not available currently.

Binding of factor VIIa to tissue factor induces alterations in gene expression in human fibroblast cells: Up-regulation of poly(A) polymerase

Tissue factor (TF) is the cellular receptor for an activated form of clotting factor VII (VIIa) and the binding of factor VII(a) to TF initiates the coagulation cascade. Sequence and structural patterns extracted from a global alignment of TF confers homology with interferon receptors of the cytokine receptor super family. Several recent studies suggested that TF could function as a genuine signal transducing receptor. However, it is unknown which biological function(s) of cells are altered upon the ligand, VIIa, binding to TF. In the present study, we examined the effect of VIIa binding to cell surface TF on cellular gene expression in fibroblasts. Differential mRNA display PCR technique was used to identify transcriptional changes in fibroblasts upon VIIa binding to TF. The display showed that VIIa binding to TF either up or down-regulated several mRNA species. The differential expression of one such transcript, VIIa-induced up-regulation, was confirmed by Northern blot analysis. Isolation of a full-length cDNA corresponding to the differentially expressed transcript revealed that VIIa-up-regulated gene was poly(A) polymerase. Northern blot analysis of various carcinomas and normal human tissues revealed an over expression of PAP in cancer tissues. Enhanced expression of PAP upon VIIa binding to tumor cell TF may potentially play an important role in tumor metastasis.

Modeling stochastic gene expression: Implications for haploinsufficiency

There is increasing recognition that stochastic processes regulate highly predictable patterns of gene expression in developing organisms, but the implications of stochastic gene expression for understanding haploinsufficiency remain largely unexplored. We have used simulations of stochastic gene expression to illustrate that gene copy number and expression deactivation rates are important variables in achieving predictable outcomes. In gene expression systems with non-zero expression deactivation rates, diploid systems had a higher probability of uninterrupted gene expression than haploid systems and were more successful at maintaining gene product above a very low threshold. Systems with relatively rapid expression deactivation rates (unstable gene expression) had more predictable responses to a gradient of inducer than systems with slow or zero expression deactivation rates (stable gene expression), and diploid systems were more predictable than haploid, with or without dosage compensation. We suggest that null mutations of a single allele in a diploid organism could decrease the probability of gene expression and present the hypothesis that some haploinsufficiency syndromes might result from an increased susceptibility to stochastic delays of gene initiation or interruptions of gene expression.

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.

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.