Identifying biological themes within lists of genes with EASE

Gene expression profiles differentiate between sterile SIRS and early sepsis

INTRODUCTION

The systemic inflammatory response syndrome (SIRS) occurs frequently in critically ill patients and presents similar clinical appearances despite diverse infectious and noninfectious etiologies. Despite similar phenotypic expression, these diverse SIRS etiologies may induce divergent genotypic expressions. We hypothesized that gene expression differences are present between sepsis and uninfected SIRS prior to the clinical appearance of sepsis.

METHODS

Critically ill uninfected SIRS patients were followed longitudinally for the development of sepsis. All patients had whole blood collected daily for gene expression analysis by Affymetrix Hg_U133 2.0 Plus microarrays. SIRS patients developing sepsis were compared with those remaining uninfected for differences in gene expression at study entry and daily for 3 days prior to conversion to sepsis. Acceptance criteria for differentially expressed genes required: >1.2 median fold change between groups and significance on univariate and multivariate analysis. Differentially expressed genes were annotated to pathways using DAVID 2.0/EASE analysis.

RESULTS

A total of 12,782 (23.4%) gene probes were differentially expressed on univariate analysis 0 to 48 hours before clinical sepsis. 626 (1.1%) probes met acceptance criteria, corresponding to 459 unique genes, 65 (14.2%) down and 395 (85.8%) up expressed. These genes annotated to 10 pathways that functionally categorized to 4 themes involving innate immunity, cytokine receptors, T cell differentiation, and protein synthesis regulation.

CONCLUSIONS

Sepsis has a unique gene expression profile that is different from uninfected inflammation and becomes apparent prior to expression of the clinical sepsis phenotype.

HIV1 protease inhibitors selectively induce inflammatory chemokine expression in primary human osteoblasts

HIV-infected patients are at increased risk of decreased bone mineral density. Several studies have implicated antiretroviral therapy as a contributor to the decreased bone mineral density seen in treated HIV-1 patients. Whilst the exact molecular mechanisms underlying decreased bone density remain to be elucidated, inflammation has been postulated to be an important pathogenomic mechanism. In this study, we have explored primary human osteoblast gene expression in response to protease inhibitors (PIs), by oligonucleotide microarray analysis. A list of dysregulated genes, correlated with the inflammatory response, increased significantly after NFV and RTV exposure. Analysis of gene and protein expression determined a selectively increase of the pro-inflammatory cytokines monocyte chemoattractant protein (MCP)-1 and interleukin-8 (IL-8) following exposure to a pharmacological concentration of NFV and RTV. These data suggested that generation of local inflammatory cascades may contribute to the development of decreased bone mineral density in highly active antiretroviral therapy (HAART)-treated HIV patients.

HIV protease inhibitors selectively induce gene expression alterations associated with reduced calcium deposition in primary human osteoblasts

HIV-infected patients are at increased risk of decreased bone mineral density. Some studies have implicated antiretroviral therapy as a contributor to the decreased bone mineral density seen in treated HIV-1 patients. In this study we explore the interactions between protease inhibitors (PI) and primary human osteoblast gene expression, highlighting a group of dysregulated genes that potentially are key factors in reducing bone formation. Runx-2 mRNA expression, calcium deposition, and alkaline phosphatase (ALP) activity decreased significantly in human osteoblast cultures after exposure to the PIs nelfinavir (NFV) and indinavir (IDV). Saquinavir (SQV), ritonavir (RTV), indinavir (IDV), or nelfinavir (NFV) exposure induced significant changes in genotypic expression as assessed by gene-chip microarray analysis. The altered genes from each group were compared to each other and a list of 8 upregulated and 13 downregulated genes only after NFV and IDV exposure was identified. This set includes TIMP-3, which has previously been demonstrated to be involved in osteoblast differentiation and extracellular matrix development processes. Silencing TIMP-3 mRNA expression using siRNA duplexes enhanced calcium deposition and ALP activity significantly, even after exposure to NFV and IDV. Our data suggest a link between reduced osteoblastic phenotype and a group of 21 altered genes following NFV and IDV treatment, and also suggest TIMP-3 may be involved in the PI-induced inhibition of osteoblast function.

Gene expression polymorphism inDrosophilapopulations

Although changes in gene expression have long been recognized as critical to evolutionary processes, the extent of natural polymorphism in gene expression has yet to be assessed, thus opening a new area of active research. We present microarray and quantitative real-time polymerase chain reaction (RT-PCR) data from Cosmopolitan and Zimbabwe morphs of Drosophila melanogaster. These morphs provide a useful model for investigations into the incipient stages of speciation because Zimbabwe females tend to preferentially mate with their own males and discriminate against Cosmopolitan males, while Cosmopolitan females mate indiscriminately. We analysed expression profiles from heads of mated and nonmated females and identified 45 candidate genes whose expression levels were associated with the behavioural morphs and were modified by mating. Genes with altered transcription levels were randomly distributed across the genome and fell into diverse categories of biological activities. Several candidate genes, such as desaturase2 and Odorant receptor 63a, were additionally subjected to quantitative RT-PCR analysis. Notably, desaturase2, which has been invoked to play a role in sexual isolation between Cosmopolitan and Zimbabwe D. melanogaster/races/strains and predicted to be translational-inactive in Cosmopolitan due to a major deletion, was found to be up-regulated in Zimbabwe and down-regulated, but still expressed, in Cosmopolitan.

Changes in skeletal muscle gene expression following clenbuterol administration

Background

Beta-adrenergic receptor agonists (BA) induce skeletal muscle hypertrophy, yet specific mechanisms that lead to this effect are not well understood. The objective of this research was to identify novel genes and physiological pathways that potentially facilitate BA induced skeletal muscle growth. The Affymetrix platform was utilized to identify gene expression changes in mouse skeletal muscle 24 hours and 10 days after administration of the BA clenbuterol.

Results

Administration of clenbuterol stimulated anabolic activity, as indicated by decreased blood urea nitrogen (BUN; P < 0.01) and increased body weight gain (P < 0.05) 24 hours or 10 days, respectively, after initiation of clenbuterol treatment. A total of 22,605 probesets were evaluated with 52 probesets defined as differentially expressed based on a false discovery rate of 10%. Differential mRNA abundance of four of these genes was validated in an independent experiment by quantitative PCR. Functional characterization of differentially expressed genes revealed several categories that participate in biological processes important to skeletal muscle growth, including regulators of transcription and translation, mediators of cell-signalling pathways, and genes involved in polyamine metabolism.

Conclusion

Global evaluation of gene expression after administration of clenbuterol identified changes in gene expression and overrepresented functional categories of genes that may regulate BA-induced muscle hypertrophy. Changes in mRNA abundance of multiple genes associated with myogenic differentiation may indicate an important effect of BA on proliferation, differentiation, and/or recruitment of satellite cells into muscle fibers to promote muscle hypertrophy. Increased mRNA abundance of genes involved in the initiation of translation suggests that increased levels of protein synthesis often associated with BA administration may result from a general up-regulation of translational initiators. Additionally, numerous other genes and physiological pathways were identified that will be important targets for further investigations of the hypertrophic effect of BA on skeletal muscle.

