TM4: a free, open-source system for microarray data management and analysis

Transcriptional program of early sporulation and stationary-phase events in Clostridium acetobutylicum

DNA microarray analysis of Clostridium acetobutylicum was used to examine the genomic-scale gene expression changes during the shift from exponential-phase growth and acidogenesis to stationary phase and solventogenesis. Self-organizing maps were used to identify novel expression patterns of functional gene classes, including aromatic and branched-chain amino acid synthesis, ribosomal proteins, cobalt and iron transporters, cobalamin biosynthesis, and lipid biosynthesis. The majority of pSOL1 megaplasmid genes (in addition to the solventogenic genes aad-ctfA-ctfB and adc) had increased expression at the onset of solventogenesis, suggesting that other megaplasmid genes may play a role in stationary-phase phenomena. Analysis of sporulation genes and comparison with published Bacillus subtilis results indicated conserved expression patterns of early sporulation genes, including spo0A, the sigF operon, and putative canonical genes of the sigma(H) and sigma(F) regulons. However, sigE expression could not be detected within 7.5 h of initial spo0A expression, consistent with the observed extended time between the appearance of clostridial forms and endospore formation. The results were compared with microarray comparisons of the wild-type strain and the nonsolventogenic, asporogenous M5 strain, which lacks the pSOL1 megaplasmid. While some results were similar, the expression of primary metabolism genes and heat shock proteins was higher in M5, suggesting a difference in metabolic regulation or a butyrate stress response in M5. The results of this microarray platform and analysis were further validated by comparing gene expression patterns to previously published Northern analyses, reporter assays, and two-dimensional protein electrophoresis data of metabolic genes (including all major solventogenesis genes), sporulation genes, heat shock proteins, and other solventogenesis-induced gene expression.

Array comparative genomic hybridisation (aCGH) analysis of premenopausal breast cancers from a nuclear fallout area and matched cases from Western New York

High-resolution array comparative genomic hybridisation (aCGH) analysis of DNA copy number aberrations (CNAs) was performed on breast carcinomas in premenopausal women from Western New York (WNY) and from Gomel, Belarus, an area exposed to fallout from the 1986 Chernobyl nuclear accident. Genomic DNA was isolated from 47 frozen tumour specimens from 42 patients and hybridised to arrays spotted with more than 3000 BAC clones. In all, 20 samples were from WNY and 27 were from Belarus. In total, 34 samples were primary tumours and 13 were lymph node metastases, including five matched pairs from Gomel. The average number of total CNAs per sample was 76 (range 35–134). We identified 152 CNAs (92 gains and 60 losses) occurring in more than 10% of the samples. The most common amplifications included gains at 8q13.2 (49%), at 1p21.1 (36%), and at 8q24.21 (36%). The most common deletions were at 1p36.22 (26%), at 17p13.2 (26%), and at 8p23.3 (23%). Belarussian tumours had more amplifications and fewer deletions than WNY breast cancers. HER2/neu negativity and younger age were also associated with a higher number of gains and fewer losses. In the five paired samples, we observed more discordant than concordant DNA changes. Unsupervised hierarchical cluster analysis revealed two distinct groups of tumours: one comprised predominantly of Belarussian carcinomas and the other largely consisting of WNY cases. In total, 50 CNAs occurred significantly more commonly in one cohort vs the other, and these included some candidate signature amplifications in the breast cancers in women exposed to significant radiation. In conclusion, our high-density aCGH study has revealed a large number of genetic aberrations in individual premenopausal breast cancer specimens, some of which had not been reported before. We identified a distinct CNA profile for carcinomas from a nuclear fallout area, suggesting a possible molecular fingerprint of radiation-associated breast cancer.

Cell cycle progression in the pericycle is not sufficient for SOLITARY ROOT/IAA14-mediated lateral root initiation in Arabidopsis thaliana

To study the mechanisms behind auxin-induced cell division, lateral root initiation was used as a model system. By means of microarray analysis, genome-wide transcriptional changes were monitored during the early steps of lateral root initiation. Inclusion of the dominant auxin signaling mutant solitary root1 (slr1) identified genes involved in lateral root initiation that act downstream of the auxin/indole-3-acetic acid (AUX/IAA) signaling pathway. Interestingly, key components of the cell cycle machinery were strongly defective in slr1, suggesting a direct link between AUX/IAA signaling and core cell cycle regulation. However, induction of the cell cycle in the mutant background by overexpression of the D-type cyclin (CYCD3;1) was able to trigger complete rounds of cell division in the pericycle that did not result in lateral root formation. Therefore, lateral root initiation can only take place when cell cycle activation is accompanied by cell fate respecification of pericycle cells. The microarray data also yielded evidence for the existence of both negative and positive feedback mechanisms that regulate auxin homeostasis and signal transduction in the pericycle, thereby fine-tuning the process of lateral root initiation.

Quality assessment of microarray experiments

OBJECTIVES

We describe a quality assurance procedure to maximize the value of oligonucleotide microarray expression profiles.

DESIGN, METHODS, AND RESULTS

Background microarray noise was 82.2+/-54.5 and 51.8+/-12.4 units, respectively, before and after enacting the program (P<0.0001). We also noted improved concordance of microarray expression fold-changes for selected genes with results of RT-PCR validation.

CONCLUSIONS

This multi-step procedure, including quantification of RNA sample degradation and detection of outlier data points, has increased data quality from our microarray facility.

Real-time expression profiling of microRNA precursors in human cancer cell lines

Our previous study described a real-time PCR method to quantify microRNA (miRNA) precursors using SYBR green detection [T. D. Schmittgen, J. Jiang, Q. Liu and L. Yang (2004) Nucleic Acids Res., 32, e43]. The present study adapted the assay to a 384-well format and expanded it to include primers to 222 human miRNA precursors. TaqMan minor groove binder probes were used to discriminate nearly identical members of the let-7 family of miRNA isoforms. The miRNA precursor expression was profiled in 32 human cell lines from lung, breast, colorectal, hematologic, prostate, pancreatic, and head and neck cancers. Some miRNA precursors were expressed at similar levels in many of the cell lines, while others were differentially expressed. Clustering analysis of the miRNA precursor expression data revealed that most of the cell lines clustered into their respective tissues from which each cell line was ostensibly derived. miRNA precursor expression by PCR paralleled the mature miRNA expression by northern blotting for most of the conditions studied. Our study provides PCR primer sequences to all of the known human miRNA precursors as of December 2004 and provides a database of the miRNA precursor expression in many commonly used human cancer cell lines.

Differential gene expression by Metarhizium anisopliae growing in root exudate and host (Manduca sexta) cuticle or hemolymph reveals mechanisms of physiological adaptation

Like many other fungal pathogens Metarhizium anisopliae is a facultative saprophyte with both soil-dwelling and insect pathogenic life-stages. In addition, as M. anisopliae traverses the cuticle and enters the hemolymph it must adapt to several different host environments. In this study, we used expressed sequence tags and cDNA microarray analyses to demonstrate that physiological adaptation by M. anisopliae to insect cuticle, insect hemolymph, bean root exudate (a model for life in the soil), and nutrient rich Sabouraud dextrose broth (SDB) involves different subsets of genes. Overall, expression patterns in cuticle and hemolymph clustered separately from expression patterns in root exudates and SDB, indicative of critical differences in transcriptional control during pathogenic and saprophytic growth. However, there were differences in gene expression between hemolymph and cuticle and these mostly involved perception mechanisms, carbon metabolism, proteolysis, cell surface properties, and synthesis of toxic metabolites. These differences suggest previously unsuspected stratagems of fungal pathogenicity that can be tested experimentally. Examples include the switch-off of cuticle-degrading proteases and a dramatic cell wall reorganization during growth in hemolymph.

