DRAGON: Database Referencing of Array Genes Online

DRAGON and DRAGON view: information annotation and visualization tools for large-scale expression data

The Database Referencing of Array Genes ONline (DRAGON) database system consists of information derived from publicly available databases including UniGene, SWISS-Prot, Pfam, and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Through a Web-accessible interface, the DRAGON Annotate tool rapidly supplies information pertaining to a range of biological characteristics of all the genes in any large-scale gene expression data set. The subsequent inclusion of this information during data analysis and visualization allows for deeper insight into gene expression patterns. The set of DRAGON View tools provides methods for the analysis and visualization of expression patterns in relation to annotated information. Instead of incorporating the standard set of clustering and graphing tools available in many large-scale expression data analysis software packages, DRAGON View has been specifically designed to allow for the analysis of expression data in relation to the biological characteristics of gene sets.

Linking microarray reporters with protein functions

Background

The analysis of microarray experiments requires accurate and up-to-date functional annotation of the microarray reporters to optimize the interpretation of the biological processes involved. Pathway visualization tools are used to connect gene expression data with existing biological pathways by using specific database identifiers that link reporters with elements in the pathways.

Results

This paper proposes a novel method that aims to improve microarray reporter annotation by BLASTing the original reporter sequences against a species-specific EMBL subset, that was derived from and crosslinked back to the highly curated UniProt database. The resulting alignments were filtered using high quality alignment criteria and further compared with the outcome of a more traditional approach, where reporter sequences were BLASTed against EnsEMBL followed by locating the corresponding protein (UniProt) entry for the high quality hits. Combining the results of both methods resulted in successful annotation of > 58% of all reporter sequences with UniProt IDs on two commercial array platforms, increasing the amount of Incyte reporters that could be coupled to Gene Ontology terms from 32.7% to 58.3% and to a local GenMAPP pathway from 9.6% to 16.7%. For Agilent, 35.3% of the total reporters are now linked towards GO nodes and 7.1% on local pathways.

Conclusion

Our methods increased the annotation quality of microarray reporter sequences and allowed us to visualize more reporters using pathway visualization tools. Even in cases where the original reporter annotation showed the correct description the new identifiers often allowed improved pathway and Gene Ontology linking. These methods are freely available at http://www.bigcat.unimaas.nl/public/publications/Gaj_Annotation/.

Microarray analysis of gliomas reveals chromosomal position-associated gene expression patterns and identifies potential immunotherapy targets

Gliomas are among the most aggressive malignant tumors and the most refractory to therapy, in part due to the propensity for malignant cells to disseminate diffusely throughout the brain. Here, we have used 27 K cDNA microarrays to investigate global gene expression changes between normal brain and high-grade glioma (glioblastoma multiforme) to try and better understand gliomagenesis and to identify new therapeutic targets. We have also included smaller groups of grade II and grade III tumors of mixed astrocytic and oligodendroglial origin as comparison. We found that the expression of hundreds of genes was significantly correlated to each group, and employed a naïve Bayesian classifier with leave-one-out cross-validation to accurately classify the samples. We developed a novel algorithm to analyze the gene expression data from the perspective of chromosomal position, and identified distinct regions of the genome that displayed coordinated expression patterns that correlated significantly to tumor grade. The regions identified corresponded to previously known genetic copy number changes in glioma (e.g. 10q23, 10q25, 7q, 7p) as well as regions not previously associated significantly with glioma (e.g. 1p13, 6p22). Furthermore, to enrich for more suitable targets for therapy, we took a bioinformatics approach and annotated our signatures with two published datasets that identified membrane/secreted genes from cytosolic genes. The resulting focused list of 31 genes included interesting novel potential targets as well as several proteins already being investigated for immunotherapy (e.g. CD44 and tenascin-C). Software for the chromosome analysis was developed and is freely available at http://base.thep.lu.se.

Basal-like phenotype is not associated with patient survival in estrogen-receptor-negative breast cancers

INTRODUCTION