Identification of a PAX-FKHR gene expression signature that defines molecular classes and determines the prognosis of alveolar rhabdomyosarcomas

Alveolar rhabdomyosarcomas (ARMS) are aggressive soft-tissue sarcomas affecting children and young adults. Most ARMS tumors express the PAX3-FKHR or PAX7-FKHR (PAX-FKHR) fusion genes resulting from the t(2;13) or t(1;13) chromosomal translocations, respectively. However, up to 25% of ARMS tumors are fusion negative, making it unclear whether ARMS represent a single disease or multiple clinical and biological entities with a common phenotype. To test to what extent PAX-FKHR determine class and behavior of ARMS, we used oligonucleotide microarray expression profiling on 139 primary rhabdomyosarcoma tumors and an in vitro model. We found that ARMS tumors expressing either PAX-FKHR gene share a common expression profile distinct from fusion-negative ARMS and from the other rhabdomyosarcoma variants. We also observed that PAX-FKHR expression above a minimum level is necessary for the detection of this expression profile. Using an ectopic PAX3-FKHR and PAX7-FKHR expression model, we identified an expression signature regulated by PAX-FKHR that is specific to PAX-FKHR-positive ARMS tumors. Data mining for functional annotations of signature genes suggested a role for PAX-FKHR in regulating ARMS proliferation and differentiation. Cox regression modeling identified a subset of genes within the PAX-FKHR expression signature that segregated ARMS patients into three risk groups with 5-year overall survival estimates of 7%, 48%, and 93%. These prognostic classes were independent of conventional clinical risk factors. Our results show that PAX-FKHR dictate a specific expression signature that helps define the molecular phenotype of PAX-FKHR-positive ARMS tumors and, because it is linked with disease outcome in ARMS patients, determine tumor behavior.

Conditional activation of RET/PTC3 and BRAFV600E in thyroid cells is associated with gene expression profiles that predict a preferential role of BRAF in extracellular matrix remodeling

Papillary thyroid cancers (PTC) are associated with nonoverlapping mutations of genes coding for mitogen-activated protein kinase signaling effectors (i.e., the TK receptors RET or NTRK and the signaling proteins RAS and BRAF). We examined the pattern of gene expression after activation of these oncoproteins in thyroid PCCL3 cells, with the goal of identifying pathways or gene subsets that may account for the phenotypic differences observed in human cancers. We hybridized cDNA from cells treated with or without doxycycline to induce expression of BRAF(V600E), RET/PTC3, or RET/PTC3 with small interfering RNA-mediated knockdown of BRAF, respectively, to slides arrayed with a rat 70-mer oligonucleotide library consisting of 27,342 oligos. Among the RET/PTC3-induced genes, 2,552 did not require BRAF as they were similarly regulated by RET/PTC3 with or without BRAF knockdown and not by expression of BRAF(V600E). Immune response and IFN-related genes were highly represented in this group. About 24% of RET/PTC3-regulated genes were BRAF dependent, as they were similarly modified by RET/PTC3 and BRAF(V600E) but not in cells expressing RET/PTC3 with knockdown of BRAF. A gene cluster coding for components of the mitochondrial electron transport chain pathway was down-regulated in this group, potentially altering regulation of cell viability. Metalloproteinases were also preferentially induced by BRAF, particularly matrix metalloproteinase 3 (MMP3), MMP9, and MMP13. Accordingly, conditional expression of BRAF was associated with markedly increased invasion into Matrigel compared with cells expressing RET/PTC3. The preferential induction of MMPs by BRAF could explain in part the more invasive behavior of thyroid cancers with BRAF mutations.

Gene expression differences in normal esophageal mucosa associated with regression and progression of mild and moderate squamous dysplasia in a high-risk Chinese population

A randomized, double-blinded, placebo-controlled 2 x 2 factorial chemoprevention trial was conducted in Linxian, China to assess the effects of selenomethionine and celecoxib on the natural history of esophageal squamous dysplasia. Results from this study indicated that asymptomatic adults with mild dysplasia were more likely to show an improvement when treated with selenomethionine compared with placebo (P = 0.02). Prompted by this finding, we examined the molecular profiles associated with regression and progression of dysplastic lesions in normal mucosa from 29 individuals, a subset of the Linxian cohort, using the Affymetrix U133A chip. Twenty differentially expressed genes were associated with regression and 129 were associated with progression when we compared the change in gene expression over time. Genes associated with immune response (n = 15), cell cycle (n = 15), metabolism (n = 15), calcium transport or calcium ion activity (n = 10), regulation of transcription (n = 9), signal transduction (n = 7), cytoskeleton and microtubules (n = 5), nucleotide processing and biosynthesis (n = 4), G-coupled signaling (n = 4), and apoptosis (n = 3) were present in the list of 149 genes. Using the Expression Analysis Systematic Explorer pathway analysis program, only the immune response pathway was significantly overrepresented among these 149 genes. Individuals whose lesions regressed seemed to have higher expression of genes associated with immune stimulation, such as antigen presentation, survival of T cells, and T-cell activation (HLA-DRA, HLA-DPA1, HLA-DBQ1, CD58, and FCER1A). In contrast, individuals whose lesions progressed had higher expression of genes involved in immune suppression and inflammation (CNR2, NFATC4, NFRKB, MBP, INHBB, CMKLR1, CRP, ORMS, SERPINA7, and SERPINA1). These data suggest that local and systemic immune responses may influence the natural history of esophageal squamous dysplasia.

Comprehensive quality control utilizing the prehybridization third-dye image leads to accurate gene expression measurements by cDNA microarrays

BACKGROUND

Gene expression profiling using microarrays has become an important genetic tool. Spotted arrays prepared in academic labs have the advantage of low cost and high design and content flexibility, but are often limited by their susceptibility to quality control (QC) issues. Previously, we have reported a novel 3-color microarray technology that enabled array fabrication QC. In this report we further investigated its advantage in spot-level data QC.

RESULTS

We found that inadequate amount of bound probes available for hybridization led to significant, gene-specific compression in ratio measurements, increased data variability, and printing pin dependent heterogeneities. The impact of such problems can be captured through the definition of quality scores, and efficiently controlled through quality-dependent filtering and normalization. We compared gene expression measurements derived using our data processing pipeline with the known input ratios of spiked in control clones, and with the measurements by quantitative real time RT-PCR. In each case, highly linear relationships (R2 > 0.94) were observed, with modest compression in the microarray measurements (correction factor < 1.17).

CONCLUSION

Our microarray analytical and technical advancements enabled a better dissection of the sources of data variability and hence a more efficient QC. With that highly accurate gene expression measurements can be achieved using the cDNA microarray technology.

The UVB-induced gene expression profile of human epidermis in vivo is different from that of cultured keratinocytes

