- Open Access
SBEAMS-Microarray: database software supporting genomic expression analyses for systems biology
BMC Bioinformatics volume 7, Article number: 286 (2006)
The biological information in genomic expression data can be understood, and computationally extracted, in the context of systems of interacting molecules. The automation of this information extraction requires high throughput management and analysis of genomic expression data, and integration of these data with other data types.
SBEAMS-Microarray, a module of the open-source Systems Biology Experiment Analysis Management System (SBEAMS), enables MIAME-compliant storage, management, analysis, and integration of high-throughput genomic expression data. It is interoperable with the Cytoscape network integration, visualization, analysis, and modeling software platform.
SBEAMS-Microarray provides end-to-end support for genomic expression analyses for network-based systems biology research.
The extraction of biological information from high throughput genomic expression data is a fundamentally network-based systems biology problem . Complex cell properties such as pathogenicity, growth control, and metabolic capabilities arise from networks of molecular interactions. Control of such cell properties involves gene activity at multiple levels, including mRNA and protein levels, and molecular modifications, localization, and interactions. The computational integration of disparate network-data types and the application of network-analysis algorithms enables the extraction of information that is not contained in individual network elements or single data types.
Systems biology research is characterized by the development and application of technologies enabling quantitative measurements that are genomic in scale. Data collected at multiple levels of gene activity are integrated, analyzed, and modeled to suggest further experiments in an iterative cycle of discovery. Microarray-based genomic expression measurement is the technology that comes closest to meeting the demands of systems biology research. Microarray-based measurement technologies are mature relative to proteomics and other genome-scale molecular interrogation methods. A major current bottleneck in systems biology is the extraction of biological information from genomic expression data integrated with other data types.
Before the development of SBEAMS-Microarray, existing microarray data-analysis software packages were evaluated based on desired attributes for the support of microarray-based systems biology research. These attributes included:
Free open-source availability
Compliance with the Minimum Information About a Microarray Experiment (MIAME) data content specification .
Flexible application of multiple microarray-data analysis methods
Strong support for large-scale datasets derived from the high-probe-density Affymetrix microarray platform
Integration of genomic expression data with other data types such as proteomics data
Interoperability with network visualization, analysis, and modeling software, for example, the Cytoscape platform .
However no single package had all of these attributes. Several packages [5–8] were developed for two-color spotted arrays, and lack support for high-probe-density Affymetrix data. Other database software projects such as BASE 2  intend to support Affymetrix data, but have not yet been released at the time of this publication. Other packages [10, 11] require proprietary software and hardware, in these cases an Oracle database and Sun server. Finally, none of the evaluated packages [5–8, 10–12] support integration of genomic expression data with other data types. Therefore, we developed the Microarray module of SBEAMS, the Systems Biology Experiment Analysis Management System . SBEAMS-Microarray combines all of the above attributes, making SBEAMS-Microarray uniquely advantageous for network-based systems biology research.
The SBEAMS-Microarray module is built using SBEAMS, a software and database framework for collecting, storing, accessing and analyzing data from different types of experimental data . SBEAMS combines a relational database management system (RDBMS) back end, a collection of tools to store, manage, and query experiment information and results in the RDBMS, a web front end for querying the database and providing integrated access to remote data sources, and an interface to other data processing and analysis programs (Figure 1).
SBEAMS, including the SBEAMS-Microarray module, uses a web-based client-server software model. Thus, the SBEAMS software runs and is updated only on a central server. Also, computationally intensive tasks are handled on the server end. On the client end, the user needs only minimal computing power, a modern web browser with Java Web Start  installed, and a network connection to the SBEAMS server. Users connect to SBEAMS via HTTPS to the web server on which SBEAMS is installed. Perl CGI scripts use the SBEAMS API to create a web interface to the back-end database. The Perl DBI module is used for database connectivity. The system is designed to support any type of RDBMS for which a DBI driver is available. SBEAMS-Microarray is known to work with Microsoft SQL Server, Sybase Adaptive Server Enterprise, and MySQL, the most popular open-source RDBMS. Work is underway to add support for PostgreSQL, another popular open-source RDBMS.
The database schema definitions are provided in a database-independent format with the SBEAMS distribution. A script that generates data definition language (DDL) commands (i.e. CREATE TABLE, etc.) for the schema for several different RDBMSs is also provided. Its use is described in the installation instructions. A schema diagram (Additional file 1) depicts the SBEAMS-Microarray database module schema, including tables containing information about samples, arrays and associated files, intensity measurements, probe set annotations, analysis runs, and the final expression results. Additional tables containing information about SBEAMS users, work groups, permissions, projects, species, etc. are part of the SBEAMS Core schema and are not reproduced here.
SBEAMS-Microarray and installation instructions can be downloaded as part of the latest release of SBEAMS at the SBEAMS-Microarray web site . Alternatively, the Subversion  version control software can be used to obtain the current development version of SBEAMS from the Subversion repository. Subversion installation allows the software to be updated easily as new additions and changes are created, as well as allowing users who continue to extend SBEAMS through software development to contribute their work. Evaluation of SBEAMS prior to installation is available through a demonstration instance of the software at the SBEAMS-Microarray web site . This site additionally provides access to two mailing lists, one for general questions and discussion, and the other for developers. SBEAMS and SBEAMS-Microarray are available under the terms of the GNU General Public License version 2.
Access to the interfaces and data in SBEAMS-Microarray is controlled through a comprehensive security model. Each user must log in to gain access to the database, and can belong to any number of work groups. A work group defines a set of privileges for each user in the work group. The software provides default administrative and user roles. Non-administrative users can create projects, and have full permissions for projects they own. By default, a user cannot view projects belonging to other users. However a user can grant varying levels of access to other users or work groups to facilitate sharing of data and analysis results. Specific result sets are shared easily by emailing hyperlinks to users that have access to a given project.
After SBEAMS-Microarray is installed, administration is handled primarily through the web interface, with a set of tools available only to users with administrative privileges. These tools, accessible via menu choices visible only to such users, provide functionality such as creating new types of arrays, modifying records for users other than themselves, and deleting incorrect records. Other administrative functions, such as loading new array annotations, are handled at the command-line. These array annotations serve as the primary source of probe information during analyses conducted through SBEAMS-Microarray, and are loaded from quarterly updates that Affymetrix produces. Command-line functions are authenticated to enforce user permissions.
Results and discussion
SBEAMS-Microarray primarily supports the high-probe-density Affymetrix platform. All current Affymetrix gene-expression microarrays are supported, and all future Affymetrix arrays should be supported as well, given the generic mechanism for loading information about new array types. Two-color microarray support exists, however it is disabled by default, given that two-color microarray quantitation and annotation formats are not standardized. Modifications could be made to the software to support additional two-color microarray data formats, allowing the two-color portion of the software to be used.
SBEAMS-Microarray functions seamlessly with the high-throughput high-probe-density Affymetrix microarray platform and GCOS system software  as shown in Figure 1. Microarray facility staff or end users enter information about their samples into the Affymetrix GCOS software. After scanning, data extraction, and initial processing, GCOS exports these data, with all of the raw data files, to an output data directory as a standard MAGE-ML file. SBEAMS-Microarray periodically scans the output data directory for new data sets, and automatically loads any complete sets into the SBEAMS database. The database stores information about each microarray and pointers to the locations of data files, which are stored in an SBEAMS-managed file tree. After automated data loading, users can access an overview of all their arrays in a web-based user interface, with the ability to view or download raw data files and data quality reports. Users can edit or add sample information and annotations to comply with the MIAME data content specifications .
SBEAMS-Microarray supports simple queries allowing quick access to data for genes of interest. Users specify a search string and select the arrays from which they want to see expression data and detection calls. The results are presented in a matrix of data with colored visual cues. MAS 5.0 signal values and detection calls  are used because normalization can be done on each microarray independently, whereas other normalization methods such as RMA  depend on the normalization group.
An advanced SQL query tool also exists, allowing more search parameters and resulting in tabular data output.
SBEAMS-Microarray incorporates widely used open-source genomics softwares, and thereby supports the flexible application of multiple microarray-data processing and analysis methods. To perform processing tasks including background correction, normalization, probe set summarization and differential expression testing, we integrated the BioConductor  open-source web interface package, webbioc, into SBEAMS-Microarray. The webbioc package implements several processing methods, including RMA , GC-RMA , VSN , MAS 5.0 , dChip , Quantile normalization  and Qspline normalization . Processing may occur on the SBEAMS server or may be submitted to a batch scheduler on a computer cluster. Email notification of completed processing jobs is available. A data-processing summary page provides several diagnostic plots to help in the identification of microarrays that failed or are inconsistent with the rest of the data set, and links to download the processed data. After processing of a specific data set, SBEAMS-Microarray supports differential-expression testing with three optional methods: simple ratio analysis, t-test, and the SAM false-discovery-rate method . SAM and t-test are available for data sets with replicate experiments, and produce test statistics for each probe set as well as a view of the probe sets with the highest scores and accompanying annotations. After a data analysis has been performed, the user can elect to load the results into the database where it can be stored and queried, or integrated with other data types.
SBEAMS-Microarray incorporates the MultiExperiment Viewer (MeV), developed at The Institute for Genomic Research. MeV provides numerous statistical tests, classification methods, and clustering algorithms  that extend the analytical capabilities of SBEAMS-Microarray. From the SBEAMS-Microarray web interface, MeV is launched using Java Web Start, so that users do not need to install MeV, and are immediately presented with their data in the MeV environment.
SBEAMS is modular in design to allow the integrated storage and access of disparate types of experiments and data, for example, microarray and proteomics experiments, molecular interaction data, and gene annotations. This integrated system is a consistent framework that combines a RDBMS back end, and a web front end providing integrated access to the data. For example, from the SBEAMS-Microarray GetExpression interface, queries can be made for gene annotations of interest or by defining threshold levels of metrics for statistical significance of expression change across one or more user-specified microarray experiments. The results can be viewed within the SBEAMS web interface, or exported in Excel, CSV, TSV, or XML formats, or accessed programmatically via HTTPS. Data sets may originate in-house or be imported from external sources. Currently, development of SBEAMS is driven mainly by the SBEAMS-Microarray and SBEAMS-Proteomics projects, with multiple modules in early stages of development.
SBEAMS-Microarray is interoperable with Cytoscape software. The results of data analyses can be loaded directly to the Cytoscape environment, launched from the SBEAMS environment. Cytoscape is an open-source bioinformatics software platform for visualizing molecular interaction networks and integrating these networks with gene expression data, proteomics data, gene annotations, and other data . A wide variety of additional functionalities are available as Cytoscape plugins. Plugins implement integrated network analyses, connection with outside databases and tools, and modeling capabilities. In SBEAMS-Microarray, query results from GetExpression can be loaded directly into Cytoscape via Java Web Start, as described above for MeV. With the data in Cytoscape, users can load molecular interaction data, annotation data, and a wide variety of other data, to generate integrated networks. Data types can be loaded directly as files, or imported from outside databases using plugins like InteractionFetcher .
Interoperability with Cytoscape enables automated data integration and subsequent network-based analyses to extract information that is not present in any one data type. For example, the Biomodules plugin  implements methods for the computational identification of groups of interacting proteins performing some collective function (modules) in integrated networks of genomic expression data, molecular interaction data, and gene annotation data. Prinz et al.  applied these methods to discover and experimentally validate molecular insights on the regulation of yeast cell differentiation from the familiar yeast form to the filamentous-invasive form.
SBEAMS-Microarray enables investigators to store, manage, analyze, and integrate genomic expression data for systems biology research projects. Investigators begin by performing a microarray experiment to answer questions about their biological system of interest. Once the primary data have been obtained, investigators log in through the SBEAMS-Microarray web interface to see their data that has been automatically loaded into the database. Before beginning analysis, it is advisable to ensure that the data are acceptable by viewing various quality control metrics and diagnostic plots provided by the software. Investigators may choose to annotate their data further, by providing greater details on the biological samples hybridized to their microarrays to aid others involved with analysis of their data and in compliance with MIAME standards . Once satisfied with their data quality, investigators may begin to gather biological information by using the querying interfaces to inspect the expression patterns of one or more genes with known (or potentially interesting) responses under their experimental conditions. After establishing confidence in their expression data, investigators use the data analysis pipeline. After applying one of several optional normalization methods, investigators have the option to launch seamlessly the MultiExperiment Viewer , which provides multiple methods to cluster and visualize the data. A second option, producing results that ultimately will lead to data integration and network visualization and analysis in Cytoscape, is differential expression analysis. Investigators choose parameters for their analysis, including a statistical method, groups of biological replicates to be compared, and thresholds for statistically significant differential expression. After analysis, result tables show genes with the greatest and most significant expression differences. Users have the option to store the results in the database. These tables of differentially expressed genes are themselves informative, but will provide further system-level insights when integrated with other biological data. From stored analyses, sets of genes with their differential expression values can be loaded directly into Cytoscape launched from SBEAMS. Expression data can be mapped to interpolated colors of nodes representing differentially expressed genes. Investigators then may begin to use Cytoscape to explore their expression data in the context of biological network information. One method is to employ the InteractionFetcher  plugin to find and integrate data on interactions among the differentially expressed genes and their products. As noted above (Data integration), the integration of these interaction data, and other data types, with the expression data produces networks that investigators analyze to find relationships that are not contained in either data type alone. Cytoscape enables integration of many data types (e.g., proteomics data, annotation data, etc.), customizable visualization of these integrated data, and computational analysis of integrated networks to extract system-level information on the questions motivating the study (e.g., refs  and ). Uses of the software are detailed in the SBEAMS-Microarray User Guide (Additional file 2).
SBEAMS is an open-source software project. Its design is intended to facilitate further development. SBEAMS allows multiple separate instances of the software to be installed on the same machine, so that software developers may have one or more developmental versions where they improve or extend SBEAMS-Microarray, without interfering with the production instance of the software. Developmental improvements are tested and eventually added to the code repository and rolled out to the production instance.
Implementation enhancements are planned. Active development is underway to support other backend database software such as PostgreSQL, so that SBEAMS – Microarray can be implemented on the local RDBMS of choice.
Support for MAGE-OM/ML database standards  is planned. Use of these standards will allow interoperability of SBEAMS-Microarray with other MAGE-compliant software and allow creation of MAGE-ML documents for submission of experimental data to public repositories such as ArrayExpress .
A major goal for SBEAMS-Microarray development is the addition of support for more types of microarrays and experimental assays. Development will be required to support the emergence of new arrays and platforms, particularly with respect to integrating results across different generations, and possibly platforms, of microarrays. Support for experiments based on genome-tiling microarrays is a priority. These microarrays enable high-throughput genome-scale investigations of alternative splicing, non-coding RNA levels, protein-DNA interaction, and comparative genomics. These assays require new data analysis methods, e.g., , which will be incorporated.
Additional tools for integrating analysis results from separate SBEAMS modules are planned. Currently SBEAMS allows for storage, access and analysis of disparate data types in their respective modules within the SBEAMS framework. External tools such as Cytoscape must be used to integrate these multiple data types. An interface for combining microarray, proteomics and interaction data within SBEAMS is currently under development.
SBEAMS-Microarray is a useful tool for both a microarray facility and its diverse user community. It is uniquely strong in its flexible incorporation of multiple data analysis methods and supporting softwares, its support of data standards, its open-source availability, and its support for data integration and network analyses. SBEAMS-Microarray is a key module in the SBEAMS database system, which has several other modules (e.g., SBEAMS-Proteomics) allowing for incorporation of disparate data types into a single framework. In combination with network-analysis tools like Cytoscape, it provides end-to-end support for systems biology research projects involving high-throughput genomic expression analysis.
Availability and requirements
Project name: SBEAMS-Microarray
Project home page: http://www.sbeams.org/Microarray
○ Application server: Linux/UNIX
○ RDBMS: Windows or Linux/UNIX
○ Apache web server
○ R 2.1.0 through 2.3.0 with Bioconductor 1.6 through 1.8
○ Microsoft SQL Server, Sybase Adaptive Server Enterprise or MySQL
○ FreeTDS 0.63
○ libgd 2.0.33
License: GNU General Public License version 2
Restrictions to use by non-academics: None
Galitski T: Molecular networks in model systems. Annu Rev Genomics Hum Genet 2004, 5: 177–187. 10.1146/annurev.genom.5.061903.180053
Ideker T, Galitski T, Hood L: A new approach to decoding life: systems biology. Annu Rev Genomics Hum Genet 2001, 2: 343–372. 10.1146/annurev.genom.2.1.343
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet 2001, 29(4):365–371. 10.1038/ng1201-365
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003, 13(11):2498–2504. 10.1101/gr.1239303
Saal LH, Troein C, Vallon-Christersson J, Gruvberger S, Borg A, Peterson C: BioArray Software Environment (BASE): a platform for comprehensive management and analysis of microarray data. Genome Biol 2002, 3(8):SOFTWARE0003. 10.1186/gb-2002-3-8-software0003
Killion PJ, Sherlock G, Iyer VR: The Longhorn Array Database (LAD): an open-source, MIAME compliant implementation of the Stanford Microarray Database (SMD). BMC Bioinformatics 2003, 4: 32. 10.1186/1471-2105-4-32
Maurer M, Molidor R, Sturn A, Hartler J, Hackl H, Stocker G, Prokesch A, Scheideler M, Trajanoski Z: MARS: microarray analysis, retrieval, and storage system. BMC Bioinformatics 2005, 6(1):101. 10.1186/1471-2105-6-101
Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, Popov D, Ryltsov A, Kostukovich E, Borisovsky I, Liu Z, Vinsavich A, Trush V, Quackenbush J: TM4: a free, open-source system for microarray data management and analysis. Biotechniques 2003, 34(2):374–378.
[http://base.thep.lu.se/] http://base.thep.lu.se/: BASE Project Site.
Ball CA, Awad IA, Demeter J, Gollub J, Hebert JM, Hernandez-Boussard T, Jin H, Matese JC, Nitzberg M, Wymore F, Zachariah ZK, Brown PO, Sherlock G: The Stanford Microarray Database accommodates additional microarray platforms and data formats. Nucleic Acids Res 2005, 33(Database issue):D580–2. 10.1093/nar/gki006
Theilhaber J, Ulyanov A, Malanthara A, Cole J, Xu D, Nahf R, Heuer M, Brockel C, Bushnell S: GECKO: a complete large-scale gene expression analysis platform. BMC Bioinformatics 2004, 5(1):195. 10.1186/1471-2105-5-195
[https://genes.med.virginia.edu] https://genes.med.virginia.edu: GEOSS Home Page.
[http://www.sbeams.org/] http://www.sbeams.org/: Systems Biology Experiment Analysis Management System.
[http://java.sun.com] http://java.sun.com: Java Technology.
[http://www.sbeams.org/Microarray] http://www.sbeams.org/Microarray: SBEAMS - Microarray.
[http://subversion.tigris.org] http://subversion.tigris.org: subversion.tigris.org.
[http://www.affymetrix.com] http://www.affymetrix.com: Affymetrix.
Hubbell E, Liu WM, Mei R: Robust estimators for expression analysis. Bioinformatics 2002, 18(12):1585–1592. 10.1093/bioinformatics/18.12.1585
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003, 4(2):249–264. 10.1093/biostatistics/4.2.249
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 2004, 5(10):R80. 10.1186/gb-2004-5-10-r80
Wu Z, Irizarry RA, Gentleman R, Murillo FM, Spencer F: A Model Based Background Adjustment for Oligonucleotide Expression Arrays. In Johns Hopkins University, Dept of Biostatistics Working Papers. Baltimore, MD ; 2004.
Huber W, von Heydebreck A, Sultmann H, Poustka A, Vingron M: Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics 2002, 18 Suppl 1: S96–104.
Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci U S A 2001, 98(1):31–36. 10.1073/pnas.011404098
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003, 19(2):185–193. 10.1093/bioinformatics/19.2.185
Workman C, Jensen LJ, Jarmer H, Berka R, Gautier L, Nielser HB, Saxild HH, Nielsen C, Brunak S, Knudsen S: A new non-linear normalization method for reducing variability in DNA microarray experiments. Genome Biol 2002, 3(9):research0048. 10.1186/gb-2002-3-9-research0048
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A 2001, 98(9):5116–5121. 10.1073/pnas.091062498
Reiss DJ, Avila-Campillo I, Thorsson V, Schwikowski B, Galitski T: Tools enabling the elucidation of molecular pathways active in human disease: application to Hepatitis C virus infection. BMC Bioinformatics 2005, 6(1):154. 10.1186/1471-2105-6-154
Prinz S, Avila-Campillo I, Aldridge C, Srinivasan A, Dimitrov K, Siegel AF, Galitski T: Control of yeast filamentous-form growth by modules in an integrated molecular network. Genome Res 2004, 14(3):380–390. 10.1101/gr.2020604
Spellman PT, Miller M, Stewart J, Troup C, Sarkans U, Chervitz S, Bernhart D, Sherlock G, Ball C, Lepage M, Swiatek M, Marks WL, Goncalves J, Markel S, Iordan D, Shojatalab M, Pizarro A, White J, Hubley R, Deutsch E, Senger M, Aronow BJ, Robinson A, Bassett D, Stoeckert CJJ, Brazma A: Design and implementation of microarray gene expression markup language (MAGE-ML). Genome Biol 2002, 3(9):RESEARCH0046. 10.1186/gb-2002-3-9-research0046
Brazma A, Parkinson H, Sarkans U, Shojatalab M, Vilo J, Abeygunawardena N, Holloway E, Kapushesky M, Kemmeren P, Lara GG, Oezcimen A, Rocca-Serra P, Sansone SA: ArrayExpress--a public repository for microarray gene expression data at the EBI. Nucleic Acids Res 2003, 31(1):68–71. 10.1093/nar/gkg091
Royce TE, Rozowsky JS, Bertone P, Samanta M, Stolc V, Weissman S, Snyder M, Gerstein M: Issues in the analysis of oligonucleotide tiling microarrays for transcript mapping. Trends Genet 2005, 21(8):466–475. 10.1016/j.tig.2005.06.007
The authors gratefully acknowledge the contributions of K Dimitrov, C Emswiler, L Hood, G Lake, S Lasky, L Mendoza, and M Whiting. This project was supported in part by NIH/NIGMS Grant P20 GM64361, Intelligent Information Systems for Systems Biology, to L Hood. T Galitski is a recipient of a Burroughs Wellcome Fund Career Award in the Biomedical Sciences.
BM contributed to software design, did software testing and debugging, and drafted the manuscript. EWD conceived the software project, designed and coded software, and drafted the manuscript. PM designed and coded software. DC coded and debugged software. MHJ coded and debugged software. TG supervised the project, participated in software design, and revised the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.