Benno Schwikowski

Integration of biological networks and gene expression data using Cytoscape

Cytoscape is a free software package for visualizing, modeling and analyzing molecular and genetic interaction networks. This protocol explains how to use Cytoscape to analyze the results of mRNA expression profiling, and other functional genomics and proteomics experiments, in the context of an interaction network obtained for genes of interest. Five major steps are described: (i) obtaining a gene or protein network, (ii) displaying the network using layout algorithms, (iii) integrating with gene expression and other functional attributes, (iv) identifying putative complexes and functional modules and (v) identifying enriched Gene Ontology annotations in the network. These steps provide a broad sample of the types of analyses performed by Cytoscape.

A network of protein–protein interactions in yeast

A global analysis of 2,709 published interactions between proteins of the yeast Saccharomyces cerevisiae has been performed, enabling the establishment of a single large network of 2,358 interactions among 1,548 proteins. Proteins of known function and cellular location tend to cluster together, with 63% of the interactions occurring between proteins with a common functional assignment and 76% occurring between proteins found in the same subcellular compartment. Possible functions can be assigned to a protein based on the known functions of its interacting partners. This approach correctly predicts a functional category for 72% of the 1,393 characterized proteins with at least one partner of known function, and has been applied to predict functions for 364 previously uncharacterized proteins.

Condition-Dependent Transcriptome Reveals High-Level Regulatory Architecture in Bacillus subtilis

Outside In Acquisition and analysis of large data sets promises to move us toward a greater understanding of the mechanisms by which biological systems are dynamically regulated to respond to external cues. Now, two papers explore the responses of a bacterium to changing nutritional conditions (see the Perspective by Chalancon et al.). Nicolas et al. (p. 1103) measured transcriptional regulation for more than 100 different conditions. Greater amounts of antisense RNA were generated than expected and appeared to be produced by alternative RNA polymerase targeting subunits called sigma factors. One transition, from malate to glucose as the primary nutrient, was studied in more detail by Buescher et al. (p. 1099) who monitored RNA abundance, promoter activity in live cells, protein abundance, and absolute concentrations of intracellular and extracellular metabolites. In this case, the bacteria responded rapidly and largely without transcriptional changes to life on malate, but only slowly adapted to use glucose, a shift that required changes in nearly half the transcription network. These data offer an initial understanding of why certain regulatory strategies may be favored during evolution of dynamic control systems. A horizontal analysis reveals the breadth of genes turned on and off as nutrients change. Bacteria adapt to environmental stimuli by adjusting their transcriptomes in a complex manner, the full potential of which has yet to be established for any individual bacterial species. Here, we report the transcriptomes of Bacillus subtilis exposed to a wide range of environmental and nutritional conditions that the organism might encounter in nature. We comprehensively mapped transcription units (TUs) and grouped 2935 promoters into regulons controlled by various RNA polymerase sigma factors, accounting for ~66% of the observed variance in transcriptional activity. This global classification of promoters and detailed description of TUs revealed that a large proportion of the detected antisense RNAs arose from potentially spurious transcription initiation by alternative sigma factors and from imperfect control of transcription termination.

Global Network Reorganization During Dynamic Adaptations of Bacillus subtilis Metabolism

Outside In Acquisition and analysis of large data sets promises to move us toward a greater understanding of the mechanisms by which biological systems are dynamically regulated to respond to external cues. Now, two papers explore the responses of a bacterium to changing nutritional conditions (see the Perspective by Chalancon et al.). Nicolas et al. (p. 1103) measured transcriptional regulation for more than 100 different conditions. Greater amounts of antisense RNA were generated than expected and appeared to be produced by alternative RNA polymerase targeting subunits called sigma factors. One transition, from malate to glucose as the primary nutrient, was studied in more detail by Buescher et al. (p. 1099) who monitored RNA abundance, promoter activity in live cells, protein abundance, and absolute concentrations of intracellular and extracellular metabolites. In this case, the bacteria responded rapidly and largely without transcriptional changes to life on malate, but only slowly adapted to use glucose, a shift that required changes in nearly half the transcription network. These data offer an initial understanding of why certain regulatory strategies may be favored during evolution of dynamic control systems. A vertical analysis reveals that a simple switch of one food for another evokes changes at many levels. Adaptation of cells to environmental changes requires dynamic interactions between metabolic and regulatory networks, but studies typically address only one or a few layers of regulation. For nutritional shifts between two preferred carbon sources of Bacillus subtilis, we combined statistical and model-based data analyses of dynamic transcript, protein, and metabolite abundances and promoter activities. Adaptation to malate was rapid and primarily controlled posttranscriptionally compared with the slow, mainly transcriptionally controlled adaptation to glucose that entailed nearly half of the known transcription regulation network. Interactions across multiple levels of regulation were involved in adaptive changes that could also be achieved by controlling single genes. Our analysis suggests that global trade-offs and evolutionary constraints provide incentives to favor complex control programs.

Alignment of LC‐MS images, with applications to biomarker discovery and protein identification