Gene expression profiles relate to SS18/SSX fusion type in synovial sarcoma

We applied 27k spotted cDNA microarray slides to assess gene expression profiles in 26 samples from 24 patients with synovial sarcomas (SS). The data were analyzed in relation to histopathologic type, cytogenetic aberrations, gene fusion type and development of distant metastases. Supervised analysis based on gene fusion type in 12 SS with SS18/SSX1 and 9 with SS18/SSX2 revealed significant differences in gene expression profiles. Among the discriminators were several genes that have previously been found to be upregulated in SS, including AXL, ZIC2, SPAG7, AGRN, FOXC1, NCAM1 and multiple metallothioneins. Histopathology and degree of cytogenetic complexity did not significantly influence expression, whereas a genetic signature that related to development of metastases could be discerned, albeit with a high false-positive rate. In conclusion, our findings demonstrate differentially expressed genes for the 2 major gene fusion variants in SS, SS18/SSX1 and SS18/SSX2, and thereby suggest that these result in different downstream effects.

List of lists-annotated (LOLA): a database for annotation and comparison of published microarray gene lists

Microarray profiling of RNA expression is a powerful tool that generates large lists of transcripts that are potentially relevant to a disease or treatment. However, because the lists of changed transcripts are embedded in figures and tables, they are typically inaccessible for search engines. Due to differences in gene nomenclatures, the lists are difficult to compare between studies. LOLA (Lists of Lists Annotated) is an internet-based database for comparing gene lists from microarray studies or other genomic-scale methods. It serves as a common platform to compare and reannotate heterogeneous gene lists from different microarray platforms or different genomic methodologies such as serial analysis of gene expression (SAGE) or proteomics. LOLA () provides researchers with a means to store, annotate, and compare gene lists produced from different studies or different analyses of the same study. It is especially useful in identifying potentially "high interest" genes which are reported as significant across multiple studies and species. Its application to the fields of stem cell, cancer, and aging research is demonstrated by comparing published papers.

Microarray analysis of gene expression following the formalin test in the infant rat☆

Injury and pain experienced by the infant results in immediate changes in pain sensitivity that last into adulthood. These long-term changes are likely initiated by altered gene expression. Here we measured how injury alters gene expression in the lumbar spinal cord early and late in the preweaning period of the rat. The expression of large numbers of genes was altered significantly at 3 days of age, when injury has long-term consequences. The functional classes of altered genes included transcription factors, cell death related and metal ion genes. The intensity of the stimulus in the 3-day-old pups induced changes in different classes of genes. Fewer changes were noted at 21 days of age. The increased expression of transcription factors and decreased expression of genes whose products are protective against cell death are hypothesized to underlie the long-term changes that are seen after injury in the neonate.

Distinct Genomic Profiles in Hereditary Breast Tumors Identified by Array-Based Comparative Genomic Hybridization

Mutations in BRCA1 and BRCA2 account for a significant proportion of hereditary breast cancers. Earlier studies have shown that inherited and sporadic tumors progress along different somatic genetic pathways and that global gene expression profiles distinguish between these groups. To determine whether genomic profiles similarly discriminate among BRCA1, BRCA2, and sporadic tumors, we established DNA copy number profiles using comparative genomic hybridization to BAC-clone microarrays providing <1 Mb resolution. Tumor DNA was obtained from BRCA1 (n = 14) and BRCA2 (n = 12) mutation carriers, as well as sporadic cases (n = 26). Overall, BRCA1 tumors had a higher frequency of copy number alterations than sporadic breast cancers (P = 0.00078). In particular, frequent losses on 4p, 4q, and 5q in BRCA1 tumors and frequent gains on 7p and 17q24 in BRCA2 tumors distinguish these from sporadic tumors. Distinct amplicons at 3q27.1-q27.3 were identified in BRCA1 tumors and at 17q23.3-q24.2 in BRCA2 tumors. A homozygous deletion on 5q12.1 was found in a BRCA1 tumor. Using a set of 169 BAC clones that detect significantly (P < 0.001) different frequencies of copy number changes in inherited and sporadic tumors, these could be discriminated into separate groups using hierarchical clustering. By comparing DNA copy number and RNA expression for genes in these regions, several candidate genes affected by up- or down-regulation were identified. Moreover, using support vector machines, we correctly classified BRCA1 and BRCA2 tumors (P < 0.0000004 and 0.00005, respectively). Further validation may prove this tumor classifier to be useful for selecting familial breast cancer cases for further mutation screening, particularly, as these data can be obtained using archival tissue.

Selective capture of Salmonella enterica serovar typhi genes expressed in macrophages that are absent from the Salmonella enterica serovar Typhimurium genome

Thirty-six Salmonella enterica serovar Typhi-specific genes, absent from the Salmonella enterica serovar Typhimurium genome, that were expressed in human macrophages were identified by selective capture of transcribed sequences. These genes are located on 15 unique loci of the serovar Typhi genome, including Salmonella pathogenicity islands (SPI-7, SPI-8, and SPI-10) and bacteriophages (ST15, ST18, and ST35).

Comparative analyses of fundamental differences in membrane transport capabilities in prokaryotes and eukaryotes

Whole-genome transporter analyses have been conducted on 141 organisms whose complete genome sequences are available. For each organism, the complete set of membrane transport systems was identified with predicted functions, and classified into protein families based on the transporter classification system. Organisms with larger genome sizes generally possessed a relatively greater number of transport systems. In prokaryotes and unicellular eukaryotes, the significant factor in the increase in transporter content with genome size was a greater diversity of transporter types. In contrast, in multicellular eukaryotes, greater number of paralogs in specific transporter families was the more important factor in the increase in transporter content with genome size. Both eukaryotic and prokaryotic intracellular pathogens and endosymbionts exhibited markedly limited transport capabilities. Hierarchical clustering of phylogenetic profiles of transporter families, derived from the presence or absence of a certain transporter family, showed that clustering patterns of organisms were correlated to both their evolutionary history and their overall physiology and lifestyles.

Membrane transporters are the cell's equivalent of delivery vehicles, garbage disposals, and communication systems—proteins that negotiate through cell membranes to deliver essential nutrients, eject waste products, and help the cell sense environmental conditions around it. Membrane transport systems play crucial roles in fundamental cellular processes of all organisms. The suite of transporters in any one organism also sheds light on its lifestyle and physiology. Up to now, analysis of membrane transporters has been limited mainly to the examination of transporter genes of individual organisms. But advances in genome sequencing have now made it possible for scientists to compare transport and other essential cellular processes across a range of organisms in all three domains of life.

Ren and Paulsen present the first comprehensive bioinformatic analysis of the predicted membrane transporter content of 141 different prokaryotic and eukaryotic organisms. The scientists developed a new computational application of the phylogenetic profiling approach to cluster together organisms that appear to have similar suites of transporters. For example, a group of obligate intracellular pathogens and endosymbionts possess only limited transporter systems in spite of the massive metabolite fluxes one would expect between the symbionts and their host. This is likely due to the relatively static nature of their intracellular environment. In contrast, a cluster of plant/soil-associated microbes encode a robust array of transporters, reflecting the organisms' versatility as well as their exposure to a wide range of different substrates in their natural environment.

Transcriptional profiling of the hyperthermophilic methanarchaeon Methanococcus jannaschii in response to lethal heat and non-lethal cold shock