In order to obtain a comprehensive picture of the molecular events regulating cutaneous photodamage of intact human epidermis, suction blister roofs obtained after a single dose of in vivo ultraviolet (UV)B exposure were used for microarray profiling. We found a changed expression of 619 genes. Half of the UVB-regulated genes had returned to pre-exposure baseline levels at 72 h, underscoring the transient character of the molecular cutaneous UVB response. Of special interest was our finding that several of the central p53 target genes remained unaffected following UVB exposure in spite of p53 protein accumulation. We next compared the in vivo expression profiles of epidermal sheets to that of cultured human epidermal keratinocytes exposed to UVB in vitro. We found 1931 genes that differed in their expression profiles between the two groups. The expression profile in intact epidemis was geared mainly towards DNA repair, whereas cultured keratinocytes responded predominantly by activating genes associated with cell-cycle arrest and apoptosis. These differences in expression profiles might reflect differences between mature differentiating keratinocytes in the suprabasal epidermal layers versus exponentially proliferating keratinocytes in cell culture. Our findings show that extreme care should be taken when extrapolating from findings based on keratinocyte cultures to changes in intact epidermis.

SOX4 expression in bladder carcinoma: clinical aspects and in vitro functional characterization

The human transcription factor SOX4 was 5-fold up-regulated in bladder tumors compared with normal tissue based on whole-genome expression profiling of 166 clinical bladder tumor samples and 27 normal urothelium samples. Using a SOX4-specific antibody, we found that the cancer cells expressed the SOX4 protein and, thus, did an evaluation of SOX4 protein expression in 2,360 bladder tumors using a tissue microarray with clinical annotation. We found a correlation (P < 0.05) between strong SOX4 expression and increased patient survival. When overexpressed in the bladder cell line HU609, SOX4 strongly impaired cell viability and promoted apoptosis. To characterize downstream target genes and SOX4-induced pathways, we used a time-course global expression study of the overexpressed SOX4. Analysis of the microarray data showed 130 novel SOX4-related genes, some involved in signal transduction (MAP2K5), angiogenesis (NRP2), and cell cycle arrest (PIK3R3) and others with unknown functions (CGI-62). Among the genes regulated by SOX4, 25 contained at least one SOX4-binding motif in the promoter sequence, suggesting a direct binding of SOX4. The gene set identified in vitro was analyzed in the clinical bladder material and a small subset of the genes showed a high correlation to SOX4 expression. The present data suggest a role of SOX4 in the bladder cancer disease.

An osteopontin splice variant induces anchorage independence in human breast cancer cells

In malignant tumors, metastasis genes are typically deregulated by aberrant expression or splicing. Osteopontin is expressed at high levels by various cancers and contributes importantly to their invasive potential. In contrast, osteopontin derived from host cells induces cellular immunity and could bolster antitumor protection by cytotoxic T lymphocytes. Here we show that breast cancer cells express multiple splice variants of osteopontin. According to RT-PCR analysis of human breast tissue specimens, the splice variant osteopontin-c is a highly specific marker for transformed cells, which is not expressed in their surrounding normal tissue. The full-length form of osteopontin aggregates in the presence of physiologic amounts of calcium and, in this state, leads to enhanced cell adhesion. Ostensibly, this effect is inhibitory for tumor cell dissemination. The shortest splice variant, osteopontin-c, does not aggregate in the presence of calcium and enhances clone formation in soft agar. According to microarray analysis, osteopontin-c induces the expression of oxidoreductases, consistent with protection from anoikis during anchorage-independent growth. These studies define a third functional domain of osteopontin, beside the C-terminal CD44-binding site and the central integrin-binding site. They also provide evidence for a bifunctional character of osteopontin, with the soluble form supporting invasiveness and the aggregated form promoting adhesion.

Inbreeding by environmental interactions affect gene expression in Drosophila melanogaster

Genomewide gene expression patterns were investigated in inbred and noninbred Drosophila melanogaster lines under benign and stressful (high temperature) environmental conditions in a highly replicated experiment using Affymetrix gene chips. We found that both heat-shock protein and metabolism genes are strongly affected by temperature stress and that genes involved in metabolism are differentially expressed in inbred compared with noninbred lines, and that this effect is accentuated after heat stress exposure. Furthermore we show that inbreeding and temperature stress cause increased between-line variance in gene expression patterns. We conclude that inbreeding and environmental stress both independently and synergistically affect gene expression patterns. Interactions between inbreeding and the environment are often observed at the phenotypic level and our results reveal some of the genes that are involved at the individual gene level. Our observation of several metabolism genes being differentially expressed in inbred lines and more so after exposure to temperature stress, together with lower fitness in the investigated inbred lines, supports the hypothesis that superiority of heterozygous individuals partly derives from increased metabolic efficiency.

Acute activation of hippocampal glucocorticoid receptors results in different waves of gene expression throughout time

Several aspects of hippocampal cell function are influenced by adrenal-secreted glucocorticoids in a delayed, genomic fashion. Previously, we used Serial Analysis of Gene Expression to identify glucocorticoid receptor (GR)-induced transcriptional changes in the hippocampus at a fixed time point. However, because changes in mRNA levels are transient and most likely precede the effects on hippocampal cell function, the aim of the current study was to assess the transcriptional changes in a broader time window by generating a time curve of GR-mediated gene expression changes. Therefore, we used rat hippocampal slices obtained from adrenalectomised rats, substituted in vivo with low corticosterone pellets, predominantly occupying the hippocampal mineralocorticoid receptors. To activate GR, slices were treated in vitro with a high (100 nM) dose of corticosterone and gene expression was profiled 1, 3 and 5 h after GR-activation. Using Affymetrix GeneChips, a striking pattern with different waves of gene expression was observed, shifting from exclusively down-regulated genes 1 h after GR-activation to both up and down regulated genes 3 h after GR-activation. After 5 h, the response was almost back to baseline. Additionally, real-time quantitative polymerase chain reaction was used for validation of a selection of responsive genes including genes involved in neurotransmission and synaptic plasticity such as the corticotropin releasing hormone receptor 1, monoamine oxidase A, LIMK1 and calmodulin 2. This permitted confirmation of GR-responsiveness of 15 out of 18 selected genes. In conclusion, direct activation of GR in hippocampal slices results in transient changes in gene expression. The pattern in which gene expression was modulated suggests that the fast genomic effects of glucocorticoids may be realised via transrepression, preceding a later wave of transactivation. Furthermore, we identified a number of interesting candidate genes which may underlie the glucocorticoid-mediated effects on hippocampal cell function.

Transcriptional profiling suggests that Barrett's metaplasia is an early intermediate stage in esophageal adenocarcinogenesis