Basal-phenotype or basal-like breast cancers are characterized by basal epithelium cytokeratin (CK5/14/17) expression, negative estrogen receptor (ER) status and distinct gene expression signature. We studied the clinical and biological features of the basal-phenotype tumors determined by immunohistochemistry (IHC) and cDNA microarrays especially within the ER-negative subgroup.

METHODS

IHC was used to evaluate the CK5/14 status of 445 stage II breast cancers. The gene expression signature of the CK5/14 immunopositive tumors was investigated within a subset (100) of the breast tumors (including 50 ER-negative tumors) with a cDNA microarray. Survival for basal-phenotype tumors as determined by CK5/14 IHC and gene expression signature was assessed.

RESULTS

From the 375 analyzable tumor specimens, 48 (13%) were immunohistochemically positive for CK5/14. We found adverse distant disease-free survival for the CK5/14-positive tumors during the first years (3 years hazard ratio (HR) 2.23, 95% confidence interval (CI) 1.17 to 4.24, p = 0.01; 5 years HR 1.80, 95% CI 1.02 to 3.15, p = 0.04) but the significance was lost at the end of the follow-up period (10 years HR 1.43, 95% CI 0.84 to 2.43, p = 0.19). Gene expression profiles of immunohistochemically determined CK5/14-positive tumors within the ER-negative tumor group implicated 1,713 differently expressed genes (p < 0.05). Hierarchical clustering analysis with the top 500 of these genes formed one basal-like and a non-basal-like cluster also within the ER-negative tumor entity. A highly concordant classification could be constructed with a published gene set (Sorlie's intrinsic gene set, concordance 90%). Both gene sets identified a basal-like cluster that included most of the CK5/14-positive tumors, but also immunohistochemically CK5/14-negative tumors. Within the ER-negative tumor entity there was no survival difference between the non-basal and basal-like tumors as identified by immunohistochemical or gene-expression-based classification.

CONCLUSION

Basal cytokeratin-positive tumors have a biologically distinct gene expression signature from other ER-negative tumors. Even if basal cytokeratin expression predicts early relapse among non-selected tumors, the clinical outcome of basal tumors is similar to non-basal ER-negative tumors. Immunohistochemically basal cytokeratin-positive tumors almost always belong to the basal-like gene expression profile, but this cluster also includes few basal cytokeratin-negative tumors.

Microarray analysis

Microarrays and related technologies have allowed investigators to ask biological questions in far greater detail than has previously been possible. Microarrays had a troubled beginning, but most of these problems resulted from the growing pains of this technology, which, like many new things, was initially more promise than delivery. Nevertheless, over the past few years, investigators have learned how to achieve optimal performance of technology, and now exciting discoveries are made using microarray-based research. Many of the advances have come from the realization that microarrays are not a magic tool but rather are like any other measurement device. Unless microarray experimentation is coupled with good experimental practices, it will not yield valid results or, worse yet, may lead to misleading results. In this chapter, we highlight some of the important steps that should be taken to successfully conduct a microarray study. These steps include a clearly stated biological question, experimental design, careful experimental conduct, complete statistical analysis, validation/verification of results, and dissemination of the data.

Cross-species comparative toxicogenomics as an aid to safety assessment

Cross-species comparative toxicogenomics has the potential for improving the understanding of the different responses of animal models to toxicants at a molecular level. This understanding could then lead to a more accurate extrapolation of the risk posed by these toxicants to humans. Cross-species comparative studies have been carried out at the genomic sequence level and using microarrays to examine changes in global mRNA profiles. However, these studies face considerable bioinformatic challenges in terms of identifying which genes are truly orthologous across species. The resources to analyse such studies, in the context of such orthologues, beg improvement. Finally, the experimental design of such studies needs to be carefully considered to make their results fully interpretable. These issues are discussed, along with the current state-of-the-art cross-species comparative toxicogenomics in this review.

An open-access long oligonucleotide microarray resource for analysis of the human and mouse transcriptomes

Two collections of oligonucleotides have been designed for preparing pangenomic human and mouse microarrays. A total of 148 993 and 121 703 oligonucleotides were designed against human and mouse transcripts. Quality scores were created in order to select 25 342 human and 24 109 mouse oligonucleotides. They correspond to: (i) a BLAST-specificity score; (ii) the number of expressed sequence tags matching each probe; (iii) the distance to the 3′ end of the target mRNA. Scores were also used to compare in silico the two microarrays with commercial microarrays. The sets described here, called RNG/MRC collections, appear at least as specific and sensitive as those from the commercial platforms. The RNG/MRC collections have now been used by an Anglo-French consortium to distribute more than 3500 microarrays to the academic community. Ad hoc identification of tissue-specific transcripts and a ∼80% correlation with hybridizations performed on Affymetrix GeneChip™ suggest that the RNG/MRC microarrays perform well. This work provides a comprehensive open resource for investigators working on human and mouse transcriptomes, as well as a generic method to generate new microarray collections in other organisms. All information related to these probes, as well as additional information about commercial microarrays have been stored in a freely-accessible database called MEDIANTE.

Comparison of gene coverage of mouse oligonucleotide microarray platforms

Background

The increasing use of DNA microarrays for genetical genomics studies generates a need for platforms with complete coverage of the genome. We have compared the effective gene coverage in the mouse genome of different commercial and noncommercial oligonucleotide microarray platforms by performing an in-house gene annotation of probes. We only used information about probes that is available from vendors and followed a process that any researcher may take to find the gene targeted by a given probe. In order to make consistent comparisons between platforms, probes in each microarray were annotated with an Entrez Gene id and the chromosomal position for each gene was obtained from the UCSC Genome Browser Database. Gene coverage was estimated as the percentage of Entrez Genes with a unique position in the UCSC Genome database that is tested by a given microarray platform.

Results

A MySQL relational database was created to store the mapping information for 25,416 mouse genes and for the probes in five microarray platforms (gene coverage level in parenthesis): Affymetrix430 2.0 (75.6%), ABI Genome Survey (81.24%), Agilent (79.33%), Codelink (78.09%), Sentrix (90.47%); and four array-ready oligosets: Sigma (47.95%), Operon v.3 (69.89%), Operon v.4 (84.03%), and MEEBO (84.03%). The differences in coverage between platforms were highly conserved across chromosomes. Differences in the number of redundant and unspecific probes were also found among arrays. The database can be queried to compare specific genomic regions using a web interface. The software used to create, update and query the database is freely available as a toolbox named ArrayGene.

Conclusion

The software developed here allows researchers to create updated custom databases by using public or proprietary information on genes for any organisms. ArrayGene allows easy comparisons of gene coverage between microarray platforms for any region of the genome. The comparison presented here reveals that the commercial microarray Sentrix, which is based on the MEEBO public oligoset, showed the best mouse genome coverage currently available. We also suggest the creation of guidelines to standardize the minimum set of information that vendors should provide to allow researchers to accurately evaluate the advantages and disadvantages of using a given platform.

Three-color cDNA microarrays with prehybridization quality control yield gene expression data comparable to that of commercial platforms

Despite their lower cost and high content flexibility, a limitation of in-house-prepared arrays has been their susceptibility to quality control (QC) issues and lack of QC standards across laboratories. Therefore, we developed a novel three-color array system that allows prehybridization QC as well as the Matarray software to facilitate acquisition of accurate gene expression data. In this study, we compared performance of our rat cDNA array to the Affymetrix RG-U34A and Agilent G4130A arrays using 2,824 UniGenes represented on all three arrays. Before data filtering, poor interplatform agreement was observed; however, after data filtering, differentially expressed UniGenes exhibited correlation coefficients of 0.91, 0.88, and 0.92 between the Affymetrix vs. Agilent, Affymetrix vs. cDNA, and Agilent vs. cDNA arrays, respectively. The Affymetrix, Agilent, and cDNA arrays agreed well with quantitative RT-PCR conducted on 42 UniGenes, yielding correlation coefficients of 0.90, 0.90, and 0.96, respectively. Each platform underestimated ratios relative to quantitative RT-PCR, possessing respective slopes of 0.86 ( R2 = 0.81), 0.65 ( R2 = 0.81), and 0.70 ( R2 = 0.92). Overall, these data show that the combination of our novel technical and analytic approaches yield an accurate platform for functional genomics that is concordant with commercial discovery arrays in terms of identifying regulated genes and pathways.

Gene expression patterns in dendritic cells infected with measles virus compared with other pathogens

Gene expression patterns supply insight into complex biological networks that provide the organization in which viruses and host cells interact. Measles virus (MV) is an important human pathogen that induces transient immunosuppression followed by life-long immunity in infected individuals. Dendritic cells (DCs) are potent antigen-presenting cells that initiate the immune response to pathogens and are postulated to play a role in MV-induced immunosuppression. To better understand the interaction of MV with DCs, we examined the gene expression changes that occur over the first 24 h after infection and compared these changes to those induced by other viral, bacterial, and fungal pathogens. There were 1,553 significantly regulated genes with nearly 60% of them down-regulated. MV-infected DCs up-regulated a core of genes associated with maturation of antigen-presenting function and migration to lymph nodes but also included genes for IFN-regulatory factors 1 and 7, 2'5' oligoadenylate synthetase, Mx, and TNF superfamily proteins 2, 7, 9, and 10 (TNF-related apoptosis-inducing ligand). MV induced genes for IFNs, ILs, chemokines, antiviral proteins, histones, and metallothioneins, many of which were also induced by influenza virus, whereas genes for protein synthesis and oxidative phosphorylation were down-regulated. Unique to MV were the induction of genes for a broad array of IFN-alphas and the failure to up-regulate dsRNA-dependent protein kinase. These results provide a modular view of common and unique DC responses after infection and suggest mechanisms by which MV may modulate the immune response.

KARMA: a web server application for comparing and annotating heterogeneous microarray platforms

We have developed a universal web server application (KARMA) that allows comparison and annotation of user-defined pairs of microarray platforms based on diverse types of genome annotation data (across different species) collected from multiple sources. The application is an effective tool for diverse microarray platforms, including arrays that are provided by (i) the Keck Microarray Resource at Yale, (ii) commercially available Affymetrix GeneChips and spotted arrays and (iii) custom arrays made by individual academics. The tool provides a web interface that allows users to input pairs of test files that represent diverse array platforms for either single or multiple species. The program dynamically identifies analogous DNA fragments spotted or synthesized on multiple microarray platforms based on the following types of information: (i) NCBI-Unigene identifiers, if the platforms being compared are within the same species or (ii) NCBI-Homologene data, if they are cross-species. The single-species comparison is implemented based on set operations: intersection, union and difference. Other forms of retrievable annotation data, including LocusLink, SwissProt and Gene Ontology (GO), are collected from multiple remote sites and stored in an integrated fashion using an Oracle database. The KARMA database, which is updated periodically, is available on line at the following URL: http://ymd.med.yale.edu/karma/cgi-bin/karma.pl.

Sperm-Induced Modification of the Oviductal Gene Expression Profile After Natural Insemination in Mice1

In mammals, the physiological interaction between spermatozoa and oviductal epithelia involves intimate and specific contact between the two cell types. Spermatozoa may undergo stringent selection processes within the female reproductive tract before they meet and fertilize oocytes. The physiological basis of the sperm selection process is largely unknown. Here we tested the hypothesis that the oviduct has a recognition system for spermatozoa that can detect the arrival of spermatozoa in the oviduct after insemination, resulting in alterations of the oviductal transcriptome. We initially performed a global screening of the oviductal transcriptome in mice 1) at the time of estrus (mating) and 2) 6 h after mating. Transcriptional alterations in the oviduct after mating were attributed to the presence of spermatozoa in the oviduct after mating and also to changes in the hormonal environment as female mice underwent the transition from estrus to diestrus. To distinguish these possibilities, female mice were then mated with T145H mutant mice, which because of spermatogenic arrest, produce seminal plasma but no spermatozoa. Focusing on two molecules that in the first experiment were upregulated after mating, it was found that adrenomedullin and prostaglandin endoperoxidase synthase 2 transcripts were upregulated in the oviducts of mice only after mating with fertile males; those mated with T145H infertile males showed significantly less response. These results indicate that it is the arrival of spermatozoa in the oviduct that activates one or more signal transduction pathways and leads to changes in the oviductal transcriptome profiles.