Temperature shock of the hyperthermophilic methanarchaeon Methanococcus jannaschii from its optimal growth temperature of 85 degrees C to 65 degrees C and 95 degrees C resulted in different transcriptional responses characteristic of both the direction of shock (heat or cold shock) and whether the shock was lethal. Specific outcomes of lethal heat shock to 95 degrees C included upregulation of genes encoding chaperones, and downregulation of genes encoding subunits of the H+ transporting ATP synthase. A gene encoding an alpha subunit of a putative prefoldin was also upregulated, which may comprise a novel element in the protein processing pathway in M. jannaschii. Very different responses were observed upon cold shock to 65 degrees C. These included upregulation of a gene encoding an RNA helicase and other genes involved in transcription and translation, and upregulation of genes coding for proteases and transport proteins. Also upregulated was a gene that codes for an 18 kDa FKBP-type PPIase, which may facilitate protein folding at low temperatures. Transcriptional profiling also revealed several hypothetical proteins that respond to temperature stress conditions.

Neurogenic potential of human mesenchymal stem cells revisited: analysis by immunostaining, time-lapse video and microarray

The possibility of generating neural cells from human bone-marrow-derived mesenchymal stem cells (hMSCs) by simple in vitro treatments is appealing both conceptually and practically. However, whether phenotypic modulations observed after chemical manipulation of such stem cells truly represent a genuine trans-lineage differentiation remains to be established. We have re-evaluated the effects of a frequently reported biochemical approach, based on treatment with butylated hydroxyanisole and dimethylsulphoxide, to bring about such phenotypic conversion by monitoring the morphological changes induced by the treatment in real time, by analysing the expression of phenotype-specific protein markers and by assessing the modulation of transcriptome. Video time-lapse microscopy showed that conversion of mesenchymal stem cells to a neuron-like morphology could be reproduced in normal primary fibroblasts as well as mimicked by addition of drugs eliciting cytoskeletal collapse and disruption of focal adhesion contacts. Analysis of markers revealed that mesenchymal stem cells constitutively expressed multi-lineage traits, including several pertaining to the neural one. However, the applied ;neural induction' protocol neither significantly modulated the expression of such markers, nor induced de novo translation of other neural-specific proteins. Similarly, global expression profiling of over 21,000 genes demonstrated that gene transcription was poorly affected. Most strikingly, we found that the set of genes whose expression was altered by the inductive treatment did not match those sets of genes differentially expressed when comparing untreated mesenchymal stem cells and immature neural tissues. Conversely, by comparing these gene expression profiles with that obtained from comparisons between the same cells and an unrelated non-neural organ, such as liver, we found that the adopted neural induction protocol was no more effective in redirecting human mesenchymal stem cells toward a neural phenotype than toward an endodermal hepatic pathway.

Automated Microarray Image Analysis Toolbox for MATLAB

The Automated Microarray Image Analysis (AMIA) Toolbox for MATLAB is a flexible, open-source, microarray image analysis tool that allows the user to customize analyses of microarray image sets. This tool provides several methods to identify and quantify spot statistics, as well as extensive diagnostic statistics and images to evaluate data quality and array processing. The open, modular nature of AMIA provides access to implementation details and encourages modification and extension of AMIA's capabilities.

AVAILABILITY

The AMIA Toolbox is freely available at http://www.pnl.gov/statistics/amia. The AMIA Toolbox requires MATLAB 6.5 (R13) (MathWorks, Inc. Natick, MA), as well as the Statistics Toolbox 4.1 and Image Processing Toolbox 4.1 for MATLAB or more recent versions.

CONTACT

amanda.white@pnl.gov

Universal mouse reference RNA derived from neonatal mice

A reproducible, transcriptionally diverse common reference RNA is required for accurate comparisons of data generated from most spotted microarray experiments in different experiments. Several methods have been proposed to make such a reference RNA, such as pooling RNAs isolated from multiple cell lines or tissues, amplifying pooled RNAs, or synthesizing RNAs or DNAs complementary to microarray features. We report an approach to prepare a large amount of mouse reference RNA from whole neonatal mice. This approach is simple, quick, reliable, reproducible, and inexpensive. The whole mouse reference RNA is highly representative when compared to two commercially available universal mouse reference RNAs isolated and pooled from multiple cell lines or organs.

Expression clustering reveals detailed co-expression patterns of functionally related proteins during B cell differentiation: a proteomic study using a combination of one-dimensional gel electrophoresis, LC-MS/MS, and stable isotope labeling by amino acids in cell culture (SILAC)

B cells play an essential role in the immune response. Upon activation they may differentiate into plasma cells that secrete specific antibodies against potentially pathogenic non-self antigens. To identify the cellular proteins that are important for efficient production of these antibodies we set out to study the B cell differentiation process at the proteome level. We performed an in-depth proteomic study to quantify dynamic relative protein expression patterns of several hundreds of proteins at five consecutive time points after lipopolysaccharide-induced activation of B lymphocytes. The proteome analysis was performed using a combination of stable isotope labeling using [13C6]leucine added to the murine B cell cultures, one-dimensional gel electrophoresis, and LC-MS/MS. In this study we identified 1,001 B cell proteins. We were able to quantify the expression levels of a quarter of all identified proteins (i.e. 234) at each of the five different time points. Nearly all proteins revealed changes in expression patterns. The quantitative dataset was further analyzed using an unbiased clustering method. Based on their expression profiles, we grouped the entire set of 234 quantified proteins into a limited number of 12 distinct clusters. Functionally related proteins showed a strong correlation in their temporal expression profiles. The quality of the quantitative data allowed us to even identify subclusters within functionally related classes of proteins such as in the endoplasmic reticulum proteins that are involved in antibody production.

Direct activation of human endothelial cells by Plasmodium falciparum-infected erythrocytes

Cytoadherence of Plasmodium falciparum-infected erythrocytes (PRBC) to endothelial cells causes severe clinical disease, presumably as a of result perfusion failure and tissue hypoxia. Cytoadherence to endothelial cells is increased by endothelial cell activation, which is believed to occur in a paracrine fashion by mediators such as tumor necrosis factor alpha (TNF-alpha) released from macrophages that initially recognize PRBC. Here we provide evidence that PRBC directly stimulate human endothelial cells in the absence of macrophages, leading to increased expression of adhesion-promoting molecules, such as intercellular adhesion molecule 1. Endothelial cell stimulation by PRBC required direct physical contact for a short time (30 to 60 min) and was correlated with parasitemia. Gene expression profiling of endothelial cells stimulated by PRBC revealed increased expression levels of chemokine and adhesion molecule genes. PRBC-stimulated endothelial cells especially showed increased expression of molecules involved in parasite adhesion but failed to express molecules promoting leukocyte adhesion, such as E-selectin and vascular cell adhesion molecule 1, even after challenge with TNF-alpha. Collectively, our data suggest that stimulation of endothelial cells by PRBC may have two effects: prevention of parasite clearance through increased cytoadherence and attenuation of leukocyte binding to endothelial cells, thereby preventing deleterious immune reactivity.

Exploring hepatic hormone actions using a compilation of gene expression profiles

Background

Microarray analysis is attractive within the field of endocrine research because regulation of gene expression is a key mechanism whereby hormones exert their actions. Knowledge discovery and testing of hypothesis based on information-rich expression profiles promise to accelerate discovery of physiologically relevant hormonal mechanisms of action. However, most studies so-far concentrate on the analysis of actions of single hormones and few examples exist that attempt to use compilation of different hormone-regulated expression profiles to gain insight into how hormone act to regulate tissue physiology. This report illustrates how a meta-analysis of multiple transcript profiles obtained from a single tissue, the liver, can be used to evaluate relevant hypothesis and discover novel mechanisms of hormonal action. We have evaluated the differential effects of Growth Hormone (GH) and estrogen in the regulation of hepatic gender differentiated gene expression as well as the involvement of sterol regulatory element-binding proteins (SREBPs) in the hepatic actions of GH and thyroid hormone.

Results