To investigate the relationship between Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC), we determined gene expression profiles of discrete pathological stages of esophageal neoplasia using a sequence-verified human cDNA microarray. Fifty one RNAs, comprising 24 normal esophagi (NE), 18 BEs, and nine EACs were hybridized to cDNA microarrays. Five statistical analyses were used for the data analysis. Genes showing significantly different expression levels among the three sample groups were identified. Genes were grouped into functional categories based on the Gene Ontology Consortium. Surprisingly, the expression pattern of BE was significantly more similar to EAC than to NE, notwithstanding the known histopathologic differences between BE and EAC. The pattern of NE was clearly distinct from that of EAC. Thirty-six genes were the most differentially modulated, according to these microarray data, in BE-associated neoplastic progression. Twelve genes were significantly differentially expressed in cancer-associated BE's plus EAC (as a single combined tissue group) vs noncancer-associated BE's. These genes represent potential biomarkers to diagnose EAC at its early stages. Our results demonstrate that molecular events at the transcriptional level in BE are remarkably similar to BE's-associated adenocarcinoma of the esophagus. This finding alarmingly implies that BE is biologically closer to cancer than to normal esophagus, and that the cancer risk of BE is perhaps higher than we had imagined. These findings suggest that changes modulated at the molecular biologic level supervene earlier than histologic changes, and that BE is an early intermediate stage in the process of EAC.

Targeted disruption of glycerol kinase gene in mice: expression analysis in liver shows alterations in network partners related to glycerol kinase activity

Glycerol kinase deficiency (GKD) is an X-linked inborn error of metabolism with metabolic and neurological crises. Liver shows the highest level of glycerol kinase (GK) activity in humans and mice. Absence of genotype-phenotype correlations in patients with GKD indicates the involvement of modifier genes, including other network partners. To understand the molecular pathogenesis of GKD, we performed microarray analysis on liver mRNA from neonatal glycerol kinase (Gyk) knockout (KO) and wild-type (WT) mice. Unsupervised learning revealed that the overall gene expression profile of the KO mice was different from that of WT. Real-time PCR confirmed the differences for selected genes. Functional gene enrichment analysis was used to find 56 increased and 37 decreased gene functional categories. PathwayAssist analysis identified changes in gene expression levels of genes involved in organic acid metabolism indicating that GK was part of the same metabolic network which correlates well with the patients with GKD having metabolic acidemia during their episodic crises. Network component analysis (NCA) showed that transcription factors sterol regulatory element-binding protein (SREBP)-1c, carbohydrate response element-binding protein (ChREBP), hepatocyte nuclear factor-4 alpha (HNF-4alpha) and peroxisome proliferative-activated receptor-alpha (PPARalpha) had increased activity in the Gyk KO mice compared with WT mice, whereas SREBP-2 was less active in the Gyk KO mice. These studies show that Gyk deletion causes alterations in expression of genes in several regulatory networks and is the first time NCA has been used to expand on microarray data from a mouse KO model of a human disease.

Searching for biomarkers of developmental toxicity with microarrays: normal eye morphogenesis in rodent embryos

Gene expression arrays reveal the potential linkage of altered gene expression with specific adverse effects leading to disease phenotypes. But how closely do microarray data reflect early physiological or pharmacological measures that predict toxic event(s)? To explore this issue, we have undertaken experiments in early mouse embryos exposed to various teratogens during neurulation stages with the aim of correlating large-scale changes in gene expression across the critical period during exposure. This study reports some of the large-scale changes in gene expression that can be detected in the optic rudiment of the developing mouse and rat embryo across the window of development during which the eye is exceedingly sensitive to teratogen-induced micro-/anophthalmia. Microarray analysis was performed on RNA from the headfold or ocular region at the optic vesicle and optic cup stages when the ocular primordium is enriched for Pax-6, a master control gene for eye morphogenesis. Statistical selection of differentially regulated genes and various clustering techniques identified groups of genes in upward or downward trajectories in the normal optic primordium during early eye development in mouse and rat species. We identified 165 genes with significant differential expression during eye development, and a smaller subset of 58 genes that showed a tight correlation between mouse-rat development. Significantly over-represented functional categories included fatty acid metabolism (up-regulated) and glycolysis (down-regulated). From studies such as these that benchmark large-scale gene expression during normal embryonic development, we may be able to identify the panel of biomarkers that best correlate with species differences and the risks for developmental toxicity.

The transforming growth factor-beta family members bone morphogenetic protein-2 and macrophage inhibitory cytokine-1 as mediators of the antiangiogenic activity of N-(4-hydroxyphenyl)retinamide

PURPOSE

Tumor growth appears to be an angiogenesis-dependent process. N-(4-hydroxyphenyl)retinamide (fenretinide; 4HPR) has been found to inhibit and/or prevent tumor growth under diverse conditions. Although 4HPR is antiangiogenic, the molecular mechanisms of this effect remain largely unknown.

EXPERIMENTAL DESIGN

Endothelial cells were treated with 4HPR in vitro to study the effects on migration, invasion, and organization, as well as gene expression by microarray and quantitative PCR studies. In vivo angiogenesis was evaluated in the Matrigel model.

RESULTS

4HPR treatment substantially modified the biological activities of endothelial cells, repressing their capacity to migrate, invade, and organize into capillary-like structures. The inhibition of invasion induced by 4HPR was also associated with decreased activities of the metalloproteases matrix metalloproteinase-2 and CD13/APN. Using oligonucleotide microarrays, we observed that bone morphogenetic protein-2 and macrophage inhibitory cytokine-1, two multifunctional cytokines of the transforming growth factor-beta family that regulate the growth, differentiation, apoptosis, and matrix accumulation of a variety of cells, are up-regulated in vitro by 4HPR. Both these molecules specifically inhibited endothelial cell growth, migration, and invasion in vitro and suppressed angiogenesis in the Matrigel plug assay in vivo. Blocking antibodies to bone morphogenetic protein-2 were able to reverse the suppressive effects of 4HPR in vitro and in vivo.

CONCLUSIONS

These data support the conclusion that 4HPR inhibits tumor growth by repression of new vessel growth and identify novel points of regulation of angiogenesis in transforming growth factor-beta family proteins.

Effects of chlorpromazine with and without UV irradiation on gene expression of HepG2 cells

Damage to DNA can trigger a variety of stress-related signals that alter the expression of genes associated with numerous biological pathways. In this study, we have used a cDNA microarray representing 1089 genes related to DNA damage and repair, cell cycle, transcription, metabolism and other toxicologically important cell functions to identify genes regulated in response to DNA damage in HepG2 cells induced by UV-activated chlorpromazine (CPZ). CPZ itself is not genotoxic but, upon UV irradiation with a non-genotoxic dose in the UVA range, it produces reactive free radical intermediates with DNA damaging properties. Genotoxicity in HepG2 cells was assessed concomitantly to gene expression profiling using the Comet assay. Kinetic studies were performed at a non-cytotoxic but clearly photogenotoxic concentration of CPZ (1.25 microg/ml) to characterize gene expression profiles at four different time points (3, 7, 15, 23 h) post short-term treatment. The results obtained from repeated experiments display a time-dependent expression pattern of up-regulated and repressed genes with distinct peaks in the number of differentially expressed genes at the 7 and 23 h time points. Most of the genes with altered expression belonged to the functional categories of cell cycle regulation and cell proliferation. A comparison with published expression profiles established in response to other genotoxic compounds showed low levels of concordance, which is most likely caused by the fact that extremely different testing conditions were used.

A microarray study of post-mortem mRNA degradation in mouse brain tissue

