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.

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.

Metallothioneins 1 and 2 Have Distinct but Overlapping Expression Patterns in Arabidopsis1

The spatial and temporal expression patterns of metallothionein (MT) isoforms MT1a and MT2a were investigated in vegetative and reproductive tissues of untreated and copper-treated Arabidopsis by in situ hybridization and by northern blotting. In control plants, MT1a mRNA was localized in leaf trichomes and in the vascular tissue in leaves, roots, flowers, and germinating embryos. In copper-treated plants, MT1a expression was also observed in the leaf mesophyll and in vascular tissue of developing siliques and seeds. In contrast, MT2a was expressed primarily in the trichomes of both untreated and copper-treated plants. In copper-treated plants, MT2a mRNA was also expressed in siliques. Northern-hybridization studies performed on developing seedlings and leaves showed temporal variations of MT1a gene expression but not of MT2a expression. The possible implications of these findings for the cellular roles of MTs in plants are discussed.

Expression Patterns Conferred by Tyrosine/Dihydroxyphenylalanine Decarboxylase Promoters from Opium Poppy Are Conserved in Transgenic Tobacco1

Opium poppy (Papaver somniferum) contains a large family of tyrosine/dihydroxyphenylalanine decarboxylase (tydc) genes involved in the biosynthesis of benzylisoquinoline alkaloids and cell wall-bound hydroxycinnamic acid amides. Eight members from two distinct gene subfamilies have been isolated, tydc1, tydc4, tydc6, tydc8, and tydc9 in one group and tydc2, tydc3, and tydc7 in the other. The tydc8 and tydc9 genes were located 3.2 kb apart on one genomic clone, suggesting that the family is clustered. Transcripts for most tydc genes were detected only in roots. Only tydc2 and tydc7 revealed expression in both roots and shoots, and TYDC3 mRNAs were the only specific transcripts detected in seedlings. TYDC1, TYDC8, and TYDC9 mRNAs, which occurred in roots, were not detected in elicitor-treated opium poppy cultures. Expression of tydc4, which contains a premature termination codon, was not detected under any conditions. Five tydc promoters were fused to the β-glucuronidase (GUS) reporter gene in a binary vector. All constructs produced transient GUS activity in microprojectile-bombarded opium poppy and tobacco (Nicotiana tabacum) cell cultures. The organ- and tissue-specific expression pattern of tydc promoter-GUS fusions in transgenic tobacco was generally parallel to that of corresponding tydc genes in opium poppy. GUS expression was most abundant in the internal phloem of shoot organs and in the stele of roots. Select tydc promoter-GUS fusions were also wound induced in transgenic tobacco, suggesting that the basic mechanisms of developmental and inducible tydc regulation are conserved across plant species.

Human semaphorins A(V) and IV reside in the 3p21.3 small cell lung cancer deletion region and demonstrate distinct expression patterns.

Semaphorins and collapsins make up a family of conserved genes that encode nerve growth cone guidance signals. We have identified two additional members of the human semaphorin family [human semaphorin A(V) and human semaphorin IV] in chromosome region 3p21.3, where several small cell lung cancer (SCLC) cell lines exhibit homozygous deletions indicative of a tumor suppressor gene. Human semaphorin A(V) has 86% amino acid homology with murine semaphorin A, whereas semaphorin IV is most closely related to murine semaphorin E, with 50% homology. These semaphorin genes are approximately 70 kb apart flanking two GTP-binding protein genes, GNAI-2 and GNAT-1. In contrast, other human semaphorin gene sequences (human semaphorin III and homologues of murine semaphorins B and C) are not located on chromosome 3. Human semaphorin A(V) is translated in vitro into a 90-kDa protein, which accumulates at the endoplasmic reticulum. The human semaphorin A(V) (3.4-kb mRNA) and IV (3.9- and 2.9-kb mRNAs) genes are expressed abundantly but differentially in a variety of human neural and nonneural tissues. Human semaphorin A(V) was expressed in only 1 out of 23 SCLCs and 7 out of 16 non-SCLCs, whereas semaphorin IV was expressed in 19 out of 23 SCLCs and 13 out of 16 non-SCLCs. Mutational analysis in semaphorin A(V) revealed mutations (germ line in one case) in 3 of 40 lung cancers. Our data suggest the need to determine the function of human semaphorins A(V) and IV in nonneural tissues and their role in the pathogenesis of lung cancer.

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.