Little similarity exists between liver transcript profiles regulated by 17-α-ethinylestradiol and those induced by the continuos infusion of bGH. On the other hand, strong correlations were found between both profiles and the female enriched transcript profile. Therefore, estrogens have feminizing effects in male rat liver which are different from those induced by GH. The similarity between bGH and T3 were limited to a small group of genes, most of which are involved in lipogenesis. An in silico promoter analysis of genes rapidly regulated by thyroid hormone predicted the activation of SREBPs by short-term treatment in vivo. It was further demonstrated that proteolytic processing of SREBP1 in the endoplasmic reticulum might contribute to the rapid actions of T3 on these genes.

Conclusion

This report illustrates how a meta-analysis of multiple transcript profiles can be used to link knowledge concerning endocrine physiology to hormonally induced changes in gene expression. We conclude that both GH and estrogen are important determinants of gender-related differences in hepatic gene expression. Rapid hepatic thyroid hormone effects affect genes involved in lipogenesis possibly through the induction of SREBP1 proteolytic processing.

Neuropeptide and calcium-binding protein gene expression profiles predict neuronal anatomical type in the juvenile rat

Neocortical neurones can be classified according to several independent criteria: morphological, physiological, and molecular expression (neuropeptides (NPs) and/or calcium-binding proteins (CaBPs)). While it has been suggested that particular NPs and CaBPs characterize certain anatomical subtypes of neurones, there is also considerable overlap in their expression, and little is known about simultaneous expression of multiple NPs and CaBPs in morphologically characterized neocortical neurones. Here we determined the gene expression profiles of calbindin (CB), parvalbumin (PV), calretinin (CR), neuropeptide Y (NPY), vasoactive intestinal peptide (VIP), somatostatin (SOM) and cholecystokinin (CCK) in 268 morphologically identified neurones located in layers 2-6 in the juvenile rat somatosensory neocortex. We used patch-clamp electrodes to label neurones with biocytin and harvest the cytoplasm to perform single-cell RT-multiplex PCR. Quality threshold clustering, an unsupervised algorithm that clustered neurones according to their entire profile of expressed genes, revealed seven distinct clusters. Surprisingly, each cluster preferentially contained one anatomical class. Artificial neural networks using softmax regression predicted anatomical types at nearly optimal statistical levels. Classification tree-splitting (CART), a simple binary neuropeptide decision tree algorithm, revealed the manner in which expression of the multiple mRNAs relates to different anatomical classes. Pruning the CART tree revealed the key predictors of anatomical class (in order of importance: SOM, PV, VIP, and NPY). We reveal here, for the first time, a strong relationship between specific combinations of NP and CaBP gene expressions and the anatomical class of neocortical neurones.

Understanding bovine trypanosomiasis and trypanotolerance: the promise of functional genomics

African bovine trypanosomiasis, caused by the protozoan parasite Trypanosoma congolense, is endemic throughout sub-Saharan Africa and is a major constraint on livestock production. A promising approach to disease control is to understand and exploit naturally evolved trypanotolerance. We describe the first attempt to investigate the transcriptional response of susceptible Boran (Bos indicus) cattle to trypanosome infection via a functional genomics approach using a bovine total leukocyte (BOTL) cDNA microarray platform. Four male Boran cattle were experimentally infected with T. congolense and peripheral blood mononuclear cells (PBMC) were collected before infection and 13, 17, 23 and 30 days post-infection (dpi). A reference experimental design was employed using a universal bovine reference RNA pool. Data were normalised to the median of a set of invariant genes (GAPDH) and BRB-Array tools was used to search for statistically significant differentially expressed genes between each time-point. Using a set of 20 microarray hybridisations, we have made a significant contribution to understand the temporal transcriptional response of bovine PBMC in vivo to a controlled trypanosome infection. The greatest changes were evident 13 dpi after parasites were first detected in the blood. Significant differences were observed in clusters of protein kinase C subunits and MHC class I/II related molecules.

Genome-wide analysis of gene expression profiles associated with cell cycle transitions in growing organs of Arabidopsis

Organ growth results from the progression of component cells through subsequent phases of proliferation and expansion before reaching maturity. We combined kinematic analysis, flowcytometry, and microarray analysis to characterize cell cycle regulation during the growth process of leaves 1 and 2 of Arabidopsis (Arabidopsis thaliana). Kinematic analysis showed that the epidermis proliferates until day 12; thereafter, cells expand until day 19 when leaves reach maturity. Flowcytometry revealed that endoreduplication occurs from the time cell division rates decline until the end of cell expansion. Analysis of 10 time points with a 6k-cDNA microarray showed that transitions between the growth stages were closely reflected in the mRNA expression data. Subsequent genome-wide microarray analysis on the three main stages allowed us to categorize known cell cycle genes into three major classes: constitutively expressed, proliferative, and inhibitory. Comparison with published expression data obtained from root zones corresponding to similar developmental stages and from synchronized cell cultures supported this categorization and enabled us to identify a high confidence set of 131 proliferation genes. Most of those had an M phase-dependent expression pattern and, in addition to many known cell cycle-related genes, there were at least 90 that were unknown or previously not associated with proliferation.

MeSHer: identifying biological concepts in microarray assays based on PubMed references and MeSH terms

MeSHer uses a simple statistical approach to identify biological concepts in the form of Medical Subject Headings (MeSH terms) obtained from the PubMed database that are significantly overrepresented within the identified gene set relative to those associated with the overall collection of genes on the underlying DNA microarray platform. As a demonstration, we apply this approach to gene lists acquired from a published study of the effects of angiotensin II (Ang II) treatment on cardiac gene expression and demonstrate that this approach can aid in the interpretation of the resulting 'significant' gene set.

AVAILABILITY

The software is available at http://www.tm4.org.

SUPPLEMENTARY INFORMATION

Results from the analysis of significant genes from the published Ang II study.

CGHAnalyzer: a stand-alone software package for cancer genome analysis using array-based DNA copy number data

SUMMARY

This synopsis provides an overview of array-based comparative genomic hybridization data display, abstraction and analysis using CGHAnalyzer, a software suite, designed specifically for this purpose. CGHAnalyzer can be used to simultaneously load copy number data from multiple platforms, query and describe large, heterogeneous datasets and export results. Additionally, CGHAnalyzer employs a host of algorithms for microarray analysis that include hierarchical clustering and class differentiation.

AVAILABILITY

CGHAnalyzer, the accompanying manual, documentation and sample data are available for download at http://acgh.afcri.upenn.edu. This is a Java-based application built in the framework of the TIGR MeV that can run on Microsoft Windows, Macintosh OSX and a variety of Unix-based platforms. It requires the installation of the free Java Runtime Environment 1.4.1 (or more recent) (http://www.java.sun.com).

Gene expression analysis of immune-mediated arrest of tumorigenesis in a transgenic mouse model of HER-2/neu-positive basal-like mammary carcinoma

We previously showed that a vaccine combining interleukin 12 and allogeneic p185(neu)-positive mammary carcinoma cells completely prevented multifocal mammary carcinogenesis in HER-2/neu transgenic mice. To identify the molecular events responsible for effective tumor prevention and to define the tumor gene expression signature, we used microarrays to analyze the expression profile of mammary tissue of untreated transgenic mice and of vaccine-treated, tumor-free mice at different time points. Mammary tissue from vaccinated mice displayed a gene expression profile different from that of untreated, tumor-bearing mice but similar to that of normal/hyperplastic mammary gland. Comparison of treated and untreated mice at 15 weeks of age revealed up-regulation of genes encoding antibodies, chemokines, gamma-interferon-induced genes and inflammatory molecules, and down-regulation of early genes induced by tumor development. The gene expression signature of HER-2/neu-transformed tumor cells showed modulation of genes promoting proliferation, angiogenesis, migration, invasion, and metastasis and inhibiting apoptosis and immune response. Meta-analysis of microarray data on human breast cancer showed that the signature of tumors arising in murine HER-2/neu transgenic model correctly classified human HER-2/neu-expressing tumors and normal breast tissue. Moreover murine and human HER-2/neu-positive tumors share the signature of basal-like breast cancers. This gene expression analysis reveals the immune events associated with prevention of tumor development and shows that HER-2/neu transgenic mice represent a good model of a poor-prognosis group of human breast tumors.