BACKGROUND

Although there is evidence that post-mortem interval (PMI) is not a major contributor to reduced overall RNA integrity, it may differentially affect a subgroup of gene transcripts that are susceptible to PMI-related degradation. This would particularly have ramifications for microarray studies that include a broad spectrum of genes.

METHOD

Brain tissue was removed from adult mice at 0, 6, 12, 18, 24, 36 and 48 h post-mortem. RNA transcript abundance was measured by hybridising RNA from the zero time point with test RNA from each PMI time point, and differential gene expression was assessed using cDNA microarrays. Sequence and ontological analyses were performed on the group of RNA transcripts showing greater than two-fold reduction.

RESULTS

Increasing PMI was associated with decreased tissue pH and increased RNA degradation as indexed by 28S/18S ribosomal RNA ratio. Approximately 12% of mRNAs detected on the arrays displayed more than a two-fold decrease in abundance by 48 h post-mortem. An analysis of nucleotide composition provided evidence that transcripts with the AUUUA motif in the 3' untranslated region (3'UTR) were more susceptible to PMI-related RNA degradation, compared to transcripts not carrying the 3'UTR AUUUA motif. Consistent with this finding, ontological analysis showed transcription factors and elements to be over-represented in the group of transcripts susceptible to degradation.

CONCLUSION

A subgroup of mammalian mRNA transcripts are particularly susceptible to PMI-related degradation, and as a group, they are more likely to carry the 3'UTR AUUUA motif. PMI should be controlled for in human and animal model post-mortem brain studies, particularly those including a broad spectrum of mRNA transcripts.

Computational systems analysis of developmental toxicity: design, development and implementation of a Birth Defects Systems Manager (BDSM)