HLA-G expression in melanoma: A way for tumor cells to escape from immunosurveillance

Considering the well established role of nonclassical HLA-G class I molecules in inhibiting natural killer (NK) cell function, the consequence of abnormal HLA-G expression in malignant cells should be the escape of tumors from immunosurveillance. To examine this hypothesis, we analyzed HLA-G expression and NK sensitivity in human malignant melanoma cells. Our analysis of three melanoma cell lines and ex vivo biopsy demonstrated that (i) IGR and M74 human melanoma cell lines exhibit a high level of HLA-G transcription with differential HLA-G isoform transcription and protein expression patterns, (ii) a higher level of HLA-G transcription ex vivo is detected in a skin melanoma metastasis biopsy compared with a healthy skin fragment from the same individual, and (iii) HLA-G protein isoforms other than membrane-bound HLA-G1 protect IGR from NK lysis. It thus appears of critical importance to consider the specific role of HLA-G expression in tumors in the design of future cancer immunotherapies.

The circadian clock controls the expression pattern of the circadian input photoreceptor, phytochrome B

Developmental and physiological responses are regulated by light throughout the entire life cycle of higher plants. To sense changes in the light environment, plants have developed various photoreceptors, including the red/far-red light-absorbing phytochromes and blue light-absorbing cryptochromes. A wide variety of physiological responses, including most light responses, also are modulated by circadian rhythms that are generated by an endogenous oscillator, the circadian clock. To provide information on local time, circadian clocks are synchronized and entrained by environmental time cues, of which light is among the most important. Light-driven entrainment of the Arabidopsis circadian clock has been shown to be mediated by phytochrome A (phyA), phytochrome B (phyB), and cryptochromes 1 and 2, thus affirming the roles of these photoreceptors as input regulators to the plant circadian clock. Here we show that the expression of PHYB∷LUC reporter genes containing the promoter and 5′ untranslated region of the tobacco NtPHYB1 or Arabidopsis AtPHYB genes fused to the luciferase (LUC) gene exhibit robust circadian oscillations in transgenic plants. We demonstrate that the abundance of PHYB RNA retains this circadian regulation and use a PHYB∷Luc fusion protein to show that the rate of PHYB synthesis is also rhythmic. The abundance of bulk PHYB protein, however, exhibits only weak circadian rhythmicity, if any. These data suggest that photoreceptor gene expression patterns may be significant in the daily regulation of plant physiology and indicate an unexpectedly intimate relationship between the components of the input pathway and the putative circadian clock mechanism in higher plants.

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.

Conservation of the expression and function of apterous orthologs in Drosophila and mammals

The Drosophila apterous (ap) gene encodes a protein of the LIM-homeodomain family. Many transcription factors of this class have been conserved during evolution; however, the functional significance of their structural conservation is generally not known. ap is best known for its fundamental role as a dorsal selector gene required for patterning and growth of the wing, but it also has other important functions required for neuronal fasciculation, fertility, and normal viability. We isolated mouse (mLhx2) and human (hLhx2) ap orthologs, and we used transgenic animals and rescue assays to investigate the conservation of the Ap protein during evolution. We found that the human protein LHX2 is able to regulate correctly ap target genes in the fly, causes the same phenotypes as Ap when ectopically produced, and most importantly rescues ap mutant phenotypes as efficiently as the fly protein. In addition, we found striking similarities in the expression patterns of the Drosophila and murine genes. Both mLhx2 and ap are expressed in the respective nerve cords, eyes, olfactory organs, brain, and limbs. These results demonstrate the conservation of Ap protein function across phyla and argue that aspects of its expression pattern have also been conserved from a common ancestor of insects and vertebrates.