DNA microarrays: from structural genomics to functional genomics. The applications of gene chips in dermatology and dermatopathology

The human genome project was successful in sequencing the entire human genome and ended earlier than expected. The vast genetic information now available will have far-reaching consequences for medicine in the twenty-first century. The knowledge gained from the mapping and sequencing of human genes on a genome-wide scale--commonly referred to as structural genomics--is prerequisite for studies that focus on the functional aspects of genes. A recently invented technique, known as gene chip, or DNA microarray, technology, allows the study of the function of thousands of genes at once, thereby opening the door to the new field of functional genomics. At its core, the DNA microarray utilizes a unique feature of DNA known as complementary hybridization. As such, it is not different from Southern (DNA) blot or northern (RNA) blot hybridizations, or the polymerase chain reaction, with the exception that it allows expression profiling of the entire human genome in a single hybridization experiment. The article highlights the principles, technology, and applications of DNA microarrays as they pertain to the field of dermatology and dermatopathology. The most important applications are the gene expression profiling of skin cancer, especially of melanoma. Other potential applications include gene expression profiling of inflammatory skin diseases, the mutational analysis of genodermatoses, and polymorphism screening, as well as drug development and chemosensitivity prediction. cDNA microarrays will shape the diagnostic approach of the dermatology and the dermatopathology of the future and may lead to new therapeutic options.

Identification by suppression subtractive hybridization of genes that are differentially expressed between near-isogenic maize lines in association with sugarcane mosaic virus resistance

The molecular mechanisms underlying the development and progression of sugarcane mosaic virus (SCMV) infection in maize are poorly understood. A study of differential expression was conducted to identify genes involved in resistance to SCMV. In this study, we combined suppression subtractive hybridization and macroarray hybridization to identify genes that are differently expressed in the near isogenic lines F7+ (SCMV resistant) and F7 (susceptible). Altogether, 302 differentially expressed genes were identified in four comparisons based on constitutively expressed and inducible genes, and on compatible and incompatible plant-virus interactions. Apart from genes related to metabolism, most of the functionally classified genes identified belonged to three pathogenesis-related categories: cell rescue, defense, cell death and ageing, signal transduction, and transcription. The latter three groups accounted for 56-66% of the genes classified. Some 19% (60 of 302) of the identified genes had previously been assigned to 29 bins distributed over all ten maize chromosomes. Among the mapped genes, 31% (18 of 58) were located within the Scmv2 and Scmv1 regions on chromosomes 3 and 6, respectively, which have been associated with resistance to SCMV. Promising candidate genes for Scmv1 have been identified, such as AA661457 (receptor-like kinase Xa21-binding protein 3). The implications of the genomic distribution of differentially expressed genes identified by this isogenic comparison are discussed in the context of breeding for resistance.

Visualization methods for statistical analysis of microarray clusters

BACKGROUND

The most common method of identifying groups of functionally related genes in microarray data is to apply a clustering algorithm. However, it is impossible to determine which clustering algorithm is most appropriate to apply, and it is difficult to verify the results of any algorithm due to the lack of a gold-standard. Appropriate data visualization tools can aid this analysis process, but existing visualization methods do not specifically address this issue.

RESULTS

We present several visualization techniques that incorporate meaningful statistics that are noise-robust for the purpose of analyzing the results of clustering algorithms on microarray data. This includes a rank-based visualization method that is more robust to noise, a difference display method to aid assessments of cluster quality and detection of outliers, and a projection of high dimensional data into a three dimensional space in order to examine relationships between clusters. Our methods are interactive and are dynamically linked together for comprehensive analysis. Further, our approach applies to both protein and gene expression microarrays, and our architecture is scalable for use on both desktop/laptop screens and large-scale display devices. This methodology is implemented in GeneVAnD (Genomic Visual ANalysis of Datasets) and is available at http://function.princeton.edu/GeneVAnD.

CONCLUSION

Incorporating relevant statistical information into data visualizations is key for analysis of large biological datasets, particularly because of high levels of noise and the lack of a gold-standard for comparisons. We developed several new visualization techniques and demonstrated their effectiveness for evaluating cluster quality and relationships between clusters.

Transcriptional co-regulation of secondary metabolism enzymes in Arabidopsis: functional and evolutionary implications

The combined knowledge of the Arabidopsis genome and transcriptome now allows to get an integrated view of the dynamics and evolution of metabolic pathways in plants. We used publicly available sets of microarray data obtained in a wide range of different stress and developmental conditions to investigate the co-expression of genes encoding enzymes of secondary metabolism pathways, in particular indoles, phenylpropanoids, and flavonoids. We performed hierarchical clustering of gene expression profiles and found that major enzymes of each pathway display a clear and robust co-expression throughout all the conditions studied. Moreover, detailed analysis evidenced that some genes display co-regulation in particular physiological conditions only, certainly reflecting their modular recruitment into stress- or developmentally regulated biosynthetic pathways. The combination of these microarray data with sequence analysis allows to draw very precise hypotheses on the function of otherwise uncharacterized genes. To illustrate this approach, we focused our analysis on secondary metabolism glycosyltransferases (UGTs), a multigenic family involved in the conjugation of small molecules to sugars like glucose. We propose that UGT74B1 and UGT74C1 may be involved in aromatic and aliphatic glucosinolates synthesis, respectively. We also suggest that UGT75C1 may function as an anthocyanin-5-O-glucosyltransferase in planta. Therefore, this data-mining approach appears very powerful for the functional prediction of unknown genes, and could be transposed to virtually any other gene family. Finally, we suggest that analysis of expression pattern divergence of duplicated genes also provides some insight into the mechanisms of metabolic pathway evolution.

Stress-induced co-expression of alternative respiratory chain components in Arabidopsis thaliana

Plant mitochondria contain non-phosphorylating bypasses of the respiratory chain, catalysed by the alternative oxidase (AOX) and alternative NADH dehydrogenases (NDH), as well as uncoupling (UCP) protein. Each of these components either circumvents or short-circuits proton translocation pathways, and each is encoded by a small gene family in Arabidopsis. Whole genome microarray experiments were performed with suspension cell cultures to examine the effects of various 3 h treatments designed to induce abiotic stress. The expression of over 60 genes encoding components of the classical, phosphorylating respiratory chain and tricarboxylic acid cycle remained largely constant when cells were subjected to a broad range of abiotic stresses, but expression of the alternative components responded differentially to the various treatments. In detailed time-course quantitative PCR analysis, specific members of both AOX and NDH gene families displayed coordinated responses to treatments. In particular, the co-expression of AOX1a and NDB2 observed under a number of treatments suggested co-regulation that may be directed by common sequence elements arranged hierarchically in the upstream promoter regions of these genes. A series of treatment sets were identified, representing the response of specific AOX and NDH genes to mitochondrial inhibition, plastid inhibition and abiotic stresses. These treatment sets emphasise the multiplicity of pathways affecting alternative electron transport components in plants.

Gene expression profiling of leukemic cell lines reveals conserved molecular signatures among subtypes with specific genetic aberrations