Birth defects and developmental disabilities remain an important public health issue worldwide. With the availability of genomic sequences from a growing number of human and model organisms and the rapid expansion of the public repositories holding large-scale gene expression datasets, a computational systems analysis of developmental toxicology can incorporate this vast digital information toward the realization of predictive models for complex disease. Here we describe the initial design, development and implementation of a Birth Defects Systems Manager (BDSM). The project was motivated by the need for a computational-bioinformatics infrastructure to manage vast digital information from functional genomics and for a new knowledge environment specifically engineered for the analysis of developmental processes and toxicities. Proof-of-concept tested BDSM using meta-analysis of gene expression data collected from different laboratories, technology platforms, and study models. The composite dataset incorporated 232 microarray comparisons of RNA samples by single or dual microarray platforms, cDNA or oligonucleotide based probes, and human or mouse sequence information. Preliminary results identified system-level features in the embryonic transcriptome as it reacted to various developmental-teratological stimuli. BDSM is open access through the worldwide web (http://systemsanalysis.louisville.edu/) and can be integrated with other bioinformatics tools and resources to advance the pace of discovery in birth defects research.

Molecular signature of retinoic acid treatment in acute promyelocytic leukemia

Acute promyelocytic leukemia (APL) is a distinct subtype of acute myeloid leukemia characterized by a block of differentiation at the promyelocytic stage. APL patients respond to pharmacological concentrations of all-trans retinoic acid (RA) and disease remission correlates with terminal differentiation of leukemic blasts. The PML/RAR oncogenic transcription factor is responsible for both the pathogenesis of APL and for its sensitivity to RA. In order to identify physiological targets of RA therapy, we analysed gene expression profiles of RA-treated APL blasts and found 1056 common target genes. Comparing these results to those obtained in RA-treated U937 cell lines revealed that transcriptional response to RA is largely dependent on the expression of PML/RAR. Several genes involved in the control of differentiation and stem cell renewal are early targets of RA regulation, and may be important effectors of RA response. Modulation of chromatin modifying genes was also observed, suggesting that specific structural changes in local chromatin domains may be required to promote RA-mediated differentiation. Computational analysis of upstream genomic regions in RA target genes revealed nonrandom distribution of transcription factor binding sites, indicating that specific transcriptional regulatory complexes may be involved in determining RA response.

Large-scale identification of proteins expressed in mouse embryonic stem cells

A protein subset expressed in the mouse embryonic stem (ES) cell line, E14-1, was characterized by mass spectrometry-based protein identification technology and data analysis. In total, 1790 proteins including 365 potential nuclear and 260 membrane proteins were identified from tryptic digests of total cell lysates. The subset contained a variety of proteins in terms of physicochemical characteristics, subcellular localization, and biological function as defined by Gene Ontology annotation groups. In addition to many housekeeping proteins found in common with other cell types, the subset contained a group of regulatory proteins that may determine unique ES cell functions. We identified 39 transcription factors including Oct-3/4, Sox-2, and undifferentiated embryonic cell transcription factor I, which are characteristic of ES cells, 88 plasma membrane proteins including cell surface markers such as CD9 and CD81, 44 potential proteinaceous ligands for cell surface receptors including growth factors, cytokines, and hormones, and 100 cell signaling molecules. The subset also contained the products of 60 ES-specific and 41 stemness genes defined previously by the DNA microarray analysis of Ramalho-Santos et al. (Ramalho-Santos et al., Science 2002, 298, 597-600), as well as a number of components characteristic of differentiated cell types such as hematopoietic and neural cells. We also identified potential post-translational modifications in a number of ES cell proteins including five Lys acetylation sites and a single phosphorylation site. To our knowledge, this study provides the largest proteomic dataset characterized to date for a single mammalian cell species, and serves as a basic catalogue of a major proteomic subset that is expressed in mouse ES cells.

Early antiretroviral therapy for simian immunodeficiency virus infection leads to mucosal CD4+ T-cell restoration and enhanced gene expression regulating mucosal repair and regeneration

Simian immunodeficiency virus (SIV) and human immunodeficiency virus (HIV) infections lead to rapid depletion of CD4(+) T cells from gut-associated lymphoid tissue (GALT). Although the administration of antiretroviral therapy (ART) has been shown to increase CD4(+) T-cell levels in the peripheral blood in both SIV and HIV infections, its efficacy in restoring intestinal mucosal CD4(+) T cells has not been well investigated. To gain insights into the molecular mechanisms of virally induced disruptions in the mucosal immune system, we have evaluated longitudinal changes in viral burden, T-cell subsets, and mucosal gene expression profiles in SIV-infected rhesus macaques in the absence or presence of ART. Our results demonstrate a dramatic suppression of mucosal viral loads and rapid reconstitution of CD4(+) T cells in GALT in animals receiving ART that were not observed in untreated SIV-infected animals. DNA microarray-based gene expression profiling indicated that CD4(+) T-cell restoration in GALT was associated with up regulation of growth factors and genes involved in repair and regeneration of the mucosal epithelium. In contrast, untreated SIV-infected animals increased expression of lymphocyte activation and inflammatory response-associated genes and did not up regulate mucosal growth and repair associated transcription. In conclusion, these data indicate that initiating ART in primary SIV infection may lead to the restoration of the mucosal immune system through reduction of inflammation and promotion of epithelial repair in the intestinal mucosa.

The in vitro and in vivo effects of re-expressing methylated von Hippel-Lindau tumor suppressor gene in clear cell renal carcinoma with 5-aza-2'-deoxycytidine

PURPOSE

Clear cell renal carcinoma (ccRCC) is strongly associated with loss of the von Hippel-Lindau (VHL) tumor suppressor gene. The VHL gene is functionally lost through hypermethylation in up to 19% of sporadic ccRCC cases. We theorized that re-expressing VHL silenced by methylation in ccRCC cells, using a hypo-methylating agent, may be an approach to treatment in patients with this type of cancer. We test the ability of two hypo-methylating agents to re-express VHL in cell culture and in mice bearing human ccRCC and evaluate the effects of re-expressed VHL in these models.

EXPERIMENTAL DESIGN

Real-time reverse transcription-PCR was used to evaluate the ability of zebularine and 5-aza-2'-deoxycytidine (5-aza-dCyd) to re-express VHL in four ccRCC cell lines with documented VHL gene silencing through hypermethylation. We evaluated if the VHL re-expressed after hypo-methylating agent treatment could recreate similar phenotypic changes in ccRCC cells observed when the VHL gene is re-expressed via transfection in cell culture and in a xenograft mouse model. Finally we evaluate global gene expression changes occurring in our cells, using microarray analysis.

RESULTS

5-Aza-dCyd was able to re-express VHL in our cell lines both in culture and in xenografted murine tumors. Well described phenotypic changes of VHL expression including decreased invasiveness into Matrigel, and decreased vascular endothelial growth factor and glucose transporter-1 expression were observed in the treated lines. VHL methylated ccRCC xenografted tumors were significantly reduced in size in mice treated with 5-aza-dCyd. Mice bearing nonmethylated but VHL-mutated tumors showed no tumor shrinkage with 5-aza-dCyd treatment.

CONCLUSION

Hypo-methylating agents may be useful in the treatment of patients having ccRCC tumors consisting of cells with methylated VHL.

Autocrine activation of an osteopontin-CD44-Rac pathway enhances invasion and transformation by H-RasV12

Activated forms of Ras family members are prevalent in many cancers where Ras mutants transduce signals essential for transformation, angiogenesis, invasion and metastasis. As a cancer progression model, we used NIH3T3 cells to explore the mechanism of Ras-induced tumorigenesis. Ras family mutants H-RasV12 and Rit79L strongly induced foci formation, while Rho family mutants RhoA-QL, Rac1-QL and Cdc42-QL were less effective. A comparison of downstream transcriptional targets of Ras and Rho family members using a 26 383 element cDNA microarray revealed that the osteopontin (OPN) gene exhibited the best correlation between magnitude of gene expression change and level of foci formation (r=0.96, P<0.001). In association with H-RasV12- and Rit79L-mediated transformation, foci secreted OPN protein and upregulated the OPN receptor CD44, suggesting the novel initiation of an aberrant OPN-CD44-Rac autocrine pathway. In support of this were the following observations. First, RGD-deficient OPN protein-binding activity was present in H-RasV12-transformed cells but not in control cells, and binding activity was inhibited by the CD44 blocking antibody. Second, foci formation, cell invasion and Rac activity were induced by H-RasV12 and inhibited by the CD44 blocking antibody. Third, foci formation by H-RasV12 was substantially reduced by a short interfering RNA (siRNA) specifically targeting OPN expression for knockdown. Fourth, H-RasV12-mediated transformation was not blocked by the GRGDS peptide, suggesting that OPN effects were not mediated by the integrins. Lastly, OPN knockdown affected the downstream expression of 160 '2nd tier' genes, and at least a subset of these genes appears to be involved in transformation. Indeed, four genes were selected for knockdown, each resulting in a disruption of foci formation and/or invasion. These results underscore the role of aberrant autocrine signaling and transcriptional networking during tumorigenesis.

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.

Proteasome gene upregulation: a possible mechanism for intestinal adaptation

BACKGROUND/PURPOSE

The mechanisms that control intestinal adaptation remain unknown. To better understand the adaptive process, microarray technology was used to analyze gene expression in a rat model of intestinal adaptation.

METHODS

Adult male Sprague-Dawley rats underwent either a massive small bowel resection (70%) with anastomosis or a sham operation with small bowel transection and reanastomosis. After 21 days, ileal mucosa RNA was extracted. Individual RNA samples (n = 5 per group) were labeled and hybridized to 10 separate RAE 230A rat GeneChips. The signal values were calculated and the 2 groups were compared using a t test with the multiple testing correction of Benjamini and Hochberg (false discovery rate of 10%). Probe sets were analyzed for overrepresented physiologic pathways using Expression Analysis Systematic Explorer (EASE).

RESULTS

Of the 15,866 probe sets on the RAE 230A GeneChip, 5437 probe sets were unexpressed and excluded. Of the remaining 10,429 probe sets, several overrepresented pathways (EASE score <0.01 after Bonferroni correction) were identified. Further analysis revealed that 13 probe sets related to proteasome degradation (an enzyme complex implicated in the regulation of cell proliferation) were significantly upregulated in the intestinal adaptation group compared to the sham group.

CONCLUSIONS

Proteasomes may play a critical role in regulating the proliferation of intestinal mucosa during intestinal adaptation.

Ras Inhibition in Glioblastoma Down-regulates Hypoxia-Inducible Factor-1α, Causing Glycolysis Shutdown and Cell Death

Active Ras and phosphatidylinositol-3-kinase–dependent pathways contribute to the malignant phenotype of glioblastoma multiformes (GBM). Here we show that the Ras inhibitor trans-farnesylthiosalicylic acid (FTS) exhibits profound antioncogenic effects in U87 GBM cells. FTS inhibited active Ras and attenuated Ras signaling to extracellular signal-regulated kinase, phosphatidylinositol-3-kinase, and Akt. Concomitantly, hypoxia-inducible factor 1α (HIF-1α) disappeared, expression of key glycolysis pathway enzymes and of other HIF-1α–regulated genes (including vascular endothelial growth factor and the Glut-1 glucose transporter) was down-regulated, and glycolysis was halted. This led to a dramatic reduction in ATP, resulting in a severe energy crisis. In addition, the expression of E2F-regulated genes was down-regulated in the FTS-treated cells. Consequently, U87 cell growth was arrested and the cells died. These results show that FTS is a potent down-regulator of HIF-1α and might therefore block invasiveness, survival, and angiogenesis in GBM.

Significance analysis of functional categories in gene expression studies: a structured permutation approach

MOTIVATION

In high-throughput genomic and proteomic experiments, investigators monitor expression across a set of experimental conditions. To gain an understanding of broader biological phenomena, researchers have until recently been limited to post hoc analyses of significant gene lists.

METHOD

We describe a general framework, significance analysis of function and expression (SAFE), for conducting valid tests of gene categories ab initio. SAFE is a two-stage, permutation-based method that can be applied to various experimental designs, accounts for the unknown correlation among genes and enables permutation-based estimation of error rates.

RESULTS

The utility and flexibility of SAFE is illustrated with a microarray dataset of human lung carcinomas and gene categories based on Gene Ontology and the Protein Family database. Significant gene categories were observed in comparisons of (1) tumor versus normal tissue, (2) multiple tumor subtypes and (3) survival times.

AVAILABILITY

Code to implement SAFE in the statistical package R is available from the authors.

SUPPLEMENTARY INFORMATION

http://www.bios.unc.edu/~fwright/SAFE.

cDNA microarray analysis identifies genes induced in common by peptide growth factors and androgen in human prostate epithelial cells

Prostate cancer cells initially require androgen for continued proliferation, but invariably become androgen independent or unresponsive and recur after treatment by androgen ablation. Exploitation of common signaling components downstream of their specific receptors (i.e., androgen receptor (AR), insulin-like growth factor 1 (IGF-1) receptor, and epidermal growth factor (EGF) receptor) could provide a mechanism by which androgen independent cells survive and proliferate. Our objective was to design and implement prostate enriched cDNA microarrays to identify genes induced in prostate epithelial cells in a similar temporal pattern by both androgen and IGF or EGF. AR positive and AR negative human prostate epithelial cells of the M12 line were exposed in parallel to DHT, EGF, or IGF for 0, 6, or 24 h. RNA extracted from each of these groups was analyzed by cDNA microarrays composed of a unique set of 6373 prostate-derived cDNA clones from the Prostate Expression Database (PEDB). We observed statistically significant changes in 20 genes induced in common after 6 and 24 h exposure to androgen or these growth factors, and validated the microarray results by RT-PCR for three or four of these genes: v-myc, isocitrate dehydrogenase, and calnexin. Androgen response element binding motifs were identified in the upstream sequence in 16 of these 20 genes. These results provide comprehensive and unique insights into potential mechanisms by which peptide growth factors provide alternate pathways to control prostate epithelial cell proliferation in malignant states.

From mice to humans: identification of commonly deregulated genes in mammary cancer via comparative SAGE studies

Genetically engineered mouse mammary cancer models have been used over the years as systems to study human breast cancer. However, much controversy exists on the utility of such models as valid equivalents to the human cancer condition. To perform an interspecies gene expression comparative study in breast cancer we used a mouse model that most closely resembles human breast carcinogenesis. This system relies on the transplant of p53 null mouse mammary epithelial cells into the cleared mammary fat pads of syngeneic hosts. Serial analysis of gene expression (SAGE) was used to obtain gene expression profiles of normal and tumor samples from this mouse mammary cancer model (>300,000 mouse mammary-specific tags). The resulting mouse data were compared with 25 of our human breast cancer SAGE libraries (>2.5 million human breast-specific tags). We observed significant similarities in the deregulation of specific genes and gene families when comparing mouse with human breast cancer SAGE data. A total of 72 transcripts were identified as commonly deregulated in both species. We observed a systematic and significant down-regulation in all of the tumors from both species of various cytokines, including CXCL1 (GRO1), LIF, interleukin 6, and CCL2. All of the mouse and most human mammary tumors also displayed decreased expression of genes known to inhibit cell proliferation, including NFKBIA (IKBalpha), GADD45B, and CDKN1A (p21); transcription-related genes such as CEBP, JUN, JUNB, and ELF1; and apoptosis-related transcripts such as IER3 and GADD34/PPP1R15A. Examples of overexpressed transcripts in tumors from both species include proliferation-related genes such as CCND1, CKS1B, and STMN1 (oncoprotein 18); and genes related to other functions such as SEPW1, SDFR1, DNCI2, and SP110. Importantly, abnormal expression of several of these genes has not been associated previously with breast cancer. The consistency of these observations was validated in independent mouse and human mammary cancer sets. This is the first interspecies comparison of mammary cancer gene expression profiles. The comparative analysis of mouse and human SAGE mammary cancer data validates this p53 null mouse tumor model as a useful system closely resembling human breast cancer development and progression. More importantly, these studies are allowing us to identify relevant biomarkers of potential use in human studies while leading to a better understanding of specific mechanisms of human breast carcinogenesis.

Local correlation of expression profiles with gene annotations—proof of concept for a general conciliatory method

MOTIVATION

Integrated analysis of expression data and gene ontology annotations is a prime example of biological data that need co-explanatory interpretation. This particular application is used to validate a new method for integrated analysis of varied biological information.

RESULTS

The proposed method consists of determining local correlation coefficients and the corresponding P-values calculated per biological entity. This measure considers the combined intensity and significance of the agreement or disagreement, between two data sources about the same biological entity. The method is applied to the integrated analysis of gene expression and annotation of two gene sets, one from yeast and other from mouse. The potential of the method to generate accurate mechanistic hypothesis is also demonstrated. Specially, negative correlation results pose a new kind of biological hypothesis. Method performance was compared with annotation enrichment methods, and optimal conditions for the superiority of local correlation results are discussed.

Complementary roles for toxicologic pathology and mathematics in toxicogenomics, with special reference to data interpretation and oscillatory dynamics

Toxicogenomics is an emerging multidisciplinary science that will profoundly impact the practice of toxicology. New generations of biologists, using evolving toxicogenomics tools, will generate massive data sets in need of interpretation. Mathematical tools are necessary to cluster and otherwise find meaningful structure in such data. The linking of this structure to gene functions and disease processes, and finally the generation of useful data interpretation remains a significant challenge. The training and background of pathologists make them ideally suited to contribute to the field of toxicogenomics, from experimental design to data interpretation. Toxicologic pathology, a discipline based on pattern recognition, requires familiarity with the dynamics of disease processes and interactions between organs, tissues, and cell populations. Optimal involvement of toxicologic pathologists in toxicogenomics requires that they communicate effectively with the many other scientists critical for the effective application of this complex discipline to societal problems. As noted by Petricoin III et al (Nature Genetics 32, 474-479, 2002), cooperation among regulators, sponsors and experts will be essential for realizing the potential of microarrays for public health. Following a brief introduction to the role of mathematics in toxicogenomics, "data interpretation" from the perspective of a pathologist is briefly discussed. Based on oscillatory behavior in the liver, the importance of an understanding of mathematics is addressed, and an approach to learning mathematics "later in life" is provided. An understanding of pathology by mathematicians involved in toxicogenomics is equally critical, as both mathematics and pathology are essential for transforming toxicogenomics data sets into useful knowledge.

Gene expression profiling of advanced ovarian cancer: characterization of a molecular signature involving fibroblast growth factor 2

Epithelial ovarian cancer (EOC) is the gynecological disease with the highest death rate. We applied an automatic class discovery procedure based on gene expression profiling to stages III-IV tumors to search for molecular signatures associated with the biological properties and progression of EOC. Using a complementary DNA microarray containing 4451 cancer-related, sequence-verified features, we identified a subset of EOC characterized by the expression of numerous genes related to the extracellular matrix (ECM) and its remodeling, along with elements of the fibroblast growth factor 2 (FGF2) signaling pathway. A total of 10 genes were validated by quantitative real-time polymerase chain reaction, and coexpression of FGF2 and fibroblast growth factor receptor 4 in tumor cells was revealed by immunohistochemistry, confirming the reliability of gene expression by cDNA microarray. Since the functional relationships among these genes clearly suggested involvement of the identified molecular signature in processes related to epithelial-stromal interactions and/or epithelial-mesenchymal cellular plasticity, we applied supervised learning analysis on ovarian-derived cell lines showing distinct cellular phenotypes in culture. This procedure enabled construction of a gene classifier able to discriminate mesenchymal-like from epithelial-like cells. Genes overexpressed in mesenchymal-like cells proved to match the FGF2 signaling and ECM molecular signature, as identified by unsupervised class discovery on advanced tumor samples. In vitro functional analysis of the cell plasticity classifier was carried out using two isogenic and immortalized cell lines derived from ovarian surface epithelium and displaying mesenchymal and epithelial morphology, respectively. The results indicated the autocrine, but not intracrine stimulation of mesenchymal conversion and cohort/scatter migration of cells by FGF2, suggesting a central role for FGF2 signaling in the maintenance of cellular plasticity of ovary-derived cells throughout the carcinogenesis process. These findings raise mechanistic hypotheses on EOC pathogenesis and progression that might provide a rational underpinning for new therapeutic modalities.

Distinctive gene expression profiles of CD34 cells from patients with myelodysplastic syndrome characterized by specific chromosomal abnormalities

Aneuploidy, especially monosomy 7 and trisomy 8, is a frequent cytogenetic abnormality in the myelodysplastic syndromes (MDSs). Patients with monosomy 7 and trisomy 8 have distinctly different clinical courses, responses to therapy, and survival probabilities. To determine disease-specific molecular characteristics, we analyzed the gene expression pattern in purified CD34 hematopoietic progenitor cells obtained from MDS patients with monosomy 7 and trisomy 8 using Affymetrix GeneChips. Two methods were employed: standard hybridization and a small-sample RNA amplification protocol for the limited amounts of RNA available from individual cases; results were comparable between these 2 techniques. Microarray data were confirmed by gene amplification and flow cytometry using individual patient samples. Genes related to hematopoietic progenitor cell proliferation and blood cell function were dysregulated in CD34 cells of both monosomy 7 and trisomy 8 MDS. In trisomy 8, up-regulated genes were primarily involved in immune and inflammatory responses, and down-regulated genes have been implicated in apoptosis inhibition. CD34 cells in monosomy 7 showed up-regulation of genes inducing leukemia transformation and tumorigenesis and apoptosis and down-regulation of genes controlling cell growth and differentiation. These results imply distinct molecular mechanisms for monosomy 7 and trisomy 8 MDS and implicate specific pathogenic pathways.

Molecular pathogenesis of uterine smooth muscle tumors from transcriptional profiling

Uterine smooth muscle tumors range from the very common benign leiomyoma to the uncommon, but frequently lethal, leiomyosarcoma. Morphological and clinical differences between these tumors are presumed to result from differences in gene expression. To test this hypothesis, RNAs from four normal uterine myometria, seven uterine leiomyomas, and nine uterine leiomyosarcomas were profiled using microarrays of oligonucleotides representing about 7,000 unique probe sets. RNAs whose levels distinguished any of the three sample types were selected by analysis of variance (ANOVA). The 153 (2.2% of the total) probe sets representing 146 unique genes with the highest test statistic selected for further analysis met minimum ratio and range thresholds between groups. Cluster analysis distinguished benign and malignant samples at the first node, and myometrium and leiomyoma were resolved in a secondary node. Downregulation of specific genes in uterine leiomyosarcoma was the most common pattern of differential gene expression selected by the three-way ANOVA. Four extrauterine leiomyosarcomas had profiles most similar to that of the uterine leiomyosarcomas. Functional analysis of the 146 genes did not reveal any strong biological theme. These genes were distributed throughout the genome, but there was slight overrepresentation of genes on 1p and 2q. These genes define a tumor signature for uterine smooth muscle neoplasia, and they suggest that the molecular pathways in leiomyoma and leiomyosarcoma are distinct.

Impaired expression and insulin-stimulated phosphorylation of Akt-2 in muscle of obese patients with atypical diabetes

Although a pharmacological dose of insulin produces a dramatic increase in phosphorylation and activity of Akt isoforms 1 and 2 in mammalian skeletal muscle, few studies have examined the effect of physiological concentrations of insulin on the phosphorylation of Akt-1 and -2 in normal and diabetic tissue. This study examined the patterns of insulin-stimulated Akt isoform phosphorylation and protein expression in muscle biopsies obtained from obese patients with atypical diabetes immediately after a hyperglycemic crisis and again after near-normoglycemic remission. In obese patients with new-onset diabetes mellitus presenting with hyperglycemic crisis (plasma glucose 30.5 +/- 4.8 mM), in vitro stimulation of vastus lateralis muscle biopsies with 100 microU/ml (0.6 nM) insulin increased insulin receptor phosphorylation threefold and Akt-1 phosphorylation on Ser(473) twofold, whereas Akt-2 phosphorylation was not stimulated. After 10-wk intensive insulin therapy that led to near-normoglycemic remission and discontinuation of insulin therapy, both Akt-2 expression and insulin-stimulated Akt-2 Ser(474) phosphorylation doubled. Hyperglycemic crisis did not affect insulin-stimulated threonine phosphorylation of either Akt-1 or Akt-2. The decreased Akt-2 expression at presentation was accompanied by reduced GLUT4 protein expression and increased expression of enzymes counterregulatory to insulin action. Thus a physiological concentration of insulin stimulated Akt-1 and Akt-2 phosphorylation in human skeletal muscle in the absence of hyperglycemia, but Akt-2 expression and stimulation appeared to be impaired in muscle of obese patients with atypical diabetes presenting with severe hyperglycemia.

DAVID: Database for Annotation, Visualization, and Integrated Discovery

BACKGROUND

Functional annotation of differentially expressed genes is a necessary and critical step in the analysis of microarray data. The distributed nature of biological knowledge frequently requires researchers to navigate through numerous web-accessible databases gathering information one gene at a time. A more judicious approach is to provide query-based access to an integrated database that disseminates biologically rich information across large datasets and displays graphic summaries of functional information.

RESULTS

Database for Annotation, Visualization, and Integrated Discovery (DAVID; http://www.david.niaid.nih.gov) addresses this need via four web-based analysis modules: 1) Annotation Tool - rapidly appends descriptive data from several public databases to lists of genes; 2) GoCharts - assigns genes to Gene Ontology functional categories based on user selected classifications and term specificity level; 3) KeggCharts - assigns genes to KEGG metabolic processes and enables users to view genes in the context of biochemical pathway maps; and 4) DomainCharts - groups genes according to PFAM conserved protein domains.

CONCLUSIONS

Analysis results and graphical displays remain dynamically linked to primary data and external data repositories, thereby furnishing in-depth as well as broad-based data coverage. The functionality provided by DAVID accelerates the analysis of genome-scale datasets by facilitating the transition from data collection to biological meaning.