iBBiG: iterative binary bi-clustering of gene sets.

MOTIVATION: Meta-analysis of genomics data seeks to identify genes associated with a biological phenotype across multiple datasets; however, merging data from different platforms by their features (genes) is challenging. Meta-analysis using functionally or biologically characterized gene sets simplifies data integration is biologically intuitive and is seen as having great potential, but is an emerging field with few established statistical methods. RESULTS: We transform gene expression profiles into binary gene set profiles by discretizing results of gene set enrichment analyses and apply a new iterative bi-clustering algorithm (iBBiG) to identify groups of gene sets that are coordinately associated with groups of phenotypes across multiple studies. iBBiG is optimized for meta-analysis of large numbers of diverse genomics data that may have unmatched samples. It does not require prior knowledge of the number or size of clusters. When applied to simulated data, it outperforms commonly used clustering methods, discovers overlapping clusters of diverse sizes and is robust in the presence of noise. We apply it to meta-analysis of breast cancer studies, where iBBiG extracted novel gene set-phenotype association that predicted tumor metastases within tumor subtypes. AVAILABILITY: Implemented in the Bioconductor package iBBiG CONTACT: aedin@jimmy.harvard.edu.

nEASE: a method for gene ontology subclassification of high-throughput gene expression data.

UNLABELLED: High-throughput technologies can identify genes whose expression profiles correlate with specific phenotypes; however, placing these genes into a biological context remains challenging. To help address this issue, we developed nested Expression Analysis Systematic Explorer (nEASE). nEASE complements traditional gene ontology enrichment approaches by determining statistically enriched gene ontology subterms within a list of genes based on co-annotation. Here, we overview an open-source software version of the nEASE algorithm. nEASE can be used either stand-alone or as part of a pathway discovery pipeline. AVAILABILITY: nEASE is implemented within the Multiple Experiment Viewer software package available at http://www.tm4.org/mev. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RNA-Seq analysis in MeV.

SUMMARY: RNA-Seq is an exciting methodology that leverages the power of high-throughput sequencing to measure RNA transcript counts at an unprecedented accuracy. However, the data generated from this process are extremely large and biologist-friendly tools with which to analyze it are sorely lacking. MultiExperiment Viewer (MeV) is a Java-based desktop application that allows advanced analysis of gene expression data through an intuitive graphical user interface. Here, we report a significant enhancement to MeV that allows analysis of RNA-Seq data with these familiar, powerful tools. We also report the addition to MeV of several RNA-Seq-specific functions, addressing the differences in analysis requirements between this data type and traditional gene expression data. These tools include automatic conversion functions from raw count data to processed RPKM or FPKM values and differential expression detection and functional annotation enrichment detection based on published methods.

Functional classification analysis of somatically mutated genes in human breast and colorectal cancers.

A recent study published by Sjoblom and colleagues [T. Sjoblom, S. Jones, L.D. Wood, D.W. Parsons, J. Lin, T.D. Barber, D. Mandelker, R.J. Leary, J. Ptak, N. Silliman, S. Szabo, P. Buckhaults, C. Farrell, P. Meeh, S.D. Markowitz, J. Willis, D. Dawson, J.K. Willson, A.F. Gazdar, J. Hartigan, L. Wu, C. Liu, G. Parmigiani, B.H. Park, K.E. Bachman, N. Papadopoulos, B. Vogelstein, K.W. Kinzler, V.E. Velculescu, The consensus coding sequences of human breast and colorectal cancers. Science 314 (2006) 268-274.] performed comprehensive sequencing of 13,023 human genes and identified mutations in genes specific to breast and colorectal tumors, providing insight into organ-specific tumor biology. Here we present a systematic analysis of the functional classifications of Sjoblom's "CAN" genes, a subset of these validated mutant genes, that identifies novel organ-specific biological themes and molecular pathways associated with disease-specific etiology. This analysis links four somatically mutated genes associated with diverse oncological types to colorectal and breast cancers through established TGF-beta1-regulated interactions, revealing mechanistic differences in these cancers and providing potential diagnostic and therapeutic targets.

Novel tricyclic glycal-based TRIB1 inducers that reprogram LDL metabolism in hepatic cells.

Increased expression of the Tribbles pseudokinase 1 gene (TRIB1) is associated with lower plasma levels of LDL cholesterol and triglycerides, higher levels of HDL cholesterol and decreased risk of coronary artery disease and myocardial infarction. We identified a class of tricyclic glycal core-based compounds that upregulate TRIB1 expression in human HepG2 cells and phenocopy the effects of genetic TRIB1 overexpression as they inhibit expression of triglyceride synthesis genes and ApoB secretion in cells. In addition to predicted effects related to downregulation of VLDL assembly and secretion these compounds also have unexpected effects as they upregulate expression of LDLR and stimulate LDL uptake. This activity profile is unique and favorably differs from profiles produced by statins or other lipoprotein targeting therapies. BRD8518, the initial lead compound from the tricyclic glycal class, exhibited stereochemically dependent activity and the potency far exceeding previously described benzofuran BRD0418. Gene expression profiling of cells treated with BRD8518 demonstrated the anticipated changes in lipid metabolic genes and revealed a broad stimulation of early response genes. Consistently, we found that BRD8518 activity is MEK1/2 dependent and the treatment of HepG2 cells with BRD8518 stimulates ERK1/2 phosphorylation. In agreement with down-regulation of genes controlling triglyceride synthesis and assembly of lipoprotein particles, the mass spectrometry analysis of cell extracts showed reduced rate of incorporation of stable isotope labeled glycerol into triglycerides in BRD8518 treated cells. Furthermore, we describe medicinal chemistry efforts that led to identification of BRD8518 analogs with enhanced potency and pharmacokinetic properties suitable for in vivo studies.

TM4 microarray software suite.

Powerful specialized software is essential for managing, quantifying, and ultimately deriving scientific insight from results of a microarray experiment. We have developed a suite of software applications, known as TM4, to support such gene expression studies. The suite consists of open-source tools for data management and reporting, image analysis, normalization and pipeline control, and data mining and visualization. An integrated MIAME-compliant MySQL database is included. This chapter describes each component of the suite and includes a sample analysis walk-through.

Therapeutic implications of GIPC1 silencing in cancer.

GIPC1 is a cytoplasmic scaffold protein that interacts with numerous receptor signaling complexes, and emerging evidence suggests that it plays a role in tumorigenesis. GIPC1 is highly expressed in a number of human malignancies, including breast, ovarian, gastric, and pancreatic cancers. Suppression of GIPC1 in human pancreatic cancer cells inhibits in vivo tumor growth in immunodeficient mice. To better understand GIPC1 function, we suppressed its expression in human breast and colorectal cancer cell lines and human mammary epithelial cells (HMECs) and assayed both gene expression and cellular phenotype. Suppression of GIPC1 promotes apoptosis in MCF-7, MDA-MD231, SKBR-3, SW480, and SW620 cells and impairs anchorage-independent colony formation of HMECs. These observations indicate GIPC1 plays an essential role in oncogenic transformation, and its expression is necessary for the survival of human breast and colorectal cancer cells. Additionally, a GIPC1 knock-down gene signature was used to interrogate publically available breast and ovarian cancer microarray datasets. This GIPC1 signature statistically correlates with a number of breast and ovarian cancer phenotypes and clinical outcomes, including patient survival. Taken together, these data indicate that GIPC1 inhibition may represent a new target for therapeutic development for the treatment of human cancers.