Hematologic malignancies are characterized by fusion genes of biological/clinical importance. Immortalized cell lines with such aberrations are today widely used to model different aspects of leukemogenesis. Using cDNA microarrays, we determined the gene expression profiles of 40 cell lines as well as of primary leukemias harboring 11q23/MLL rearrangements, t(1;19)[TCF3/PBX1], t(12;21)[ETV6/RUNX1], t(8;21)[RUNX1/CBFA2T1], t(8;14)[IGH@/MYC], t(8;14)[TRA@/MYC], t(9;22)[BCR/ABL1], t(10;11)[PICALM/MLLT10], t(15;17)[PML/RARA], or inv(16)[CBFB/MYH11]. Unsupervised classification revealed that hematopoietic cell lines of diverse origin, but with the same primary genetic changes, segregated together, suggesting that pathogenetically important regulatory networks remain conserved despite numerous passages. Moreover, primary leukemias cosegregated with cell lines carrying identical genetic rearrangements, further supporting that critical regulatory pathways remain intact in hematopoietic cell lines. Transcriptional signatures correlating with clinical subtypes/primary genetic changes were identified and annotated based on their biological/molecular properties and chromosomal localization. Furthermore, the expression profile of tyrosine kinase-encoding genes was investigated, identifying several differentially expressed members, segregating with primary genetic changes, which may be targeted with tyrosine kinase inhibitors. The identified conserved signatures are likely to reflect regulatory networks of importance for the transforming abilities of the primary genetic changes and offer important pathogenetic insights as well as a number of targets for future rational drug design.

Independence and reproducibility across microarray platforms

Microarrays have been widely used for the analysis of gene expression, but the issue of reproducibility across platforms has yet to be fully resolved. To address this apparent problem, we compared gene expression between two microarray platforms: the short oligonucleotide Affymetrix Mouse Genome 430 2.0 GeneChip and a spotted cDNA array using a mouse model of angiotensin II-induced hypertension. RNA extracted from treated mice was analyzed using Affymetrix and cDNA platforms and then by quantitative RT-PCR (qRT-PCR) for validation of specific genes. For the 11,710 genes present on both arrays, we assessed the relative impact of experimental treatment and platform on measured expression and found that biological treatment had a far greater impact on measured expression than did platform for more than 90% of genes, a result validated by qRT-PCR. In the small number of cases in which platforms yielded discrepant results, qRT-PCR generally did not confirm either set of data, suggesting that sequence-specific effects may make expression predictions difficult to make using any technique.

BarleyBase--an expression profiling database for plant genomics

BarleyBase (BB) (www.barleybase.org) is an online database for plant microarrays with integrated tools for data visualization and statistical analysis. BB houses raw and normalized expression data from the two publicly available Affymetrix genome arrays, Barley1 and Arabidopsis ATH1 with plans to include the new Affymetrix 61K wheat, maize, soybean and rice arrays, as they become available. BB contains a broad set of query and display options at all data levels, ranging from experiments to individual hybridizations to probe sets down to individual probes. Users can perform cross-experiment queries on probe sets based on observed expression profiles and/or based on known biological information. Probe set queries are integrated with visualization and analysis tools such as the R statistical toolbox, data filters and a large variety of plot types. Controlled vocabularies for gene and plant ontologies, as well as interconnecting links to physical or genetic map and other genomic data in PlantGDB, Gramene and GrainGenes, allow users to perform EST alignments and gene function prediction using Barley1 exemplar sequences, thus, enhancing cross-species comparison.

ArrayExpress--a public repository for microarray gene expression data at the EBI

ArrayExpress is a public repository for microarray data that supports the MIAME (Minimum Information About a Microarray Experiment) requirements and stores well-annotated raw and normalized data. As of November 2004, ArrayExpress contains data from approximately 12,000 hybridizations covering 35 species. Data can be submitted online or directly from local databases or LIMS in a standard format, and password-protected access to prepublication data is provided for reviewers and authors. The data can be retrieved by accession number or queried by various parameters such as species, author and array platform. A facility to query experiments by gene and sample properties is provided for a growing subset of curated data that is loaded in to the ArrayExpress data warehouse. Data can be visualized and analysed using Expression Profiler, the integrated data analysis tool. ArrayExpress is available at http://www.ebi.ac.uk/arrayexpress.

Laser capture microdissection for the analysis of gene expression during embryogenesis of Arabidopsis

It is during embryogenesis that the body plan of the developing plant is established. Analysis of gene expression during embryogenesis has been limited due to the technical difficulty of accessing the developing embryo. Here we demonstrate that laser capture microdissection can be applied to the analysis of embryogenesis. We show how this technique can be used in concert with DNA microarray for the large-scale analysis of gene expression in apical and basal domains of the globular-stage and heart-stage embryo, respectively, when critical events of polarity, symmetry and biochemical differentiation are established. This high resolution spatial analysis shows that up to approximately 65% of the genome is expressed in the developing embryo, and that differential expression of a number of gene classes can be detected. We discuss the validity of this approach for the functional analysis of both published and previously uncharacterized essential genes.

Microarray analyses reveal strong influence of DNA copy number alterations on the transcriptional patterns in pancreatic cancer: implications for the interpretation of genomic amplifications

DNA copy number alterations are believed to play a major role in the development and progression of human neoplasms. Although most of these genomic imbalances have been associated with dysregulation of individual genes, their large-scale transcriptional consequences remain unclear. Pancreatic carcinomas frequently display gene copy number variation of entire chromosomes as well as of chromosomal subregions. These changes range from homozygous deletions to high-level amplifications and are believed to constitute key genetic alterations in the cellular transformation of this tumor type. To investigate the transcriptional consequences of the most drastic genomic changes, that is, genomic amplifications, and to analyse the genome-wide transcriptional effects of DNA copy number changes, we performed expression profiling of 29 pancreatic carcinoma cell lines and compared the results with matching genomic profiling data. We show that a strong association between DNA copy numbers and mRNA expression levels is present in pancreatic cancer, and demonstrate that as much as 60% of the genes within highly amplified genomic regions display associated overexpression. Consequently, we identified 67 recurrently overexpressed genes located in seven precisely mapped commonly amplified regions. The presented findings indicate that more than one putative target gene may be of importance in most pancreatic cancer amplicons.

A robust method for the amplification of RNA in the sense orientation

BACKGROUND

Small quantities of RNA (1-4 microg total RNA) available from biological samples frequently require a single round of amplification prior to analysis, but current amplification strategies have limitations that may restrict their usefulness in downstream genomic applications. The Eberwine amplification method has been extensively validated but is limited by its ability to produce only antisense RNA. Alternatives lack extensive validation and are often confounded by problems with bias or yield attributable to their greater biological and technical complexity.

RESULTS

To overcome these limitations, we have developed a straightforward and robust protocol for amplification of RNA in the sense orientation. This protocol is based upon Eberwine's method but incorporates elements of more recent amplification techniques while avoiding their complexities. Our technique yields greater than 100-fold amplification, generates long transcript, and produces mRNA that is well suited for use with microarray applications. Microarrays performed with RNA amplified using this protocol demonstrate minimal amplification bias and high reproducibility.

CONCLUSION

The protocol we describe here is readily adaptable for the production of sense or antisense, labeled or unlabeled RNA from intact or partially-degraded prokaryotic or eukaryotic total RNA. The method outperforms several commercial RNA amplification kits and can be used in conjunction with a variety of microarray platforms, such as cDNA arrays, oligonucleotide arrays, and Affymetrix GeneChip arrays.

Issues for consideration in the analysis of microarray data in behavioural studies

Microarrays are one of several technologies that allow for measurement the expression of thousands of genes simultaneously. This technological advance provides a challenge for the analysis of these data. In this review we discuss these analytical issues from the initial quality control to normalization, differential expression, clustering and finally functional pathway analysis. We focus on Affymetrix data but many of the issues are the same for other array platforms.

Global transcription profiling reveals comprehensive insights into hypoxic response in Arabidopsis