Differential expression of dystrophin isoforms and utrophin during dibutyryl-cAMP-induced morphological differentiation of rat brain astrocytes

We have identified isoforms of dystrophin and utrophin, a dystrophin homologue, expressed in astrocytes and examined their expression patterns during dibutyryl-cAMP (dBcAMP)-induced morphological differentiation of astrocytes. Immunoblot and immunocytochemical analyses showed that full-length-type dystrophin (427 kDa), utrophin (395 kDa), and Dp71 (75 kDa), a small-type dystrophin isoform, were coexpressed in cultured nondifferentiated rat brain astrocytes and were found to be located in the cell membrane. During morphological differentiation of the astrocytes induced by 1 mM dBcAMP, the amount of Dp71 markedly increased, whereas that of dystrophin and utrophin decreased. Northern blot analyses revealed that dBcAMP regulates the mRNA levels of Dp71 and dystrophin but not that of utrophin. dBcAMP slightly increased the amount of the β-dystroglycan responsible for anchoring dystrophin isoforms and utrophin to the cell membrane. Immunocytochemical analyses showed that most utrophin was observed in the cytoplasmic area during astrocyte differentiation, whereas Dp71 was found along the cell membrane of the differentiated astrocytes. These findings suggest that most of the dystrophin/utrophin-dystroglycan complex on cell membrane in cultured astrocytes was replaced by the Dp71-dystroglycan complex during morphological differentiation. The cell biological roles of Dp71 are discussed.

Expression of homeobox genes shows chelicerate arthropods retain their deutocerebral segment

Expression patterns of six homeobox containing genes in a model chelicerate, the oribatid mite Archegozetes longisetosus, were examined to establish homology of chelicerate and insect head segments and to investigate claims that the chelicerate deutocerebral segment has been reduced or lost. engrailed (en) expression, which has been used to demonstrate the presence of segments in insects, fails to demonstrate a reduced deutocerebral segment. Expression patterns of the chelicerate homologs of the Drosophila genes Antennapedia (Antp), Sex combs reduced (Scr), Deformed (Dfd), proboscipedia (pb), and orthodenticle (otd) confirm direct correspondence of head segments. The chelicerate deutocerebral segment has not been reduced or lost. We make further inferences concerning the evolution of heads and Hox genes in arthropods.

BodyMap incorporated PCR-based expression profiling data and a gene ranking system

BodyMap is a human and mouse gene expression database that is based on site-directed 3′-expressed sequence tags generated at Osaka University. To date, it contains more than 300 000 tag sequences from 64 human and 39 mouse tissues. For the recent release, the precise anatomical expression patterns for more than half of the human gene entries were generated by introduced amplified fragment length polymorphism (iAFLP), which is a PCR-based high-throughput expression profiling method. The iAFLP data incorporated into BodyMap describe the relative contents of more than 12 000 transcripts across 30 tissue RNAs. In addition, a newly developed gene ranking system helps users obtain lists of genes that have desired expression patterns according to their significance. BodyMap supports complete transfer of unique data sets and provides analysis that is accessible through the WWW at http://bodymap.ims.u-tokyo.ac.jp.

Delineating developmental and metabolic pathways in vivo by expression profiling using the RIKEN set of 18,816 full-length enriched mouse cDNA arrays

We have systematically characterized gene expression patterns in 49 adult and embryonic mouse tissues by using cDNA microarrays with 18,816 mouse cDNAs. Cluster analysis defined sets of genes that were expressed ubiquitously or in similar groups of tissues such as digestive organs and muscle. Clustering of expression profiles was observed in embryonic brain, postnatal cerebellum, and adult olfactory bulb, reflecting similarities in neurogenesis and remodeling. Finally, clustering genes coding for known enzymes into 78 metabolic pathways revealed a surprising coordination of expression within each pathway among different tissues. On the other hand, a more detailed examination of glycolysis revealed tissue-specific differences in profiles of key regulatory enzymes. Thus, by surveying global gene expression by using microarrays with a large number of elements, we provide insights into the commonality and diversity of pathways responsible for the development and maintenance of the mammalian body plan.