Gene expression alterations over large chromosomal regions in cancers include multiple genes unrelated to malignant progression

In solid tumors, the relationship between DNA copy number and global expression over large chromosomal regions has not been systematically explored. We used a 12,626-gene expression array analysis of head and neck squamous cell carcinoma and normal oral mucosa and annotated gene expression levels to specific chromosomal loci. Expression alterations correlated with reported data using comparative genomic hybridization. When genes with significant differences in expression between normal and malignant lesions, as defined by significance analysis of microarrays (SAM), were compared to nonsignificant genes, similar chromosomal patterns of alteration in expression were noted. Individual tumors underwent microsatellite analysis and chi(2) analysis of expression at 3p and 22q. Significant 3p underexpression and 22q overexpression were found in all primary tumors with 3p and 22q allelic imbalance, respectively, whereas no tumor without allelic imbalance on these chromosomal arms demonstrated expression differences. Loss and gain of chromosomal material in solid cancers can alter gene expression over large chromosomal regions, including multiple genes unrelated to malignant progression.

Gene expression profile analysis of 4-phenylbutyrate treatment of IB3-1 bronchial epithelial cell line demonstrates a major influence on heat-shock proteins

Most individuals with cystic fibrosis (CF) carry one or two mutations that result in a maturation defect of the full-length CFTR protein. The ΔF508 mutation results in a mutant protein that is degraded by the proteosome instead of progressing to the apical membrane where it functions as a cAMP-regulated chloride channel. 4-Phenylbutyrate (PBA) modulates heat-shock protein expression and promotes trafficking of ΔF508, thus permitting maturation and membrane insertion. The goal of this study was to gain insight into the genetic mechanism of PBA action through a large-scale analysis of gene expression. The Affymetrix genome-spanning U133 microarray set was used to compare mRNA expression levels in untreated IB3-1 cell line cultures with cultures treated with 1 mM PBA for 12 and 24 h. The most notable changes in mRNA levels were transient elevations in heat-shock proteins. The majority of genes downregulated throughout the application period were functionally associated with control of gene expression. Another set of genes increased in expression starting at 24 h, suggesting these are downstream effects of altered gene expression initiated by PBA. More than one-third of the genes in this late expressing set were identified as having potential significance in understanding the pathology of CF. Our results demonstrate the usefulness of gene expression profile analysis in understanding the consequences of PBA treatment and provide insights in how this drug exerts its effect on the trafficking of CFTR.