Plants have evolved adaptation mechanisms to sense oxygen deficiency in their environments and make coordinated physiological and structural adjustments to enhance their hypoxic tolerance. To gain insight into how plants respond to low-oxygen stress, gene expression profiling using whole-genome DNA amplicon microarrays was carried out at seven time points over 24 h, in wild-type and transgenic P(SAG12):ipt Arabidopsis (Arabidopsis thaliana) plants under normoxic and hypoxic conditions. Transcript levels of genes involved in glycolysis and fermentation pathways, ethylene synthesis and perception, calcium signaling, nitrogen utilization, trehalose metabolism, and alkaloid synthesis were significantly altered in response to oxygen limitation. Analysis based on gene ontology assignments suggested a significant down-regulation of genes whose functions are associated with cell walls, nucleosome structures, water channels, and ion transporters and a significant up-regulation of genes involved in transcriptional regulation, protein kinase activity, and auxin responses under conditions of oxygen shortage. Promoter analysis on a cluster of up-regulated genes revealed a significant overrepresentation of the AtMYB2-binding motif (GT motif), a sugar response element-like motif, and a G-box-related sequence, and also identified several putative anaerobic response elements. Finally, quantitative real-time polymerase chain reactions using 29 selected genes independently verified the microarray results. This study represents one of the most comprehensive analyses conducted to date investigating hypoxia-responsive transcriptional networks in plants.

Background correction for cDNA microarray images using the TV+L1 model

MOTIVATION

Background correction is an important preprocess in cDNA microarray data analysis. A variety of methods have been used for this purpose. However, many kinds of backgrounds, especially inhomogeneous ones, cannot be estimated correctly using any of the existing methods. In this paper, we propose the use of the TV+L1 model, which minimizes the total variation (TV) of the image subject to an L1-fidelity term, to correct background bias. We demonstrate its advantages over the existing methods by both analytically discussing its properties and numerically comparing it with morphological opening.

RESULTS

Experimental results on both synthetic data and real microarray images demonstrate that the TV+L1 model gives the restored intensity that is closer to the true data than morphological opening. As a result, this method can serve an important role in the preprocessing of cDNA microarray data.

Pathway-specific profiling identifies the NF-kappa B-dependent tumor necrosis factor alpha-regulated genes in epidermal keratinocytes

Identification of tumor necrosis factor alpha (TNF alpha) as the key agent in inflammatory disorders led to new therapies specifically targeting TNF alpha and avoiding many side effects of earlier anti-inflammatory drugs. However, because of the wide spectrum of systems affected by TNF alpha, drugs targeting TNF alpha have a potential risk of delaying wound healing, secondary infections, and cancer. Indeed, increased risks of tuberculosis and carcinogenesis have been reported as side effects after anti-TNF alpha therapy. TNF alpha regulates many processes (e.g. immune response, cell cycle, and apoptosis) through several signal transduction pathways that convey the TNF alpha signals to the nucleus. Hypothesizing that specific TNF alpha-dependent pathways control specific processes and that inhibition of a specific pathway may yield even more precisely targeted therapies, we used oligonucleotide microarrays and parthenolide, an NF-kappa B-specific inhibitor, to identify the NF-kappa B-dependent set of the TNF alpha-regulated genes in human epidermal keratinocytes. Expression of approximately 40% of all TNF alpha-regulated genes depends on NF-kappa B; 17% are regulated early (1-4 h post-treatment), and 23% are regulated late (24-48 h). Cytokines and apoptosis-related and cornification proteins belong to the "early" NF-kappa B-dependent group, and antigen presentation proteins belong to the "late" group, whereas most cell cycle, RNA-processing, and metabolic enzymes are not NF-kappa B-dependent. Therefore, inflammation, immunomodulation, apoptosis, and differentiation are on the NF-kappa B pathway, and cell cycle, metabolism, and RNA processing are not. Most early genes contain consensus NF-kappaB binding sites in their promoter DNA and are, presumably, directly regulated by NF-kappa B, except, curiously, the cornification markers. Using siRNA silencing, we identified cFLIP/CFLAR as an essential NF-kappa B-dependent antiapoptotic gene. The results confirm our hypothesis, suggesting that inhibiting a specific TNF alpha-dependent signaling pathway may inhibit a specific TNF alpha-regulated process, leaving others unaffected. This could lead to more specific anti-inflammatory agents that are both more effective and safer.

Analytical trapping: extraction of time-independent structures from time-dependent crystallographic data

All chemical and biological reactions involve atomic motion, embodied in dynamic structural changes. Identifying these changes is the goal of time-resolved crystallography. The "raw" output of a time-resolved macromolecular crystallography experiment is the time-dependent set of difference electron density maps that span the desired time range and display the time-dependent changes in density (and underlying structure) as the reaction progresses. The goal is to interpret such data in terms of a small number of crystallographically refinable, time-independent structures, each associated with a reaction intermediate; to establish the pathways and rate coefficients by which the intermediates interconvert; and thus to establish a chemical kinetic mechanism. We review briefly the various strategies that may be used to achieve this goal and concentrate on two promising advances: singular value decomposition and cluster analysis. The strategies are illustrated by using data on the photocycle of the bacterial blue light photoreceptor, photoactive yellow protein.

Murine macrophage transcriptional responses to Bacillus anthracis infection and intoxication

Interactions between Bacillus anthracis and host macrophages represent critical early events in anthrax pathogenesis, but their details are not clearly understood. Here we report the first genomewide characterization of the transcriptional changes within macrophages infected with B. anthracis and the identification of several hundred host genes that were differentially expressed during this intracellular stage of infection. These loci included both genes that are known to be regulated differentially in response to many other bacterial pathogens and those that appear to be differentially regulated in response to B. anthracis but not other bacterial species that have been tested. These data provide a transcriptional basis for a variety of physiological changes observed during infection, including the induction of apoptosis caused by the infecting bacteria. The expression patterns underlying B. anthracis-induced apoptosis led us to test further the importance of one very highly induced macrophage gene, that for ornithine decarboxylase. Our data show that this enzyme plays an important and previously unrecognized role in suppressing apoptosis in B. anthracis-infected cells. We have also characterized the transcriptional response to anthrax lethal toxin in activated macrophages and found that, following toxin treatment, many of the host inflammatory response pathways are dampened. These data provide insights into B. anthracis pathogenesis as well as potential leads for the development of new diagnostic and therapeutic options.

Gene expression profile in response to Xanthomonas axonopodis pv. manihotis infection in cassava using a cDNA microarray

A cassava cDNA microarray based on a large cassava EST database was constructed and used to study the incompatible interaction between cassava and Xanthomonas axonopodis pv. manihotis (Xam) strain CIO151. For microarray construction, 5700 clones from the cassava unigene set were amplified by polymerase chain reaction (PCR) and printed on glass slides. Microarray hybridization was performed using cDNA from cassava plants (resistant variety MBra685) collected at 12, 24, 48 h and 7 and 15 days post-infection as treatment and cDNA from mock-inoculated plants as control. A total of 199 genes were found to be differentially expressed (126 up-regulated and 73 down-regulated). A greater proportion of differentially-expressed genes was observed at 7 days after inoculation. Expression profiling and cluster analyses indicate that, in response to inoculation with Xam, cassava induces dozens of genes, including principally those involved in oxidative burst, protein degradation and pathogenesis-related (PR) genes. In contrast, genes encoding proteins that are involved in photosynthesis and metabolism were down regulated. In addition, various other genes encoding proteins with unknown function or showing no similarity to other proteins were also induced. Quantitative real time PCR experiments confirmed the reliability of our microarray data. In addition we showed that some genes are induced more rapidly in the resistant than in the susceptible cultivar.

Microarray analysis confirms the specificity of a Chlamydomonas reinhardtii chloroplast RNA stability mutant