LC‐MS‐based approaches have gained considerable interest for the analysis of complex peptide or protein mixtures, due to their potential for full automation and high sampling rates. Advances in resolution and accuracy of modern mass spectrometers allow new analytical LC‐MS‐based applications, such as biomarker discovery and cross‐sample protein identification. Many of these applications compare multiple LC‐MS experiments, each of which can be represented as a 2‐D image. In this article, we survey current approaches to LC‐MS image alignment. LC‐MS image alignment corrects for experimental variations in the chromatography and represents a computational key technology for the comparison of LC‐MS experiments. It is a required processing step for its two major applications: biomarker discovery and protein identification. Along with descriptions of the computational analysis approaches, we discuss their relative merits and potential pitfalls.

GOlorize: a Cytoscape plug-in for network visualization with Gene Ontology-based layout and coloring

UNLABELLED We have implemented a graph layout algorithm that exposes Gene Ontology (GO) class structure on the network nodes. It can be used in conjunction with BiNGO plug-in to Cytoscape, which finds the GO categories over-represented in a given network. Our plug-in, named GOlorize, first highlights the class members with category-specific color-coding and then constructs an enhanced visualization of the network using a class-directed layout algorithm. AVAILABILITY http://www.cytoscape.org/plugins2.php. SUPPLEMENTARY INFORMATION Installation instructions and tutorial at http://www.cytoscape.org/plugins/GOlorize/GOlorizeUserGuide.pdf.

Universal DNA tag systems: a combinatorial design scheme

Custom-designed DNA arrays offer the possibility of simultaneously monitoring thousands of hybridization reactions These arrays show great potential for many medical and scientific applications such as polymorphism analysis and genotyping. Relatively high costs are associated with the need to specifically design and synthesize problem specific arrays. Recently, an alternative approach was suggested that utilizes fixed, universal arrays. This approach presents an interesting design problem—the arrays should contain as many probes as possible, while minimizing experimental errors caused by cross-hybridization. We use a simple thermodynamic model to cast this design problem in a formal mathematical framework. Employing new combinatorial ideas, we derive an efficient construction for the design problem, and prove that our construction is near-optimal.

Algorithms for phylogenetic footprinting.

Phylogenetic footprinting is a technique that identifies regulatory elements by finding unusually well conserved regions in a set of orthologous noncoding DNA sequences from multiple species. We introduce a new motif-finding problem, the Substring Parsimony Problem, which is a formalization of the ideas behind phylogenetic footprinting, and we present an exact dynamic programming algorithm to solve it. We then present a number of algorithmic optimizations that allow our program to run quickly on most biologically interesting datasets. We show how to handle data sets in which only an unknown subset of the sequences contains the regulatory element. Finally, we describe how to empirically assess the statistical significance of the motifs found. Each technique is implemented and successfully identifies a number of known binding sites, as well as several highly conserved but uncharacterized regions. The program is available at http://bio.cs.washington.edu/software.html.

Functional Analysis via Standardized Whole-Blood Stimulation Systems Defines the Boundaries of a Healthy Immune Response to Complex Stimuli

Standardization of immunophenotyping procedures has become a high priority. We have developed a suite of whole-blood, syringe-based assay systems that can be used to reproducibly assess induced innate or adaptive immune responses. By eliminating preanalytical errors associated with immune monitoring, we have defined the protein signatures induced by (1) medically relevant bacteria, fungi, and viruses; (2) agonists specific for defined host sensors; (3) clinically employed cytokines; and (4) activators of T cell immunity. Our results provide an initial assessment of healthy donor reference values for induced cytokines and chemokines and we report the failure to release interleukin-1α as a common immunological phenotype. The observed naturally occurring variation of the immune response may help to explain differential susceptibility to disease or response to therapeutic intervention. The implementation of a general solution for assessment of functional immune responses will help support harmonization of clinical studies and data sharing.

Assessing Bias in Experiment Design for Large Scale Mass Spectrometry-based Quantitative Proteomics

Mass spectrometry-based proteomics holds great promise as a discovery tool for biomarker candidates in the early detection of diseases. Recently much emphasis has been placed upon producing highly reliable data for quantitative profiling for which highly reproducible methodologies are indispensable. The main problems that affect experimental reproducibility stem from variations introduced by sample collection, preparation, and storage protocols and LC-MS settings and conditions. On the basis of a formally precise and quantitative definition of similarity between LC-MS experiments, we have developed Chaorder, a fully automatic software tool that can assess experimental reproducibility of sets of large scale LC-MS experiments. By visualizing the similarity relationships within a set of experiments, this tool can form the basis of systematic quality control and thus help assess the comparability of mass spectrometry data over time, across different laboratories, and between instruments. Applying Chaorder to data from multiple laboratories and a range of instruments, experimental protocols, and sample complexities revealed biases introduced by the sample processing steps, experimental protocols, and instrument choices. Moreover we show that reducing bias by correcting for just a few steps, for example randomizing the run order, does not provide much gain in statistical power for biomarker discovery.