Microarray analysis of tumor necrosis factor alpha induced gene expression in U373 human glioblastoma cells

BACKGROUND

Tumor necrosis factor alpha (TNF) is able to induce a variety of biological responses in the nervous system including inflammation and neuroprotection. Human astrocytoma cells U373 have been widely used as a model for inflammatory cytokine actions in the nervous system. Here we used cDNA microarrays to analyze the time course of the transcriptional response from 1 h up to 12 h post TNF treatment in comparison to untreated U373 cells. TNF activated strongly the NF-kappaB transcriptional pathway and is linked to other pathways via the NF-kappaB target genes JUNB and IRF-1. Part of the TNF-induced gene expression could be inhibited by pharmacological inhibition of NF-kappaB with pyrrolidine-dithiocarbamate (PDTC). NF-kappaB comprises a family of transcription factors which are involved in the inducible expression of genes regulating neuronal survival, inflammatory response, cancer and innate immunity.

RESULTS

In this study we show that numerous genes responded to TNF (> 880 from 7500 tested) with a more than two-fold induction rate. Several novel TNF-responsive genes (about 60% of the genes regulated by a factor > or = 3) were detected. A comparison of our TNF-induced gene expression profiles of U373, with profiles from 3T3 and Hela cells revealed a striking cell-type specificity. SCYA2 (MCP-1, CCL2, MCAF) was induced in U373 cells in a sustained manner and at the highest level of all analyzed genes. MCP-1 protein expression, as monitored with immunofluorescence and ELISA, correlated exactly with microarray data. Based on these data and on evidence from literature we suggest a model for the potential neurodegenerative effect of NF-kappaB in astroglia: Activation of NF-kappaB via TNF results in a strongly increased production of MCP-1. This leads to a exacerbation of neurodegeneration in stoke or Multiple Sclerosis, presumably via infiltration of macrophages.