The expression of chloroplast and mitochondrial genes depends on nucleus-encoded proteins, some of which control processing, stability, and/or translation of organellar RNAs. To test the specificity of one such RNA stability factor, we used two known Chlamydomonas reinhardtii nonphotosynthetic mutants carrying mutations in the Mcd1 nuclear gene (mcd1-1 and mcd1-2). We previously reported that these mutants fail to accumulate the chloroplast petD mRNA and its product, subunit IV of the cytochrome b6/f complex, which is essential for photosynthesis. Such mutants are generally presumed to be gene specific but are not tested rigorously. Here, we have used microarray analysis to assess changes in chloroplast, mitochondrial, and nuclear RNAs, and since few other RNAs were significantly altered in these mutants, conclude that Mcd1 is indeed specifically required for petD mRNA accumulation. In addition, a new unlinked nuclear mutation was discovered in mcd1-2, which greatly reduced chloroplast atpA mRNA accumulation. Genetic analyses showed failure to complement mda1-ncc1, where atpA-containing transcripts are similarly affected (D. Drapier, J. Girard-Bascou, D.B. Stern, F.-A. Wollman [2002] Plant J 31: 687-697), and we have named this putative new allele mda1-2. We conclude that DNA microarrays are efficient and useful for characterizing the specificity of organellar RNA accumulation mutants.

SP100 inhibits ETS1 activity in primary endothelial cells

SP100 was first identified as a nuclear autoimmune antigen and is a constituent of the nuclear body. SP100 interacts with the ETS1 transcription factor, and we have previously shown that SP100 reduces ETS1-DNA binding and inhibits ETS1 transcriptional activity on the MMP1 and uPA promoters. We now demonstrate that SP100 expression is upregulated by interferons, which have been shown to be antiangiogenic, in primary endothelial cells. As ETS1 is functionally important in promoting angiogenesis, we tested the hypothesis that ETS1 activity is negatively modulated by SP100 in endothelial cells. SP100 directly antagonizes ETS1-mediated morphological changes in human umbilical vein endothelial cell (HUVEC) network formation and reduces HUVEC migration and invasion. To further understand the functional relationship between ETS1 and SP100, cDNA microarray analysis was utilized to assess reprogramming of gene expression by ETS1 and SP100. A subset of the differentially regulated genes, including heat-shock proteins (HSPs) H11, HSPA1L, HSPA6, HSPA8, HSPE1 and AXIN1, BRCA1, CD14, CTGF (connective tissue growth factor), GABRE (gamma-aminobutyric acid A receptor epsilon), ICAM1, SNAI1, SRD5A1 (steroid-5-alpha-reductase 1) and THY1, were validated by real-time PCR and a majority showed reciprocal expression in response to ETS1 and SP100. Interestingly, genes that are negatively regulated by ETS1 and upregulated by SP100 have antimigratory or antiangiogenic properties. Collectively, these data indicate that SP100 negatively modulates ETS1-dependent downstream biological processes.

Genetic regulation of endotoxin-induced airway disease

To identify novel genes regulating the biologic response to lipopolysaccharide (LPS), we used a combination of quantitative trait locus (QTL) analysis and microarray-based gene expression studies of C57BL/6J x DBA/2J(BXD) F2 and recombinant inbred (RI) mice. A QTL affecting pulmonary TNF-alpha production was identified on chromosome 2, and a region affecting both polymorphonuclear leukocyte recruitment and TNF-alpha levels was identified on chromosome 11. Microarray analyses of unchallenged and LPS-challenged BXD RI strains identified approximately 500 genes whose expression was significantly changed by inhalation of LPS. Of these genes, 28 reside within the chromosomal regions identified by the QTL analyses, implicating these genes as high priority candidates for functional studies. Additional high priority candidate genes were identified based on their differential expression in mice having high and low responses to LPS. Functional studies of these genes are expected to reveal important molecular mechanisms regulating the magnitude of biologic responses to LPS.

Microarray and real-time RT-PCR analyses of differential human gene expression patterns induced by severe acute respiratory syndrome (SARS) coronavirus infection of Vero cells

Vero E6 African green monkey kidney cells are highly susceptible to infection with the newly emerging severe acute respiratory syndrome coronavirus (SARS-CoV), and they are permissive for rapid viral replication, with resultant cytopathic effects. We employed cDNA microarray analysis to characterize the cellular transcriptional responses of homologous human genes at 12 h post-infection. Seventy mRNA transcripts belonging to various functional classes exhibited significant alterations in gene expression. There was considerable induction of heat shock proteins that are crucial to the immune response mechanism. Modified levels of several transcripts involved in pro-inflammatory and anti-inflammatory processes exemplified the balance between opposing forces during SARS pathogenesis. Other genes displaying altered transcription included those associated with host translation, cellular metabolism, cell cycle, signal transduction, transcriptional regulation, protein trafficking, protein modulators, and cytoskeletal proteins. Alterations in the levels of several novel transcripts encoding hypothetical proteins and expressed sequence tags were also identified. In addition, transcription of apoptosis-related genes DENN and hIAP1 was upregulated in contrast to FAIM. Elevated Mx1 expression signified a strong host response to mediate antiviral resistance. Also expressed in infected cells was the C-terminal alternative splice variant of the p53 tumor suppressor gene encoding a modified truncated protein that can influence the activity of wild-type p53. We observed the interplay between various mechanisms to favor virus multiplication before full-blown apoptosis and the triggering of several pathways in host cells in an attempt to eliminate the pathogen. Microarray analysis identifies the critical host–pathogen interactions during SARS-CoV infection and provides new insights into the pathophysiology of SARS.

Gene Expression Profiles of Bovine Trophoblastic Cell Line (BT-1) Analyzed by a Custom cDNA Microarray

Gene expression of bovine trophoblast cell line (BT-1) was analyzed with a custom utero-placenta complementary DNA microarray. Expression comparison with in vivo tissues of trophoblast derivation was performed to investigate characteristics of the expression in BT-1. BT-1 is a cell line established without feeder cells using trophoblast cells that are separated from bovine blastocysts. The bovine in vivo tissues of the trophoblast derivation were collected on day 17 through 56 of gestation. Of 1,773 genes, 933 exhibited an expression difference exceeding two-fold between BT-1 and in vivo tissues. These genes were analyzed by the k-means clustering method and were distributed into six clusters. Some genes, such as placental lactogen, interferon-tau and pregnancy-associated glycoproteins, all known as trophoblast-specific genes, were detected in BT-1 as well as in vivo tissues throughout the experiment period. These trophoblast-specific genes and octamer-binding transcription factor-4, known as a marker for undifferentiation of cells in mice, were detected by RT-PCR in both BT-1 and in vivo trophoblast tissues. The overall gene expression profile in BT-1 suggests that this cell line contains trophoblast-specific characteristics and is similar to trophoblast cells around the implantation period.

Transcriptional divergence of the duplicated oxidative stress-responsive genes in the Arabidopsis genome

Previous studies have indicated that Arabidopsis thaliana experienced a genome-wide duplication event shortly before its divergence from Brassica followed by extensive chromosomal rearrangements and deletions. While a large number of the duplicated genes have significantly diverged or lost their sister genes, we found 4222 pairs that are still highly conserved, and as a result had similar functional assignments during the annotation of the genome sequence. Using whole-genome DNA microarrays, we identified 906 duplicated gene pairs in which at least one member exhibited a significant response to oxidative stress. Among these, only 117 pairs were up- or down-regulated in both pairs and many of these exhibited dissimilar patterns of expression. Examination of the expression patterns of PAL1 and PAL2, ACD1 and ACD2, genes coding for two Hsp20s, various P450s, and electron transfer flavoproteins suggests Arabidopsis evolved a number of distinct oxidative stress response mechanisms using similar gene sets following the duplication of its genome.