CONCLUSIONS

The vast majority of genes regulated more than 3-fold were previously not linked to tumor necrosis factor alpha as a search in published literature revealed. Striking co-regulation for several functional groups such as proteasome and ribosomal proteins were detected.

A design and statistical perspective on microarray gene expression studies in nutrition: the need for playful creativity and scientific hard-mindedness

OBJECTIVES

Our purpose is to highlight some of the past and potential future uses of microarray in nutrition research, while also commenting on some aspects of the design conduct and analysis of microarray data that will leave to improved data quality.

METHODS

In this review article we outline some of the aspects of microarray experimentation that must be considered before and during these experiments. These topics include: identification of the experiment's objective (hypothesis), the experimental design, sample size, statistical analysis, data verification, data handling, and experimental interpretation.

RESULTS

In order to illustrate the principles we outline in this article we use the methods to layout the design of a microarray experiment to study one aspect of the observation that a diet high in soy is associated with lower rates of breast cancer.

CONCLUSIONS

Microarrays are a very powerful tool for studying virtually every nutrition-related disease and trait and can provide valuable insights that are not obtainable with other techniques. However, unless nutrition researchers conduct their studies with scientific hard-mindedness, the studies will be of lower power at least if not completely misleading.

Global up-regulation of chromosome 21 gene expression in the developing Down syndrome brain

Down syndrome (DS) results from complete or partial triplication of human chromosome 21. It is assumed that the neurological and other symptoms are caused by the overexpression of genes on chromosome 21, but this hypothesis has not yet been assessed on a chromosome-wide basis. Here we show that expression of genes localized to chromosome 21 is globally up-regulated in human fetal trisomy 21 cases, both in cerebral cortex extracts and in astrocytic cell lines cultured from cerebral cortex. This abnormal regulation of gene expression is specific to chromosome 21. Our data describe transcriptional changes that are specific to many genes assigned to chromosome 21 and do not directly measure the clinical phenotype of DS. However, it is possible that these gene expression changes ultimately relate to the phenotypic variability of DS.

Integrative biology and the developing limb bud

The identification or selective construction of mutations within genes has allowed researchers to explore the downstream effects of gene disruption. Although these approaches have been successful, a limitation in our assessment of the consequences of conditional changes, and thereby our understanding of roles or function of genes, limits the degree to which we examine the effects of our manipulations. It is also clear that linear associations are incorrect models for describing development, and newer methods now give us an opportunity to practice an integrative biology. In our attempts to explore the consequences of Hoxa13 disruption in mice and humans, it has become clear that a better understanding of the consequences of gene alteration may be achievable by taking a broader approach with a long-term view. Fundamental questions regarding Hox gene function in vertebrates, including those related to the number of target genes; the degree of overlap of target gene regulation among paralogs; the magnitude of modulation exerted; and the identity of genes that are activated versus repressed need to be explored if a more thorough mechanistic understanding is to be achieved. To begin to address these questions, we undertook a comprehensive analysis of the expression of genes within developing limb buds of mice, and here we present some of our preliminary results. Our efforts will further (1) the exploration of the broader genetic relationships of expressed genes, (2) the determination of parallels or variations in target usage for a given gene in different tissues and between different organisms, (3) the evaluation of limb patterning mechanisms in other animal model systems, and (4) the exploration of gene expression hierarchies regulated by HOX proteins in developmental systems.

Altered patterns of cellular gene expression in dermal microvascular endothelial cells infected with Kaposi's sarcoma-associated herpesvirus

Kaposi's sarcoma (KS)-associated herpesvirus (KSHV; also called human herpesvirus 8) is believed to be the etiologic agent of Kaposi's sarcoma, multicentric Castleman's disease, and AIDS-associated primary effusion lymphoma. KSHV infection of human dermal microvascular endothelial cells (DMVEC) in culture results in the conversion of cobblestone-shaped cells to spindle-shaped cells, a characteristic morphological feature of cells in KS lesions. All spindle-shaped cells in KSHV-infected DMVEC cultures express the latency-associated nuclear protein LANA1, and a subfraction of these cells undergo spontaneous lytic cycle induction that can be enhanced by tetradecanoyl phorbol acetate (TPA) treatment. To study the cellular response to infection by KSHV, we used two different gene array screening systems to examine the expression profile of either 2,350 or 9,180 human genes in infected compared to uninfected DMVEC cultures in both the presence and absence of TPA. In both cases, between 1.4 and 2.5% of the genes tested were found to be significantly upregulated or downregulated. Further analysis by both standard and real-time reverse transcription-PCR procedures directly confirmed these results for 14 of the most highly upregulated and 13 of the most highly downregulated genes out of a total of 37 that were selected for testing. These included strong upregulation of interferon-responsive genes such as interferon response factor 7 (IRF7) and myxovirus resistance protein R1, plus upregulation of exodus 2 beta-chemokine, RDC1 alpha-chemokine receptor, and transforming growth factor beta3, together with strong downregulation of cell adhesion factors alpha(4)-integrin and fibronectin plus downregulation of bone morphogenesis protein 4, matrix metalloproteinase 2, endothelial plasminogen activator inhibitor 1, connective tissue growth factor, and interleukin-8. Significant dysregulation of several other cytokine-related genes or receptors, as well as endothelial cell and macrophage markers, and various other genes associated with angiogenesis or transformation was also detected. Western immunoblot and immunohistochemical analyses confirmed that the cellular IRF7 protein levels were strongly upregulated during the early lytic cycle both in KSHV-infected DMVEC and in the body cavity-based lymphoma BCBL1 PEL cell line.

Microarray analysis of differential gene expression in lead-exposed astrocytes

The toxic metal lead is a widespread environmental health hazard that can adversely affect human health. In an effort to better understand the cellular and molecular consequences of lead exposure, we have employed cDNA microarrays to analyze the effects of acute lead exposure on large-scale gene expression patterns in immortalized rat astrocytes. Our studies identified many genes previously reported to be differentially regulated by lead exposure. Additionally, we have identified novel putative targets of lead-mediated toxicity, including members of the family of calcium/phospholipid binding annexins, the angiogenesis-inducing thrombospondins, collagens, and tRNA synthetases. We demonstrate the ability to distinguish lead-exposed samples from control or sodium samples solely on the basis of large-scale gene expression patterns using two complementary clustering methods. We have confirmed the altered expression of candidate genes and their encoded proteins by RT-PCR and Western blotting, respectively. Finally, we show that the calcium-dependent phospholipid binding protein annexin A5, initially identified as a differentially regulated gene by our microarray analysis, is directly bound and activated by nanomolar concentrations of lead. We conclude that microarray technology is an effective tool for the identification of lead-induced patterns of gene expression and molecular targets of lead.

Gene expression profiling in postmortem Rett Syndrome brain: differential gene expression and patient classification

The identification of mutations in the transcriptional repressor methyl-CpG-binding protein 2 (MECP2) gene in Rett Syndrome (RTT) suggests that an inappropriate release of transcriptional silencing may give rise to RTT neuropathology. Despite this progress, the molecular basis of RTT neuropathogenesis remains unclear. Using multiple cDNA microarray technologies, subtractive hybridization, and conventional biochemistry, we generated comprehensive gene expression profiles of postmortem brain tissue from RTT patients and matched controls. Many glial transcripts involved in known neuropathological mechanisms were found to have increased expression in RTT brain, while decreases were observed in the expression of multiple neuron-specific mRNAs. Dramatic and consistent decreases in transcripts encoding presynaptic markers indicated a specific deficit in presynaptic development. Employing multiple clustering algorithms, it was possible to accurately segregate RTT from control brain tissue samples based solely on gene expression profile. Although previously achieved in cancers, our results constitute the first report of human disease classification using gene expression profiling in a complex tissue source such as brain.